Publications

2024

Protein turnover is critical for proteostasis, but turnover quantification is challenging, and even in well-studied E. coli, proteome-wide measurements remain scarce. Here, we quantify the turnover rates of ~3200 E. coli proteins under 13 conditions by combining heavy isotope labeling with complement reporter ion quantification and find that cytoplasmic proteins are recycled when nitrogen is limited. We use knockout experiments to assign substrates to the known cytoplasmic ATP-dependent proteases. Surprisingly, none of these proteases are responsible for the observed cytoplasmic protein degradation in nitrogen limitation, suggesting that a major proteolysis pathway in E. coli remains to be discovered. Lastly, we show that protein degradation rates are generally independent of cell division rates. Thus, we present broadly applicable technology for protein turnover measurements and provide a rich resource for protein half-lives and protease substrates in E. coli, complementary to genomics data, that will allow researchers to study the control of proteostasis.

Phase separation of biomolecules can facilitate their spatiotemporally regulated self-assembly within living cells. Due to the selective yet dynamic exchange of biomolecules across condensate interfaces, condensates can function as reactive hubs by concentrating enzymatic components for faster kinetics. The principles governing this dynamic exchange between condensate phases, however, are poorly understood. In this work, we systematically investigate the influence of client-sticker interactions on the exchange dynamics of protein molecules across condensate interfaces. We show that increasing affinity between a model protein scaffold and its client molecules causes the exchange of protein chains between the dilute and dense phases to slow down and that beyond a threshold interaction strength, this slowdown in exchange becomes substantial. Investigating the impact of interaction symmetry, we found that chain exchange dynamics are also considerably slower when client molecules interact equally with different sticky residues in the protein. The slowdown of exchange is due to a sequestration effect, by which there are fewer unbound stickers available at the interface to which dilute phase chains may attach. These findings highlight the fundamental connection between client-scaffold interaction networks and condensate exchange dynamics.

Phages-viruses that infect bacteria-have evolved over billions of years to overcome bacterial defenses. Temperate phage, upon infection, can "choose" between two pathways: lysis-in which the phage create multiple new phage particles, which are then liberated by cell lysis, and lysogeny-where the phage's genetic material is added to the bacterial DNA and transmitted to the bacterial progeny. It was recently discovered that some phages can read information from the environment related to the density of bacteria or the number of nearby infection attempts. Such information may help phage make the right choice between the two pathways. Here, we develop a theoretical model that allows an infecting phage to change its strategy (i.e. the ratio of lysis to lysogeny) depending on an outside signal, and we find the optimal strategy that maximizes phage proliferation. While phages that exploit extra information naturally win in competition against phages with a fixed strategy, there may be costs to information, e.g. as the necessary extra genes may affect the growth rate of a lysogen or the burst size of new phage for the lysis pathway. Surprisingly, even when phages pay a large price for information, they can still maintain an advantage over phages that lack this information, indicating the high benefit of intelligence gathering in phage-bacteria warfare.

When phage infect their bacterial hosts, they may either lyse the cell and generate a burst of new phage, or lysogenize the bacterium, incorporating the phage genome into it. Phage lysis/lysogeny strategies are assumed to be highly optimized, with the optimal tradeoff depending on environmental conditions. However, in nature, phage of radically different lysis/lysogeny strategies coexist in the same environment, preying on the same bacteria. How can phage preying on the same bacteria coexist if one is more optimal than the other? Here, we address this conundrum within a modeling framework, simulating the population dynamics of communities of phage and their lysogens. We find that coexistence between phage of different lysis/lysogeny strategies is a natural outcome of chaotic population dynamics that arise within sufficiently diverse communities, which ensure no phage is able to absolutely dominate its competitors. Our results further suggest a bet-hedging mechanism at the level of the phage pan-genome, wherein obligate lytic (virulent) strains typically outcompete temperate strains, but also more readily fluctuate to extinction within a local community.

2023

