Publications

2025

It was recently discovered that eukaryotic cells are host to a multiplicity of biomolecular condensates. These condensates typically contain protein and/or RNA components with intrinsically disordered regions (IDRs). While IDRs have been proposed and demonstrated to play many roles in condensate biology, we suggest here an additional crucial role of IDRs, which is to exclude unwanted “intruders” from condensates. This exclusion effect arises from the large conformational entropy of IDRs; i.e., there is a high free-energy cost to occupying space that would otherwise be available to the IDRs. By combining polymer theory with sticker-spacer simulations, we show that the relevant insertion free energy increases with the concentration of IDRs in the condensate as well as with intruder size, attaining a linear scaling with surface area for large intruders. We find that at realistic IDR concentrations, particles as small as the size of an average protein (4 nm in diameter) can be 97% excluded from condensates. To overcome this entropic barrier, molecules must interact favorably with condensate components to be recruited as clients. Application of the developed size-exclusion theory to biological condensates suggests that condensate IDRs may play a generic exclusionary role across organisms and types of condensates. Published by the American Physical Society 2025

In microbial communities, various cell types often coexist by occupying distinct spatial domains. What determines the shape of the interface between such domains—which, in turn, influences the interactions between cells and overall community function? Here, we address this question by developing a continuum model of a 2D spatially structured microbial community with two distinct cell types. We find that, depending on the balance of the different cell proliferation rates and substrate friction coefficients, the interface between domains is either stable and smooth or unstable and develops fingerlike protrusions. We establish quantitative principles describing when these different interfacial behaviors arise and find good agreement with both the results of previous experimental reports as well as new experiments performed here. Our work, thus, helps to provide a biophysical basis for understanding the interfacial morphodynamics of proliferating microbial communities as well as a broader range of proliferating active systems. Published by the American Physical Society 2025

Germinal centres are specialized microenvironments where B cells undergo affinity maturation. B cells expressing antibodies whose affinity is improved by somatic hypermutation are selected for expansion by limiting numbers of T follicular helper cells. Cell division is accompanied by mutation of the immunoglobulin genes, at what is believed to be a fixed rate of around 1 × 10 per base pair per cell division. As mutagenesis is random, the probability of acquiring deleterious mutations outweighs the probability of acquiring affinity-enhancing mutations. This effect might be heightened, and even become counterproductive, in B cells that express high-affinity antibodies and undergo the greatest number of cell divisions. Here we experimentally examine a theoretical model that explains how affinity maturation could be optimized by varying the rate of somatic hypermutation such that cells that express higher-affinity antibodies divide more but mutate less per division. Data obtained from mice immunized with SARS-CoV-2 vaccines or a model antigen align with the theoretical model and show that cells producing high-affinity antibodies shorten the G0/G1 phases of the cell cycle and reduce their mutation rates. We propose that these mechanisms safeguard high-affinity B cell lineages and enhance the outcomes of antibody affinity maturation.

Bacteria commonly exist in multicellular, surface-attached communities called biofilms. Biofilms are central to ecology, medicine, and industry. The Vibrio cholerae pathogen forms biofilms from single founder cells that, via cell division, mature into three-dimensional structures with distinct, yet reproducible, regional architectures. To define mechanisms underlying biofilm developmental transitions, we establish a single-molecule fluorescence in situ hybridization (smFISH) approach that enables accurate quantitation of spatiotemporal gene-expression patterns in biofilms at cell-scale resolution. smFISH analyses of V. cholerae biofilm regulatory and structural genes demonstrate that, as biofilms mature, overall matrix gene expression decreases, and simultaneously, a pattern emerges in which matrix gene expression becomes largely confined to peripheral biofilm cells. Both quorum sensing and c-di-GMP-signaling are required to generate the proper temporal pattern of matrix gene expression. Quorum sensing signaling is uniform across the biofilm, and thus, c-di-GMP-signaling alone sets the regional matrix gene expression pattern. The smFISH strategy provides insight into mechanisms conferring particular fates to individual biofilm cells.

Fluid flows are dominant features of many bacterial environments, and flow can often impact bacterial behaviors in unexpected ways. For example, the most common type of cardiovascular infection is heart valve colonization by gram-positive bacteria like and (endocarditis). This behavior is counterintuitive because heart valves experience high shear rates that would naively be expected to reduce colonization. To determine whether these bacteria preferentially colonize higher shear rate environments, we developed a microfluidic system to quantify the effect of flow conditions on the colonization of and . We find that the preferential colonization in high flow of both species is not specific to heart valves and can be found in simple configurations lacking any host factors. This behavior enables bacteria that are outcompeted in low flow to dominate in high flow. Surprisingly, experimental and computational studies reveal that the two species achieve this behavior via distinct mechanisms. grows in cell clusters and produces a dispersal signal whose transport is affected by shear rate. Meanwhile, grows in linear chains whose mechanical properties result in less dispersal in the presence of higher shear force. In addition to establishing two divergent mechanisms by which these bacteria each preferentially colonize high-flow environments, our findings highlight the importance of understanding bacterial behaviors at the level of collective interactions among cells. These results suggest that distinct multicellular nanocolony morphologies have previously unappreciated costs and benefits in different environments, like those introduced by fluid flow.

