Models for synthetic biology
© Kaznessis; licensee BioMed Central Ltd. 2007
Received: 18 October 2007
Accepted: 06 November 2007
Published: 06 November 2007
Synthetic biological engineering is emerging from biology as a distinct discipline based on quantification. The technologies propelling synthetic biology are not new, nor is the concept of designing novel biological molecules. What is new is the emphasis on system behavior.
The objective is the design and construction of new biological devices and systems to deliver useful applications. Numerous synthetic gene circuits have been created in the past decade, including bistable switches, oscillators, and logic gates, and possible applications abound, including biofuels, detectors for biochemical and chemical weapons, disease diagnosis, and gene therapies.
More than fifty years after the discovery of the molecular structure of DNA, molecular biology is mature enough for real quantification that is useful for biological engineering applications, similar to the revolution in modeling in chemistry in the 1950s. With the excitement that synthetic biology is generating, the engineering and biological science communities appear remarkably willing to cross disciplinary boundaries toward a common goal.
Synthetic biological engineering is emerging from biology as a distinct discipline based on quantification [1–5]. The objective is the design and construction of new biological devices and systems to deliver useful applications. Numerous synthetic gene circuits have been created in the past decade, including bistable switches, oscillators, and logic gates [[1–5], and references therein], and possible applications abound, ranging from biofuels, to detectors for biochemical and chemical weapons, to disease diagnosis, to gene therapies.
Certainly, the technologies propelling synthetic biology are not new, nor is the concept of designing novel biological molecules [6, 7]. What is perhaps new is the emphasis on system behavior, designing DNA sequences with synthetic phenotypes exhibiting prescribed dynamic responses.
Despite the initial successes of synthetic designs [1–5], the paradigm of biological sciences as descriptive disciplines may not rapidly assist in rationally engineering novel gene networks, despite the increasing volume of components that can be used in constructing synthetic networks. Genome projects identify the components of gene networks in biological organisms, gene after gene, and DNA microarray experiments discover the network connections. Yet, the static pictures of networks these experiments provide cannot adequately explain biomolecular phenomena or enable rational engineering of dynamic gene expression regulation. In other words, as an engineering discipline, synthetic biology cannot rely on endless trial and error methods driven by verbal description of biomolecular interaction networks.
The challenge facing the scientific and engineering communities is then to reduce the enormous volume and complexity of biological data into concise theoretical formulations with predictive ability, ultimately associating synthetic DNA sequences to dynamic phenotypes. The paradigm is not new either: In the 1940s and 1950s chemistry was a well matured discipline for pioneers like Neil Amundson, Byron Bird and Rutherford Aris to develop mathematical models that captured the enormous complexity of chemical processes in a way useful for chemical engineering applications [8–10]. Quantitative models of chemical processes led to the establishment of the chemical engineering discipline and the emergence of a strong chemical/petroleum industry. Although arguments can be made about the detrimental role of this industry on the environment, there can be no doubt of the overall positive effects on human life.
But what types of models are appropriate for synthetic biology? Because of the large number of participating species and the complexity of their interactions, only detailed modeling can allow the investigation of dynamic gene expression in a way fit for analysis and design. Designs can be detailed at the molecular level with dynamic models of all the biomolecular interactions involved in transcription, translation, regulation, transport and induction. We contrast this to a posteriori modeling of synthetic networks. For example in their seminal 2000 paper , Gardner and co-workers developed a very elegant model that captures and explains the observed dynamic behavior of the bistable switch and provides additional insight in the biological mechanism. This formalism may abide well with Occam's razor, but cannot guide the choice of specific DNA sequences and their regulatory relations to achieve a bistable switch. More specifically, it will be challenging to use reduced models to choose, for example, between lactose, arabinose or tetracycline operators, or any one of dozens of their mutant variants, for building a new, different bistable switch.
In engineering, descriptive models that are succinct and lucid are appreciated, but the ones used will be at the level of design degrees of freedom. For example, Bernoulli's equation can explain the aerodynamic lift of an airplane, but modern aircraft design is based on simulations that include all the components of flight in detail. Turning to synthetic biology, model-driven rational engineering of synthetic gene networks is possible at two levels:
First, the level of network topologies, where biomolecules control the concentration of other biomolecules, e.g. DNA binding proteins regulate the expression of specific genes by either activation or repression. By combining simple regulatory interactions, such as negative and positive feedback and feed forward loops, one may create more intricate networks that precisely control the production of protein molecules, such as bistable switches, oscillators, and filters. In the laboratory, these networks can be created using existing libraries of regulatory proteins and their corresponding operator sites. The now classical example is the aforementioned bistable switch Gardner and co-workers built : they connected two regulatory proteins repressing one another and this resulted in a bistable switch they could control. Another is the repressilator of Elowitz and Leibler : three regulatory proteins repressing one another in a sequential loop resulted in oscillating concentration profiles.
Secondly, the level of molecular components, which describes the kinetics and strengths of biomolecular interactions within the system. Indeed, the dynamical behavior of the system is a complex function of the kinetic interactions of the components. By altering the characteristics of the components, such as DNA-binding proteins and their corresponding DNA sites, one can modify the system's dynamical behavior without modifying the network topology. In the laboratory, the DNA sequences that yield the desired characteristics of each component can be engineered to achieve the desired protein-protein, protein-RNA, or protein-DNA binding constants and enzymatic activities. For example, Alon and co-workers  showed how simple mutations on the DNA sequence of the lactose operon can result in widely different phenotypic behaviors.
Ultimately, the large number of variants (interaction topologies and strengths) for these two types of design degrees of freedom requires sophisticated computational modeling, since the cost of experimentally changing these components and the kinetics of their interactions can quickly become prohibitive. Computer simulations enable exhaustive searches of different network connectivities and molecular thermodynamic/kinetic parameters, greatly advancing the development of design principles that seek to simplify the complicated behavior of the network into a brief, usable framework.
with k1 and k-1 the association and dissociation kinetic constants, respectively. If we considered the cell as a well-stirred reactor we could calculate the behavior of the network using a set of ordinary differential equations, which determine concentration changes as prescribed by kinetic laws. However, the underlying assumption of such continuous-deterministic models, that the number of molecules approaches the thermodynamic limit (i.e. that the volume of the system is infinite), can be invalid for biological systems, since for some components (DNA for example) there are only a few copies available.
More than fifty years after the discovery of the molecular structure of DNA, molecular biology is mature enough for quantification useful for biological engineering applications, similar to chemistry in the 1950s. With the excitement synthetic biology is generating, the engineering and biological science communities appear remarkably willing to cross disciplinary boundaries toward this common goal.
- Andrianantoandro E, Basu S, Karig DK, Weiss R: Synthetic biology: new engineering rules for an emerging discipline. Mol Syst Biol. 2006, 2: 2006.0028- 10.1038/msb4100073PubMed CentralPubMedView ArticleGoogle Scholar
- Church GM: From systems biology to synthetic biology. Mol Syst Biol. 2005, 1: 2005.0032-PubMed CentralPubMedGoogle Scholar
- Keasling J: The promise of synthetic biology. Bridge Natl Acad Eng. 2005, 35: 18-21.Google Scholar
- Alon U: Biological networks: the tinkerer as an engineer. Science. 2003, 301: 1866-1867. 10.1126/science.1089072PubMedView ArticleGoogle Scholar
- Kobayashi H, Kaern M, Araki M, Chung K, Gardner TS, Cantor CR, Collins JJ: Programmable cells: interfacing natural and engineered gene networks. Proc Natl Acad Sci USA. 2004, 101: 8414-8419. 10.1073/pnas.0402940101PubMed CentralPubMedView ArticleGoogle Scholar
- Cohen SN, Chang ACY, Boyer H, Helling RB: Construction of biologically functional bacterial plasmids in vitro. Proc Natl Acad Sci USA. 1973, 70: 3240-3244. 10.1073/pnas.70.11.3240PubMed CentralPubMedView ArticleGoogle Scholar
- Eisen H, Brachet P, Pereira da Silva L, Jacob F: Regulation of repressor expression in λ. Proc Natl Acad Sci USA. 1970, 66: 855-862. 10.1073/pnas.66.3.855PubMed CentralPubMedView ArticleGoogle Scholar
- Acrivos A, Amundson NR: Applications of Matrix Mathematics to Chemical Engineering Problems. Ind Eng Chem. 1955, 47 (8): 1533-1541. 10.1021/ie50548a027.View ArticleGoogle Scholar
- Vectors, Tensors and the Basic Equations of Fluid Mechanics. 1962, Rutherford Aris, Prentice-Hall
- Fredrickson A, Bird RB: Non-Newtonian Flow in Annuli. Ind Eng Chem. 1958, 50 (3): 347-352. 10.1021/ie50579a035.View ArticleGoogle Scholar
- Gardner TS, Cantor CR, Collins JJ: Construction of a genetic toggle switch in Escherichia coli. Nature. 2000, 403: 339-342. 10.1038/35002131PubMedView ArticleGoogle Scholar
- Elowitz MB, Leibler S: A synthetic oscillatory network of transcriptional regulators. Nature. 2000, 403: 335-338. 10.1038/35002125PubMedView ArticleGoogle Scholar
- Mayo AE, Setty Y, Shavit S, Zaslaver A, Alon U: Plasticity of the cis-regulatory input function of a gene. PLoS Biol. 2006, 4: e45- 10.1371/journal.pbio.0040045PubMed CentralPubMedView ArticleGoogle Scholar
- McQuarrie DA: Stochastic Approach to Chemical Kinetics. Journal of Applied Probability. 1967, 4 (3): 413-478. 10.2307/3212214.View ArticleGoogle Scholar
- Grabert H, Hanggi P, Oppenheim I: Fluctuations in Reversible Chemical Reactions. Physica. 1983, l17A: 300-316.View ArticleGoogle Scholar
- Gillespie DT: A general method for numerically simulating the stochastic time evolution of coupled chemical reactions. Journal of Computational Physics. 1976, 22: 403-434. 10.1016/0021-9991(76)90041-3.View ArticleGoogle Scholar
- Gibson MA, Bruck J: Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels. J Phys Chem. 2000, 104: 1876-1889.View ArticleGoogle Scholar
- Puchalka J, Kierzek AM: Bridging the gap between stochastic and deterministic regimes in the kinetic simulations of the biochemical reaction networks. Biophys J. 2004, 1357-1372.Google Scholar
- Rathinam M, Petzold LR, Cao Y: Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method. J Chem Phys. 2005, 119 (24): 12784-12794. 10.1063/1.1627296.View ArticleGoogle Scholar
- Chatterjee A, Mayawala K, Edwards JS, Vlachos DG: Time accelerated Monte Carlo simulations of biological networks using the binomial t-leap method. Bioinformatics. 2005, 21 (9): 2136-2137. 10.1093/bioinformatics/bti308PubMedView ArticleGoogle Scholar
- Haseltine EL, Rawlings JB: Approximate simulation of coupled fast and slow reactions for stochastic chemical kinetics. J Chem Phys. 2002, 117 (15): 6959-6969. 10.1063/1.1505860.View ArticleGoogle Scholar
- Salis H, Kaznessis Y: Accurate Hybrid Stochastic Simulation of a System of Coupled Chemical or Biochemical Reactions. J Chem Phys. 2005, 122: 054103-10.1063/1.1835951.View ArticleGoogle Scholar
- Salis H, Kaznessis YN: An equation-free probabilistic steady-state approximation: dynamic application to the stochastic simulation of biochemical reaction networks. J Chem Phys. 2005, 123 (21): 214106- 10.1063/1.2131050PubMedView ArticleGoogle Scholar
- Kaznessis Y: Multi-Scale Models for Gene Network Engineering. Chemical Engineering Science. 2006, 61 (3): 940-953. 10.1016/j.ces.2005.06.033.View ArticleGoogle Scholar
- Salis H, Kaznessis Y: Computer-aided design of modular protein devices: Boolean AND gene activation. Phys Biol. 2006, 3 (4): 295-310. 10.1088/1478-3975/3/4/007PubMedView ArticleGoogle Scholar
- Sotiropoulos V, Kaznessis Y: Synthetic tetracycline-inducible regulatory networks: computer-aided design of dynamic phenotypes. BMC Systems Biology. 2007, 1: 7- 10.1186/1752-0509-1-7PubMed CentralPubMedView ArticleGoogle Scholar
- Kiehl TR, Mattheyses RM, Simmons MK: Hybrid simulation of cellular behavior. Bioinformatics. 2004, 316-322.Google Scholar
- De Jong H: Modeling and simulation of genetic regulatory systems: A literature review. J Comput Biol. 2002, 9 (1): 67-103. 10.1089/10665270252833208PubMedView ArticleGoogle Scholar
- Salis H, Sotiropoulos V, Kaznessis YN: Multiscale Hy3S: hybrid stochastic simulation for supercomputers. BMC Bioinformatics. 2006, 24 (7): 93-10.1186/1471-2105-7-93.View ArticleGoogle Scholar