Cell signaling and gene transcription occur at faster time scales compared to cellular death, division, and evolution. Bridging these multiscale events in a model is computationally challenging. We introduce a framework for the systematic development of multiscale cell population models. Using message passing interface (MPI) parallelism, the framework creates a population model from a single-cell biochemical network model. It launches parallel simulations on a single-cell model and treats each stand-alone parallel process as a cell object. MPI mediates cell-to-cell and cell-to-environment communications in a server-client fashion. In the framework, model-specific higher level rules link the intracellular molecular events to cellular functions, such as death, division, or phenotype change. Cell death is implemented by terminating a parallel process, while cell division is carried out by creating a new process (daughter cell) from an existing one (mother cell). We first demonstrate these capabilities by creating two simple example models. In one model, we consider a relatively simple scenario where cells can evolve independently. In the other model, we consider interdependency among the cells, where cellular communication determines their collective behavior and evolution under a temporally evolving growth condition. We then demonstrate the framework’s capability by simulating a full-scale model of bacterial quorum sensing, where the dynamics of a population of bacterial cells is dictated by the intercellular communications in a time-evolving growth environment.
more »
« less
Influence of confinement on the spreading of bacterial populations
The spreading of bacterial populations is central to processes in agriculture, the environment, and medicine. However, existing models of spreading typically focus on cells in unconfined settings—despite the fact that many bacteria inhabit complex and crowded environments, such as soils, sediments, and biological tissues/gels, in which solid obstacles confine the cells and thereby strongly regulate population spreading. Here, we develop an extended version of the classic Keller-Segel model of bacterial spreading via motility that also incorporates cellular growth and division, and explicitly considers the influence of confinement in promoting both cell-solid and cell-cell collisions. Numerical simulations of this extended model demonstrate how confinement fundamentally alters the dynamics and morphology of spreading bacterial populations, in good agreement with recent experimental results. In particular, with increasing confinement, we find that cell-cell collisions increasingly hinder the initial formation and the long-time propagation speed of chemotactic pulses. Moreover, also with increasing confinement, we find that cellular growth and division plays an increasingly dominant role in driving population spreading—eventually leading to a transition from chemotactic spreading to growth-driven spreading via a slower, jammed front. This work thus provides a theoretical foundation for further investigations of the influence of confinement on bacterial spreading. More broadly, these results help to provide a framework to predict and control the dynamics of bacterial populations in complex and crowded environments.
more »
« less
- PAR ID:
- 10394412
- Editor(s):
- Wodarz, Dominik
- Date Published:
- Journal Name:
- PLOS Computational Biology
- Volume:
- 18
- Issue:
- 5
- ISSN:
- 1553-7358
- Page Range / eLocation ID:
- e1010063
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Csikász-Nagy, Attila (Ed.)Microbial populations show striking diversity in cell growth morphology and lifecycle; however, our understanding of how these factors influence the growth rate of cell populations remains limited. We use theory and simulations to predict the impact of asymmetric cell division, cell size regulation and single-cell stochasticity on the population growth rate. Our model predicts that coarse-grained noise in the single-cell growth rate λ decreases the population growth rate, as previously seen for symmetrically dividing cells. However, for a given noise in λ we find that dividing asymmetrically can enhance the population growth rate for cells with strong size control (between a “sizer” and an “adder”). To reconcile this finding with the abundance of symmetrically dividing organisms in nature, we propose that additional constraints on cell growth and division must be present which are not included in our model, and we explore the effects of selected extensions thereof. Further, we find that within our model, epigenetically inherited generation times may arise due to size control in asymmetrically dividing cells, providing a possible explanation for recent experimental observations in budding yeast. Taken together, our findings provide insight into the complex effects generated by non-canonical growth morphologies.more » « less
-
Efficient allocation of energy resources to key physiological functions allows living organisms to grow and thrive in diverse environments and adapt to a wide range of perturbations. To quantitatively understand how unicellular organisms utilize their energy resources in response to changes in growth environment, we introduce a theory of dynamic energy allocation that describes cellular growth dynamics by partitioning metabolizable energy into key physiological functions: growth, division, cell shape regulation, energy storage and loss through dissipation. By optimizing the energy flux for growth, we develop the equations governing the time evolution of cell morphology and growth rate in diverse environments. The resulting model accurately captures experimentally observed dependencies of bacterial cell size on growth rate, superlinear scaling of metabolic rate with cell size and predicts nutrient-dependent trade-offs between energy expended for growth, division and shape maintenance. By calibrating model parameters with experimental data for the model organismEscherichia coli, our model describes bacterial growth control in dynamic conditions, particularly during nutrient shifts and osmotic shocks. Integrating both the mechanical properties of the cell and underlying biochemical regulation, our model predicts the driving factors behind a wide range of observed morphological and growth phenomena with minimal added complexity.more » « less
-
I. ABSTRACT Bacteriophage (phage) infect, lyse, and propagate within bacterial populations. However, physiological changes in bacterial cell state can protect against infection even within genetically susceptible populations. One such example is the generation of endospores byBacillusand its relatives, characterized by a reversible state of reduced metabolic activity that protects cells against stressors including desiccation, energy limitation, antibiotics, and infection by phage. Here we tested how sporulation at the cellular scale impacts phage dynamics at population scales when propagating amongstB. subtilisin spatially structured environments. Initially, we found that plaques resulting from infection and lysis were approximately 3-fold smaller on lawns of sporulating wild-type bacteria vs. non-sporulating bacteria. Notably, plaque size was reduced due to an early termination of expanding phage plaques rather than the reduction of plaque growth speed. Microscopic imaging of the plaques revealed ‘sporulation rings’, i.e., spores enriched around plaque edges relative to phage-free regions. We developed a series of mathematical models of phage, bacteria, spore, and small molecules that recapitulate plaque dynamics and identify a putative mechanism: sporulation rings arise in response to lytic activity. In aggregate, sporulation rings inhibit phage from accessing susceptible cells even when sufficient resources are available for further infection and lysis. Together, our findings identify how dormancy can self-limit phage infections at population scales, opening new avenues to explore the entangled fates of phages and their bacterial hosts in environmental and therapeutic contexts.more » « less
-
Abstract The inside of a cell is highly crowded with proteins and other biomolecules. How proteins express their specific functions together with many off-target proteins in crowded cellular environments is largely unknown. Here, we investigate an inhibitor binding with c-Src kinase using atomistic molecular dynamics (MD) simulations in dilute as well as crowded protein solution. The populations of the inhibitor, 4-amino-5-(4-methylphenyl)−7-(t-butyl)pyrazolo[3,4-d]pyrimidine (PP1), in bulk solution and on the surface of c-Src kinase are reduced as the concentration of crowder bovine serum albumins (BSAs) increases. This observation is consistent with the reduced PP1 inhibitor efficacy in experimental c-Src kinase assays in addition with BSAs. The crowded environment changes the major binding pathway of PP1 toward c-Src kinase compared to that in dilute solution. This change is explained based on the population shift mechanism of local conformations near the inhibitor binding site in c-Src kinase.more » « less
An official website of the United States government