Colonies of the social bacterium Myxococcus xanthus go through a morphological transition from a thin colony of cells to three-dimensional droplet-like fruiting bodies as a strategy to survive starvation. The biological pathways that control the decision to form a fruiting body have been studied extensively. However, the mechanical events that trigger the creation of multiple cell layers and give rise to droplet formation remain poorly understood. By measuring cell orientation, velocity, polarity, and force with cell-scale resolution, we reveal a stochastic local polar order in addition to the more obvious nematic order. Average cell velocity and active force at topological defects agree with predictions from active nematic theory, but their fluctuations are substantially larger than the mean due to polar active forces generated by the self-propelled rod-shaped cells. We find that M. xanthus cells adjust their reversal frequency to tune the magnitude of this local polar order, which in turn controls the mechanical stresses and triggers layer formation in the colonies.

Many bacteria live in polymeric fluids, such as mucus, environmental polysaccharides, and extracellular polymers in biofilms. However, laboratory studies typically focus on cells in polymer-free fluids. Here, we show that interactions with polymers shape a fundamental feature of bacterial life-how they proliferate in space in multicellular colonies. Using experiments, we find that when polymer is sufficiently concentrated, cells generically and reversibly form large serpentine "cables" as they proliferate. By combining experiments with biophysical theory and simulations, we demonstrate that this distinctive form of colony morphogenesis arises from an interplay between polymer-induced entropic attraction between neighboring cells and their hindered ability to diffusely separate from each other in a viscous polymer solution. Our work thus reveals a pivotal role of polymers in sculpting proliferating bacterial colonies, with implications for how they interact with hosts and with the natural environment, and uncovers quantitative principles governing colony morphogenesis in such complex environments.

2024

Bacteria have evolved many defenses against invading viruses (phage). Despite the many bacterial defenses and phage counterdefenses, in most environments, bacteria and phage coexist, with neither driving the other to extinction. How is coexistence realized in the context of the bacteria/phage arms race, and how are immune repertoire sizes determined in conditions of coexistence? Here we develop a simple mathematical model to consider the evolutionary and ecological dynamics of competing bacteria and phage with different immune/counterimmune repertoires. We find an ecologically stable fixed point exhibiting coexistence, in agreement with the experimental observation that each individual bacterium typically carries multiple defense systems, though fewer than the maximum number possible. However, in simulations, the populations typically remain dynamic, exhibiting chaotic fluctuations around this fixed point. These dynamics enable coexistence even when phage (predator) strains outnumber bacteria (prey) strains. We obtain quantitative predictions for the mean, amplitude, and timescale of these dynamics. Our results provide a framework for understanding the evolutionary and ecological dynamics of the bacteria/phage arms race and demonstrate how bacteria/phage coexistence can stably arise from the coevolution of bacterial defense systems and phage counterdefense systems.

Season length and its associated variables can influence the expression of social behaviours, including the occurrence of eusociality in insects. Eusociality can vary widely across environmental gradients, both within and between different species. Numerous theoretical models have been developed to examine the life history traits that underlie the emergence and maintenance of eusociality, yet the impact of seasonality on this process is largely uncharacterized. Here, we present a theoretical model that incorporates season length and offspring development time into a single, individual-focused model to examine how these factors can shape the costs and benefits of social living. We find that longer season lengths and faster brood development times are sufficient to favour the emergence and maintenance of a social strategy, while shorter seasons favour a solitary one. We also identify a range of season lengths where social and solitary strategies can coexist. Moreover, our theoretical predictions are well matched to the natural history and behaviour of two flexibly eusocial bee species, suggesting that our model can make realistic predictions about the evolution of different social strategies. Broadly, this work reveals the crucial role that environmental conditions can have in shaping social behaviour and its evolution and it underscores the need for further models that explicitly incorporate such variation to study the evolutionary trajectories of eusociality.

A hallmark of biomolecular condensates formed via liquid-liquid phase separation is that they dynamically exchange material with their surroundings, and this process can be crucial to condensate function. Intuitively, the rate of exchange can be limited by the flux from the dilute phase or by the mixing speed in the dense phase. Surprisingly, a recent experiment suggests that exchange can also be limited by the dynamics at the droplet interface, implying the existence of an 'interface resistance'. Here, we first derive an analytical expression for the timescale of condensate material exchange, which clearly conveys the physical factors controlling exchange dynamics. We then utilize sticker-spacer polymer models to show that interface resistance can arise when incident molecules transiently touch the interface without entering the dense phase, i.e., the molecules 'bounce' from the interface. Our work provides insight into condensate exchange dynamics, with implications for both natural and synthetic systems.

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.