Biomolecular condensates are membraneless organelles formed via phase separation of macromolecules, typically consisting of bond-forming “stickers” connected by flexible “linkers.” Linkers have diverse roles, such as occupying space and facilitating interactions. To understand how linker length relative to other lengths affects condensation, we focus on the pyrenoid, which enhances photosynthesis in green algae. Specifically, we apply coarse-grained simulations and analytical theory to the pyrenoid proteins of Chlamydomonas reinhardtii: the rigid holoenzyme Rubisco and its flexible partner EPYC1. Remarkably, halving EPYC1 linker lengths decreases critical concentrations by tenfold. We attribute this difference to the molecular “fit,” i.e., the number of stickers of EPYC1 that can bind to a single Rubisco given the constraint of EPYC1 linker length. We find an inverse relationship between molecular fit and the tendency of EPYC1 and Rubisco to phase separate. Moreover, by computationally varying Rubisco sticker locations we discover that the naturally occurring sticker locations yield the poorest fit for all EPYC1 linker lengths; thus the natural locations optimize phase separation. Surprisingly, shorter linkers mediate a transition to a gas of rods as Rubisco stickers approach the poles. These findings illustrate how intrinsically disordered proteins affect phase separation through the interplay of molecular length scales.

Bacterial populations are highly adaptive. They can respond to stress and survive in shifting environments. How the behaviours of individual bacteria vary during stress, however, is poorly understood. To identify and characterize rare bacterial subpopulations, technologies for single-cell transcriptional profiling have been developed. Existing approaches show some degree of limitation, for example, in terms of number of cells or transcripts that can be profiled. Due in part to these limitations, few conditions have been studied with these tools. Here we develop massively-parallel, multiplexed, microbial sequencing (M3-seq)-a single-cell RNA-sequencing platform for bacteria that pairs combinatorial cell indexing with post hoc rRNA depletion. We show that M3-seq can profile bacterial cells from different species under a range of conditions in single experiments. We then apply M3-seq to hundreds of thousands of cells, revealing rare populations and insights into bet-hedging associated with stress responses and characterizing phage infection.

The fascinating patterns of collective motion created by autonomously driven particles have fuelled active-matter research for over two decades. So far, theoretical active-matter research has often focused on systems with a fixed number of particles. This constraint imposes strict limitations on what behaviours can and cannot emerge. However, a hallmark of life is the breaking of local cell number conservation by replication and death. Birth and death processes must be taken into account, for example, to predict the growth and evolution of a microbial biofilm, the expansion of a tumour, or the development from a fertilized egg into an embryo and beyond. In this Perspective, we argue that unique features emerge in these systems because proliferation represents a distinct form of activity: not only do the proliferating entities consume and dissipate energy, they also inject biomass and degrees of freedom capable of further self-proliferation, leading to myriad dynamic scenarios. Despite this complexity, a growing number of studies document common collective phenomena in various proliferating soft-matter systems. This generality leads us to propose proliferation as another direction of active-matter physics, worthy of a dedicated search for new dynamical universality classes. Conceptual challenges abound, from identifying control parameters and understanding large fluctuations and nonlinear feedback mechanisms to exploring the dynamics and limits of information flow in self-replicating systems. We believe that, by extending the rich conceptual framework developed for conventional active matter to proliferating active matter, researchers can have a profound impact on quantitative biology and reveal fascinating emergent physics along the way.

It has recently become appreciated that cells self-organize their interiors through the formation of biomolecular condensates. These condensates, typically formed through liquid-liquid phase separation of proteins, nucleic acids, and other biopolymers, exhibit reversible assembly/disassembly in response to changing conditions. Condensates play many functional roles, aiding in biochemical reactions, signal transduction, and sequestration of certain components. Ultimately, these functions depend on the physical properties of condensates, which are encoded in the microscopic features of the constituent biomolecules. In general, the mapping from microscopic features to macroscopic properties is complex, but it is known that near a critical point, macroscopic properties follow power laws with only a small number of parameters, making it easier to identify underlying principles. How far does this critical region extend for biomolecular condensates and what principles govern condensate properties in the critical regime? Using coarse-grained molecular-dynamics simulations of a representative class of biomolecular condensates, we found that the critical regime can be wide enough to cover the full physiological range of temperatures. Within this critical regime, we identified that polymer sequence influences surface tension predominately via shifting the critical temperature. Finally, we show that condensate surface tension over a wide range of temperatures can be calculated from the critical temperature and a single measurement of the interface width.

