Bridging time scales in cellular decision making with a stochastic bistable switch
© Waldherr et al; licensee BioMed Central Ltd. 2010
Received: 1 March 2010
Accepted: 9 August 2010
Published: 9 August 2010
Cellular transformations which involve a significant phenotypical change of the cell's state use bistable biochemical switches as underlying decision systems. Some of these transformations act over a very long time scale on the cell population level, up to the entire lifespan of the organism.
In this work, we aim at linking cellular decisions taking place on a time scale of years to decades with the biochemical dynamics in signal transduction and gene regulation, occuring on a time scale of minutes to hours. We show that a stochastic bistable switch forms a viable biochemical mechanism to implement decision processes on long time scales. As a case study, the mechanism is applied to model the initiation of follicle growth in mammalian ovaries, where the physiological time scale of follicle pool depletion is on the order of the organism's lifespan. We construct a simple mathematical model for this process based on experimental evidence for the involved genetic mechanisms.
Despite the underlying stochasticity, the proposed mechanism turns out to yield reliable behavior in large populations of cells subject to the considered decision process. Our model explains how the physiological time constant may emerge from the intrinsic stochasticity of the underlying gene regulatory network. Apart from ovarian follicles, the proposed mechanism may also be of relevance for other physiological systems where cells take binary decisions over a long time scale.
The dynamics of biological systems span a wide range of temporal and spatial scales. The interactions among dynamical properties on different scales govern the overall behavior of the biological system, and thus form an important area of computational research in biology . A particularly interesting question in this context is how the behavior on a slow time scale emerges mechanistically from the dynamics on fast time scales. For example, how do cell population dynamics in tissues, which may evolve on a time scale of months, years or even decades, originate from the dynamics of the underlying gene regulatory networks, with a time scale of just minutes to hours?
In this work, we aim at bridging the time scale from gene regulation to cellular transformation processes on the tissue or cell population level. We specifically consider cellular transformation processes based on a bistable biochemical switch. Such switches have two distinct stable stationary states, and the cell initiates a transformation when the switch changes from one stable state to the other one. Bistable switches have previously been used to model a large number of cellular transformation events, such as progression through cell cycle arrest in the maturation of Xenopus oocytes [2, 3] or initiation of programmed cell death  and cellular differentiation  in higher organisms. Most models for these systems are constructed as deterministic models, and thus an external stimulus is required to induce changes in the switch's state. In addition, stochastic models for biochemical switches within a variety of biological processes have been formulated, for example the lac operon in E. coli[6, 7], the genetic toggle switch , or a generic phosphorylation/dephosporylation cycle . The typical questions that have been adressed by stochastic switch models are for example the steady state probability distribution of the different possible states of the switch , or the residence times in these states . In the previously proposed stochastic models of bistable biochemical switches, cells are able to switch forth and back between the possible qualitative states of the switch. While this is appropriate if the switch serves to choose a cellular state based on environmental conditions, such as for example in the galactose utilization network in yeast , this feature should not be held up for transformation processes. In transformation processes, subsequent mechanisms, which are not included in the model description, are in place to ensure irreversibility once the switch changed its qualitative state from the initial condition. The most obvious example for such mechanisms is cell death, where the model of the biochemical switch does not hold anymore once the cell transitions to the "dead" state.
In this work, we consider irreversible transformation processes based on a stochastic switch model, which apparently do not require any external stimulus to be initiated, where the transition is based only on stochastic fluctuations. Despite the stochasticity, we see in this paper that the dynamics of the switch still follow reliable temporal characteristics. Reliable thereby means that in a large population of cells, the number of cells that have already initiated the transformation can be described deterministically with high accuracy. We propose a generic transformation process, where a phenotypical change in the state of a cell is initiated as soon as a bistable biochemical switch changes its internal state. In previous studies, random switching caused by internal fluctuations is usually attributed to pathological events . In the mechanism proposed here, random switching has a regular physiological function.
A striking example for the kind of transformation processes we aim to describe is involved in mammalian oocyte maturation. In mammalian females, all or almost all of the oocytes that will ovulate through the organism's life-span are already present at birth or shortly thereafter as a population of so-called primordial follicles. Throughout the organism's reproductive life, follicles undergo the primordial to primary transition, which marks the start of a development process that will eventually lead to either ovulation or removal of the oocyte through atresia [12, 13]. In this way, there is a steady supply of mature follicles for ovulation, while the pool of primordial follicles is gradually depleted. The mechanisms through which the follicle transition is initiated are largely unknown, although a number of ovarian factors that may be relevant have been identified experimentally [14–16]. Importantly, the transition seems to be regulated locally in the ovary, and not through the endocrine system . An astonishing observation in this process is that in one follicle, the transition may occur already a few months after generation of the primordial follicle pool, while another follicle may stay several decades (for organisms with a sufficiently long lifespan) in the resting stage before growth is initiated. From the medical side, a misregulation of this process is implicated in premature ovarian failure due to follicle depletion, which is a major reason for infertility in human females. By way of a case study, we apply the proposed transformation mechanism to the problem of growth initiation in ovarian follicles. Including also cell-cell interactions supported by experimental evidence, we obtain a physiologically plausible model for this process, showing very good agreement with human clinical data on a time scale of several decades.
Deterministic model of a bistable switch
Nominal parameter values for the bistable switch model (1).
with suitable parameter L. With model parameters as given in Table 1, a suitable choice which we will use in this work is L = 55.
Stochastic model of a bistable switch
for , where the reaction propensities v i , i = 1,..., 5, are the same expressions as in the deterministic model (1).
In the stochastic description, we can compute the probabilities that the switch is in any of its two macrostates directly from a solution of the CME. Define p off (t) and p on (t) as the probabilities that the switch is off and on, respectively. Given a solution of the CME, these can be computed by summing up the probabilities that the system is in the corresponding microstates, i.e. , and equivalently for p off (t).
A transformation process modelled with a stochastic switch
Cellular transformation processes are often based on a bistable biochemical or genetic switch. In the initial state of the cell, the switch would be in the off state. Switching to the on state implies a significant change in the amount of an involved signaling molecule, e.g. a transcription factor. If the on state is maintained for some time, this change would result in a larger phenotypical change of the cell, e.g. through significant changes in gene expression. The mechanisms that induce this change are not included in the stochastic switch model, but from a signaling perspective downstream of it.
Most transformation processes rely on specific external stimuli, and the cell will initiate the transformation upon encountering the required stimulus. There are however examples where such a stimulus is not strictly required, and this is the case that we are dealing with in this paper. Moreover, we will focus on the behavior of cell populations, studying the problem how the temporal dynamics of the transformation process evolve in a pool of many cells.
The basic mechanism that actually triggers the bistable switch in our model without an external stimulus are the intrinsic fluctuations of concentrations in any biochemical reaction network, that are due to the stochastic nature of chemical reactions. As a rare event, these fluctuations may become so large that the microstate of the system crosses the separatrix between the domains of attraction in the deterministic system. As a consequence, the microstate around the other stable equilibrium point will become strongly attractive, and the switch will change its macrostate to on with a high probability. In this paper, we assume that the transformation is irreversible, which fits well to the process of follicle growth initiation. Also other processes such as programmed cell death are irreversible.
The described transformation process is easily modelled as a continuous-time Markov process. If the switch is in the macrostate off, then we directly use the microstates and transition probabilities of the underlying biochemical reaction network to model the transformation process. To account for the irreversibility of the transformation, all microstates are collapsed to one state of the Markov process, labeled with "on" in Figure 2B, which is an absorbing state. The transitions of other microstates to the absorbing state are governed by the propensity functions for the corresponding transitions in the underlying biochemical network. The resulting state space for the Markov process model of the transformation process is shown in Figure 2B.
where A ∈ ℝn × nis the state transition matrix. The matrix A can be computed directly from the values of the reaction propensity functions in each microstate . The differential equation (5) can be solved using standard tools for numerical integration. For the results described in this paper, we used the ode15s function in MATLAB (The MathWorks, Natick, MA) to obtain a numerical solution of (5).
Results and Discussion
A hypothetical mechanism for oocyte maturation
In this section, we suggest a biochemical mechanism that offers a molecular explanation for the large depletion times of several decades in the human oocyte pool. The model is based on experimental evidence obtained in a very informative series of studies by Skinner and colleagues (see  for a review), where the influence of several ovarian factors on the primordial to primary transition as well as some interactions between them have been studied. Because a positive feedback loop is necessary for a bistable switch , we have specifically searched for such an interconnection.
Primary ovarian follicles are composed by three main cell types: a single oocyte as the main component, and granulosa and theca cells surrounding the oocyte . Experimental evidence suggests a positive feedback circuit involving two ovarian factors that are relevant in the primordial to primary transition: the factor KIT ligand (KITL) is produced by granulosa cells and stimulates both the oocyte and surrounding theca cells to promote follicle development. Moreover, KITL stimulates the production of both keratinocyte growth factor (KGF) and hepatocyte growth factor (HGF) in the surrounding theca cells. KGF and HGF themself stimulate the production of KITL in the granulosa cells, thus providing a positive feedback loop . Moreover, the oocyte of primordial and developing follicles produces basic fibroblast growth factor (bFGF), which acts on surrounding granulosa cells and has been shown to increase the expression of KITL .
The stochastic switch generates reliable long-term behavior
where c1 > 0, 0 <λ1 < Re (λ i ) for i = 2, ..., n, and the c i (t) are polynomials in t. The mathematical derivation of (6) is provided in the appendix.
The probability distribution P (N r (t) = N) for the population size in the transformation process considered in this paper is shown in Figure 5B as a function of both cell number N and time t. The number of initial cells N0 = 106 was chosen from the reported range of ovarian follicles, 7 · 105 to 2 · 106 in human females at birth . For each point in time, the distribution has a very sharp peak, which indicates that the average value E[N r (t)] is a reliable prediction for the number of cells that have already undergone the transformation at a given time.
The probability density function for the depletion time is computed by taking the derivative of the cumulative probability distribution function (9). The resulting probability density function for the depletion time in follicle growth initiation is shown in Figure 5C. From the density function, the expected value and the standard deviation are obtained as E[T d ] = 59.1 years and years, respectively.
The expected value for T d can also be computed by solving . Using (6), it can thus be approximated by , where λ1 is the dominant decay rate of the process.
Next, we compare the computed statistical characteristics of the follicle depletion process to medical data. Explicit follicle counts are only sparsely available. However, the available pieces of data indicate that fluctuations in actual follicle numbers are larger than predicted by our model . Concerning the depletion time, a recent medical study suggests an average age of 51.1 years for the onset of menopause, with a standard deviation of 3.8 years . Our model predicts a depletion time of T d = 59.1 years, which is reasonably close to the experimentally observed depletion time. However, the standard deviation of 0.27 years in our model is significantly less than observed from medical data. In summary, even though our model is based on a highly stochastic process, the analysis reveals that it leads to much more reliable temporal characteristics than observed in the real system. This indicates that stochastic effects alone may not be sufficient to explain the heterogeneity observed in the follicle depletion process.
Expected depletion times (years) in the single cell model (5).
1.2 · 105
2.0 · 107
9.6 · 103
Increased robustness by interactions on the population level
In the last section, we have characterized the properties of the transformation process based on a bistable switch, with the depletion time of a pool of cells subject to the transformation as characteristic output of the model. We have shown that the proposed model produces reliable depletion times, in the sense of a small standard deviation, for fixed values of the biochemical parameters. However, we have also observed that the average depletion time in the basic model is quite sensitive to variations in the biochemical parameters. Clearly, this large sensitivity is not acceptable in a model that should be a meaningful representation of the primordial to primary follicle transition. In this section, we propose an additional mechanism that reduces the sensitivity of the average depletion time significantly.
Nominal parameter values for the population model (11).
8.2 · 104
Expected depletion times (years) in the population model (11).
In this paper, we deal with a fundamental question in the development of multicellular organisms: How can biochemical reactions and genetic mechanism acting on the scale of minutes in individual cells generate dynamics with characteristic times of years to decades on the tissue level? As a possible mechanism to achieve this, we propose a generic transformation process based on a bistable stochastic switch. From the underlying genetic interactions and biochemical reactions, the process can be modelled as a continuous-time Markov process. We show that the proposed stochastic mechanism generates reliable long-term behavior on the population level, with cells undergoing the transformation with an exponentially decaying rate. Thereby, the decay rate is equal to the dominant eigenvalue of the transition matrix describing the underlying biochemical network. Due to bistability of the considered switch, this dominant eigenvalue corresponds to very slow dynamics, thus leading to the very long timescale as observed in the simulations. We pose the hypothesis that a biological instance of this mechanism is present in the development of ovarian follicles. To describe this process, we constructed a simple model of a bistable switch in the primordial to primary transition for ovarian follicles. The model is based on experimentally determined factors and their interactions in the different cell types constituting the ovarian follicles. Although it is not assured that a bistable switch in ovarian follicles will indeed be based on the factors that we have used here, the basic mechanism would work equivalently well with other factors.
Despite its simplicity, our model explains well how the long-term characteristics of follicle development may reliably be generated by biochemical reactions occurring on much shorter time scales. Keeping the model simple serves two purposes: first, it shows that the dynamics of follicle growth initiation can be generated by a quite simple mechanism. Clearly, additional pathways and regulatory feedback interactions that we have not included in this model can be expected to be present in the system. These may serve to increase robustness of the network, or to provide additional inputs to control the transition rate, e.g. for the endocrine system. Second, the simplicity of the model allows us to solve the chemical master equation for the network numerically, and thus to obtain a good quantitative description of the model behavior.
As a possible shortcoming of the basic model on the single cell level, we observe an unrealistic large sensitivity of the follicle depletion time with respect to parameter variations. By adding the experimentally supported inhibition of follicle growth initiation by later-stage growing follicles, the sensitivity of the depletion time could be reduced significantly. Apart from the inhibition included in the model, other interactions among individual follicles seem to play a role in the primordial to primary transition . We envision that the inclusion of more regulatory interactions may further decrease the sensitivity of the depletion time with respect to parameter variations to a physiologically plausible level.
Appendix: Computation of the transition probability
In this section, we prove that the probability that a given cell has undergone the considered transformation process is given by p on (t) as in (6). The proof is based on considering the solution of the underlying CME (5).
where A rev ∈ ℝ(n - 1) × (n - 1)describes the interactions among the non-absorbing states, and a abs ∈ ℝ1 × (n - 1)describes the transition propensities to the absorbing state.
for i = 1, ..., n - 1 i.e. A rev is diagonally dominant. Thus, Gersgorin's theorem  asserts that all eigenvalues of A rev have a non-positive real part. Even more, since a abs is non-zero, (13) holds with a strict inequality for at least one i. Thus, by Theorem 10.7.2 in , all eigenvalues of A rev have negative real part. By the properties of the considered biochemical network, A rev is irreducible, and its off-diagonal elements are non-negative. From Corollary 4.3.2 in , it follows that A rev has an eigenvalue λ1 ∈ ℝ with algebraic multiplicity 1 and a strictly positive corresponding eigenvector v1 such that Re λ <λ1 for all λ ≠ λ1 in the spectrum of A rev .
where and .
SW and FA acknowledge support by the German Research Foundation (DFG) through the Cluster of Excellence in Simulation Technology (EXC 310) at the University of Stuttgart.
- Martins ML, Ferreira SC, Vilela MJ: Multiscale models for biological systems. Curr Opin Colloid Interface Sci. 2010, 15: 18-23. 10.1016/j.cocis.2009.04.004.View ArticleGoogle Scholar
- Ferrell JE, Xiong W: Bistability in cell signaling: How to make continuous processes discontinuous, and reversible processes irreversible. Chaos. 2001, 11: 227-236. 10.1063/1.1349894View ArticlePubMedGoogle Scholar
- Ferrell JE, Machleder EM: The biochemical basis of an all-or-none cell fate switch in Xenopus oocytes. Science. 1998, 280 (5365): 895-898. 10.1126/science.280.5365.895View ArticlePubMedGoogle Scholar
- Eissing T, Conzelmann H, Gilles ED, Allgöwer F, Bullinger E, Scheurich P: Bistability Analyses of a Caspase Activation Model for Receptor-induced Apoptosis. J Biol Chem. 2004, 279 (35): 36892-97. 10.1074/jbc.M404893200View ArticlePubMedGoogle Scholar
- Chickarmane V, Enver T, Peterson C: Computational modeling of the hematopoietic erythroid-myeloid switch reveals insights into cooperativity, priming, and irreversibility. PLoS Comput Biol. 2009, 5: e1000268- 10.1371/journal.pcbi.1000268PubMed CentralView ArticlePubMedGoogle Scholar
- Mettetal JT, Muzzey D, Pedraza JM, Ozbudak EM, van Oudenaarden A: Predicting stochastic gene expression dynamics in single cells. Proc Natl Acad Sci. 2006, 103: 7304-9. 10.1073/pnas.0509874103PubMed CentralView ArticlePubMedGoogle Scholar
- Kaufmann BB, Yang Q, Mettetal JT, van Oudenaarden A: Heritable Stochastic Switching Revealed by Single-Cell Genealogy. PLoS Biology. 2007, 5 (9): e239- 10.1371/journal.pbio.0050239PubMed CentralView ArticlePubMedGoogle Scholar
- Tian T, Burrage K: Stochastic models for regulatory networks of the genetic toggle switch. Proc Natl Acad Sci. 2006, 103 (22): 8372-8377. 10.1073/pnas.0507818103PubMed CentralView ArticlePubMedGoogle Scholar
- Krishnamurthy S, Smith E, Krakauer D, Fontana W: The stochastic behavior of a molecular switching circuit with feedback. Biol Direct. 2007, 2: 13- 10.1186/1745-6150-2-13PubMed CentralView ArticlePubMedGoogle Scholar
- Song C, Phenix H, Abedi V, Scott M, Ingalls BP, Kaern M, Perkins TJ: Estimating the stochastic bifurcation structure of cellular networks. PLoS Comput Biol. 2010, 6 (3): e1000699- 10.1371/journal.pcbi.1000699PubMed CentralView ArticlePubMedGoogle Scholar
- Isaacs FJ, Hasty J, Cantor CR, Collins JJ: Prediction and measurement of an autoregulatory genetic module. Proc Natl Acad Sci. 2003, 100 (13): 7714-7719. 10.1073/pnas.1332628100PubMed CentralView ArticlePubMedGoogle Scholar
- Fortune JE, Cushman RA, Wahl CM, Kito S: The primordial to primary follicle transition. Mol Cell Endocrinol. 2000, 163 (1-2): 53-60. 10.1016/S0303-7207(99)00240-3View ArticlePubMedGoogle Scholar
- Skinner MK: Regulation of primordial follicle assembly and development. Hum Reprod Update. 2005, 11 (5): 461-471. 10.1093/humupd/dmi020View ArticlePubMedGoogle Scholar
- Parrott JA, Skinner MK: Kit-ligand/stem cell factor induces primordial follicle development and initiates folliculogenesis. Endocrinol. 1999, 140 (9): 4262-4271. 10.1210/en.140.9.4262.Google Scholar
- Castrillon DH, Miao L, Kollipara R, Horner JW, DePinho RA: Suppression of ovarian follicle activation in mice by the transcription factor Foxo3a. Science. 2003, 301 (5630): 215-218. 10.1126/science.1086336View ArticlePubMedGoogle Scholar
- Nilsson EE, Skinner MK: Kit ligand and basic fibroblast growth factor interactions in the induction of ovarian primordial to primary follicle transition. Mol Cell Endocrinol. 2004, 214 (1-2): 19-25. 10.1016/j.mce.2003.12.001View ArticlePubMedGoogle Scholar
- Braw-Tal R: The initiation of follicle growth: the oocyte or the somatic cells?. Mol Cell Endocrinol. 2002, 187 (1-2): 11-18. 10.1016/S0303-7207(01)00699-2View ArticlePubMedGoogle Scholar
- Ross J: mRNA stability in mammalian cells. Microbiol Rev. 1995, 59 (3): 423-450.PubMed CentralPubMedGoogle Scholar
- Hargrove JL, Hulsey MG, Beale EG: The kinetics of mammalian gene expression. Bioessays. 1991, 13 (12): 667-674. 10.1002/bies.950131209View ArticlePubMedGoogle Scholar
- Raj A, Peskin CS, Tranchina D, Vargas DY, Tyagi S: Stochastic mRNA Synthesis in Mammalian Cells. PLoS Biol. 2006, 4 (10): e309- 10.1371/journal.pbio.0040309PubMed CentralView ArticlePubMedGoogle Scholar
- Warren L, Bryder D, Weissman IL, Quake SR: Transcription factor profiling in individual hematopoietic progenitors by digital RT-PCR. Proc Natl Acad Sci. 2006, 103 (47): 17807-17812. 10.1073/pnas.0608512103PubMed CentralView ArticlePubMedGoogle Scholar
- Eissing T, Waldherr S, Allgöwer F, Scheurich P, Bullinger E: Steady state and (bi-)stability evaluation of simple protease signalling networks. BioSystems. 2007, 90: 591-601. 10.1016/j.biosystems.2007.01.003View ArticlePubMedGoogle Scholar
- Cherry JL, Adler FR: How to make a Biological Switch. J Theor Biol. 2000, 203 (2): 117-133. 10.1006/jtbi.2000.1068View ArticlePubMedGoogle Scholar
- van Kampen NG: Stochastic processes in physics and chemistry. 1981, North-Holland AmsterdamGoogle Scholar
- Gillespie DT: A rigorous derivation of the chemical master equation. Physica A: Statist Theor Phys. 1992, 188 (1-3): 404-425. 10.1016/0378-4371(92)90283-V.View ArticleGoogle Scholar
- Rao CV, Arkin AP: Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm. J Chem Phys. 2003, 118: 4999-5010. 10.1063/1.1545446.View ArticleGoogle Scholar
- Nuño J, Tarazona P: Lifetimes of small catalytic networks. Bull Math Biol. 1994, 56 (5): 875-898.View ArticleGoogle Scholar
- Munsky B, Khammash M: The finite state projection algorithm for the solution of the chemical master equation. J Chem Phys. 2006, 124 (4): 044104- 10.1063/1.2145882View ArticlePubMedGoogle Scholar
- Kaufman M, Soule C, Thomas R: A new necessary condition on interaction graphs for multistationarity. J Theor Biol. 2007, 248 (4): 675-685. 10.1016/j.jtbi.2007.06.016View ArticlePubMedGoogle Scholar
- Parrott JA, Skinner MK: Thecal cell-granulosa cell interactions involve a positive feedback loop among keratinocyte growth factor, hepatocyte growth factor, and Kit ligand during ovarian follicular development. Endocrinol. 1998, 139 (5): 2240-2245. 10.1210/en.139.5.2240.Google Scholar
- Yeh J, Adashi EY: The ovarian life cycle. Reproductive Endocrinology. Edited by: Yen SSC, Jaffe RB. 1999, 153-190. Saunders PhiladelphiaGoogle Scholar
- Faddy MJ, Gosden RG: A model conforming the decline in follicle numbers to the age of menopause in women. Hum Reprod. 1996, 11 (7): 1484-1486.View ArticlePubMedGoogle Scholar
- Gougeon A, Ecochard R, Thalabard JC: Age-related changes of the population of human ovarian follicles: increase in the disappearance rate of non-growing and early-growing follicles in aging women. Biol Reprod. 1994, 50 (3): 653-663. 10.1095/biolreprod50.3.653View ArticlePubMedGoogle Scholar
- de Bruin JP, Bovenhuis H, van Noord PA, Pearson PL, van Arendonk JA, te Velde ER, Kuurman WW, Dorland M: The role of genetic factors in age at natural menopause. Hum Reprod. 2001, 16 (9): 2014-2018. 10.1093/humrep/16.9.2014View ArticlePubMedGoogle Scholar
- Nilsson E, Rogers N, Skinner MK: Actions of anti-Mullerian hormone on the ovarian transcriptome to inhibit primordial to primary follicle transition. Reproduction. 2007, 134 (2): 209-221. 10.1530/REP-07-0119View ArticlePubMedGoogle Scholar
- Visser JA, Themmen APN: Anti-Müllerian hormone and folliculogenesis. Mol Cell Endocrinol. 2005, 234 (1-2): 81-86. 10.1016/j.mce.2004.09.008View ArticlePubMedGoogle Scholar
- Faddy MJ, Gosden RG: A mathematical model of follicle dynamics in the human ovary. Hum Reprod. 1995, 10 (4): 770-775.PubMedGoogle Scholar
- Silva-Buttkus PD, Marcelli G, Franks S, Stark J, Hardy K: Inferring biological mechanisms from spatial analysis: Prediction of a local inhibitor in the ovary. Proc Natl Acad Sci. 2009, 106: 456-461. 10.1073/pnas.0810012106PubMed CentralView ArticlePubMedGoogle Scholar
- Lancaster P, Tismenetsky M: The Theory of Matrices. 1985, San Diego: Academic PressGoogle Scholar
- Smith HL: Monotone Dynamical Systems. An Introduction to the Theory of Competitive and Cooperative Systems, Volume 41 of. 1995, Mathematical Surveys and Monographs. Providence, Rhode Island: American Mathematical SocietyGoogle Scholar
- Gillespie DT: Exact stochastic simulation of coupled chemical reactions. J Phys Chem. 1977, 81 (25): 2340-2361. 10.1021/j100540a008.View ArticleGoogle Scholar
- Ramsey S, Orrell D, Bolouri H: Dizzy: Stochastic Simulation of Large-Scale Genetic Regulatory Networks. J Bioinform Comput Biol. 2005, 3 (2): 415-436. 10.1142/S0219720005001132View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.