A computational method for the investigation of multistable systems and its application to genetic switches
- Miriam Leon^{1},
- Mae L. Woods^{1},
- Alex J. H. Fedorec^{1} and
- Chris P. Barnes^{1, 2}Email authorView ORCID ID profile
https://doi.org/10.1186/s12918-016-0375-z
© The Author(s) 2016
Received: 3 August 2016
Accepted: 13 November 2016
Published: 7 December 2016
Abstract
Background
Genetic switches exhibit multistability, form the basis of epigenetic memory, and are found in natural decision making systems, such as cell fate determination in developmental pathways. Synthetic genetic switches can be used for recording the presence of different environmental signals, for changing phenotype using synthetic inputs and as building blocks for higher-level sequential logic circuits. Understanding how multistable switches can be constructed and how they function within larger biological systems is therefore key to synthetic biology.
Results
Here we present a new computational tool, called StabilityFinder, that takes advantage of sequential Monte Carlo methods to identify regions of parameter space capable of producing multistable behaviour, while handling uncertainty in biochemical rate constants and initial conditions. The algorithm works by clustering trajectories in phase space, and iteratively minimizing a distance metric. Here we examine a collection of models of genetic switches, ranging from the deterministic Gardner toggle switch to stochastic models containing different positive feedback connections. We uncover the design principles behind making bistable, tristable and quadristable switches, and find that rate of gene expression is a key parameter. We demonstrate the ability of the framework to examine more complex systems and examine the design principles of a three gene switch. Our framework allows us to relax the assumptions that are often used in genetic switch models and we show that more complex abstractions are still capable of multistable behaviour.
Conclusions
Our results suggest many ways in which genetic switches can be enhanced and offer designs for the construction of novel switches. Our analysis also highlights subtle changes in correlation of experimentally tunable parameters that can lead to bifurcations in deterministic and stochastic systems. Overall we demonstrate that StabilityFinder will be a valuable tool in the future design and construction of novel gene networks.
Keywords
Genetic switches Sequential Monte Carlo Design of genetic circuitsBackground
Synthetic biology has seen the development of many simple gene circuits such as switches [1–6], oscillators [7–9] and pulse generators [10]. Larger systems have been constructed [11], but the leap from building low-level circuits to assembling them into complex networks is still a major challenge [12, 13]. Efforts to do so are plagued by circuit crosstalk, retroactivity, chassis loading effects, and cellular noise, which can render synthetic networks non-functional in vivo [14, 15]. Although standardization and better part design can partially lower this barrier [12, 16–19], design processes that enable the informed selection of appropriate parts are crucial [11, 20, 21].
One of the foundational constructs in synthetic biology is the genetic toggle switch. The toggle switch consists of a set of transcription factors that mutually repress each other [1, 22–24]. Genetic switches play a major role in binary cell fate decisions such as stem cell differentiation, as they are capable of exhibiting bistable behaviour, which gives rise to the existence of two distinct phenotypic states. This allows populations of cells to maintain a heterogeneous response to environmental cues and can increase fitness by bet-hedging [25]. Switches are powerful building blocks; they underlie electronics and logic systems, and have great potential in synthetic biology. The genetic toggle switch has been used for a number of applications including the construction of a synthetic genetic clock [22], the regulation of mammalian gene expression [2, 5], the development of a predictable genetic timer [26], and the formation of biofilms in response to engineered stimuli [27].
The stability of the toggle switch has been investigated extensively in the literature, but the conclusions drawn vary according to model abstraction. Numerous studies have concluded that cooperativity is a necessary condition for bistability to arise [1, 28–31]. However, Lipshtat et al. found that stochastic effects can give rise to bistability even without cooperativity [32]. In another study, Ma et al. found that stochastic fluctuations can stabilize the unstable steady state in the deterministic system, giving rise to tristability [33]. In addition, Biancalani et al. identified multiplicative noise as the source of bistability in the stochastic case [34]. As is clear from the above, there is yet to exist a consensus on the stability a switch is capable of, and the most appropriate method of modelling it. Most of these studies assumed the quasi-steady state approximation (QSSA) [35], which cannot always be assumed to hold in vivo [36].
In terms of system design, extensions of the basic toggle switch motif, including additional positive feedback mechanisms, have been investigated [37, 38], and optimization methods have been used to identify topologies and parameter values for bistable and tristable genetic switches [39–42]. For stochastic switch design, control theoretic approaches [43], and simulation-based frameworks [44], have been developed. However, none of these existing approaches can be be applied to reasonably sized models, under the assumption of deterministic and stochastic dynamics, and identify regions of parameter space for which switching occurs, which we argue is critical in designing systems under considerable uncertainty.
Methods
StabilityFinder is based on a statistical inference method that combines approximate Bayesian computation (ABC) with sequential Monte Carlo [46]. This simulation-based method uses an iterative process to arrive at a distribution of parameter values that can give rise to observed data or a desired system behaviour [44]. ABC methods are used for inferring the posterior distribution in cases where the likelihood is intractable or is too computationally expensive to evaluate. Instead of computing the likelihood, ABC methods simulate the data and then compare the simulated and observed data through a distance function [46]. Given the prior distribution π(θ) we can approximate the posterior distribution, π (θ|x)∝f(x|θ)π (θ), where f(x|θ) is the likelihood of a parameter, θ, given the data, x. There are a number of different variations of the ABC algorithm depending on how the approximate posterior distribution is sampled.
The simplest ABC algorithm is the ABC rejection sampler [47]. In this method, parameters are sampled from the prior and data simulated through the data generating model. For each simulated data set, x ^{∗}, a distance from that of the desired behaviour is calculated, ρ(x ^{∗},y), and if greater than a threshold, ε, the sample is rejected, otherwise it is accepted.
The main disadvantage of this method is that if the prior distribution is very different from the posterior, the acceptance rate is very low [46]. An alternative method is ABC Markov Chain Monte Carlo (MCMC) [48]. The disadvantage of this method is that if it gets stuck in an area of low probability it can be very slow to converge [49]. The method used here is based on sequential Monte Carlo, which avoids both issues faced by the rejection and MCMC methods. It propagates the prior through a series of intermediate distributions in order to arrive at an approximation of the posterior. The tolerance, ε, for the distance of the simulated data to the desired data is made smaller at each iteration. When ε is sufficiently small, the result will approximate the posterior distribution [46].
To investigate the multistable behaviour of systems, a number of extensions to existing approaches are required. For a given set of parameter values, sample points are taken across initial conditions using latin hypercube sampling [50], and the ensemble system simulated in time until steady state. The distance function in ABC is replaced by a distance on the desired stability of the simulated model. To do this we cluster the steady state coordinates using K-means clustering [51] and use the Gap statistic to determine the number of clusters [52]. At each iteration, the number of steady states is determined by the number of clusters in phase space. A particle is accepted only if the number of clusters present is within an acceptable distance from the threshold ε. The algorithm is summarised below.
Results and discussion
The Gardner switch under deterministic and stochastic dynamics
where u is the concentration of repressor 1, v the concentration of repressor 2, α _{1} and α _{2} denote the effective rates of synthesis of repressors 1 and 2 respectively, and β and γ are the cooperativity of repression of promoter 1 and of repressor 2 respectively. Gardner et al. studied the deterministic case and concluded that there are two conditions for bistability for this model; that α _{1} and α _{2} are balanced and that β, γ >1 [1]. In order to test StabilityFinder we used it to find the posterior distribution for which this model exhibits bistable behaviour. We therefore set the desired behaviour to two (stable) steady states, and using a wide range of values as priors as shown in the Additional file 1, we used StabilityFinder to find the parameter values necessary for bistability to occur. The posterior distribution calculated by StabilityFinder for the Gardner deterministic case is shown in Fig. 2 c. These results agree with the results reported by Gardner et al. [1]. For this switch to be bistable α _{1} and α _{2} must be balanced while β and γ must both be >1, as can be seen in the marginal distributions of β and γ in Fig. 2 c.
We next applied StabilityFinder to the case of the Gardner switch under stochastic dynamics using the same priors as the deterministic case, and again searched the parameter space for bistable behaviour. The posterior is shown in Fig. 2 d. We can see that the conditions on the parameters required for bistability in the deterministic case generally still stand in the stochastic case. There appears to be slightly looser requirements on the parameters of the stochastic model (wider marginal distributions), which is expected due to the nature of clustering deterministic steady states versus stochastic steady states. The Gap statistic is used in the case of the stochastic case, as it is capable of dealing with noisier data whereas a simpler and faster algorithm is used for clustering the deterministic solutions (see Additional file 1). These results demonstrate that StabilityFinder can be used to find the parameter values that produce a desired stability and allow us to confidently apply the methodology to more complex models.
Repressor degradation rates are key for achieving bistablity
and g _{ I } represents the production rate, k _{ I } the degradation rate, n _{ I } the Hill coefficient, x _{ I } the Hill threshold concentration and λ _{ I } the fold change of the transcription rates, and I∈{x y,y x,x x,y y} (see Additional file 1 for the details of all models used).
The addition of symmetric and asymmetric positive autoregulation
It is known that the addition of positive autoregulation to the classical toggle switch can induce tristability [37, 38, 58]. Here we investigate the interplay of positive autoregulation on the values of the other parameters in the model. We extended the analysis presented in Lu et al. by including the switch with single positive autoregulation (model LU-SP), where a single positive autoregulatory feedback is present on one of the genes. This system topology has also been constructed previously [23, 59]. The advantage of using StabilityFinder over traditional bifurcation analysis is that the full parameter space is explored rather than solving the system for a single set of parameters. This allows us to deduce model properties that could not otherwise be identified. Robustness to parameter fluctuations can be explored, as well as parameter correlations and restrictions on the values they can take while still producing the desired behaviour.
We find that the switch with single positive autoregulation is capable of bistable behaviour as seen in Fig. 3 b, but this is only possible when the strength of the promoter under positive autoregulation, g _{ x }, is small. There appear to be no such constraints on the strength of the original, unmodified, promoter, g _{ y }. We also find that unlike the LU-CS and LU-DP models, the steady states of the bistable LU-SP are not symmetric. The levels of Y are around zero and always lower than the levels of X. The levels of both are lower than those found in the LU-CS and LU-DP models.
Upon examination of the LU-DP model, we also find that tristability in the switch is relatively robust, as this phenotype is found across a large range of parameter values, with no parameters strongly constrained (see Additional file 1) but the two parameters for gene expression, g _{ x } and g _{ y } tend to be small compared to the priors. Two types of tristable behaviour are identified, one where the third steady state is at (0,0) and and one where the third steady state has non-zero values. This result agrees with previous work where it was found that a switch can exhibit two kinds of tristability, one in which the third steady state is high (III _{ H }) and one in which it is low (III _{ L }) [37].
Design principles for a switch capable of two, three and four steady states
Using StabilityFinder we obtained posterior distributions for the bistable, tristable and quadristable phenotypes (Fig. 4). Upon examination of the these distributions, we observe that a subset of the parameter values are different for the three behaviours, although the differences are surprisingly subtle. We find differences in the univariate distribution of the parameters for gene expression, g _{ x }, as highlighted in Fig. 4 c, box 1. This parameter must be small for four steady states to occur but there are no such restraints for a bistable or a tristable switch. Furthermore, parameter x _{ xx } (the dissociation constant for autoregulation) must be small for tristable and quadristable behaviour to be achieved, but there are no such restraints for a bistable switch, as can be seen in Fig. 4 c, box 2.
Additionally, we find a difference in the bivariate distributions in the posterior. Most notably, we find that parameters x _{ xx } and g _{ x } are tightly constrained in the tristable and quadristable cases, where both parameters are required to be small, but less so in the bistable case (Fig. 4 c, box 3). Another notable difference is between parameters x _{ xx } and n _{ xx } shown in Fig. 4 c, box 5, where they are constrained in the tristable and quadristable cases but not the bistable case. Interestingly, we also find parameter correlations conserved between the three behaviours, as seen in Fig. 4 c, box 4, where parameters l _{ xx } and g _{ x }, (positive autoregulation and gene expression) are negatively correlated in both cases. This highlights the importance of treating unknown parameters as distributions rather than fixed values when studying complex system behaviour. These ensemble based methods are capable of uncovering not only the ranges and values required for certain behaviour, but also the correlations between parameters, which would be missed by optimisation methods.
Bistability and tristability using more realistic abstractions
In order to study the switch system in the most realistic way, we avoid using the quasi-steady state approximation (QSSA) that is often used in modelling the toggle switch. Under the assumption of mass action kinetics, the two-equation system becomes a system of 14 equations and 9 parameters in the classical switch case. These are shown in the Additional file 1. In the cellular environment stochastic effects can be non-negligible and should therefore be taken into account when trying to elucidate the behaviours that a system is capable of. Therefore, we model these switches using stochastic dynamics.
Design principles of bistable and tristable switches
CS-MA | DP-MA | |||
---|---|---|---|---|
Bistable | Tristable | Bistable | Tristable | |
Dimerisation | High | Low | High | Low |
Protein degradation | - | - | - | Low |
Dimer degradation | Low | - | Low | - |
Achieving higher multistability
Conclusions
We have developed an algorithm that can identify the parameter regions necessary for a model to achieve a given number of stable steady states. The novelty in our framework over existing methodology is that complex models can be analyzed assuming both deterministic and stochastic dynamics. We have shown that the algorithm can be used to infer the parameter ranges that give rise to specific behaviour in various models. We uncovered the design principles that make a bistable, a tristable and a quadristable switch. We also found that a three-node switch is capable of hexastability. Importantly, we removed assumptions made to simplify the switch models and showed that they are still capable of bistable and tristable behaviour.
Although we only examined models containing combined transcription and translation, our approach could be applied to any models of switching behaviour, including more detailed kinetic models and more complex multistable switches that exist in natural biological systems, such as developmental pathways. We also limited our framework to the objective behaviour of a given number of stable steady states. However, this approach is extremely flexible, and could be extended to find systems with a given switching rate, or systems robust to a particular set of perturbations, both of which could be of great importance for building more complex genetic circuits.
One limitation of our approach is that we cannot rule out a specific behaviour; it is always possible that some part of parameter space remains unexplored, or because the search space must be predefined, interesting regions are not included in the search. In the Bayesian sense, this predefined space is the prior distribution for the parameters that give rise to the stability under investigation. In principle, once our knowledge of these biochemical rate constants grows, we can incorporate these data into the prior regions for exploration. Another limitation is that of scalability. Our framework can currently be applied to small and medium size gene networks since the computational time is exponential in size, whereas optimization methods are more scalable [39–41]. This is a manifestation of a general tradeoff between finding an optimal value and exploring a parameter space. However, we argue that for current and relevant problems in synthetic biology, this computational burden is acceptable.
Approaches based on parameter space exploration are indispensible tools for providing understanding of general system properties and guiding more detailed experimental and theoretical studies. They will also be key for the design and construction of synthetic gene networks. By selecting standardized parts accordingly, such as promoters, RBS sequences and other untranslated regions [18, 63–65], in vivo systems can be matched to parameter regions with a high probability of function.
More generally our results highlight that changing the level of abstraction, in addition to the modification of the feedback structure and parameter values, can significantly alter the qualitative behaviour of a system model. These results advocate the need for a programme of experimental work, combined with systems modelling, to understand the rules of thumb for abstraction in model based design of synthetic biological systems.
Declarations
Acknowledgements
All the authors would like to acknowledge that the work presented here made use of the Emerald High Performance Computing facility made available by the Centre for Innovation. The Centre is formed by the universities of Oxford, Southampton, Bristol, and University College London in partnership with the STFC Rutherford-Appleton Laboratory.
Funding
CPB and MLW acknowledge funding from the Wellcome Trust through a Research Career Development Fellowship (097319/Z/11/Z). ML and AJHF acknowledge funding through the UCL Impact Award scheme and the UCL CoMPLEX doctoral training centre respectively.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its Additional file 1. The StabilityFinder package and the models run in this analysis can be found in the public GitHub repository ucl-cssb (https://github.com/ucl-cssb).
Authors’ contributions
The project was conceived by CPB. ML developed the StabilityFinder code, developed the switch models, ran the Monte Carlo analyses and interpreted the results. MLW, AJHF and CPB performed additional analysis of the results. ML and CPB wrote the manuscript. All authors read and approved the final manuscript.
Competing interests
The authors declare that they have no competing interests.
Consent for publication
Not applicable
Ethics approval and consent to participate
Not applicable
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
References
- Gardner TS, Cantor CR, Collins JJ. Construction of a genetic toggle switch in Escherichia coli. Nature. 2000; 403(6767):339–42.View ArticlePubMedGoogle Scholar
- Kramer BP, Viretta AU, Daoud-El-Baba M, Aubel D, Weber W, Fussenegger M. An engineered epigenetic transgene switch in mammalian cells. Nat Biotechnol. 2004; 22(7):867–70.View ArticlePubMedGoogle Scholar
- Isaacs FJ, Hasty J, Cantor CR, Collins JJ. Proc Nat Acad Sci USA. 2003; 100(13):7714–9.Google Scholar
- Ham TS, Lee SK, Keasling JD, Arkin AP. Design and Construction of a Double Inversion Recombination Switch for Heritable Sequential Genetic Memory. PLoS ONE. 2008; 3(7):2815.View ArticleGoogle Scholar
- Deans TL, Cantor CR, Collins JJ. A Tunable Genetic Switch Based on RNAi and Repressor Proteins for Regulating Gene Expression in Mammalian Cells. Cell. 2007; 130(2):363–72.View ArticlePubMedGoogle Scholar
- Friedland AE, Lu TK, Wang X, Shi D, Church G, Collins JJ. Synthetic gene networks that count. Science. 2009; 324(5931):1199–202.View ArticlePubMedPubMed CentralGoogle Scholar
- Stricker J, Cookson S, Bennett MR, Mather WH, Tsimring LS, Hasty J. A fast, robust and tunable synthetic gene oscillator. Nature. 2008; 456(7221):516–9.View ArticlePubMedGoogle Scholar
- Fung E, Wong WW, Suen JK, Bulter T, Lee SG, Liao JC. A synthetic gene–metabolic oscillator. Nature. 2005; 435(7038):118–22.View ArticlePubMedGoogle Scholar
- Tigges M, Marquez-Lago TT, Stelling J, Fussenegger M. A tunable synthetic mammalian oscillator. Nature. 2009; 457(7227):309–12.View ArticlePubMedGoogle Scholar
- Basu S, Mehreja R, Thiberge S, Chen MT, Weiss R. Proc Nat Acad Sci USA. 2004; 101(17):6355–360.Google Scholar
- Nielsen AA, Der BS, Shin J, Vaidyanathan P, Paralanov V, Strychalski EA, Ross D, Densmore D, Voigt CA. Genetic circuit design automation. Science. 2016; 352(6281):aac7341. doi:http://dx.doi.org/10.1126/science.aac7341.
- Lu TK, Khalil AS, Collins JJ. Next-generation synthetic gene networks. Nat Biotechnol. 2009; 27(12):1139–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Cardinale S, Arkin AP. Contextualizing context for synthetic biology Ű identifying causes of failure of synthetic biological systems. Biotechnol J.2012.doi:http://dx.doi.org/10.1002/biot.201200085.
- Del Vecchio D. Modularity, context-dependence, and insulation in engineered biological circuits. Trends Biotechnol. 2015; 33(2):111–9.View ArticlePubMedGoogle Scholar
- Ceroni F, Algar R, Stan GB, Ellis T. Quantifying cellular capacity identifies gene expression designs with reduced burden. Nat Methods. 2015; 12(5):415–8. doi:http://dx.doi.org/10.1038/nmeth.3339.
- Shetty RP, Endy D, Knight TF. Engineering BioBrick vectors from BioBrick parts. J Biol Eng. 2008; 2:5–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Galdzicki M, Rodriguez C, Chandran D, Sauro HM, Gennari JH. Standard biological parts knowledgebase. PLoS ONE. 2011; 6(2):17005.View ArticleGoogle Scholar
- Mutalik VK, Guimaraes JC, Cambray G, Lam C, Christoffersen MJJ, Mai Q-AA, Tran AB, Paull M, Keasling JD, Arkin AP, Endy D. Precise and reliable gene expression via standard transcription and translation initiation elements. Nat Methods. 2013; 10(4):354–60. doi:http://dx.doi.org/10.1038/nmeth.2404.
- Nielsen AA, Segall-Shapiro TH, Voigt CA. Advances in genetic circuit design: novel biochemistries, deep part mining, and precision gene expression. Curr Opinion Chem Biol. 2013; 17(6):878–92. doi:http://dx.doi.org/10.1016/j.cbpa.2013.10.003.
- Beal J, Weiss R, Densmore D, Adler A, Appleton E, Babb J, Bhatia S, Davidsohn N, Haddock T, Loyall J, Schantz R, Vasilev V, Yaman F. An end-to-end workflow for engineering of biological networks from high-level specifications. ACS Synthetic Biol. 2012; 1(8):317–31. doi:http://dx.doi.org/10.1021/sb300030d.
- Yaman F, Bhatia S, Adler A, Densmore D, Beal J. Automated selection of synthetic biology parts for genetic regulatory networks. ACS Synthetic Biol. 2012; 1(8):332–44.View ArticleGoogle Scholar
- Atkinson MR, Savageau MA, Myers JT, Ninfa AJ. Development of genetic circuitry exhibiting toggle switch or oscillatory behavior in Escherichia coli. Cell. 2003; 113(5):597–607.View ArticlePubMedGoogle Scholar
- Lou C, Liu X, Ni M, Huang Y, Huang Q, Huang L, Jiang L, Lu D, Wang M, Liu C, Chen D, Chen C, Chen X, Yang L, Ma H, Chen J, Ouyang Q. Synthesizing a novel genetic sequential logic circuit: a push-on push-off switch. Mol Syst Biol. 2010; 6. doi:http://dx.doi.org/10.1038/msb.2010.2.
- Litcofsky KD, Afeyan RB, Krom RJ, Khalil AS, Collins JJ. Iterative plug-and-play methodology for constructing and modifying synthetic gene networks. Nat Methods. 2012; 9(11):1077–80.View ArticlePubMedPubMed CentralGoogle Scholar
- Veening JW, Smits WK, Kuipers OP. Bistability, epigenetics, and bet-hedging in bacteria. Microbiology. 2008; 62:193–210.View ArticleGoogle Scholar
- Ellis T, Wang X, Collins JJ. Diversity-based, model-guided construction of synthetic gene networks with predicted functions. Nat Biotechnol. 2009; 27(5):465–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Kobayashi H, Kaern M, Araki M, Chung K, Gardner TS, Cantor CR, Collins JJ. Programmable cells: interfacing natural and engineered gene networks. Proc Nat Acad Sci USA. 2004; 101(22):8414–419.View ArticlePubMedPubMed CentralGoogle Scholar
- Cherry JL, Adler FR. How to make a biological switch,. J Theor Biol. 2000; 203(2):117–33.View ArticlePubMedGoogle Scholar
- Warren PB, ten Wolde PR. Enhancement of the Stability of Genetic Switches by Overlapping Upstream Regulatory Domains. Phys Rev Lett. 2004; 92(12):128101.View ArticlePubMedGoogle Scholar
- Walczak AM, Onuchic JN, Wolynes PG. Absolute rate theories of epigenetic stability. Proc Nat Acad Sci USA. 2005; 102(52):18926–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Warren PB, ten Wolde PR. Chemical models of genetic toggle switches. J Phys Chem B. 2005; 109(14):6812–23.View ArticlePubMedGoogle Scholar
- Lipshtat A, Loinger A, Balaban NQ, Biham O. Genetic toggle switch without cooperative binding. Phys Rev Lett. 2006; 96(18):188101.View ArticlePubMedGoogle Scholar
- Ma R, Wang J, Hou Z, Liu H. Small-number effects: a third stable state in a genetic bistable toggle switch. Phys Rev Lett. 2012; 109(24):248107.View ArticlePubMedGoogle Scholar
- Biancalani T, Assaf M. Genetic Toggle Switch in the Absence of Cooperative Binding: Exact Results. Phys Rev Lett. 2015; 115:208101.View ArticlePubMedGoogle Scholar
- Loinger A, Lipshtat A, Balaban NQ, Biham O. Stochastic simulations of genetic switch systems. Phys Rev E Stat Nonlin Soft Matter Phys. 2007; 75(2 Pt 1):021904.View ArticlePubMedGoogle Scholar
- Pedersen MG, Bersani AM, Bersani E. Quasi steady-state approximations in complex intracellular signal transduction networks – a word of caution. J Math Chem. 2007; 43(4):1318–44.View ArticleGoogle Scholar
- Guantes R, Poyatos JF. Multistable decision switches for flexible control of epigenetic differentiation. PLoS Comput Biol. 2008; 4(11):1000235.View ArticleGoogle Scholar
- Lu M, Onuchic J, Ben-Jacob E. Construction of an Effective Landscape for Multistate Genetic Switches. Phys Rev Lett. 2014; 113(7):078102.View ArticlePubMedGoogle Scholar
- Dasika MS, Maranas CD. Optcircuit: an optimization based method for computational design of genetic circuits. BMC Syst Biol. 2008; 2(1):1.View ArticleGoogle Scholar
- Otero-Muras I, Banga JR. Multicriteria global optimization for biocircuit design. BMC Syst Biol. 2014; 8:113.View ArticlePubMedPubMed CentralGoogle Scholar
- Otero-Muras I, Banga JR. Exploring design principles of gene regulatory networks via pareto optimality. IFAC-PapersOnLine. 2016; 49(7):809–14. doi:http://dx.doi.org/10.1016/j.ifacol.2016.07.289. 11th {IFAC} Symposium on Dynamics and Control of Process SystemsIncluding Biosystems DYCOPS-CAB 2016Trondheim, Norway, 6–8 June 2016.
- Rodrigo G, Carrera J, Jaramillo A. Computational design of synthetic regulatory networks from a genetic library to characterize the designability of dynamical behaviors. Nucleic Acids Res. 2011; 39(20):138–8.View ArticleGoogle Scholar
- Baetica AA, Yuan Y, Gonçalves JM, Murray RM. A stochastic framework for the design of transient and steady state behavior of biochemical reaction networks. In: 54th IEEE Conference on Decision and Control, CDC 2015, Osaka, Japan, December 15–18, 2015: 2015. p. 3199–205. doi:http://dx.doi.org/10.1109/CDC.2015.7402699. http://dx.doi.org/10.1109/CDC.2015.7402699.
- Barnes CP, Silk D, Sheng X, Stumpf MPH. Proc Nat Acad Sci USA. 2011; 108(37):15190–15195.Google Scholar
- Del Moral P, Doucet A, Jasra A. Sequential monte carlo samplers. J R Stat Soc Series B (Statistical Methodology). 2006; 68(3):411–36.View ArticleGoogle Scholar
- Toni T, Welch D, Strelkowa N, Ipsen A, Stumpf MPH. Approximate Bayesian computation scheme for parameter inference and model selection in dynamical systems. J R Soc Interface / R Soc. 2009; 6(31):187–202.View ArticleGoogle Scholar
- Pritchard JK, Seielstad MT, Perez-Lezaun A, Feldman MW. Population growth of human Y chromosomes: A study of Y chromosome microsatellites. Mol Biol Evol. 1999; 16(12):1791–8.View ArticlePubMedGoogle Scholar
- Marjoram P, Molitor J, Plagnol V, Tavare S. Markov chain Monte Carlo without likelihoods. Proc Nat Acad Sci USA. 2003; 100(26):15324–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Sisson SA, Fan Y, Tanaka MM. Sequential Monte Carlo without likelihoods. Proc Nat Acad Sci USA. 2007; 104(6):1760–1765.View ArticlePubMedPubMed CentralGoogle Scholar
- McKay MD, Beckman RJ, Conover WJ. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics. 2000; 42(1):55–61.View ArticleGoogle Scholar
- Lloyd SP. Least squares quantization in PCM. IEEE Trans Inf Theory. 1982; 28(2):129–137.View ArticleGoogle Scholar
- Tibshirani R, Walther G, Hastie T. Estimating the number of clusters in a data set via the gap statistic. J R Stat Soc B. 2001; 63:411–23.View ArticleGoogle Scholar
- Hucka M, Finney A, Sauro HM, Bolouri H, Doyle JC, Kitano H, Arkin AP, Bornstein BJ, Bray D, Cornish-Bowden A, Cuellar AA, Dronov S, Gilles ED, Ginkel M, Gor V, Goryanin II, Hedley WJ, Hodgman TC, Hofmeyr JH, Hunter PJ, Juty NS, Kasberger JL, Kremling A, Kummer U, Le Novère N, Loew LM, Lucio D, Mendes P, Minch E, Mjolsness ED, Nakayama Y, Nelson MR, Nielsen PF, Sakurada T, Schaff JC, Shapiro BE, Shimizu TS, Spence HD, Stelling J, Takahashi K, Tomita M, Wagner J, Wang J. SBML Forum. Bioinformatics. 2003; 19(4):524–31.View ArticlePubMedGoogle Scholar
- Bornstein BJ, Keating SM, Jouraku A, Hucka M. LibSBML: An API Library for SBML. Bioinformatics. 2008; 24(6):880–1. doi:http://dx.doi.org/10.1093/bioinformatics/btn051.
- Kirk DB, Hwu W-mW. Programming Massively Parallel Processors. A Hands-on Approach. Burlington: Morgan Kaufmann; 2010.Google Scholar
- Wong WW, Tsai TY, Liao JC. Single-cell zeroth-order protein degradation enhances the robustness of synthetic oscillator. Mol Syst Biol. 2007; 3:130.View ArticlePubMedPubMed CentralGoogle Scholar
- Woods ML, Leon M, Perez-Carrasco R, Barnes CP. A Statistical Approach Reveals Designs for the Most Robust Stochastic Gene Oscillators. ACS Synthetic Biol. 2016; 5(6):459–70.View ArticleGoogle Scholar
- Lu M, Jolly MK, Gomoto R, Huang B, Onuchic J, Ben-Jacob E. Tristability in cancer-associated microRNA-TF chimera toggle switch. J Phys Chem B. 2013; 117(42):13164–74.View ArticlePubMedGoogle Scholar
- Huang D, Holtz WJ, Maharbiz MM. A genetic bistable switch utilizing nonlinear protein degradation. J Biol Eng. 2012; 6(1):1–13. doi:http://dx.doi.org/10.1186/1754-1611-6-9.
- Clewley R. Hybrid models and biological model reduction with PyDSTool. PLoS Comput Biol. 2012; 8(8):1002628.View ArticleGoogle Scholar
- Ghaffarizadeh A, Flann NS, Podgorski GJ. Multistable switches and their role in cellular differentiation networks. BMC Bioinformatics. 2014; 15 Suppl 7:7.View ArticleGoogle Scholar
- Cinquin O, Demongeot J. High-dimensional switches and the modelling of cellular differentiation. J Theor Biol. 2005; 233(3):391–411.View ArticlePubMedGoogle Scholar
- Canton B, Labno A, Endy D. Refinement and standardization of synthetic biological parts and devices. Nat Biotechnol. 2008; 26(7):787–93.View ArticlePubMedGoogle Scholar
- Kelly JR, Rubin AJ, Davis JH, Ajo-Franklin CM, Cumbers J, Czar MJ, de Mora K, Glieberman AL, Monie DD, Endy D. Measuring the activity of BioBrick promoters using an in vivo reference standard. J Biol Eng. 2009; 3(1):4.View ArticlePubMedPubMed CentralGoogle Scholar
- Salis HM, Mirsky EA, Voigt CA. Automated design of synthetic ribosome binding sites to control protein expression. Nat Biotechnol. 2009; 27(10):946–50.View ArticlePubMedPubMed CentralGoogle Scholar