Quorum-Sensing Synchronization of Synthetic Toggle Switches: A Design based on Monotone Dynamical Systems Theory

Synthetic constructs in biotechnology, biocomputing, and modern gene therapy interventions are often based on plasmids or transfected circuits which implement some form of “on-off” switch. For example, the expression of a protein used for therapeutic purposes might be triggered by the recognition of a specific combination of inducers (e.g., antigens), and memory of this event should be maintained across a cell population until a specific stimulus commands a coordinated shut-off. The robustness of such a design is hampered by molecular (“intrinsic”) or environmental (“extrinsic”) noise, which may lead to spontaneous changes of state in a subset of the population and is reflected in the bimodality of protein expression, as measured for example using flow cytometry. In this context, a “majority-vote” correction circuit, which brings deviant cells back into the desired state, is highly desirable, and quorum-sensing has been suggested as a way for cells to broadcast their states to the population as a whole so as to facilitate consensus. In this paper, we propose what we believe is the first such design that has mathematically guaranteed properties of stability and auto-correction. Our approach is guided by concepts and theory from the field of “monotone” dynamical systems developed by M. Hirsch, H. Smith, and others. We benchmark our design by comparing it to an existing design which has been the subject of experimental and theoretical studies, illustrating its superiority in stability and self-correction of synchronization errors. Our stability analysis, based on dynamical systems theory, guarantees global convergence to steady states, ruling out unpredictable (“chaotic”) behaviors and even sustained oscillations. These results are valid no matter what are the values of parameters, and are based only on the wiring diagram. The theory is complemented by extensive computational bifurcation analysis, performed for a biochemically-detailed and biologically-relevant model that we developed. Another novel feature of our approach is that our theorems on stability of steady states for homogeneous or mixed populations are valid independently of the number of cells in the population, and depend only on the relative proportions of each type of state. While monotone systems theory has been used previously for systems biology analysis, the current work illustrates its power for synthetic biology design, and thus has wider significance well beyond the application to the important problem of coordination of toggle switches. Author Summary For the last decade, outstanding progress has been made, and considerable practical experience has accumulated, in the construction of elementary genetic circuits that perform various tasks, such as memory storage and logical operations, in response to both exogenous and endogenous stimuli. Using modern molecular “plug-and-play” technologies, various (re-)programmable cellular populations can be engineered in a routine manner, and they can be combined into more complex cellular systems. Among all engineered synthetic circuits, a toggle, a robust bistable switch leading to a binary response dynamics, is the simplest basic synthetic biology device, analogous to the “flip-flop” or latch in electronic design, and it plays a key role in biotechnology, biocomputing, and proposed gene therapies. However, despite many remarkable properties of the existing toggle designs, they must be tightly controlled in order to avoid spontaneous switching between different expression states (loss of long-term memory) or even the breakdown of stability through the generation of stable oscillations. To address this concrete challenge, we have developed a new design for quorum-sensing synthetic toggles, based on monotone dynamical systems theory. Our design is endowed with strong theoretical guarantees that completely exclude unpredictable chaotic behaviors, as well as undesired stable oscillations, and leads to robust consensus states.

[1]  J. A. Kuznecov Elements of applied bifurcation theory , 1998 .

[2]  Ido Golding,et al.  Single-molecule studies of repressor-DNA interactions show long-range interactions. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[3]  J. Collins,et al.  Construction of a genetic toggle switch in Escherichia coli , 2000, Nature.

[4]  Kresimir Josic,et al.  Engineered temperature compensation in a synthetic genetic clock , 2014, Proceedings of the National Academy of Sciences.

[5]  J. Keener,et al.  A mathematical model for quorum sensing in Pseudomonas aeruginosa , 2001, Bulletin of mathematical biology.

[6]  Peter Lancaster,et al.  The theory of matrices , 1969 .

[7]  B. Bassler,et al.  Quorum sensing in bacteria. , 2001, Annual review of microbiology.

[8]  Guy Katriel,et al.  Synchronization of oscillators coupled through an environment , 2008, 0804.3734.

[9]  Yiannis N Kaznessis,et al.  Model-driven designs of an oscillating gene network. , 2005, Biophysical journal.

[10]  M M Domach,et al.  Computer model for glucose‐limited growth of a single cell of Escherichia coli B/r‐A , 1984, Biotechnology and bioengineering.

[11]  W. Hillen,et al.  Lack of a 5' non-coding region in Tn1721 encoded tetR mRNA is associated with a low efficiency of translation and a short half-life in Escherichia coli. , 1991, Nucleic acids research.

[12]  P. Seed,et al.  Regulation of las and rhl quorum sensing in Pseudomonas aeruginosa , 1997, Journal of bacteriology.

[13]  Willy Govaerts,et al.  MATCONT: A MATLAB package for numerical bifurcation analysis of ODEs , 2003, TOMS.