Phase separation of biomolecules into condensates has emerged as a mechanism for intracellular organization and affects many intracellular processes, including reaction pathways through the clustering of enzymes and pathway intermediates. Precise and rapid spatiotemporal control of reactions by condensates requires tuning of their sizes. However, the physical processes that govern the distribution of condensate sizes remain unclear. Here we show that both native and synthetic condensates display an exponential size distribution, which is captured by Monte Carlo simulations of fast nucleation followed by coalescence. In contrast, pathological aggregates exhibit a power-law size distribution. These distinct behaviours reflect the relative importance of nucleation and coalescence kinetics. We demonstrate this by utilizing a combination of synthetic and native condensates to probe the underlying physical mechanisms determining condensate size. The appearance of exponential distributions for abrupt nucleation versus power-law distributions under continuous nucleation may reflect a general principle that determines condensate size distributions.

While most studies of biomolecular phase separation have focused on the condensed phase, relatively little is known about the dilute phase. Theory suggests that stable complexes form in the dilute phase of two-component phase-separating systems, impacting phase separation; however, these complexes have not been interrogated experimentally. We show that such complexes indeed exist, using an in vitro reconstitution system of a phase-separated organelle, the algal pyrenoid, consisting of purified proteins Rubisco and EPYC1. Applying fluorescence correlation spectroscopy (FCS) to measure diffusion coefficients, we found that complexes form in the dilute phase with or without condensates present. The majority of these complexes contain exactly one Rubisco molecule. Additionally, we developed a simple analytical model which recapitulates experimental findings and provides molecular insights into the dilute phase organization. Thus, our results demonstrate the existence of protein complexes in the dilute phase, which could play important roles in the stability, dynamics, and regulation of condensates.

2022

How do growing bacterial colonies get their shapes? While colony morphogenesis is well studied in two dimensions, many bacteria grow as large colonies in three-dimensional (3D) environments, such as gels and tissues in the body or subsurface soils and sediments. Here, we describe the morphodynamics of large colonies of bacteria growing in three dimensions. Using experiments in transparent 3D granular hydrogel matrices, we show that dense colonies of four different species of bacteria generically become morphologically unstable and roughen as they consume nutrients and grow beyond a critical size-eventually adopting a characteristic branched, broccoli-like morphology independent of variations in the cell type and environmental conditions. This behavior reflects a key difference between two-dimensional (2D) and 3D colonies; while a 2D colony may access the nutrients needed for growth from the third dimension, a 3D colony inevitably becomes nutrient limited in its interior, driving a transition to unstable growth at its surface. We elucidate the onset of the instability using linear stability analysis and numerical simulations of a continuum model that treats the colony as an "active fluid" whose dynamics are driven by nutrient-dependent cellular growth. We find that when all dimensions of the colony substantially exceed the nutrient penetration length, nutrient-limited growth drives a 3D morphological instability that recapitulates essential features of the experimental observations. Our work thus provides a framework to predict and control the organization of growing colonies-as well as other forms of growing active matter, such as tumors and engineered living materials-in 3D environments.

Bacteria orchestrate collective behaviors and accomplish feats that would be unsuccessful if carried out by a lone bacterium. Processes undertaken by groups of bacteria include bioluminescence, biofilm formation, virulence factor production, and release of public goods that are shared by the community. Collective behaviors are controlled by signal transduction networks that integrate sensory information and transduce the information internally. Here, we discuss network features and mechanisms that, even in the face of dramatically changing environments, drive precise execution of bacterial group behaviors. We focus on representative quorum-sensing and second-messenger cyclic dimeric GMP (c-di-GMP) signal relays. We highlight ligand specificity versus sensitivity, how small-molecule ligands drive discrimination of kin versus nonkin, signal integration mechanisms, single-input sensory systems versus coincidence detectors, and tuning of input-output dynamics via feedback regulation. We summarize how different features of signal transduction systems allow groups of bacteria to successfully interpret and collectively react to dynamically changing environments.

