1
|
Helenek C, Krzysztoń R, Petreczky J, Wan Y, Cabral M, Coraci D, Balázsi G. Synthetic gene circuit evolution: Insights and opportunities at the mid-scale. Cell Chem Biol 2024; 31:1447-1459. [PMID: 38925113 PMCID: PMC11330362 DOI: 10.1016/j.chembiol.2024.05.018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/12/2024] [Revised: 05/07/2024] [Accepted: 05/30/2024] [Indexed: 06/28/2024]
Abstract
Directed evolution focuses on optimizing single genetic components for predefined engineering goals by artificial mutagenesis and selection. In contrast, experimental evolution studies the adaptation of entire genomes in serially propagated cell populations, to provide an experimental basis for evolutionary theory. There is a relatively unexplored gap at the middle ground between these two techniques, to evolve in vivo entire synthetic gene circuits with nontrivial dynamic function instead of single parts or whole genomes. We discuss the requirements for such mid-scale evolution, with hypothetical examples for evolving synthetic gene circuits by appropriate selection and targeted shuffling of a seed set of genetic components in vivo. Implementing similar methods should aid the rapid generation, functionalization, and optimization of synthetic gene circuits in various organisms and environments, accelerating both the development of biomedical and technological applications and the understanding of principles guiding regulatory network evolution.
Collapse
Affiliation(s)
- Christopher Helenek
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA; Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY 11794, USA
| | - Rafał Krzysztoń
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA
| | - Julia Petreczky
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA; Department of Chemistry, Stony Brook University, Stony Brook, NY 11794, USA
| | - Yiming Wan
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA
| | - Mariana Cabral
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA; Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY 11794, USA
| | - Damiano Coraci
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA
| | - Gábor Balázsi
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY 11794, USA; Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY 11794, USA; Stony Brook Cancer Center, Stony Brook University, Stony Brook, NY 11794, USA.
| |
Collapse
|
2
|
Kohram M, Sanderson AE, Loui A, Thompson PV, Vashistha H, Shomar A, Oltvai ZN, Salman H. Nonlethal deleterious mutation-induced stress accelerates bacterial aging. Proc Natl Acad Sci U S A 2024; 121:e2316271121. [PMID: 38709929 PMCID: PMC11098108 DOI: 10.1073/pnas.2316271121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2023] [Accepted: 03/29/2024] [Indexed: 05/08/2024] Open
Abstract
Random mutagenesis, including when it leads to loss of gene function, is a key mechanism enabling microorganisms' long-term adaptation to new environments. However, loss-of-function mutations are often deleterious, triggering, in turn, cellular stress and complex homeostatic stress responses, called "allostasis," to promote cell survival. Here, we characterize the differential impacts of 65 nonlethal, deleterious single-gene deletions on Escherichia coli growth in three different growth environments. Further assessments of select mutants, namely, those bearing single adenosine triphosphate (ATP) synthase subunit deletions, reveal that mutants display reorganized transcriptome profiles that reflect both the environment and the specific gene deletion. We also find that ATP synthase α-subunit deleted (ΔatpA) cells exhibit elevated metabolic rates while having slower growth compared to wild-type (wt) E. coli cells. At the single-cell level, compared to wt cells, individual ΔatpA cells display near normal proliferation profiles but enter a postreplicative state earlier and exhibit a distinct senescence phenotype. These results highlight the complex interplay between genomic diversity, adaptation, and stress response and uncover an "aging cost" to individual bacterial cells for maintaining population-level resilience to environmental and genetic stress; they also suggest potential bacteriostatic antibiotic targets and -as select human genetic diseases display highly similar phenotypes, - a bacterial origin of some human diseases.
Collapse
Affiliation(s)
- Maryam Kohram
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA15260
| | - Amy E. Sanderson
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA15260
| | - Alicia Loui
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA15260
| | | | - Harsh Vashistha
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA15260
| | - Aseel Shomar
- Department of Chemical Engineering, Technion–Israel Institute of Technology, Haifa32000, Israel
| | - Zoltán N. Oltvai
- Department of Pathology, University of Pittsburgh, Pittsburgh, PA15260
- Department of Computational and Systems Biology, University of Pittsburgh, Pittsburgh, PA15260
- Department of Pathology and Laboratory Medicine, University of Rochester, Rochester, NY14627
| | - Hanna Salman
- Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA15260
| |
Collapse
|
3
|
Stevanovic M, Teuber Carvalho JP, Bittihn P, Schultz D. Dynamical model of antibiotic responses linking expression of resistance genes to metabolism explains emergence of heterogeneity during drug exposures. Phys Biol 2024; 21:036002. [PMID: 38412523 PMCID: PMC10988634 DOI: 10.1088/1478-3975/ad2d64] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Revised: 01/25/2024] [Accepted: 02/27/2024] [Indexed: 02/29/2024]
Abstract
Antibiotic responses in bacteria are highly dynamic and heterogeneous, with sudden exposure of bacterial colonies to high drug doses resulting in the coexistence of recovered and arrested cells. The dynamics of the response is determined by regulatory circuits controlling the expression of resistance genes, which are in turn modulated by the drug's action on cell growth and metabolism. Despite advances in understanding gene regulation at the molecular level, we still lack a framework to describe how feedback mechanisms resulting from the interdependence between expression of resistance and cell metabolism can amplify naturally occurring noise and create heterogeneity at the population level. To understand how this interplay affects cell survival upon exposure, we constructed a mathematical model of the dynamics of antibiotic responses that links metabolism and regulation of gene expression, based on the tetracycline resistancetetoperon inE. coli. We use this model to interpret measurements of growth and expression of resistance in microfluidic experiments, both in single cells and in biofilms. We also implemented a stochastic model of the drug response, to show that exposure to high drug levels results in large variations of recovery times and heterogeneity at the population level. We show that stochasticity is important to determine how nutrient quality affects cell survival during exposure to high drug concentrations. A quantitative description of how microbes respond to antibiotics in dynamical environments is crucial to understand population-level behaviors such as biofilms and pathogenesis.
Collapse
Affiliation(s)
- Mirjana Stevanovic
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States of America
| | - João Pedro Teuber Carvalho
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States of America
| | - Philip Bittihn
- Max Planck Institute for Dynamics and Self-Organization, Göttingen, Germany
- Institute for the Dynamics of Complex Systems, University of Göttingen, Göttingen, Germany
| | - Daniel Schultz
- Department of Microbiology and Immunology, Geisel School of Medicine at Dartmouth, Hanover, NH, United States of America
| |
Collapse
|
4
|
Glass DS, Bren A, Vaisbourd E, Mayo A, Alon U. A synthetic differentiation circuit in Escherichia coli for suppressing mutant takeover. Cell 2024; 187:931-944.e12. [PMID: 38320549 PMCID: PMC10882425 DOI: 10.1016/j.cell.2024.01.024] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Revised: 11/27/2023] [Accepted: 01/16/2024] [Indexed: 02/08/2024]
Abstract
Differentiation is crucial for multicellularity. However, it is inherently susceptible to mutant cells that fail to differentiate. These mutants outcompete normal cells by excessive self-renewal. It remains unclear what mechanisms can resist such mutant expansion. Here, we demonstrate a solution by engineering a synthetic differentiation circuit in Escherichia coli that selects against these mutants via a biphasic fitness strategy. The circuit provides tunable production of synthetic analogs of stem, progenitor, and differentiated cells. It resists mutations by coupling differentiation to the production of an essential enzyme, thereby disadvantaging non-differentiating mutants. The circuit selected for and maintained a positive differentiation rate in long-term evolution. Surprisingly, this rate remained constant across vast changes in growth conditions. We found that transit-amplifying cells (fast-growing progenitors) underlie this environmental robustness. Our results provide insight into the stability of differentiation and demonstrate a powerful method for engineering evolutionarily stable multicellular consortia.
Collapse
Affiliation(s)
- David S Glass
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel.
| | - Anat Bren
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Elizabeth Vaisbourd
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Avi Mayo
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel
| | - Uri Alon
- Department of Molecular Cell Biology, Weizmann Institute of Science, Rehovot 76100, Israel.
| |
Collapse
|
5
|
Ruess J, Ballif G, Aditya C. Stochastic chemical kinetics of cell fate decision systems: From single cells to populations and back. J Chem Phys 2023; 159:184103. [PMID: 37937934 DOI: 10.1063/5.0160529] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2023] [Accepted: 10/14/2023] [Indexed: 11/09/2023] Open
Abstract
Stochastic chemical kinetics is a widely used formalism for studying stochasticity of chemical reactions inside single cells. Experimental studies of reaction networks are generally performed with cells that are part of a growing population, yet the population context is rarely taken into account when models are developed. Models that neglect the population context lose their validity whenever the studied system influences traits of cells that can be selected in the population, a property that naturally arises in the complex interplay between single-cell and population dynamics of cell fate decision systems. Here, we represent such systems as absorbing continuous-time Markov chains. We show that conditioning on non-absorption allows one to derive a modified master equation that tracks the time evolution of the expected population composition within a growing population. This allows us to derive consistent population dynamics models from a specification of the single-cell process. We use this approach to classify cell fate decision systems into two types that lead to different characteristic phases in emerging population dynamics. Subsequently, we deploy the gained insights to experimentally study a recurrent problem in biology: how to link plasmid copy number fluctuations and plasmid loss events inside single cells to growth of cell populations in dynamically changing environments.
Collapse
Affiliation(s)
- Jakob Ruess
- Inria Saclay, 91120 Palaiseau, France
- Institut Pasteur, Université Paris Cité, 75015 Paris, France
| | | | - Chetan Aditya
- Institut Pasteur, Université Paris Cité, 75015 Paris, France
- Department of Molecular Biology, Princeton University, Princeton, New Jersey 08544, USA
| |
Collapse
|
6
|
Torres A, Cockerell S, Phillips M, Balázsi G, Ghosh K. MaxCal can infer models from coupled stochastic trajectories of gene expression and cell division. Biophys J 2023; 122:2623-2635. [PMID: 37218129 PMCID: PMC10397576 DOI: 10.1016/j.bpj.2023.05.017] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2023] [Revised: 05/03/2023] [Accepted: 05/18/2023] [Indexed: 05/24/2023] Open
Abstract
Gene expression is inherently noisy due to small numbers of proteins and nucleic acids inside a cell. Likewise, cell division is stochastic, particularly when tracking at the level of a single cell. The two can be coupled when gene expression affects the rate of cell division. Single-cell time-lapse experiments can measure both fluctuations by simultaneously recording protein levels inside a cell and its stochastic division. These information-rich noisy trajectory data sets can be harnessed to learn about the underlying molecular and cellular details that are often not known a priori. A critical question is: How can we infer a model given data where fluctuations at two levels-gene expression and cell division-are intricately convoluted? We show the principle of maximum caliber (MaxCal)-integrated within a Bayesian framework-can be used to infer several cellular and molecular details (division rates, protein production, and degradation rates) from these coupled stochastic trajectories (CSTs). We demonstrate this proof of concept using synthetic data generated from a known model. An additional challenge in data analysis is that trajectories are often not in protein numbers, but in noisy fluorescence that depends on protein number in a probabilistic manner. We again show that MaxCal can infer important molecular and cellular rates even when data are in fluorescence, another example of CST with three confounding factors-gene expression noise, cell division noise, and fluorescence distortion-all coupled. Our approach will provide guidance to build models in synthetic biology experiments as well as general biological systems where examples of CSTs are abundant.
Collapse
Affiliation(s)
- Andrew Torres
- Department of Physics and Astronomy, University of Denver, Denver, Colorado
| | - Spencer Cockerell
- Department of Physics and Astronomy, University of Denver, Denver, Colorado
| | - Michael Phillips
- Department of Physics and Astronomy, University of Denver, Denver, Colorado
| | - Gábor Balázsi
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York
| | - Kingshuk Ghosh
- Molecular and Cellular Biophysics, University of Denver, Denver, Colorado; Department of Physics and Astronomy, University of Denver, Denver, Colorado.
| |
Collapse
|
7
|
Wan Y, Cohen J, Szenk M, Farquhar KS, Coraci D, Krzysztoń R, Azukas J, Van Nest N, Smashnov A, Chern YJ, De Martino D, Nguyen LC, Bien H, Bravo-Cordero JJ, Chan CH, Rosner MR, Balázsi G. Nonmonotone invasion landscape by noise-aware control of metastasis activator levels. Nat Chem Biol 2023; 19:887-899. [PMID: 37231268 PMCID: PMC10299915 DOI: 10.1038/s41589-023-01344-z] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Accepted: 04/18/2023] [Indexed: 05/27/2023]
Abstract
A major pharmacological assumption is that lowering disease-promoting protein levels is generally beneficial. For example, inhibiting metastasis activator BACH1 is proposed to decrease cancer metastases. Testing such assumptions requires approaches to measure disease phenotypes while precisely adjusting disease-promoting protein levels. Here we developed a two-step strategy to integrate protein-level tuning, noise-aware synthetic gene circuits into a well-defined human genomic safe harbor locus. Unexpectedly, engineered MDA-MB-231 metastatic human breast cancer cells become more, then less and then more invasive as we tune BACH1 levels up, irrespective of the native BACH1. BACH1 expression shifts in invading cells, and expression of BACH1's transcriptional targets confirm BACH1's nonmonotone phenotypic and regulatory effects. Thus, chemical inhibition of BACH1 could have unwanted effects on invasion. Additionally, BACH1's expression variability aids invasion at high BACH1 expression. Overall, precisely engineered, noise-aware protein-level control is necessary and important to unravel disease effects of genes to improve clinical drug efficacy.
Collapse
Affiliation(s)
- Yiming Wan
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Joseph Cohen
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Mariola Szenk
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Kevin S Farquhar
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
- Genetics and Epigenetics Graduate Program, The University of Texas MD Anderson Cancer Center, UT Health Graduate School of Biomedical Sciences, Houston, TX, USA
| | - Damiano Coraci
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Rafał Krzysztoń
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Joshua Azukas
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
| | - Nicholas Van Nest
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Alex Smashnov
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Yi-Jye Chern
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY, USA
- Department of Molecular and Cellular Biology, Roswell Park Cancer Institute, Buffalo, NY, USA
| | - Daniela De Martino
- Department of Medicine, Division of Hematology and Oncology, Tisch Cancer Institute, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Long Chi Nguyen
- Ben May Department for Cancer Research, University of Chicago, Chicago, IL, USA
| | - Harold Bien
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA
| | - Jose Javier Bravo-Cordero
- Department of Medicine, Division of Hematology and Oncology, Tisch Cancer Institute, Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Chia-Hsin Chan
- Department of Pharmacological Sciences, Stony Brook University, Stony Brook, NY, USA
- Department of Molecular and Cellular Biology, Roswell Park Cancer Institute, Buffalo, NY, USA
| | - Marsha Rich Rosner
- Ben May Department for Cancer Research, University of Chicago, Chicago, IL, USA
| | - Gábor Balázsi
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, USA.
- Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA.
- Stony Brook Cancer Center, Stony Brook University, Stony Brook, NY, USA.
| |
Collapse
|
8
|
Kong LW, Shi W, Tian XJ, Lai YC. Effects of growth feedback on gene circuits: A dynamical understanding. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.06.543915. [PMID: 37333159 PMCID: PMC10274713 DOI: 10.1101/2023.06.06.543915] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/20/2023]
Abstract
The successful integration of engineered gene circuits into host cells remains a significant challenge in synthetic biology due to circuit-host interactions, such as growth feedback, where the circuit influences cell growth and vice versa. Understanding the dynamics of circuit failures and identifying topologies resilient to growth feedback are crucial for both fundamental and applied research. Utilizing transcriptional regulation circuits with adaptation as a paradigm, we systematically study 435 distinct topological structures and uncover six categories of failures. Three dynamical mechanisms of circuit failures are identified: continuous deformation of the response curve, strengthened or induced oscillations, and sudden switching to coexisting attractors. Our extensive computations also uncover a scaling law between a circuit robustness measure and the strength of growth feedback. Despite the negative effects of growth feedback on the majority of circuit topologies, we identify a few circuits that maintain optimal performance as designed, a feature important for applications.
Collapse
|
9
|
Barua A, Hatzikirou H. Cell Decision Making through the Lens of Bayesian Learning. ENTROPY (BASEL, SWITZERLAND) 2023; 25:e25040609. [PMID: 37190396 PMCID: PMC10137733 DOI: 10.3390/e25040609] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 03/20/2023] [Accepted: 03/29/2023] [Indexed: 05/17/2023]
Abstract
Cell decision making refers to the process by which cells gather information from their local microenvironment and regulate their internal states to create appropriate responses. Microenvironmental cell sensing plays a key role in this process. Our hypothesis is that cell decision-making regulation is dictated by Bayesian learning. In this article, we explore the implications of this hypothesis for internal state temporal evolution. By using a timescale separation between internal and external variables on the mesoscopic scale, we derive a hierarchical Fokker-Planck equation for cell-microenvironment dynamics. By combining this with the Bayesian learning hypothesis, we find that changes in microenvironmental entropy dominate the cell state probability distribution. Finally, we use these ideas to understand how cell sensing impacts cell decision making. Notably, our formalism allows us to understand cell state dynamics even without exact biochemical information about cell sensing processes by considering a few key parameters.
Collapse
Affiliation(s)
- Arnab Barua
- Departement de Biochimie, Université de Montréal, Montréal, QC H3T 1C5, Canada
- Centre Robert-Cedergren en Bio-Informatique et Génomique, Université de Montréal, Montréal, QC H3C 3J7, Canada
| | - Haralampos Hatzikirou
- Center for Information Services and High Performance Computing, Technische Univesität Dresden, 01062 Dresden, Germany
- Mathematics Department, Khalifa University, Abu Dhabi P.O. Box 127788, United Arab Emirates
| |
Collapse
|
10
|
Pillai M, Hojel E, Jolly MK, Goyal Y. Unraveling non-genetic heterogeneity in cancer with dynamical models and computational tools. NATURE COMPUTATIONAL SCIENCE 2023; 3:301-313. [PMID: 38177938 DOI: 10.1038/s43588-023-00427-0] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Accepted: 03/03/2023] [Indexed: 01/06/2024]
Abstract
Individual cells within an otherwise genetically homogenous population constantly undergo fluctuations in their molecular state, giving rise to non-genetic heterogeneity. Such diversity is being increasingly implicated in cancer therapy resistance and metastasis. Identifying the origins of non-genetic heterogeneity is therefore crucial for making clinical breakthroughs. We discuss with examples how dynamical models and computational tools have provided critical multiscale insights into the nature and consequences of non-genetic heterogeneity in cancer. We demonstrate how mechanistic modeling has been pivotal in establishing key concepts underlying non-genetic diversity at various biological scales, from population dynamics to gene regulatory networks. We discuss advances in single-cell longitudinal profiling techniques to reveal patterns of non-genetic heterogeneity, highlighting the ongoing efforts and challenges in statistical frameworks to robustly interpret such multimodal datasets. Moving forward, we stress the need for data-driven statistical and mechanistically motivated dynamical frameworks to come together to develop predictive cancer models and inform therapeutic strategies.
Collapse
Affiliation(s)
- Maalavika Pillai
- Department of Cell and Developmental Biology, Feinberg School of Medicine, Northwestern University, Chicago, IL, USA
- Center for Synthetic Biology, Northwestern University, Chicago, IL, USA
- Robert H. Lurie Comprehensive Cancer Center, Northwestern University Feinberg School of Medicine, Chicago, IL, USA
- Centre for BioSystems Science and Engineering, Indian Institute of Science, Bangalore, India
| | - Emilia Hojel
- Department of Cell and Developmental Biology, Feinberg School of Medicine, Northwestern University, Chicago, IL, USA
- Center for Synthetic Biology, Northwestern University, Chicago, IL, USA
- Department of Biomedical Engineering, Northwestern University McCormick School of Engineering, Evanston, IL, USA
| | - Mohit Kumar Jolly
- Centre for BioSystems Science and Engineering, Indian Institute of Science, Bangalore, India.
| | - Yogesh Goyal
- Department of Cell and Developmental Biology, Feinberg School of Medicine, Northwestern University, Chicago, IL, USA.
- Center for Synthetic Biology, Northwestern University, Chicago, IL, USA.
- Robert H. Lurie Comprehensive Cancer Center, Northwestern University Feinberg School of Medicine, Chicago, IL, USA.
- Department of Biomedical Engineering, Northwestern University McCormick School of Engineering, Evanston, IL, USA.
| |
Collapse
|
11
|
Wang T, Weiss A, Aqeel A, Wu F, Lopatkin AJ, David LA, You L. Horizontal gene transfer enables programmable gene stability in synthetic microbiota. Nat Chem Biol 2022; 18:1245-1252. [PMID: 36050493 PMCID: PMC10018779 DOI: 10.1038/s41589-022-01114-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2021] [Accepted: 07/15/2022] [Indexed: 11/09/2022]
Abstract
The functions of many microbial communities exhibit remarkable stability despite fluctuations in the compositions of these communities. To date, a mechanistic understanding of this function-composition decoupling is lacking. Statistical mechanisms have been commonly hypothesized to explain such decoupling. Here, we proposed that dynamic mechanisms, mediated by horizontal gene transfer (HGT), also enable the independence of functions from the compositions of microbial communities. We combined theoretical analysis with numerical simulations to illustrate that HGT rates can determine the stability of gene abundance in microbial communities. We further validated these predictions using engineered microbial consortia of different complexities transferring one or more than a dozen clinically isolated plasmids, as well as through the reanalysis of data from the literature. Our results demonstrate a generalizable strategy to program the gene stability of microbial communities.
Collapse
Affiliation(s)
- Teng Wang
- Department of Biomedical Engineering, Duke University, Durham, NC, USA
| | - Andrea Weiss
- Department of Biomedical Engineering, Duke University, Durham, NC, USA
| | - Ammara Aqeel
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA
| | - Feilun Wu
- Department of Biomedical Engineering, Duke University, Durham, NC, USA
| | - Allison J Lopatkin
- Department of Chemical Engineering, University of Rochester, Rochester, NY, USA
| | - Lawrence A David
- Department of Biomedical Engineering, Duke University, Durham, NC, USA
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA
- Center for Genomic and Computational Biology, Duke University, Durham, NC, USA
| | - Lingchong You
- Department of Biomedical Engineering, Duke University, Durham, NC, USA.
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA.
- Center for Genomic and Computational Biology, Duke University, Durham, NC, USA.
| |
Collapse
|
12
|
Guthrie J, Charlebois D. Non-genetic resistance facilitates survival while hindering the evolution of drug resistance due to intraspecific competition. Phys Biol 2022; 19. [PMID: 35998624 DOI: 10.1088/1478-3975/ac8c17] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2022] [Accepted: 08/23/2022] [Indexed: 11/11/2022]
Abstract
Rising rates of resistance to antimicrobial drugs threaten the effective treatment of infections across the globe. Drug resistance has been established to emerge from non-genetic mechanisms as well as from genetic mechanisms. However, it is still unclear how non-genetic resistance affects the evolution of genetic drug resistance. We develop deterministic and stochastic population models that incorporate resource competition to quantitatively investigate the transition from non-genetic to genetic resistance during the exposure to static and cidal drugs. We find that non-genetic resistance facilitates the survival of cell populations during drug treatment while hindering the development of genetic resistance due to competition between the non-genetically and genetically resistant subpopulations. Non-genetic resistance in the presence of subpopulation competition increases the fixation times of drug resistance mutations, while increasing the probability of mutation before population extinction during cidal drug treatment. Intense intraspecific competition during drug treatment leads to extinction of susceptible and non-genetically resistant subpopulations. Alternating between drug and no drug conditions results in oscillatory population dynamics, increased resistance mutation fixation timescales, and reduced population survival. These findings advance our fundamental understanding of the evolution of resistance and may guide novel treatment strategies for patients with drug-resistant infections.
Collapse
Affiliation(s)
- Joshua Guthrie
- Department of Physics, University of Alberta, 11455 Saskatchewan Drive NW, Edmonton, Alberta, T6G 2E1, CANADA
| | - Daniel Charlebois
- Departments of Physics and Biological Sciences, University of Alberta, 11455 Saskatchewan Drive NW, Edmonton, Alberta, T6G 2E1, CANADA
| |
Collapse
|
13
|
Lunz D, Bonnans JF, Ruess J. Revisiting moment-closure methods with heterogeneous multiscale population models. Math Biosci 2022; 350:108866. [PMID: 35753520 DOI: 10.1016/j.mbs.2022.108866] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Revised: 04/10/2022] [Accepted: 06/08/2022] [Indexed: 11/29/2022]
Abstract
Stochastic chemical kinetics at the single-cell level give rise to heterogeneous populations of cells even when all individuals are genetically identical. This heterogeneity can lead to nonuniform behaviour within populations, including different growth characteristics, cell-fate dynamics, and response to stimuli. Ultimately, these diverse behaviours lead to intricate population dynamics that are inherently multiscale: the population composition evolves based on population-level processes that interact with stochastically distributed single-cell states. Therefore, descriptions that account for this heterogeneity are essential to accurately model and control chemical processes. However, for real-world systems such models are computationally expensive to simulate, which can make optimisation problems, such as optimal control or parameter inference, prohibitively challenging. Here, we consider a class of multiscale population models that incorporate population-level mechanisms while remaining faithful to the underlying stochasticity at the single-cell level and the interplay between these two scales. To address the complexity, we study an order-reduction approximations based on the distribution moments. Since previous moment-closure work has focused on the single-cell kinetics, extending these techniques to populations models prompts us to revisit old observations as well as tackle new challenges. In this extended multiscale context, we encounter the previously established observation that the simplest closure techniques can lead to non-physical system trajectories. Despite their poor performance in some systems, we provide an example where these simple closures outperform more sophisticated closure methods in accurately, efficiently, and robustly solving the problem of optimal control of bioproduction in a microbial consortium model.
Collapse
Affiliation(s)
- Davin Lunz
- Inria Paris, 2 rue Simone Iff, 75012 Paris, France; Institut Pasteur, 28 rue du Docteur Roux, 75015 Paris, France.
| | - J Frédéric Bonnans
- Université Paris-Saclay, CNRS, CentraleSupélec, Inria, Laboratory of signals and systems, 91190, Gif-sur-Yvette, France
| | - Jakob Ruess
- Inria Paris, 2 rue Simone Iff, 75012 Paris, France; Institut Pasteur, 28 rue du Docteur Roux, 75015 Paris, France
| |
Collapse
|
14
|
Aditya C, Bertaux F, Batt G, Ruess J. Using single-cell models to predict the functionality of synthetic circuits at the population scale. Proc Natl Acad Sci U S A 2022; 119:e2114438119. [PMID: 35271387 PMCID: PMC8931247 DOI: 10.1073/pnas.2114438119] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2021] [Accepted: 02/07/2022] [Indexed: 12/16/2022] Open
Abstract
SignificanceAt the single-cell level, biochemical processes are inherently stochastic. For many natural systems, the resulting cell-to-cell variability is exploited by microbial populations. In synthetic biology, however, the interplay of cell-to-cell variability and population processes such as selection or growth often leads to circuits not functioning as predicted by simple models. Here we show how multiscale stochastic kinetic models that simultaneously track single-cell and population processes can be obtained based on an augmentation of the chemical master equation. These models enable us to quantitatively predict complex population dynamics of a yeast optogenetic differentiation system from a specification of the circuit's components and to demonstrate how cell-to-cell variability can be exploited to purposefully create unintuitive circuit functionality.
Collapse
Affiliation(s)
- Chetan Aditya
- Inria Paris, Inria, 75012 Paris, France
- Department of Computational Biology, Institut Pasteur, 75015 Paris, France
- Université de Paris, 75006 Paris, France
| | - François Bertaux
- Inria Paris, Inria, 75012 Paris, France
- Department of Computational Biology, Institut Pasteur, 75015 Paris, France
| | - Gregory Batt
- Inria Paris, Inria, 75012 Paris, France
- Department of Computational Biology, Institut Pasteur, 75015 Paris, France
| | - Jakob Ruess
- Inria Paris, Inria, 75012 Paris, France
- Department of Computational Biology, Institut Pasteur, 75015 Paris, France
| |
Collapse
|
15
|
Bose I. Tipping the Balance: A Criticality Perspective. ENTROPY 2022; 24:e24030405. [PMID: 35327916 PMCID: PMC8947304 DOI: 10.3390/e24030405] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 03/10/2022] [Accepted: 03/12/2022] [Indexed: 01/02/2023]
Abstract
Cell populations are often characterised by phenotypic heterogeneity in the form of two distinct subpopulations. We consider a model of tumour cells consisting of two subpopulations: non-cancer promoting (NCP) and cancer-promoting (CP). Under steady state conditions, the model has similarities with a well-known model of population genetics which exhibits a purely noise-induced transition from unimodality to bimodality at a critical value of the noise intensity σ2. The noise is associated with the parameter λ representing the system-environment coupling. In the case of the tumour model, λ has a natural interpretation in terms of the tissue microenvironment which has considerable influence on the phenotypic composition of the tumour. Oncogenic transformations give rise to considerable fluctuations in the parameter. We compute the λ−σ2 phase diagram in a stochastic setting, drawing analogies between bifurcations and phase transitions. In the region of bimodality, a transition from a state of balance to a state of dominance, in terms of the competing subpopulations, occurs at λ = 0. Away from this point, the NCP (CP) subpopulation becomes dominant as λ changes towards positive (negative) values. The variance of the steady state probability density function as well as two entropic measures provide characteristic signatures at the transition point.
Collapse
Affiliation(s)
- Indrani Bose
- Department of Physics, Bose Institute, 93/1, A. P. C. Road, Kolkata 700009, India
| |
Collapse
|
16
|
Traets JJ, van der Burght SN, Rademakers S, Jansen G, van Zon JS. Mechanism of life-long maintenance of neuron identity despite molecular fluctuations. eLife 2021; 10:66955. [PMID: 34908528 PMCID: PMC8735970 DOI: 10.7554/elife.66955] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2021] [Accepted: 12/14/2021] [Indexed: 11/13/2022] Open
Abstract
Cell fate is maintained over long timescales, yet molecular fluctuations can lead to spontaneous loss of this differentiated state. Our simulations identified a possible mechanism that explains life-long maintenance of ASE neuron fate in Caenorhabditis elegans by the terminal selector transcription factor CHE-1. Here, fluctuations in CHE-1 level are buffered by the reservoir of CHE-1 bound at its target promoters, which ensures continued che-1 expression by preferentially binding the che-1 promoter. We provide experimental evidence for this mechanism by showing that che-1 expression was resilient to induced transient CHE-1 depletion, while both expression of CHE-1 targets and ASE function were lost. We identified a 130 bp che-1 promoter fragment responsible for this resilience, with deletion of a homeodomain binding site in this fragment causing stochastic loss of ASE identity long after its determination. Because network architectures that support this mechanism are highly conserved in cell differentiation, it may explain stable cell fate maintenance in many systems.
Collapse
Affiliation(s)
| | | | | | - Gert Jansen
- Department of Cell Biology, Erasmus MC, Rotterdam, Netherlands
| | - Jeroen S van Zon
- Quantitative Developmental Biology, AMOLF, Amsterdam, Netherlands
| |
Collapse
|
17
|
Potentiating antibiotic efficacy via perturbation of non-essential gene expression. Commun Biol 2021; 4:1267. [PMID: 34741116 PMCID: PMC8571399 DOI: 10.1038/s42003-021-02783-x] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2020] [Accepted: 09/21/2021] [Indexed: 02/06/2023] Open
Abstract
Proliferation of multidrug-resistant (MDR) bacteria poses a threat to human health, requiring new strategies. Here we propose using fitness neutral gene expression perturbations to potentiate antibiotics. We systematically explored 270 gene knockout-antibiotic combinations in Escherichia coli, identifying 90 synergistic interactions. Identified gene targets were subsequently tested for antibiotic synergy on the transcriptomic level via multiplexed CRISPR-dCas9 and showed successful sensitization of E. coli without a separate fitness cost. These fitness neutral gene perturbations worked as co-therapies in reducing a Salmonella enterica intracellular infection in HeLa. Finally, these results informed the design of four antisense peptide nucleic acid (PNA) co-therapies, csgD, fnr, recA and acrA, against four MDR, clinically isolated bacteria. PNA combined with sub-minimal inhibitory concentrations of trimethoprim against two isolates of Klebsiella pneumoniae and E. coli showed three cases of re-sensitization with minimal fitness impacts. Our results highlight a promising approach for extending the utility of current antibiotics.
Collapse
|
18
|
Akiyama T, Kim M. Stochastic response of bacterial cells to antibiotics: its mechanisms and implications for population and evolutionary dynamics. Curr Opin Microbiol 2021; 63:104-108. [PMID: 34325154 DOI: 10.1016/j.mib.2021.07.002] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2021] [Revised: 06/22/2021] [Accepted: 07/01/2021] [Indexed: 11/20/2022]
Abstract
The effectiveness of antibiotics against bacterial infections has been declining due to the emergence of resistance. Precisely understanding the response of bacteria to antibiotics is critical to maximizing antibiotic-induced bacterial eradication while minimizing the emergence of antibiotic resistance. Cell-to-cell heterogeneity in antibiotic susceptibility is observed across various bacterial species for a wide range of antibiotics. Heterogeneity in antibiotic susceptibility is not always due to the genetic differences. Rather, it can be caused by non-genetic mechanisms such as stochastic gene expression and biased partitioning upon cell division. Heterogeneous susceptibility leads to the stochastic growth and death of individual cells and stochastic fluctuations in population size. These fluctuations have important implications for the eradication of bacterial populations and the emergence of genotypic resistance.
Collapse
Affiliation(s)
- Tatsuya Akiyama
- Department of Physics, Emory University, Atlanta, GA, 30322, USA; Graduate Division of Biological and Biomedical Sciences, Emory University, Atlanta, GA, 30322, USA
| | - Minsu Kim
- Department of Physics, Emory University, Atlanta, GA, 30322, USA; Graduate Division of Biological and Biomedical Sciences, Emory University, Atlanta, GA, 30322, USA; Emory Antibiotic Resistance Center, Emory University, Atlanta, GA, 30322, USA.
| |
Collapse
|
19
|
Barua A, Beygi A, Hatzikirou H. Close to Optimal Cell Sensing Ensures the Robustness of Tissue Differentiation Process: The Avian Photoreceptor Mosaic Case. ENTROPY 2021; 23:e23070867. [PMID: 34356408 PMCID: PMC8303396 DOI: 10.3390/e23070867] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Revised: 07/05/2021] [Accepted: 07/05/2021] [Indexed: 12/22/2022]
Abstract
The way that progenitor cell fate decisions and the associated environmental sensing are regulated to ensure the robustness of the spatial and temporal order in which cells are generated towards a fully differentiating tissue still remains elusive. Here, we investigate how cells regulate their sensing intensity and radius to guarantee the required thermodynamic robustness of a differentiated tissue. In particular, we are interested in finding the conditions where dedifferentiation at cell level is possible (microscopic reversibility), but tissue maintains its spatial order and differentiation integrity (macroscopic irreversibility). In order to tackle this, we exploit the recently postulated Least microEnvironmental Uncertainty Principle (LEUP) to develop a theory of stochastic thermodynamics for cell differentiation. To assess the predictive and explanatory power of our theory, we challenge it against the avian photoreceptor mosaic data. By calibrating a single parameter, the LEUP can predict the cone color spatial distribution in the avian retina and, at the same time, suggest that such a spatial pattern is associated with quasi-optimal cell sensing. By means of the stochastic thermodynamics formalism, we find out that thermodynamic robustness of differentiated tissues depends on cell metabolism and cell sensing properties. In turn, we calculate the limits of the cell sensing radius that ensure the robustness of differentiated tissue spatial order. Finally, we further constrain our model predictions to the avian photoreceptor mosaic.
Collapse
Affiliation(s)
- Arnab Barua
- Centre for Information Services and High Performance Computing, Technische Universität Dresden, Nöthnitzer Straße 46, 01062 Dresden, Germany; (A.B.); (A.B.)
| | - Alireza Beygi
- Centre for Information Services and High Performance Computing, Technische Universität Dresden, Nöthnitzer Straße 46, 01062 Dresden, Germany; (A.B.); (A.B.)
| | - Haralampos Hatzikirou
- Centre for Information Services and High Performance Computing, Technische Universität Dresden, Nöthnitzer Straße 46, 01062 Dresden, Germany; (A.B.); (A.B.)
- Mathematics Department, Khalifa University, Abu Dhabi P.O. Box 127788, United Arab Emirates
- Correspondence:
| |
Collapse
|
20
|
Farquhar KS, Rasouli Koohi S, Charlebois DA. Does transcriptional heterogeneity facilitate the development of genetic drug resistance? Bioessays 2021; 43:e2100043. [PMID: 34160842 DOI: 10.1002/bies.202100043] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Revised: 05/30/2021] [Accepted: 06/07/2021] [Indexed: 12/24/2022]
Abstract
Non-genetic forms of antimicrobial (drug) resistance can result from cell-to-cell variability that is not encoded in the genetic material. Data from recent studies also suggest that non-genetic mechanisms can facilitate the development of genetic drug resistance. We speculate on how the interplay between non-genetic and genetic mechanisms may affect microbial adaptation and evolution during drug treatment. We argue that cellular heterogeneity arising from fluctuations in gene expression, epigenetic modifications, as well as genetic changes contribute to drug resistance at different timescales, and that the interplay between these mechanisms enhance pathogen resistance. Accordingly, developing a better understanding of the role of non-genetic mechanisms in drug resistance and how they interact with genetic mechanisms will enhance our ability to combat antimicrobial resistance. Also see the video abstract here: https://youtu.be/aefGpdh-bgU.
Collapse
Affiliation(s)
| | - Samira Rasouli Koohi
- Department of Physics, University of Alberta, Edmonton, Alberta, T6G-2E1, Canada
| | - Daniel A Charlebois
- Department of Physics, University of Alberta, Edmonton, Alberta, T6G-2E1, Canada.,Department of Biological Sciences, University of Alberta, Edmonton, Alberta, Canada
| |
Collapse
|
21
|
Abstract
Temperature is an important determinant of bacterial growth. While the dependence of bacterial growth on different temperatures has been well studied for many bacterial species, prediction of bacterial growth rate for dynamic temperature changes is relatively unclear. Here, the authors address this issue using a combination of experimental measurements of the growth, at the resolution of 5 min, of Escherichia coli and mathematical models. They measure growth curves at different temperatures and estimate model parameters to predict bacterial growth profiles subject to dynamic temperature changes. They compared these predicted growth profiles for various step‐like temperature changes with experimental measurements using the coefficient of determination and mean square error and based on this comparison, ranked the different growth models, finding that the generalised logistic growth model gave the smallest error. They note that as the maximum specific growth increases the duration of this growth predominantly decreases. These results provide a basis to compute the dependence of the growth rate parameter in biomolecular circuits on dynamic temperatures and may be useful for designing biomolecular circuits that are robust to temperature.
Collapse
Affiliation(s)
- Abhishek Dey
- Department of Electrical Engineering, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Venkat Bokka
- Department of Electrical Engineering, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India
| | - Shaunak Sen
- Department of Electrical Engineering, Indian Institute of Technology Delhi, Hauz Khas, New Delhi 110016, India.
| |
Collapse
|
22
|
Castle SD, Grierson CS, Gorochowski TE. Towards an engineering theory of evolution. Nat Commun 2021; 12:3326. [PMID: 34099656 PMCID: PMC8185075 DOI: 10.1038/s41467-021-23573-3] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2021] [Accepted: 05/04/2021] [Indexed: 02/07/2023] Open
Abstract
Biological technologies are fundamentally unlike any other because biology evolves. Bioengineering therefore requires novel design methodologies with evolution at their core. Knowledge about evolution is currently applied to the design of biosystems ad hoc. Unless we have an engineering theory of evolution, we will neither be able to meet evolution's potential as an engineering tool, nor understand or limit its unintended consequences for our biological designs. Here, we propose the evotype as a helpful concept for engineering the evolutionary potential of biosystems, or other self-adaptive technologies, potentially beyond the realm of biology.
Collapse
Affiliation(s)
- Simeon D Castle
- School of Biological Sciences, University of Bristol, Bristol, UK
| | - Claire S Grierson
- School of Biological Sciences, University of Bristol, Bristol, UK
- BrisSynBio, University of Bristol, Bristol, UK
| | - Thomas E Gorochowski
- School of Biological Sciences, University of Bristol, Bristol, UK.
- BrisSynBio, University of Bristol, Bristol, UK.
| |
Collapse
|
23
|
Grandel NE, Reyes Gamas K, Bennett MR. Control of synthetic microbial consortia in time, space, and composition. Trends Microbiol 2021; 29:1095-1105. [PMID: 33966922 DOI: 10.1016/j.tim.2021.04.001] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Revised: 04/02/2021] [Accepted: 04/07/2021] [Indexed: 02/07/2023]
Abstract
While synthetic microbial systems are becoming increasingly complicated, single-strain systems cannot match the complexity of their multicellular counterparts. Such complexity, however, is much more difficult to control. Recent advances have increased our ability to control temporal, spatial, and community compositional organization, including modular adhesive systems, strain growth relationships, and asymmetric cell division. While these systems generally work independently, combining them into unified systems has proven difficult. Once such unification is proven successful we will unlock a new frontier of synthetic biology and open the door to the creation of synthetic biological systems with true multicellularity.
Collapse
Affiliation(s)
- Nicolas E Grandel
- Graduate Program in Systems, Synthetic, and Physical Biology, Rice University, Houston, TX, USA
| | - Kiara Reyes Gamas
- Graduate Program in Systems, Synthetic, and Physical Biology, Rice University, Houston, TX, USA
| | - Matthew R Bennett
- Department of Biosciences, Rice University, Houston, TX, USA; Department of Bioengineering, Rice University, Houston, TX, USA.
| |
Collapse
|
24
|
Capp J. Interplay between genetic, epigenetic, and gene expression variability: Considering complexity in evolvability. Evol Appl 2021; 14:893-901. [PMID: 33897810 PMCID: PMC8061278 DOI: 10.1111/eva.13204] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2020] [Revised: 01/27/2021] [Accepted: 01/30/2021] [Indexed: 12/11/2022] Open
Abstract
Genetic variability, epigenetic variability, and gene expression variability (noise) are generally considered independently in their relationship with phenotypic variation. However, they appear to be intrinsically interconnected and influence it in combination. The study of the interplay between genetic and epigenetic variability has the longest history. This article rather considers the introduction of gene expression variability in its relationships with the two others and reviews for the first time experimental evidences over the four relationships connected to gene expression noise. They show how introducing this third source of variability complicates the way of thinking evolvability and the emergence of biological novelty. Finally, cancer cells are proposed to be an ideal model to decipher the dynamic interplay between genetic, epigenetic, and gene expression variability when one of them is either experimentally increased or therapeutically targeted. This interplay is also discussed in an evolutionary perspective in the context of cancer cell drug resistance.
Collapse
Affiliation(s)
- Jean‐Pascal Capp
- Toulouse Biotechnology InstituteINSACNRSINRAEUniversity of ToulouseToulouseFrance
| |
Collapse
|
25
|
Critical Comparison of MaxCal and Other Stochastic Modeling Approaches in Analysis of Gene Networks. ENTROPY 2021; 23:e23030357. [PMID: 33802879 PMCID: PMC8002683 DOI: 10.3390/e23030357] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Revised: 03/09/2021] [Accepted: 03/10/2021] [Indexed: 11/24/2022]
Abstract
Learning the underlying details of a gene network with feedback is critical in designing new synthetic circuits. Yet, quantitative characterization of these circuits remains limited. This is due to the fact that experiments can only measure partial information from which the details of the circuit must be inferred. One potentially useful avenue is to harness hidden information from single-cell stochastic gene expression time trajectories measured for long periods of time—recorded at frequent intervals—over multiple cells. This raises the feasibility vs. accuracy dilemma while deciding between different models of mining these stochastic trajectories. We demonstrate that inference based on the Maximum Caliber (MaxCal) principle is the method of choice by critically evaluating its computational efficiency and accuracy against two other typical modeling approaches: (i) a detailed model (DM) with explicit consideration of multiple molecules including protein-promoter interaction, and (ii) a coarse-grain model (CGM) using Hill type functions to model feedback. MaxCal provides a reasonably accurate model while being significantly more computationally efficient than DM and CGM. Furthermore, MaxCal requires minimal assumptions since it is a top-down approach and allows systematic model improvement by including constraints of higher order, in contrast to traditional bottom-up approaches that require more parameters or ad hoc assumptions. Thus, based on efficiency, accuracy, and ability to build minimal models, we propose MaxCal as a superior alternative to traditional approaches (DM, CGM) when inferring underlying details of gene circuits with feedback from limited data.
Collapse
|
26
|
Yin H, Liu S, Wen X. Optimal transition paths of phenotypic switching in a non-Markovian self-regulation gene expression. Phys Rev E 2021; 103:022409. [PMID: 33736096 DOI: 10.1103/physreve.103.022409] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2020] [Accepted: 01/06/2021] [Indexed: 11/07/2022]
Abstract
Gene expression is a complex biochemical process involving multiple reaction steps, creating molecular memory because the probability of waiting time between consecutive reaction steps no longer follows exponential distributions. What effect the molecular memory has on metastable states in gene expression remains not fully understood. Here, we study transition paths of switching between bistable states for a non-Markovian model of gene expression equipped with a self-regulation. Employing the large deviation theory for this model, we analyze the optimal transition paths of switching between bistable states in gene expression, interestingly finding that dynamic behaviors in gene expression along the optimal transition paths significantly depend on the molecular memory. Moreover, we discover that the molecular memory can prolong the time of switching between bistable states in gene expression along the optimal transition paths. Our results imply that the molecular memory may be an unneglectable factor to affect switching between metastable states in gene expression.
Collapse
Affiliation(s)
- Hongwei Yin
- School of Mathematics and Statistics, Xuzhou University of Technology, Xuzhou 221111, People's Republic of China
- School of Science, Nanchang University, Nanchang 330031, People's Republic of China
| | - Shuqin Liu
- School of Science, Nanchang University, Nanchang 330031, People's Republic of China
| | - Xiaoqing Wen
- School of Mathematics and Statistics, Xuzhou University of Technology, Xuzhou 221111, People's Republic of China
- School of Science, Nanchang University, Nanchang 330031, People's Republic of China
| |
Collapse
|
27
|
Kuwahara H, Gao X. Stable maintenance of hidden switches as a strategy to increase the gene expression stability. NATURE COMPUTATIONAL SCIENCE 2021; 1:62-70. [PMID: 38217152 DOI: 10.1038/s43588-020-00001-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/30/2020] [Accepted: 11/02/2020] [Indexed: 01/15/2024]
Abstract
In response to severe genetic and environmental perturbations, wild-type organisms can express hidden alternative phenotypes adaptive to such adverse conditions. While our theoretical understanding of the population-level fitness advantage and evolution of phenotypic switching under variable environments has grown, the mechanism by which these organisms maintain phenotypic switching capabilities under static environments remains to be elucidated. Here, using computational simulations, we analyzed the evolution of gene circuits under natural selection and found that different strategies evolved to increase the gene expression stability near the optimum level. In a population comprising bistable individuals, a strategy of maintaining bistability and raising the potential barrier separating the bistable regimes was consistently taken. Our results serve as evidence that hidden bistable switches can be stably maintained during environmental stasis-an essential property enabling the timely release of adaptive alternatives with small genetic changes in the event of substantial perturbations.
Collapse
Affiliation(s)
- Hiroyuki Kuwahara
- Computational Bioscience Research Center (CBRC), Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia
| | - Xin Gao
- Computational Bioscience Research Center (CBRC), Computer, Electrical and Mathematical Sciences and Engineering (CEMSE) Division, King Abdullah University of Science and Technology (KAUST), Thuwal, Saudi Arabia.
| |
Collapse
|
28
|
Balázsi G. Discovering evolutionary hidden treasures. NATURE COMPUTATIONAL SCIENCE 2021; 1:18-19. [PMID: 38217151 DOI: 10.1038/s43588-020-00004-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/15/2024]
Affiliation(s)
- Gábor Balázsi
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, USA.
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, USA.
| |
Collapse
|
29
|
Roach TNF. Use and Abuse of Entropy in Biology: A Case for Caliber. ENTROPY 2020; 22:e22121335. [PMID: 33266519 PMCID: PMC7760317 DOI: 10.3390/e22121335] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 11/18/2020] [Accepted: 11/19/2020] [Indexed: 11/16/2022]
Abstract
Here, I discuss entropy and its use as a tool in fields of biology such as bioenergetics, ecology, and evolutionary biology. Statistical entropy concepts including Shannon's diversity, configurational entropy, and informational entropy are discussed in connection to their use in describing the diversity, heterogeneity, and spatial patterning of biological systems. The use of entropy as a measure of biological complexity is also discussed, and I explore the extension of thermodynamic entropy principles to open, nonequilibrium systems operating in finite time. I conclude with suggestions for use of caliber, a metric similar to entropy but for time-dependent trajectories rather than static distributions, and propose the complementary notion of path information.
Collapse
Affiliation(s)
- Ty N F Roach
- Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, Kāne'ohe, HI 96744, USA
| |
Collapse
|
30
|
Weistuch C, Agozzino L, Mujica-Parodi LR, Dill KA. Inferring a network from dynamical signals at its nodes. PLoS Comput Biol 2020; 16:e1008435. [PMID: 33253160 PMCID: PMC7728228 DOI: 10.1371/journal.pcbi.1008435] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2020] [Revised: 12/10/2020] [Accepted: 10/12/2020] [Indexed: 12/26/2022] Open
Abstract
We give an approximate solution to the difficult inverse problem of inferring the topology of an unknown network from given time-dependent signals at the nodes. For example, we measure signals from individual neurons in the brain, and infer how they are inter-connected. We use Maximum Caliber as an inference principle. The combinatorial challenge of high-dimensional data is handled using two different approximations to the pairwise couplings. We show two proofs of principle: in a nonlinear genetic toggle switch circuit, and in a toy neural network.
Collapse
Affiliation(s)
- Corey Weistuch
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York, USA
- Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, New York, USA
| | - Luca Agozzino
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York, USA
| | - Lilianne R. Mujica-Parodi
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York, USA
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York, USA
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, New York, USA
- Program in Neuroscience, Stony Brook University, Stony Brook, New York, USA
- Athinoula A. Martinos Center for Biomedical Imaging, Massachusetts General Hospital and Harvard Medical School, Boston, Massachusetts, USA
| | - Ken A. Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York, USA
- Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York, USA
- Department of Chemistry, Stony Brook University, Stony Brook, New York, USA
| |
Collapse
|
31
|
Guinn MT, Wan Y, Levovitz S, Yang D, Rosner MR, Balázsi G. Observation and Control of Gene Expression Noise: Barrier Crossing Analogies Between Drug Resistance and Metastasis. Front Genet 2020; 11:586726. [PMID: 33193723 PMCID: PMC7662081 DOI: 10.3389/fgene.2020.586726] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 09/29/2020] [Indexed: 12/12/2022] Open
Affiliation(s)
- Michael Tyler Guinn
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States.,Stony Brook Medical Scientist Training Program, Stony Brook, NY, United States
| | - Yiming Wan
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States
| | - Sarah Levovitz
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States
| | - Dongbo Yang
- Ben May Department for Cancer Research, The University of Chicago, Chicago, IL, United States
| | - Marsha R Rosner
- Ben May Department for Cancer Research, The University of Chicago, Chicago, IL, United States
| | - Gábor Balázsi
- Biomedical Engineering Department, Stony Brook University, Stony Brook, NY, United States.,Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States
| |
Collapse
|
32
|
Farquhar KS, Flohr H, Charlebois DA. Advancing Antimicrobial Resistance Research Through Quantitative Modeling and Synthetic Biology. Front Bioeng Biotechnol 2020; 8:583415. [PMID: 33072732 PMCID: PMC7530828 DOI: 10.3389/fbioe.2020.583415] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2020] [Accepted: 09/02/2020] [Indexed: 11/13/2022] Open
Abstract
Antimicrobial resistance (AMR) is an emerging global health crisis that is undermining advances in modern medicine and, if unmitigated, threatens to kill 10 million people per year worldwide by 2050. Research over the last decade has demonstrated that the differences between genetically identical cells in the same environment can lead to drug resistance. Fluctuations in gene expression, modulated by gene regulatory networks, can lead to non-genetic heterogeneity that results in the fractional killing of microbial populations causing drug therapies to fail; this non-genetic drug resistance can enhance the probability of acquiring genetic drug resistance mutations. Mathematical models of gene networks can elucidate general principles underlying drug resistance, predict the evolution of resistance, and guide drug resistance experiments in the laboratory. Cells genetically engineered to carry synthetic gene networks regulating drug resistance genes allow for controlled, quantitative experiments on the role of non-genetic heterogeneity in the development of drug resistance. In this perspective article, we emphasize the contributions that mathematical, computational, and synthetic gene network models play in advancing our understanding of AMR to discover effective therapies against drug-resistant infections.
Collapse
Affiliation(s)
| | - Harold Flohr
- Department of Physics, University of Alberta, Edmonton, AB, Canada
| | | |
Collapse
|
33
|
Abstract
Cells adapt to changing environments. Perturb a cell and it returns to a point of homeostasis. Perturb a population and it evolves toward a fitness peak. We review quantitative models of the forces of adaptation and their visualizations on landscapes. While some adaptations result from single mutations or few-gene effects, others are more cooperative, more delocalized in the genome, and more universal and physical. For example, homeostasis and evolution depend on protein folding and aggregation, energy and protein production, protein diffusion, molecular motor speeds and efficiencies, and protein expression levels. Models provide a way to learn about the fitness of cells and cell populations by making and testing hypotheses.
Collapse
Affiliation(s)
- Luca Agozzino
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, USA; .,Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, USA
| | - Gábor Balázsi
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, USA; .,Department of Biomedical Engineering, Stony Brook University, Stony Brook, New York 11794, USA
| | - Jin Wang
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, USA; .,Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, USA.,Department of Chemistry, Stony Brook University, Stony Brook, New York 11790, USA
| | - Ken A Dill
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, USA; .,Department of Physics and Astronomy, Stony Brook University, Stony Brook, New York 11794, USA.,Department of Chemistry, Stony Brook University, Stony Brook, New York 11790, USA
| |
Collapse
|
34
|
Lopatkin AJ, Collins JJ. Predictive biology: modelling, understanding and harnessing microbial complexity. Nat Rev Microbiol 2020; 18:507-520. [DOI: 10.1038/s41579-020-0372-5] [Citation(s) in RCA: 50] [Impact Index Per Article: 12.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/15/2020] [Indexed: 12/11/2022]
|
35
|
Jiang Y, AkhavanAghdam Z, Li Y, Zid BM, Hao N. A protein kinase A-regulated network encodes short- and long-lived cellular memories. Sci Signal 2020; 13:eaay3585. [PMID: 32430291 PMCID: PMC7302112 DOI: 10.1126/scisignal.aay3585] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Cells can store memories of prior experiences to modulate their responses to subsequent stresses, as seen for the protein kinase A (PKA)-mediated general stress response in yeast, which is required for resistance against future stressful conditions. Using microfluidics and time-lapse microscopy, we quantitatively analyzed how the cellular memory of stress adaptation is encoded in single yeast cells. We found that cellular memory was biphasic. Short-lived memory was mediated by trehalose synthase and trehalose metabolism. Long-lived memory was mediated by PKA-regulated stress-responsive transcription factors and cytoplasmic messenger ribonucleoprotein granules. Short- and long-lived memory could be selectively induced by different priming input dynamics. Computational modeling revealed how the PKA-mediated regulatory network could encode previous stimuli into memories with distinct dynamics. This biphasic memory-encoding scheme might represent a general strategy to prepare for future challenges in rapidly changing environments.
Collapse
Affiliation(s)
- Yanfei Jiang
- Section of Molecular Biology, Division of Biological Sciences, University of California San Diego, La Jolla, CA 92093, USA
| | - Zohreh AkhavanAghdam
- Section of Molecular Biology, Division of Biological Sciences, University of California San Diego, La Jolla, CA 92093, USA
| | - Yutian Li
- Section of Molecular Biology, Division of Biological Sciences, University of California San Diego, La Jolla, CA 92093, USA
| | - Brian M Zid
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, CA 92093, USA
| | - Nan Hao
- Section of Molecular Biology, Division of Biological Sciences, University of California San Diego, La Jolla, CA 92093, USA.
| |
Collapse
|
36
|
Abstract
Ever since Clausius in 1865 and Boltzmann in 1877, the concepts of entropy and of its maximization have been the foundations for predicting how material equilibria derive from microscopic properties. But, despite much work, there has been no equally satisfactory general variational principle for nonequilibrium situations. However, in 1980, a new avenue was opened by E.T. Jaynes and by Shore and Johnson. We review here maximum caliber, which is a maximum-entropy-like principle that can infer distributions of flows over pathways, given dynamical constraints. This approach is providing new insights, particularly into few-particle complex systems, such as gene circuits, protein conformational reaction coordinates, network traffic, bird flocking, cell motility, and neuronal firing.
Collapse
Affiliation(s)
- Kingshuk Ghosh
- Department of Physics and Astronomy, University of Denver, Denver, Colorado 80209, USA
| | - Purushottam D. Dixit
- Department of Systems Biology, Columbia University, New York, NY 10032, USA,Department of Physics, University of Florida, Gainesville, Florida 32611, USA
| | - Luca Agozzino
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, USA
| | - Ken A. Dill
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
37
|
Zhang R, Li J, Melendez-Alvarez J, Chen X, Sochor P, Goetz H, Zhang Q, Ding T, Wang X, Tian XJ. Topology-dependent interference of synthetic gene circuit function by growth feedback. Nat Chem Biol 2020; 16:695-701. [PMID: 32251409 PMCID: PMC7246135 DOI: 10.1038/s41589-020-0509-x] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2019] [Accepted: 02/28/2020] [Indexed: 11/21/2022]
Abstract
Growth-mediated feedback between synthetic gene circuits and host organisms leads to diverse emerged behaviors, including growth bistability and enhanced ultrasensitivity. However, the range of possible impacts of growth feedback on gene circuits remains underexplored. Here, we mathematically and experimentally demonstrated that growth feedback affects the functions of memory circuits in a network topology-dependent way. Specifically, the memory of the self-activation switch is quickly lost due to the growth-mediated dilution of the circuit products. Decoupling of growth feedback reveals its memory, manifested by its hysteresis property across a broad range of inducer concentration. On the contrary, the toggle switch is more refractory to growth-mediated dilution and can retrieve its memory after the fast-growth phase. The underlying principle lies in the different dependence of active and repressive regulations in these circuits on the growth-mediated dilution. Our results unveil the topology-dependent mechanism on how growth-mediated feedback influences the behaviors of gene circuits.
Collapse
Affiliation(s)
- Rong Zhang
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA
| | - Jiao Li
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA.,Department of Food Science and Nutrition, Zhejiang University, Hangzhou, Zhejiang, China
| | - Juan Melendez-Alvarez
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA
| | - Xingwen Chen
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA
| | - Patrick Sochor
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA
| | - Hanah Goetz
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA
| | - Qi Zhang
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA
| | - Tian Ding
- Department of Food Science and Nutrition, Zhejiang University, Hangzhou, Zhejiang, China
| | - Xiao Wang
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA.
| | - Xiao-Jun Tian
- School of Biological and Health Systems Engineering, Arizona State University, Tempe, AZ, USA.
| |
Collapse
|
38
|
Abstract
Natural or synthetic genetic modules can lose their function over long-term evolution if the function is costly. How populations can evolve to restore such broken function is poorly understood. To test the reversibility of evolutionary breakdown, we use yeast cell populations with a chromosomally integrated synthetic gene circuit. In previous evolution experiments the gene circuit lost its costly function through various mutations. By exposing such mutant populations to conditions where regaining gene circuit function would be beneficial, we find adaptation scenarios with or without repairing lost gene circuit function. These results are important for drug resistance or future synthetic biology applications where evolutionary loss and regain of function play a significant role. Evolutionary reversibility—the ability to regain a lost function—is an important problem both in evolutionary and synthetic biology, where repairing natural or synthetic systems broken by evolutionary processes may be valuable. Here, we use a synthetic positive-feedback (PF) gene circuit integrated into haploid Saccharomyces cerevisiae cells to test if the population can restore lost PF function. In previous evolution experiments, mutations in a gene eliminated the fitness costs of PF activation. Since PF activation also provides drug resistance, exposing such compromised or broken mutants to both drug and inducer should create selection pressure to regain drug resistance and possibly PF function. Indeed, evolving 7 PF mutant strains in the presence of drug revealed 3 adaptation scenarios through genomic, PF-external mutations that elevate PF basal expression, possibly by affecting transcription, translation, degradation, and other fundamental cellular processes. Nonfunctional mutants gained drug resistance without ever developing high expression, while quasifunctional and dysfunctional PF mutants developed high expression nongenetically, which then diminished, although more slowly for dysfunctional mutants where revertant clones arose. These results highlight how intracellular context, such as the growth rate, can affect regulatory network dynamics and evolutionary dynamics, which has important consequences for understanding the evolution of drug resistance and developing future synthetic biology applications.
Collapse
|
39
|
Farquhar KS, Charlebois DA, Szenk M, Cohen J, Nevozhay D, Balázsi G. Role of network-mediated stochasticity in mammalian drug resistance. Nat Commun 2019; 10:2766. [PMID: 31235692 PMCID: PMC6591227 DOI: 10.1038/s41467-019-10330-w] [Citation(s) in RCA: 63] [Impact Index Per Article: 12.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Accepted: 05/03/2019] [Indexed: 11/11/2022] Open
Abstract
A major challenge in biology is that genetically identical cells in the same environment can display gene expression stochasticity (noise), which contributes to bet-hedging, drug tolerance, and cell-fate switching. The magnitude and timescales of stochastic fluctuations can depend on the gene regulatory network. Currently, it is unclear how gene expression noise of specific networks impacts the evolution of drug resistance in mammalian cells. Answering this question requires adjusting network noise independently from mean expression. Here, we develop positive and negative feedback-based synthetic gene circuits to decouple noise from the mean for Puromycin resistance gene expression in Chinese Hamster Ovary cells. In low Puromycin concentrations, the high-noise, positive-feedback network delays long-term adaptation, whereas it facilitates adaptation under high Puromycin concentration. Accordingly, the low-noise, negative-feedback circuit can maintain resistance by acquiring mutations while the positive-feedback circuit remains mutation-free and regains drug sensitivity. These findings may have profound implications for chemotherapeutic inefficiency and cancer relapse. The role of gene expression noise in the evolution of drug resistance in mammalian cells is unclear. Here, by uncoupling noise from mean expression of a drug resistance gene in CHO cells the authors show that noisy expression aids adaptation to high drug levels, but delays it at low drug levels.
Collapse
Affiliation(s)
- Kevin S Farquhar
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA.,Genetics and Epigenetics Graduate Program, The University of Texas MD Anderson Cancer Center, UT Health Graduate School of Biomedical Sciences, Houston, TX, 77030, USA
| | - Daniel A Charlebois
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA.,Department of Physics, University of Alberta, Edmonton, AB, 4-181 CCIS, T6G-2E1, Canada
| | - Mariola Szenk
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA.,Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, 11794, USA
| | - Joseph Cohen
- Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, 11794, USA
| | - Dmitry Nevozhay
- School of Biomedicine, Far Eastern Federal University, 8 Sukhanova Street, Vladivostok, 690950, Russia.,Department of Systems Biology, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA
| | - Gábor Balázsi
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, 11794, USA. .,Genetics and Epigenetics Graduate Program, The University of Texas MD Anderson Cancer Center, UT Health Graduate School of Biomedical Sciences, Houston, TX, 77030, USA. .,Department of Biomedical Engineering, Stony Brook University, Stony Brook, NY, 11794, USA. .,Department of Systems Biology, The University of Texas MD Anderson Cancer Center, Houston, TX, 77030, USA.
| |
Collapse
|
40
|
Roy A, Klumpp S. Simulating Genetic Circuits in Bacterial Populations with Growth Heterogeneity. Biophys J 2019; 114:484-492. [PMID: 29401445 DOI: 10.1016/j.bpj.2017.11.3745] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2017] [Revised: 10/31/2017] [Accepted: 11/17/2017] [Indexed: 01/21/2023] Open
Abstract
We computationally study genetic circuits in bacterial populations with heterogeneities in the growth rate. To that end, we present a stochastic simulation method for gene circuits in populations of cells and propose an efficient implementation that we call the "Next Family Method". Within this approach, we implement different population setups, specifically Chemostat-type growth and growth in an ideal Mother Machine and show that the population structure and its statistics are different for the different setups whenever there is growth heterogeneity. Such dependence on the population setup is demonstrated, in the case of bistable systems with different growth rates in the stable states, to have distinctive signatures on quantities including the distributions of protein concentration and growth rates, and hysteresis curves. Applying this method to a bistable antibiotic resistance circuit, we find that as a result of the different statistics in different population setups, the estimated minimal inhibitory concentration of the antibiotic becomes dependent on the population setup in which it is measured.
Collapse
Affiliation(s)
- Anjan Roy
- Max Planck Institute of Colloids and Interfaces, Potsdam, Germany; Institute for Nonlinear Dynamics, University of Göttingen, Göttingen, Germany
| | - Stefan Klumpp
- Max Planck Institute of Colloids and Interfaces, Potsdam, Germany; Institute for Nonlinear Dynamics, University of Göttingen, Göttingen, Germany.
| |
Collapse
|
41
|
Firman T, Amgalan A, Ghosh K. Maximum Caliber Can Build and Infer Models of Oscillation in a Three-Gene Feedback Network. J Phys Chem B 2019; 123:343-355. [PMID: 30507199 DOI: 10.1021/acs.jpcb.8b07465] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Single-cell protein expression time trajectories provide rich temporal data quantifying cellular variability and its role in dictating fitness. However, theoretical models to analyze and fully extract information from these measurements remain limited for three reasons: (i) gene expression profiles are noisy, rendering models of averages inapplicable, (ii) experiments typically measure only a few protein species while leaving other molecular actors-necessary to build traditional bottom-up models-unnoticed, and (iii) measured data are in fluorescence, not particle number. We recently addressed these challenges in an alternate top-down approach using the principle of Maximum Caliber (MaxCal) to model genetic switches with one and two protein species. In the present work we address scalability and broader applicability of MaxCal by extending to a three-gene (A, B, C) feedback network that exhibits oscillation, commonly known as the repressilator. We test MaxCal's inferential power by using synthetic data of noisy protein number time traces-serving as a proxy for experimental data-generated from a known underlying model. We notice that the minimal MaxCal model-accounting for production, degradation, and only one type of symmetric coupling between all three species-reasonably infers several underlying features of the circuit such as the effective production rate, degradation rate, frequency of oscillation, and protein number distribution. Next, we build models of higher complexity including different levels of coupling between A, B, and C and rigorously assess their relative performance. While the minimal model (with four parameters) performs remarkably well, we note that the most complex model (with six parameters) allowing all possible forms of crosstalk between A, B, and C slightly improves prediction of rates, but avoids ad hoc assumption of all the other models. It is also the model of choice based on Bayesian information criteria. We further analyzed time trajectories in arbitrary fluorescence (using synthetic trajectories) to mimic realistic data. We conclude that even with a three-protein system including both fluorescence noise and intrinsic gene expression fluctuations, MaxCal can faithfully infer underlying details of the network, opening future directions to model other network motifs with many species.
Collapse
|
42
|
Charlebois DA, Balázsi G. Modeling cell population dynamics. In Silico Biol 2019; 13:21-39. [PMID: 30562900 PMCID: PMC6598210 DOI: 10.3233/isb-180470] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2018] [Revised: 09/13/2018] [Accepted: 10/16/2018] [Indexed: 12/27/2022]
Abstract
Quantitative modeling is quickly becoming an integral part of biology, due to the ability of mathematical models and computer simulations to generate insights and predict the behavior of living systems. Single-cell models can be incapable or misleading for inferring population dynamics, as they do not consider the interactions between cells via metabolites or physical contact, nor do they consider competition for limited resources such as nutrients or space. Here we examine methods that are commonly used to model and simulate cell populations. First, we cover simple models where analytic solutions are available, and then move on to more complex scenarios where computational methods are required. Overall, we present a summary of mathematical models used to describe cell population dynamics, which may aid future model development and highlights the importance of population modeling in biology.
Collapse
Affiliation(s)
- Daniel A. Charlebois
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, NY, USA
- Department of Physics, University of Alberta, Edmonton, AB, Canada
| | - Gábor Balázsi
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, NY, USA
- Department of Biomedical Engineering, Stony Brook University, NY, USA
| |
Collapse
|
43
|
Abstract
Quantitative modeling is quickly becoming an integral part of biology, due to the ability of mathematical models and computer simulations to generate insights and predict the behavior of living systems. Single-cell models can be incapable or misleading for inferring population dynamics, as they do not consider the interactions between cells via metabolites or physical contact, nor do they consider competition for limited resources such as nutrients or space. Here we examine methods that are commonly used to model and simulate cell populations. First, we cover simple models where analytic solutions are available, and then move on to more complex scenarios where computational methods are required. Overall, we present a summary of mathematical models used to describe cell population dynamics, which may aid future model development and highlights the importance of population modeling in biology.
Collapse
Affiliation(s)
- Daniel A Charlebois
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, NY, USA.,Department of Physics, University of Alberta, Edmonton, AB, Canada
| | - Gábor Balázsi
- The Louis and Beatrice Laufer Center for Physical and Quantitative Biology, Stony Brook University, NY, USA.,Department of Biomedical Engineering, Stony Brook University, NY, USA
| |
Collapse
|
44
|
Firman T, Amgalan A, Ghosh K. Maximum Caliber Can Build and Infer Models of Oscillation in a Three-Gene Feedback Network. J Phys Chem A 2018. [DOI: 10.1021/acs.jpca.8b07465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
45
|
Meredith HR, Andreani V, Ma HR, Lopatkin AJ, Lee AJ, Anderson DJ, Batt G, You L. Applying ecological resistance and resilience to dissect bacterial antibiotic responses. SCIENCE ADVANCES 2018; 4:eaau1873. [PMID: 30525104 PMCID: PMC6281428 DOI: 10.1126/sciadv.aau1873] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/15/2018] [Accepted: 11/07/2018] [Indexed: 05/14/2023]
Abstract
An essential property of microbial communities is the ability to survive a disturbance. Survival can be achieved through resistance, the ability to absorb effects of a disturbance without a notable change, or resilience, the ability to recover after being perturbed by a disturbance. These concepts have long been applied to the analysis of ecological systems, although their interpretations are often subject to debate. Here, we show that this framework readily lends itself to the dissection of the bacterial response to antibiotic treatment, where both terms can be unambiguously defined. The ability to tolerate the antibiotic treatment in the short term corresponds to resistance, which primarily depends on traits associated with individual cells. In contrast, the ability to recover after being perturbed by an antibiotic corresponds to resilience, which primarily depends on traits associated with the population. This framework effectively reveals the phenotypic signatures of bacterial pathogens expressing extended-spectrum β-lactamases (ESBLs) when treated by a β-lactam antibiotic. Our analysis has implications for optimizing treatment of these pathogens using a combination of a β-lactam and a β-lactamase (Bla) inhibitor. In particular, our results underscore the need to dynamically optimize combination treatments based on the quantitative features of the bacterial response to the antibiotic or the Bla inhibitor.
Collapse
Affiliation(s)
| | - Virgile Andreani
- Inria Saclay–Île-de-France, Palaiseau, France
- Institut Pasteur, Paris, France
| | - Helena R. Ma
- Department of Biomedical Engineering, Duke University, Durham, NC, USA
| | | | - Anna J. Lee
- Department of Biomedical Engineering, Duke University, Durham, NC, USA
| | - Deverick J. Anderson
- Division of Infectious Diseases, Department of Medicine, Duke University School of Medicine, Durham, NC, USA
- Duke Center for Antimicrobial Stewardship and Infection Prevention, Duke University School of Medicine, Durham, NC, USA
| | - Gregory Batt
- Inria Saclay–Île-de-France, Palaiseau, France
- Institut Pasteur, Paris, France
| | - Lingchong You
- Department of Biomedical Engineering, Duke University, Durham, NC, USA
- Center for Genomic and Computational Biology, Duke University, Durham, NC, USA
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA
| |
Collapse
|
46
|
Abstract
Fluctuating environments such as changes in ambient temperature represent a fundamental challenge to life. Cells must protect gene networks that protect them from such stresses, making it difficult to understand how temperature affects gene network function in general. Here, we focus on single genes and small synthetic network modules to reveal four key effects of nonoptimal temperatures at different biological scales: (i) a cell fate choice between arrest and resistance, (ii) slower growth rates, (iii) Arrhenius reaction rates, and (iv) protein structure changes. We develop a multiscale computational modeling framework that captures and predicts all of these effects. These findings promote our understanding of how temperature affects living systems and enables more robust cellular engineering for real-world applications. Most organisms must cope with temperature changes. This involves genes and gene networks both as subjects and agents of cellular protection, creating difficulties in understanding. Here, we study how heating and cooling affect expression of single genes and synthetic gene circuits in Saccharomyces cerevisiae. We discovered that nonoptimal temperatures induce a cell fate choice between stress resistance and growth arrest. This creates dramatic gene expression bimodality in isogenic cell populations, as arrest abolishes gene expression. Multiscale models incorporating population dynamics, temperature-dependent growth rates, and Arrhenius scaling of reaction rates captured the effects of cooling, but not those of heating in resistant cells. Molecular-dynamics simulations revealed how heating alters the conformational dynamics of the TetR repressor, fully explaining the experimental observations. Overall, nonoptimal temperatures induce a cell fate decision and corrupt gene and gene network function in computationally predictable ways, which may aid future applications of engineered microbes in nonstandard temperatures.
Collapse
|
47
|
Otoupal PB, Cordell WT, Bachu V, Sitton MJ, Chatterjee A. Multiplexed deactivated CRISPR-Cas9 gene expression perturbations deter bacterial adaptation by inducing negative epistasis. Commun Biol 2018; 1:129. [PMID: 30272008 PMCID: PMC6123780 DOI: 10.1038/s42003-018-0135-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2018] [Accepted: 08/08/2018] [Indexed: 12/21/2022] Open
Abstract
The ever-increasing threat of multi-drug resistant bacteria, a shrinking antibiotic pipeline, and the innate ability of microorganisms to adapt necessitates long-term strategies to slow the evolution of antibiotic resistance. Here we develop an approach, dubbed Controlled Hindrance of Adaptation of OrganismS or CHAOS, involving induction of epistasis between gene perturbations to deter adaption. We construct a combinatorial library of multiplexed, deactivated CRISPR-Cas9 devices to systematically perturb gene expression in Escherichia coli. While individual perturbations improved fitness during antibiotic exposure, multiplexed perturbations caused large fitness loss in a significant epistatic fashion. Strains exhibiting epistasis adapted significantly more slowly over three to fourteen days, and loss in adaptive potential was shown to be sustainable. Finally, we show that multiplexed peptide nucleic acids increase the antibiotic susceptibility of clinically isolated Carbapenem-resistant E. coli in an epistatic fashion. Together, these results suggest a new therapeutic strategy for restricting the evolution of antibiotic resistance. Peter Otoupal et al. present CHAOS, an approach for preventing the development of antibiotic resistance in E. coli through CRISPR-Cas9-based perturbation of gene expression. They show that multiplexed perturbations decrease fitness of clinically-isolated Carbapenem-resistant E. coli upon antibiotic exposure.
Collapse
Affiliation(s)
- Peter B Otoupal
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, CO, 80303, USA
| | - William T Cordell
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, CO, 80303, USA
| | - Vismaya Bachu
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, CO, 80303, USA
| | - Madeleine J Sitton
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, CO, 80303, USA
| | - Anushree Chatterjee
- Department of Chemical and Biological Engineering, University of Colorado at Boulder, Boulder, CO, 80303, USA. .,BioFrontiers Institute, University of Colorado at Boulder, Boulder, CO, 80303, USA.
| |
Collapse
|
48
|
Bandyopadhyay A, Wang H, Ray JCJ. Lineage space and the propensity of bacterial cells to undergo growth transitions. PLoS Comput Biol 2018; 14:e1006380. [PMID: 30133447 PMCID: PMC6122811 DOI: 10.1371/journal.pcbi.1006380] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2018] [Revised: 09/04/2018] [Accepted: 07/19/2018] [Indexed: 11/18/2022] Open
Abstract
The molecular makeup of the offspring of a dividing cell gradually becomes phenotypically decorrelated from the parent cell by noise and regulatory mechanisms that amplify phenotypic heterogeneity. Such regulatory mechanisms form networks that contain thresholds between phenotypes. Populations of cells can be poised near the threshold so that a subset of the population probabilistically undergoes the phenotypic transition. We sought to characterize the diversity of bacterial populations around a growth-modulating threshold via analysis of the effect of non-genetic inheritance, similar to conditions that create antibiotic-tolerant persister cells and other examples of bet hedging. Using simulations and experimental lineage data in Escherichia coli, we present evidence that regulation of growth amplifies the dependence of growth arrest on cellular lineage, causing clusters of related cells undergo growth arrest in certain conditions. Our simulations predict that lineage correlations and the sensitivity of growth to changes in toxin levels coincide in a critical regime. Below the critical regime, the sizes of related growth arrested clusters are distributed exponentially, while in the critical regime clusters sizes are more likely to become large. Furthermore, phenotypic diversity can be nearly as high as possible near the critical regime, but for most parameter values it falls far below the theoretical limit. We conclude that lineage information is indispensable for understanding regulation of cellular growth.
Collapse
Affiliation(s)
- Arnab Bandyopadhyay
- Center for Computational Biology, Department of Molecular Biosciences, University of Kansas, Lawrence, KS United States of America
| | - Huijing Wang
- Center for Computational Biology, Department of Molecular Biosciences, University of Kansas, Lawrence, KS United States of America
| | - J. Christian J. Ray
- Center for Computational Biology, Department of Molecular Biosciences, University of Kansas, Lawrence, KS United States of America
- * E-mail:
| |
Collapse
|
49
|
Scott LH, Mathews JC, Flematti GR, Filipovska A, Rackham O. An Artificial Yeast Genetic Circuit Enables Deep Mutational Scanning of an Antimicrobial Resistance Protein. ACS Synth Biol 2018; 7:1907-1917. [PMID: 29979580 DOI: 10.1021/acssynbio.8b00121] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Understanding the molecular mechanisms underlying antibiotic resistance requires concerted efforts in enzymology and medicinal chemistry. Here we describe a new synthetic biology approach to antibiotic development, where the presence of tetracycline antibiotics is linked to a life-death selection in Saccharomyces cerevisiae. This artificial genetic circuit allowed the deep mutational scanning of the tetracycline inactivating enzyme TetX, revealing key functional residues. We used both positive and negative selections to confirm the importance of different residues for TetX activity, and profiled activity hotspots for different tetracyclines to reveal substrate-specific activity determinants. We found that precise positioning of FAD and hydrophobic shielding of the tetracycline are critical for enzymatic inactivation of doxycycline. However, positioning of FAD is suboptimal in the case of anhydrotetracycline, potentially explaining its comparatively poor degradation and potential as an inhibitor for this family of enzymes. By combining artificial genetic circuits whose function can be modulated by antimicrobial resistance determinants, we establish a framework to select for the next generation of antibiotics.
Collapse
Affiliation(s)
- Louis H. Scott
- Harry Perkins Institute of Medical Research and Centre for Medical Research, The University of Western Australia, Nedlands 6009, Australia
| | - James C. Mathews
- Harry Perkins Institute of Medical Research and Centre for Medical Research, The University of Western Australia, Nedlands 6009, Australia
| | - Gavin R. Flematti
- School of Molecular Sciences, The University of Western Australia, Crawley 6009, Australia
| | - Aleksandra Filipovska
- Harry Perkins Institute of Medical Research and Centre for Medical Research, The University of Western Australia, Nedlands 6009, Australia
- School of Molecular Sciences, The University of Western Australia, Crawley 6009, Australia
| | - Oliver Rackham
- Harry Perkins Institute of Medical Research and Centre for Medical Research, The University of Western Australia, Nedlands 6009, Australia
- School of Molecular Sciences, The University of Western Australia, Crawley 6009, Australia
| |
Collapse
|
50
|
Abstract
Systems biology seeks to combine experiments with computation to predict biological behaviors. However, despite tremendous data and knowledge, biological models make less-accurate predictions compared with other fields. By analyzing single-cell, single-molecule measurements of mRNA during yeast stress response, we explore why and how the shapes of experimental distributions control prediction accuracy. We show how asymmetric data distributions with long tails cause standard modeling approaches to yield excellent fits but make meaningless predictions. We show how these biases arise from the violation of fundamental assumptions in standard modeling approaches. We demonstrate how advanced computational tools solve this dilemma and achieve predictive understanding of spatiotemporal mechanisms of transcription control including RNA polymerase initiation and elongation and mRNA accumulation, transport, and decay. Despite substantial experimental and computational efforts, mechanistic modeling remains more predictive in engineering than in systems biology. The reason for this discrepancy is not fully understood. One might argue that the randomness and complexity of biological systems are the main barriers to predictive understanding, but these issues are not unique to biology. Instead, we hypothesize that the specific shapes of rare single-molecule event distributions produce substantial yet overlooked challenges for biological models. We demonstrate why modern statistical tools to disentangle complexity and stochasticity, which assume normally distributed fluctuations or enormous datasets, do not apply to the discrete, positive, and nonsymmetric distributions that characterize mRNA fluctuations in single cells. As an example, we integrate single-molecule measurements and advanced computational analyses to explore mitogen-activated protein kinase induction of multiple stress response genes. Through systematic analyses of different metrics to compare the same model to the same data, we elucidate why standard modeling approaches yield nonpredictive models for single-cell gene regulation. We further explain how advanced tools recover precise, reproducible, and predictive understanding of transcription regulation mechanisms, including gene activation, polymerase initiation, elongation, mRNA accumulation, spatial transport, and decay.
Collapse
|