[14]  G Ambika,et al.  Synchronized states in chaotic systems coupled indirectly through a dynamic environment. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  Priscilla E. M. Purnick,et al.  The second wave of synthetic biology: from modules to systems , 2009, Nature Reviews Molecular Cell Biology.

[16]  K. P. Hadeler,et al.  Quasimonotone systems and convergence to equilibrium in a population genetic model , 1983 .

[17]  Jiajun Zhang,et al.  Noise-induced switches in network systems of the genetic toggle switch , 2007, BMC Systems Biology.

[18]  Kim Sneppen,et al.  Quorum sensing regulation in Aeromonas hydrophila. , 2010, Journal of molecular biology.

[19]  H. Sauro Enzyme Kinetics for Systems Biology , 2012 .

[20]  M L Shuler,et al.  Ribosomal protein limitations in Escherichia coli under conditions of high translational activity , 1994, Biotechnology and bioengineering.

[21]  F. Wisniewski-Dyé,et al.  The regulatory locus cinRI in Rhizobium leguminosarum controls a network of quorum‐sensing loci , 2000, Molecular microbiology.

[22]  J. Collins,et al.  Programmable cells: interfacing natural and engineered gene networks. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[23]  J. Hasty,et al.  Noise-based switches and amplifiers for gene expression. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[24]  K. Postle,et al.  Nucleotide sequence of the repressor gene of the TN10 tetracycline resistance determinant. , 1984, Nucleic acids research.

[25]  C Cruz,et al.  Genetic studies of the lac repressor. XIV. Analysis of 4000 altered Escherichia coli lac repressors reveals essential and non-essential residues, as well as "spacers" which do not require a specific sequence. , 1994, Journal of molecular biology.

[26]  E. Greenberg,et al.  Acyl-homoserine lactone quorum sensing in gram-negative bacteria: a signaling mechanism involved in associations with higher organisms. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[27]  T. Baker,et al.  Overlapping recognition determinants within the ssrA degradation tag allow modulation of proteolysis , 2001, Proceedings of the National Academy of Sciences of the United States of America.

[28]  M. Elowitz,et al.  A synthetic oscillatory network of transcriptional regulators , 2000, Nature.

[29]  Martin Krupa,et al.  Bifurcations of relative equilibria , 1990 .

[30]  M. Hirsch Systems of Differential Equations that are Competitive or Cooperative II: Convergence Almost Everywhere , 1985 .

[31]  Robert T Sauer,et al.  Altered specificity of a AAA+ protease. , 2007, Molecular cell.

[32]  Raquel Tobes,et al.  The TetR Family of Transcriptional Repressors , 2005, Microbiology and Molecular Biology Reviews.

[33]  J E Bailey,et al.  Mathematical Modeling and Analysis in Biochemical Engineering: Past Accomplishments and Future Opportunities , 1998, Biotechnology progress.

[34]  Howard J. Li,et al.  Rapid and tunable post-translational coupling of genetic circuits , 2014, Nature.

[35]  L. Ju,et al.  Degradation and synthesis kinetics of quorum-sensing autoinducer in Pseudomonas aeruginosa cultivation. , 2005, Journal of biotechnology.

[36]  Lian-Hui Zhang,et al.  Quenching quorum-sensing-dependent bacterial infection by an N-acyl homoserine lactonase , 2001, Nature.

[37]  Eduardo D. Sontag,et al.  Mathematical Control Theory: Deterministic Finite Dimensional Systems , 1990 .

[38]  Lingchong You,et al.  Optimal tuning of bacterial sensing potential , 2009, Molecular systems biology.

[39]  Stephen Smale,et al.  THE DYNAMICAL SYSTEMS APPROACH TO DIFFERENTIAL EQUATIONS , 2007 .

[40]  Michael L. Shuler,et al.  A MATHEMATICAL MODEL FOR THE GROWTH OF A SINGLE BACTERIAL CELL * , 1979 .

[41]  G. Yarrington Molecular Cell Biology , 1987, The Yale Journal of Biology and Medicine.

[42]  M L Shuler,et al.  Incorporating genome-wide DNA sequence information into a dynamic whole-cell model of Escherichia coli: application to DNA replication. , 2008, IET systems biology.

[43]  Arkady B. Khodursky,et al.  Global analysis of mRNA decay and abundance in Escherichia coli at single-gene resolution using two-color fluorescent DNA microarrays , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[44]  Craig McAnulla,et al.  Quorum-sensing-regulated transcriptional initiation of plasmid transfer and replication genes in Rhizobium leguminosarum biovar viciae. , 2007, Microbiology.

[45]  Chris Cosner,et al.  Book Review: Monotone dynamical systems: An introduction to the theory of competitive and cooperative systems , 1996 .

[46]  Bonnie L Bassler,et al.  Observing bacteria through the lens of social evolution , 2008, Journal of biology.

[47]  F. Harary On the notion of balance of a signed graph. , 1953 .