Multivalent associative proteins with strong complementary interactions play a crucial role in phase separation of intracellular liquid condensates. We study the internal dynamics of such "bond-network" condensates comprising two complementary proteins via scaling analysis and molecular dynamics. We find that when stoichiometry is balanced, relaxation slows down dramatically due to a scarcity of alternative binding partners following bond breakage. This microscopic slow-down strongly affects the bulk diffusivity, viscosity, and mixing, which provides a means to experimentally test this prediction.

Cells can achieve internal organization by exploiting liquid-liquid phase separation to form biomolecular condensates. Here we focus on the surface properties of condensates composed of two multivalent associative polymers held together by one-to-one "sticker" bonds. Using coarse-grained molecular-dynamics simulations, we study the influence of component stoichiometry on condensate surface properties. We find that unequal stoichiometry results in enrichment of the majority species at the interface and a sharp reduction of surface tension. To relate these two effects, we show that the reduction in surface tension scales linearly with the excess concentration of free binding sites at the interface. Our results imply that each excess free site contributes an approximately fixed additional energy and entropy to the interface, with the latter dominating such that enrichment of free majority sites lowers the surface tension. Our work provides insight into novel physical mechanisms by which cells can regulate condensate surface properties.

Bacterial biofilms are multicellular communities that collectively overcome environmental threats and clinical treatments. To regulate the biofilm lifecycle, bacteria commonly transduce sensory information via the second messenger molecule cyclic diguanylate (c-di-GMP). Using experimental and modeling approaches, we quantitatively capture c-di-GMP signal transmission via the bifunctional polyamine receptor NspS-MbaA, from ligand binding to output, in the pathogen Vibrio cholerae. Upon binding of norspermidine or spermidine, NspS-MbaA synthesizes or degrades c-di-GMP, respectively, which, in turn, drives alterations specifically to biofilm gene expression. A long-standing question is how output specificity is achieved via c-di-GMP, a diffusible molecule that regulates dozens of effectors. We show that NspS-MbaA signals locally to specific effectors, sensitizing V. cholerae to polyamines. However, local signaling is not required for specificity, as changes to global cytoplasmic c-di-GMP levels can selectively regulate biofilm genes. This work establishes the input-output dynamics underlying c-di-GMP signaling, which could be useful for developing bacterial manipulation strategies.

Cross-feeding, the exchange of nutrients between organisms, is ubiquitous in microbial communities. Despite its importance in natural and engineered microbial systems, our understanding of how inter-species cross-feeding arises is incomplete, with existing theories limited to specific scenarios. Here, we introduce a novel theory for the emergence of such cross-feeding, which we term noise-averaging cooperation (NAC). NAC is based on the idea that, due to their small size, bacteria are prone to noisy regulation of metabolism which limits their growth rate. To compensate, related bacteria can share metabolites with each other to 'average out' noise and improve their collective growth. According to the Black Queen Hypothesis, this metabolite sharing among kin, a form of 'leakage', then allows for the evolution of metabolic interdependencies among species including de novo speciation via gene deletions. We first characterize NAC in a simple ecological model of cell metabolism, showing that metabolite leakage can in principle substantially increase growth rate in a community context. Next, we develop a generalized framework for estimating the potential benefits of NAC among real bacteria. Using single-cell protein abundance data, we predict that bacteria suffer from substantial noise-driven growth inefficiencies, and may therefore benefit from NAC. We then discuss potential evolutionary pathways for the emergence of NAC. Finally, we review existing evidence for NAC and outline potential experimental approaches to detect NAC in microbial communities.

Many eukaryotic photosynthetic organisms enhance their carbon uptake by supplying concentrated CO to the CO-fixing enzyme Rubisco in an organelle called the pyrenoid. Ongoing efforts seek to engineer this pyrenoid-based CO-concentrating mechanism (PCCM) into crops to increase yields. Here we develop a computational model for a PCCM on the basis of the postulated mechanism in the green alga Chlamydomonas reinhardtii. Our model recapitulates all Chlamydomonas PCCM-deficient mutant phenotypes and yields general biophysical principles underlying the PCCM. We show that an effective and energetically efficient PCCM requires a physical barrier to reduce pyrenoid CO leakage, as well as proper enzyme localization to reduce futile cycling between CO and HCO. Importantly, our model demonstrates the feasibility of a purely passive CO uptake strategy at air-level CO, while active HCO uptake proves advantageous at lower CO levels. We propose a four-step engineering path to increase the rate of CO fixation in the plant chloroplast up to threefold at a theoretical cost of only 1.3 ATP per CO fixed, thereby offering a framework to guide the engineering of a PCCM into land plants.