[48]  M L Shuler,et al.  Prediction of effects of amino acid supplementation on growth of E. coli B/r , 1991, Biotechnology and bioengineering.

[49]  H. Bujard,et al.  Independent and tight regulation of transcriptional units in Escherichia coli via the LacR/O, the TetR/O and AraC/I1-I2 regulatory elements. , 1997, Nucleic acids research.

[50]  Frank Allgöwer,et al.  Network-level dynamics of diffusively coupled cells , 2012, 2012 IEEE 51st IEEE Conference on Decision and Control (CDC).

[51]  M. Shuler,et al.  A structured, segregated model for genetically modified Escherichia coli cells and its use for prediction of plasmid stability , 1990, Biotechnology and bioengineering.

[52]  E. Seibert,et al.  Fundamentals of enzyme kinetics. , 2014, Methods in molecular biology.

[53]  Y. Kuznetsov,et al.  Continuation techniques and interactive software for bifurcation analysis of ODEs and iterated maps: physics , 1993 .

[54]  Sandeep Krishna,et al.  The effect of LacI autoregulation on the performance of the lactose utilization system in Escherichia coli , 2013, Nucleic acids research.

[55]  R. Weiss,et al.  Ultrasensitivity and noise propagation in a synthetic transcriptional cascade. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[56]  Willy Govaerts,et al.  Numerical Continuation of Bifurcations of Limit Cycles in MATLAB , 2005, SIAM J. Sci. Comput..

[57]  M. Elowitz,et al.  Modeling a synthetic multicellular clock: repressilators coupled by quorum sensing. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[58]  Mauricio Barahona,et al.  Tuning the dials of Synthetic Biology , 2013, Microbiology.

[59]  Michael L. Shuler,et al.  Computer models of bacterial cells: from generalized coarsegrained to genome-specific modular models , 2006 .

[60]  Nancy Kopell,et al.  Synchrony in a Population of Hysteresis-Based Genetic Oscillators , 2004, SIAM J. Appl. Math..

[61]  D. Ohman,et al.  Synthesis of multiple exoproducts in Pseudomonas aeruginosa is under the control of RhlR-RhlI, another set of regulators in strain PAO1 with homology to the autoinducer-responsive LuxR-LuxI family , 1995, Journal of bacteriology.

[62]  M. Golubitsky,et al.  Singularities and groups in bifurcation theory , 1985 .

[63]  S. Smale On the differential equations of species in competition , 1976, Journal of mathematical biology.

[64]  M. Shuler,et al.  Structured model of genetic control via the lac promoter in Escherichia coli , 1994, Biotechnology and bioengineering.

[65]  Akiko Takaya,et al.  Negative Regulation of Quorum-Sensing Systems in Pseudomonas aeruginosa by ATP-Dependent Lon Protease , 2008, Journal of bacteriology.

[66]  Paul C. Bressloff,et al.  Stochastic Processes in Cell Biology , 2014, Interdisciplinary Applied Mathematics.

[67]  Walter Fast,et al.  The quorum-quenching lactonase from Bacillus thuringiensis is a metalloprotein. , 2005, Biochemistry.

[68]  G. Church,et al.  Global RNA half-life analysis in Escherichia coli reveals positional patterns of transcript degradation. , 2003, Genome research.

[69]  Marcus J. Tindall,et al.  A mathematical model of the sterol regulatory element binding protein 2 cholesterol biosynthesis pathway , 2014, Journal of theoretical biology.

[70]  E. Greenberg,et al.  Acyl homoserine-lactone quorum-sensing signal generation. , 1999, Proceedings of the National Academy of Sciences of the United States of America.

[71]  M L Shuler,et al.  Single-cell models: promise and limitations. , 1999, Journal of biotechnology.

[72]  E. Shnol,et al.  The synchronization of oscillators which interact via a medium , 1987 .

[73]  H. Bolouri Computational Modeling of Gene Regulatory Networks - A Primer , 2008 .

[74]  Eduardo D. Sontag,et al.  Monotone and near-monotone biochemical networks , 2007, Systems and Synthetic Biology.

[75]  David Angeli,et al.  Monotone control systems , 2003, IEEE Trans. Autom. Control..

[76]  J. Garey,et al.  The evolution of bacterial LuxI and LuxR quorum sensing regulators. , 2001, Microbiology.

[77]  Leonid V. Kalachev,et al.  The Boundary Function Method for Singular Perturbation Problems , 1995 .

[78]  F. Neidhardt,et al.  Physiology of the bacterial cell : a molecular approach , 1990 .

[79]  M L Shuler,et al.  A modular minimal cell model: purine and pyrimidine transport and metabolism. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[80]  J. Reiser,et al.  Autoinducer-mediated regulation of rhamnolipid biosurfactant synthesis in Pseudomonas aeruginosa. , 1995, Proceedings of the National Academy of Sciences of the United States of America.

[81]  M. Bennett,et al.  A fast, robust, and tunable synthetic gene oscillator , 2008, Nature.