Bacteria require membrane fission for both cell division and endospore formation. In Bacillus subtilis, sporulation initiates with an asymmetric division that generates a large mother cell and a smaller forespore that contains only a quarter of its genome. As the mother cell membranes engulf the forespore, a DNA translocase pumps the rest of the chromosome into the small forespore compartment, inflating it due to increased turgor. When the engulfing membrane undergoes fission, the forespore is released into the mother cell cytoplasm. The B. subtilis protein FisB catalyzes membrane fission during sporulation, but the molecular basis is unclear. Here, we show that forespore inflation and FisB accumulation are both required for an efficient membrane fission. Forespore inflation leads to higher membrane tension in the engulfment membrane than in the mother cell membrane, causing the membrane to flow through the neck connecting the two membrane compartments. Thus, the mother cell supplies some of the membrane required for the growth of the membranes surrounding the forespore. The oligomerization of FisB at the membrane neck slows the equilibration of membrane tension by impeding the membrane flow. This leads to a further increase in the tension of the engulfment membrane, promoting its fission through lysis. Collectively, our data indicate that DNA translocation has a previously unappreciated second function in energizing the FisB-mediated membrane fission under energy-limited conditions.

How well mRNA transcript levels represent protein abundances has been a controversial issue. Particularly across different environments, correlations between mRNA and protein exhibit remarkable variability from gene to gene. Translational regulation is likely to be one of the key factors contributing to mismatches between mRNA level and protein abundance in bacteria. Here, we quantified genome-wide transcriptome and relative translation efficiency (RTE) under 12 different conditions in Escherichia coli. By quantifying the mRNA-RTE correlation both across genes and across conditions, we uncovered a diversity of gene-specific translational regulations, cooperating with transcriptional regulations, in response to carbon (C), nitrogen (N), and phosphate (P) limitations. Intriguingly, we found that many genes regulating translation are themselves subject to translational regulation, suggesting possible feedbacks. Furthermore, a random forest model suggests that codon usage partially predicts a gene's cross-condition variability in translation efficiency; such cross-condition variability tends to be an inherent quality of a gene, independent of the specific nutrient limitations. These findings broaden the understanding of translational regulation under different environments and provide novel strategies for the control of translation in synthetic biology. In addition, our data offers a resource for future multi-omics studies.

How do growing bacterial colonies get their shapes? While colony morphogenesis is well studied in two dimensions, many bacteria grow as large colonies in three-dimensional (3D) environments, such as gels and tissues in the body or subsurface soils and sediments. Here, we describe the morphodynamics of large colonies of bacteria growing in three dimensions. Using experiments in transparent 3D granular hydrogel matrices, we show that dense colonies of four different species of bacteria generically become morphologically unstable and roughen as they consume nutrients and grow beyond a critical size-eventually adopting a characteristic branched, broccoli-like morphology independent of variations in the cell type and environmental conditions. This behavior reflects a key difference between two-dimensional (2D) and 3D colonies; while a 2D colony may access the nutrients needed for growth from the third dimension, a 3D colony inevitably becomes nutrient limited in its interior, driving a transition to unstable growth at its surface. We elucidate the onset of the instability using linear stability analysis and numerical simulations of a continuum model that treats the colony as an "active fluid" whose dynamics are driven by nutrient-dependent cellular growth. We find that when all dimensions of the colony substantially exceed the nutrient penetration length, nutrient-limited growth drives a 3D morphological instability that recapitulates essential features of the experimental observations. Our work thus provides a framework to predict and control the organization of growing colonies-as well as other forms of growing active matter, such as tumors and engineered living materials-in 3D environments.