Skip to main content
  • Research article
  • Open access
  • Published:

The smallest chemical reaction system with bistability



Bistability underlies basic biological phenomena, such as cell division, differentiation, cancer onset, and apoptosis. So far biologists identified two necessary conditions for bistability: positive feedback and ultrasensitivity.


Biological systems are based upon elementary mono- and bimolecular chemical reactions. In order to definitely clarify all necessary conditions for bistability we here present the corresponding minimal system. According to our definition, it contains the minimal number of (i) reactants, (ii) reactions, and (iii) terms in the corresponding ordinary differential equations (decreasing importance from i-iii). The minimal bistable system contains two reactants and four irreversible reactions (three bimolecular, one monomolecular).

We discuss the roles of the reactions with respect to the necessary conditions for bistability: two reactions comprise the positive feedback loop, a third reaction filters out small stimuli thus enabling a stable 'off' state, and the fourth reaction prevents explosions. We argue that prevention of explosion is a third general necessary condition for bistability, which is so far lacking discussion in the literature.

Moreover, in addition to proving that in two-component systems three steady states are necessary for bistability (five for tristability, etc.), we also present a simple general method to design such systems: one just needs one production and three different degradation mechanisms (one production, five degradations for tristability, etc.). This helps modelling multistable systems and it is important for corresponding synthetic biology projects.


The presented minimal bistable system finally clarifies the often discussed question for the necessary conditions for bistability. The three necessary conditions are: positive feedback, a mechanism to filter out small stimuli and a mechanism to prevent explosions. This is important for modelling bistability with simple systems and for synthetically designing new bistable systems. Our simple model system is also well suited for corresponding teaching purposes.


Bistability is key for understanding basic phenomena of cellular functioning, such as decision-making processes in cell cycle progression, cell differentiation, and apoptosis [1]. It is also involved in loss of cellular homeostasis associated with early events in cancer onset [2] and in prion diseases [3]. A recent review discussed different bistability phenomena in bacteria, such as different phenotypes in clonal populations being important for the origin of new species [4].

Bistable switches are typically enabled by positive feedback loops in signal transduction networks. Here a sufficiently strong (external) signal switches on a self-amplifying process leading to expression of the corresponding target genes. Due to the corresponding hysteresis effect this process can retain its activity without a persistent signal. Such switches are therefore called 'decision-making'. One often discussed example is the restriction point control for the regulation of G1-S transition of the mammalian cell cycle, where recently a detailed small ordinary differential equation (ODE) model was presented [5]. The G2-M transition was also described as a toggle-switch [6]. Oocyte maturation is an example for the involvement of a bistable system in cell differentiation [7, 8]. Biochemical switches have also been found in nutrient utilization in bacteria [9], mating response in yeast [10], and synaptic memory processing [11]. Interestingly, just quite small parts of the large signal transduction systems (for instance a single layer of a MAPK cascade) can already induce bistability. This was demonstrated for the epidermal growth factor receptor system [12] and other kinase phosphatase systems [13, 14].

Given its outstanding biological importance, it is clear that bistable switches also attracted the attention of theoretical biologists. A frequently discussed problem is the necessary (and/or sufficient) condition for bistability. The central result goes back to the work of Clarke [15] and Thomas [16]: autonomous differential systems can possess multiple steady states only under the presence of positive feedback loops [17, 18]. It was also argued that the feedback needs 'some type of non-linearity' or 'ultrasensitivity' for inducing bistability [19]. Different types of ultrasensitivity have been discussed [20], but formulated in the most general manner that means the system needs some mechanism for filtering out small stimuli to enable a stable 'off' state [2022]. Feinberg's chemical reaction network theory (CRNT) even gives necessary and sufficient conditions for bistability, by restricting to special mass-action kinetic (MAK) systems [23]. In a recent application example a single layer of a MAPK cascade was studied and the region in parameter space being relevant for bistability was analytically described [24]. However, for larger systems CRNT leads to cumbersome calculations, but the applicability could recently be improved by just studying important subnetworks [25] which are based on the concept of elementary flux modes [26].

Another approach for identifying necessary structural conditions for any dynamic behaviour is the identification of the corresponding minimal systems. Such systems have the advantage of being "simple enough to understand at an intuitive level" [19] and are well suited for different basic studies. For instance, the Lotka-Volterra system [27, 28], the Higgins-Selkov-oscillator [2931], and the "Brusselator" [32] have been studied extensively. Some years ago we identified the smallest chemical system with Hopf bifurcation [33]. Minimal MAK systems are summarized in Table 1. Recently, different 'smallest' or 'minimal' bistable systems for cell polarity [34] and G protein signalling [35] have been presented, as well as 'the smallest multistationary mass-preserving chemical reaction network' [36]. However, these systems are still too large to represent a minimal bistable system according to definition (1). Many different reaction topologies with 3 and 4 molecules have also been analysed computationally for the possibility of bistability [37]. Although this type of bistability detection may miss some bistable systems, the authors found nevertheless many topologies with switching behaviour (10% of tested configurations). The identified 'minimal' system contained 3 variables (5 reactants, 2 conservation relations) and 6 reactions. The bistable one-dimensional Schloegl system [38] contains trimolecular reactions and can therefore not represent a realistic elementary chemical system. Elementary chemical reactions are at most bimolecular.

Table 1 Distinguished minimal MAK systems

Here we present and discuss the smallest bistable chemical reaction system. Application of our previously presented Instability Causing Structure Analysis (ICSA [39]) leads to additional insight into system functioning.


The smallest bistable chemical reaction system

We define the smallest chemical system (contains only mono- and bimolecular reactions, reversible reactions are considered as two irreversible ones) by the following criteria in decreasing order of importance:


According to this definition, the following bistable system is unique (Methods section contains the proof for this statement).


Assuming spatially homogeneous conditions, the system can be described by the two-component mass-action kinetic ODE system (S is incorporated into k1):


Due to its simplicity, the mathematical analysis of the system is simple as well. The system has two elementary flux modes [26], following directly from the two nullspace vectors of the corresponding stoichiometric matrix


One mode uses reactions 1-3, and the other reactions 1,2, and 4. Bistability can, of course, only arise if all reactions are active. Introducing dimensionless quantities (x/cx, y/cy, k1/(k2c) → k1, k3/k2k3, k4/(k2c) → k4, tk2ct), we set k2 = 1, without restriction of generality (k1, k3, k4 > 0). The system has three steady states: , , with the discriminant D = k1 - 4k3k4. A saddle-node bifurcation occurs at D = 0, the three steady states are real if D > 0. The second and third steady state is always positive.

Generally, in two-component systems a steady state is locally stable if the trace tr and determinant det of the Jacobian at this point are negative and positive, respectively (node if 4det <tr2, focus otherwise). If the corresponding determinant is negative, the steady state is a saddle point. It can be seen from the Jacobian that its trace is always negative (phase flow of system (2) is confined to the positive part of the phase space). This excludes Hopf bifurcations (arising at tr = 0) and it means the system is dissipative, i.e. phase-space contracting all over the phase space (trace = two-dimensional Lyapunov exponent). The determinants of the Jacobian at the three steady states read k1k4, , and , respectively. Therefore, the first and third steady state are always locally stable, the second is always locally unstable, a saddle-point (determinant of Jacobian at second steady state always negative, cf. point 3 in Methods). Simple calculation shows further that 4det <tr2 at the first and third steady state, so these are always stable nodes.

For k1 = 8, k3 = 1, k4 = 1.5 the second and third steady state are and , respectively. Figure 1 shows a corresponding signal-response curve, also called bifurcation diagram [6]. The signal is the concentration of the constant outer substance S (assuming for the bimolecular rate constant k bi = 1, the concentration S is identical to the apparent rate constant k1 = k bi S), the response is the steady state concentration of an internal reactant, here X. The saddle-node bifurcation occurs at S = 3/4. Beyond that point the system has two stable steady states. It follows from that (for fixed other parameters), so the real dynamic behaviour of the system is that of a toggle-switch (Figure 1): for sufficiently large k1 small fluctuations in the concentrations would drive the system to the positive steady state (the 'on' state).

Figure 1
figure 1

Signal-response curve (bifurcation diagram) of system (2) for the parameters k1 = 8, k2 = 1, k3 = 1, k4 = 1.5. Solid lines indicate locally stable steady states, the dashed line locally unstable steady states. The inset shows the signal-response curve if an additional small constant influx into X (here 0.6) is assumed (enabling a positive 'off' state, leaving the 'on' state and bifurcation point nearly unchanged). This is the classical toggle switch (terminology of Tyson et al. (6), others use the term toggle switch to describe a double negative (i.e. positive) feedback loop (4)) picture enabling the hysteresis cycle: starting with low values and increasing the signal continuously increases the response, until the saddle-node bifurcation at about S = 1.7 is reached. Further increase of the signal leads to a sudden jump of the response to the upper steady state. Decreasing the signal now leads to a continuous decrease of the response, the systems stays in the upper steady state until the left bifurcation point is reached where the response jumps back to the lower steady state.

Figure 2 shows rate curves of system (2). It can be seen that the three crossings of production and degradation rate (i.e. the three steady states) are due to the different contributions of the three degradation terms. This implies a simple general procedure for designing bi- or multistable systems: a bistable system can be created with one function for production and three different functions for degradation, e.g. a linear, a quadratic, and a cubic one as in our simple example system. Accordingly, tristable systems require 5 different functions to enable 5 crossings (three stable and two unstable steady states, cf. point 3 in Methods), and so forth for more steady states. This observation helps constructing minimal and/or realistic models of more complicated multistable systems. It can also be a starting point for the design of real bistable systems, for instance in synthetic biology. Note that all enzyme kinetic rate laws can be modelled with polynomial ODEs [40].

Figure 2
figure 2

Rate curves [6] of system (2) for the parameters k1 = 8, k2 = 1, k3 = 1, k4 = 1.5. The thick solid line is the rate of the removal of reactant X (sum of the negative terms in ) and the thick dashed line the rate of production (positive term in ). The three crossings indicate the three steady states . The thin lines show the contributions of the three degradation terms separately: quadratic term k2x2 dashed, the effectively cubic term k3 xy solid, and the linear term k4 x dotdashed. The inset shows a zoomed version for x < 2.1.

The Instability Causing Structure Analysis (ICSA) of system (2)

Recently we presented a new method for topological network analysis of dynamical systems, the Instability Causing Structure Analysis (ICSA [39]). Standard stoichiometric network analyses (such as elementary flux mode calculations [26]) are based on the assumption of steady states and lead to linear constraints in flux space. ICSA is a nonlinear network analysis. It is based on the assumption of locally stable steady states. The additional demand for local stability yields additional nonlinear constraints for the steady state flux space [39].

ICSA needs no knowledge about kinetic details. It can be applied to (bio)chemical systems where just the stoichiometric matrix is known (or even just the signs of its elements) or to signal transduction/gene regulatory networks represented by interaction graphs [18] (also called incidence graph [21] or causal influence graph [41]). ICSA leads to additional insight into system functioning by identifying all contained feedback loops and all the corresponding instability causing structures (ICS). An ICS is either a single feedback loop or a special combination of feedback loops [39]. ICSA yields a necessary condition for local instability of any steady state of the system: if there is no ICS all potential steady states are locally stable. For two-dimensional systems this also implies that the system has just one steady state (cf. point 3 in Methods).

We apply ICSA for additional analysis of system (2). The stoichiometric matrix S is given in (2c). Multiplication of S with the reaction velocity substrates ector (contains the substrates for each reaction) (v1 (y) v2 (x) v3 (x, y) v4 (x))Tleads to . Differentiation yields the general Jacobian JG (39) of system (2):


where the indices x and y denote the corresponding partial differentiation. The off-diagonal elements represent the fundamental activating and inhibiting interactions in the system: the positive v2xin JG 21shows that x activates y by the second reaction, and equivalently for the two terms in JG 12: v1y→ y activates x by the first reaction, - v3y→ y inhibits x by the third reaction. Figure 3 shows the corresponding interaction (incidence) graph summarizing these interactions. Interaction graphs can often be found in the biological literature and corresponding databases (KEGG [42]; BIOBASE [43, 44]; Dynamic Signaling Maps ICSA [39] was developed for structural analyses of (bio)chemical systems (KEGG [42]; BRENDA [45]) AND such interaction graphs.

Figure 3
figure 3

Interaction graph of system (2). It follows directly from the off-diagonal elements of the general Jacobian (3). The positive feedback loop is the only instability causing structure (ICS) in the system, allowing for a locally unstable steady state (presupposition for bistability).

The system contains one positive and one negative feedback loop. The positive loop is the necessary structural condition for bistability [18]. The second often discussed ingredient for bistability is any mechanism for filtering out small stimuli to make the off-state stable. Ferrell and Xiong [20] discussed different such mechanisms, as ultrasensitivity and back reaction saturation. Our system has the simplest mechanism to stabilize the off-state: the monomolecular efflux reaction 4. Analysis shows that without this reaction the second (unstable) steady state merges with the zero off-state making it (weakly) unstable. In fact, Figure 2 (especially the inset) shows that it is the linear degradation term that makes the overall degradation rate higher then production for low X concentrations.

System (2) is the minimal chemical reaction system with bistability. Therefore, any ingredient is essential. That also means, without the negative loop (without reaction 3 being catalyzed by Y) this system cannot be bistable. In fact, the negative feedback prevents explosion of the system: without reaction 3 the system has just one locally stable off-state and one unstable positive steady state. Figure 2 shows that without the cubic term there are just two crossings of the rate curves.

Summarizing, system (2) contains the three different necessary conditions for bistability: (i) a positive feedback loop, (ii) a mechanism for filtering out small stimuli, and (iii) a mechanism for preventing explosion. Interestingly, the third condition lacks discussion in the literature [19, 20, 22]. Moreover, these three conditions can be related to the three degradation terms in the balance equation of X and so to the last three of the four reactions in (2). The first reaction represents the production of the system, it is the only input. The second reaction closes the positive feedback cycle. Reactions three and four are degradation reactions. The third is the main degradation for higher concentrations, it prevents explosion. The fourth reaction filters out the small stimuli.

The first step in the general ICSA procedure is the analysis of feedback cycles resulting from the off-diagonal terms of the general Jacobian. The second step is the identification of the topological structures that actually cause instability. As mentioned above, standard local stability analysis of two-component systems needs consideration of the Jacobian's trace and determinant det. In ICSA we study the general Jacobian. It follows from (3) that the trace is always negative and det = v2xv3y+ v1y(v3x+ v4x- v2x). The only negative term in det comprises reactions 1 and 2, so this positive feedback cycle is the only ICS of the system. The negative feedback cycle cannot cause any instability here.


The minimal bistable chemical reaction system (2) is based on definition (1). This definition is based on a chemical/physical point of view, but other definitions might also be possible. The definition for the smallest chemical reaction system with Hopf bifurcation [33], for instance, was more mathematically motivated: here minimal number of quadratic terms in the ODEs had a higher priority than minimal number of reactions. So far it has not been clear whether any at most bimolecular 3-variable system with only four (irreversible) reactions and Hopf bifurcation exists. Meanwhile we found a corresponding system (, , ) which has a supercritical Hopf bifurcation (e.g. at k1 = k2 = 1, ). This might be the minimal chemical reaction system with sustained oscillations according to definition (1).

Interestingly, system (2) is similar to a previously presented "minimal reaction network" [12] modelling activation of the epidermal growth factor receptor (EGFR). Our equation resembles the balance equation of the phosphorylated receptor tyrosine kinase (this superfamily contains EGFR) and is similar to the differential equation for active protein tyrosine phosphatase.

Bistable systems play important roles also beyond biology. They are usually depicted by the mechanical example of a ball rolling into two different valley basins. Bistable chemical systems, in particular, have been studied extensively to analyse relaxation kinetics [46], non-equilibrium thermodynamics [47], stochastic resonance [48], as well as climate change [49].

Positive feedback is clearly associated with bi/multistability. Negative feedback, in contrast, is often discussed in the context of oscillations [6], and we recently conjectured that this is indeed a necessary condition for (sustained) oscillations [39]. However, this contradicts statements as "sustained oscillations can occur in models based on positive or negative feedback" [50]. Obviously, different understandings of feedback exist. We suggest to use the general Jacobian J G (for system (2) it is given in (3)) for a simple definition of feedback: if the J G terms close any cycles, then feedback exists (e.g. the positive and negative feedback cycles of system (2)), otherwise not. Analysis of J G guarantees a unique identification of all in a system contained feedback cycles. Goldbeter [50] mentioned different examples where the oscillations should be based on a positive feedback, such as glycolytic, Ca2+, and cAMP oscillations. However, a more detailed analysis of these systems shows that a negative feedback (according to our definition) is always contained (results unpublished). An example is the simplest model for glycolytic oscillations, the Higgins-Selkov oscillator ([2931], Table 1): the corresponding general Jacobian is . Its trace and determinant are v2y- v2x- v3yand v2xv3y, respectively. Obviously, v2yin JG 22is the only instability causing term (a positive feedback), so the second reaction is the only ICS in the system. However, inspection of the off-diagonal elements of J G reveals a negative feedback as well: the larger x, the larger becomes y, but the larger y, the smaller becomes x. The same is realised in the bistable system (2), where the positive feedback is the only ICS and another negative feedback is contained. These examples show how the analysis of the general Jacobian helps clarifying the discussion of feedback loops.

We have shown that a mechanism for preventing explosions is a third necessary condition for bistability (complementing the previously discussed two other conditions positive feedback and filtering out of small stimuli). In system (2) this is achieved by a negative feedback. Other bistable systems contain negative feedbacks as well (e.g. ERK pathway [2]), so we hypothesize that this is indeed a typical feature of bistable systems. Interestingly, also oscillating systems typically contain (besides the necessary negative feedback) a positive feedback (for better tunable frequency, evolvability and robustness [51]). Thus, oscillating and bistable systems are practically based on the same set of feedback cycles.


Bi/multistability and oscillations are the two most important dynamic phenomena in biology. Limit cycle oscillations are associated with biological clocks and cell signalling [52], and spatial oscillations with proper cell division [53]. The fundamental importance of bistability is discussed in the introduction. Some years ago we presented the smallest chemical system with limit cycles [33]. Here we have derived the smallest chemical system with bistability (2).

Minimal systems are well suited for basic studies and for teaching purposes. This explains the great success of, for instance, the Lotka-Volterra [27, 28] and the Higgins-Selkov system [30]. We have demonstrated for the minimal chemical reaction system with Hopf bifurcation [33] that it is accessible for detailed mathematic-analytical examination [54] and a good example system for thermodynamic considerations [55]. We hope that also the minimal chemical system with bistability will serve for such purposes in the future.


System (2) is the smallest bistable chemical reaction system according to definition (1) - proof (inductive proof systematically considering all possibilities):

  1. 1.

    The Schloegl system (Table 1) is the smallest bistable one-variable (1d) system: 1d systems need an unstable steady state to separate the attractor regions of two stable steady states, so we need at least three steady states to realize a bistable system. The simplest function (which is realizable as a MAK system) f(x) with three zeros is the cubic polynomial. To realize a stable 'on' state the sign of the cubic term needs to be "-" (). For three different non-negative steady states we also need a positive quadratic and a negative linear term: f(x) = -a x3 +b x2-c x+d (a, b, c>0) possesses two positive extrema as can simply be seen considering f'(x) = 0. Thus, the minimal bistable 1d system reads . The minimal corresponding MAK system is the Schloegl system shown in Table 1 (a reversible monomolecular efflux reaction (d>0) allows for two positive stable steady states). It follows that a chemically realistic bistable system with only mono- and bimolecular reactions needs at least two variables.

  2. 2.

    Using our general quasi-steady-state-approximation procedure [56], we have previously shown that any irreversible trimolecular reaction can be understood as limit case of a reversible bimolecular reaction and another irreversible bimolecular reaction by introducing one additional intermediate [57]. Transforming the Schloegl system accordingly proves that at most bimolecular bistable 2d systems exist, i.e. the number of variable reactants in the minimal bistable chemical system is fixed to two (cf. definition (1)).

  3. 3.

    The lemma of the index sum [58] states that the sum of indices of all steady states within a two-dimensional confined set (closed region in phase space where all trajectories point inwards [59]) equals one. The index values of a node, a focus, and a saddle are +1, +1, and -1, respectively (for stable and unstable nodes and foci [58]). Our minimal system should therefore contain one unstable (saddle-point) and two stable steady states (we are only considering non-exploding systems, such that a confined set could simply be constructed, trajectories point inward at the boundary of the positive orthant anyway). To get three steady states we need at least a cubic steady state equation f(x) = 0, i.e. an x3 term (higher order polynomials would require more bimolecular reactions). In 2d MAK systems this can only be realized by one ODE with an xy term and the other with x2 and y terms. Inserting the corresponding steady state expression y = x2..., into the other ODE's xy term gives the cubic term. A direct x3 term is forbidden in bimolecular systems. The symmetric case y2 and x needs no extra consideration. These terms already correspond to at least three reactions. One can show that such three reactions are not sufficient to give a bistable system (the explicit proof is not necessary, because it turns out indirectly from the following analysis). So the minimal bistable system has at least four reactions.

  4. 4.

    To realize a -x3 term in the corresponding steady state equation, one either needs the terms y and x2 with different signs and -xy in the other ODE, or y and x2 with the same signs and +xy in the other ODE. But the same signs variant cannot work: is not chemical, would contain at least one trimolecular reaction, is not chemical, and is also impossible. The latter applies because: (i) the equation must also contain a term -xy (the corresponding reaction cannot be bimolecular otherwise: a term +xy in implies the bimolecular reaction X+Y->2X, implying the term -xy in ), (ii) detailed analysis of the system shows that the cubic term in the steady state equation always vanishes, independent of any other added mono- and/or bimolecular reactions. So the y and x2 terms must have different signs in one ODE.

  5. 5.

    There are two variants for the different-sign-case: 1. 2. (note that the x2 term has to appear also in the equation to be interpretable as bimolecular reaction). We show that no bistable system with only four reactions exists that corresponds to the first case: to realize a positive quadratic term in the cubic polynomial of the steady state equation (cf. discussion concerning the Schloegl system in 1.) we have three options: (i) , (ii) , (iii) . The second case could only be realized by 4 different reactions and would need a fifth reaction to yield a linear term in the cubic polynomial (cf. discussion concerning the Schloegl system in 1.). The same holds for the third case, even if another x2 term (still realizable by four reactions) would be added in the equation. The first case is realizable by three reactions. However, it is easy to see that a linear term in the steady state equation needs at least two additional reactions.

  6. 6.

    However, the second system can yield a bistable at most bimolecular system with only 4 reactions: The system is realizable by three reactions, we add y as the positive term in (still three reactions) and a fourth reaction, a simple efflux of x to ensure a linear term in the cubic polynomial (to yield three steady states): . Analysis shows that the reaction from y to x needs to be bimolecular (reaction 1 in (2a)) to enable three nonnegative steady states. In fact, this is the only input into the system. Note that this system is quite unique, there is no really different bistable chemical system with only four reactions: the positive term in has to be y to enable a positive quadratic term in the steady state equation . System (2) is unique with respect to definition (1). The very similar system where reaction 2XX + Y is replaced with 2X → 2Y is mathematically equivalent (factor 2 can be incorporated into the kinetic constant). □

The slightly modified system replacing reaction X + YY + P with X + YP is bistable as well (it has also just four reactions, but is not as minimal as system (2) for the third criterion of definition (1). Modifying further by replacing 2XX + Y with 2XY and S + Y → 2X with S + Y → 3X gives another bistable system resembling a system which can be derived from the Schloegl-system using our previously discussed general transformation rules [57].


  1. Eissing T, Conzelmann H, Gilles ED, Allgoewer F, Bullinger E, Scheurich P: Bistability analyses of a caspase activation model for receptor-induced apoptosis. J Biol Chem. 2004, 279: 36892-36897. 10.1074/jbc.M404893200

    Article  CAS  PubMed  Google Scholar 

  2. Kim D, Rath O, Kolch W, Cho K-H: A hidden oncogenic positive feedback loop caused by crosstalk between Wnt and ERK pathways. Oncogene. 2007, 26: 4571-4579. 10.1038/sj.onc.1210230

    Article  CAS  PubMed  Google Scholar 

  3. Kellershohn N, Laurent M: Prion diseases: dynamics of the infection and properties of the bistable transition. Biophys J. 2001, 81: 2517-2529. 10.1016/S0006-3495(01)75897-3

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  4. Veening J-W, Smiths WK, Kuipers OP: Bistability, epigenetics, and bet-hedging in bacteria. Annual Rev Microbiol. 2008, 62: 193-210. 10.1146/annurev.micro.62.081307.163002.

    Article  CAS  Google Scholar 

  5. Yao G, Lee TJ, Mori S, Nevins JR, You L: A bistable Rb-E2F switch underlies the restriction point. Nature Cell Biol. 2008, 10: 476-482. 10.1038/ncb1711

    Article  CAS  PubMed  Google Scholar 

  6. Tyson JJ, Chen KC, Novak B: Sniffers, buzzers, toggles and blinkers: dynamics of regulatory and signalling pathways in the cell. Curr Opin Cell Biol. 2003, 15: 221-231. 10.1016/S0955-0674(03)00017-6

    Article  CAS  PubMed  Google Scholar 

  7. Ferrell JE, Machleder EM: The biochemical basis of an all-or-none cell fate switch in Xenopus oocytes. Science. 1998, 280: 895-898. 10.1126/science.280.5365.895

    Article  CAS  PubMed  Google Scholar 

  8. Xiong W, Ferrell JE: A positive-feedback-based bistable 'memory module' that governs a cell fate decision. Nature. 2003, 426: 460-465. 10.1038/nature02089

    Article  CAS  PubMed  Google Scholar 

  9. Ozbudak EM, Thattai M, Lim HN, Shraiman BI, Van Oudenaarden A: Multistability in the lactose utilization network of Escherichia coli. Nature. 2004, 427: 737-740. 10.1038/nature02298

    Article  CAS  PubMed  Google Scholar 

  10. Paliwal S, Iglesias PA, Campbell K, Hilioti Z, Groisman A, Levchenko A: MAPK-mediated bimodal gene expression and adaptive gradient sensing in yeast. Nature. 2007, 446: 46-51. 10.1038/nature05561

    Article  CAS  PubMed  Google Scholar 

  11. Miller P, Zhabotinsky AM, Lisman JE, Wang X-J: The stability of a stochastic CaMKII switch: dependence on the number of enzyme molecules and protein turnover. PLoS Biol. 2005, 3: e107- 10.1371/journal.pbio.0030107

    Article  PubMed Central  PubMed  Google Scholar 

  12. Reynolds AR, Tischer C, Verveer PJ, Rocks O, Bastiaens PIH: EGFR activation coupled to inhibition of tyrosine phosphatases causes lateral signal propagation. Nature Cell Biol. 2003, 5: 447-453. 10.1038/ncb981

    Article  CAS  PubMed  Google Scholar 

  13. Bhalla US, Ram PT, Iyengar R: MAP kinase phosphotase as a locus of flexibility in a mitogen-activated protein kinase signalling network. Science. 2002, 297: 1018-1023. 10.1126/science.1068873

    Article  CAS  PubMed  Google Scholar 

  14. Markevich NI, Hoek JB, Kholodenko BN: Signalling switches and bistability arising from multisite phosphorylation in protein kinase cascades. J Cell Biol. 2004, 164: 353-359. 10.1083/jcb.200308060

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Clarke BL: Stability of complex reaction networks. Adv Chem Phys. 1980, 43: 1-216. full_text. full_text

    CAS  Google Scholar 

  16. Thomas R: The role of feedback circuits: Positive feedbacl circuits are a necessary condition for positive real eigenvalues of the Jacobian matrix. Ber Bunsenges Phys Chem. 1994, 98: 1148-1151.

    Article  Google Scholar 

  17. Cinquin O, Demongeot J: Positive and negative feedback: Striking a balance between necessary antagonists. J Theor Biol. 2002, 216: 229-241. 10.1006/jtbi.2002.2544

    Article  CAS  PubMed  Google Scholar 

  18. Soulé C: Graphic requirements for multistationarity. ComPlexUs. 2003, 1: 123-133. 10.1159/000076100.

    Article  Google Scholar 

  19. Ferrell JE: Self-perpetuating states in signal transduction: positive feedback, double-negative feedback and bistability. Curr Opin Chem Biol. 2002, 6: 140-148.

    Article  Google Scholar 

  20. Ferrell JE, Xiong W: Bistability in cell signalling: How to make continuous processes discontinuous, and reversible processes irreversible. Chaos. 2001, 11: 227-236. 10.1063/1.1349894

    Article  CAS  PubMed  Google Scholar 

  21. Angeli D, Ferrell JE, Sontag ED: Detection of multistability, bifurcations, and hysteresis in a large class of biological positive-feedback systems. Proc Natl Acad Sci USA. 2004, 101: 1822-1827. 10.1073/pnas.0308265100

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Eissing T, Waldherr S, Allgoewer 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.003

    Article  CAS  PubMed  Google Scholar 

  23. Craciun G, Tang Y, Feinberg M: Understanding bistability in complex enzyme-driven reaction networks. Proc Natl Acad Sci USA. 2006, 103: 8697-8702. 10.1073/pnas.0602767103

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Conradi C, Flockerzi D, Raisch J: Multistationarity in the activation of a MAPK: Parametrizing the relevant region in parameter space. Math Biosci. 2008, 211: 105-131. 10.1016/j.mbs.2007.10.004

    Article  CAS  PubMed  Google Scholar 

  25. Conradi C, Flockerzi D, Raisch J, Stelling J: Subnetwork analysis reveals dynamic features of complex (bio)chemical networks. Proc Natl Acad Sci USA. 2007, 104: 19175-19180. 10.1073/pnas.0705731104

    Article  PubMed Central  PubMed  Google Scholar 

  26. Schuster S, von Kamp A, Pachkov M: Understanding the roadmap of metabolism by pathway analysis. Methods Mol Biol. 2007, 358: 199-226. full_text

    Article  CAS  PubMed  Google Scholar 

  27. Lotka AJ: Undamped oscillations derived from the law of mass action. J Am Chem Soc. 1920, 42: 1595-1599. 10.1021/ja01453a010.

    Article  CAS  Google Scholar 

  28. Volterra V: Fluctuations in the abundance of a species considered mathematically. Nature. 1926, 118: 558-560. 10.1038/118558a0.

    Article  Google Scholar 

  29. Higgins J: The theory of oscillating reactions. Ind Eng Chem. 1967, 59: 18-62. 10.1021/ie50689a006.

    Article  CAS  Google Scholar 

  30. Selkov EE: Self-oscillations in glycolysis. 1. A simple kinetic model. Eur J Biochem. 1968, 4: 79-86. 10.1111/j.1432-1033.1968.tb00175.x

    Article  CAS  Google Scholar 

  31. Schnakenberg J: Simple chemical reaction systems with limit cycle behaviour. J Theor Biol. 1979, 81: 389-400. 10.1016/0022-5193(79)90042-0

    Article  CAS  PubMed  Google Scholar 

  32. Prigogine I, Lefever R: Symmetry breaking instabilities in dissipative systems. II. J Chem Phys. 1968, 48: 1695-1700. 10.1063/1.1668896.

    Article  Google Scholar 

  33. Wilhelm T, Heinrich R: Smallest chemical reaction system with Hopf bifurcation. J Math Chem. 1995, 17: 1-14. 10.1007/BF01165134.

    Article  CAS  Google Scholar 

  34. Mori Y, Jilkine A, Edelstein-Keshet L: Wave-pinning and cell polarity from a bistable reaction-diffusion system. Biophys J. 2008, 94: 3684-3697. 10.1529/biophysj.107.120824

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Csercik D, Hangos KM, Nagy GM: A simple reaction kinetic model of rapid (G protein dependent) and slow (β-Arrestin dependent) transmission. J Theor Biol. 2008, 255: 119-128. 10.1016/j.jtbi.2008.07.032

    Article  Google Scholar 

  36. Shiu A: The smallest multistationary mass-preserving chemical reaction network. Lecture Notes Comp Sci. 2008, 5147: 172-184. full_text. full_text

    Article  Google Scholar 

  37. Ramakrishnan N, Bhalla US: Memory switches in chemical reaction space. PloS Comput Biol. 2008, 4: e1000122- 10.1371/journal.pcbi.1000122

    Article  PubMed Central  PubMed  Google Scholar 

  38. Schloegl F: Chemical reaction models for non-equilibrium phase transitions. Z Physik. 1972, 253: 147-161. 10.1007/BF01379769.

    Article  Google Scholar 

  39. Wilhelm T: Analysis of structures causing instabilities. Phys Rev E. 2007, 76: 011911-10.1103/PhysRevE.76.011911.

    Article  Google Scholar 

  40. Wilhelm T, Hoffmann-Klipp E, Heinrich R: An evolutionary approach to enzyme kinetics: optimization of ordered mechanisms. Bull Math Biol. 1994, 56: 65-106.

    Article  CAS  Google Scholar 

  41. Klamt S, Saez-Rodriguez J, Lindquist JA, Simeoni L, Gilles ED: A methodology for the structural and functional analysis of signalling and regulatory networks. BMC Bioinformatics. 2006, 7: 56- 10.1186/1471-2105-7-56

    Article  PubMed Central  PubMed  Google Scholar 

  42. Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita KF, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M: From genomics to chemical genomics: new developments in KEGG. Nucleic Acids Res. 2006, 34: D354-357. 10.1093/nar/gkj102

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  43. Krull M, Pistor S, Voss N, Kel A, Reuter I, Kronenberg D, Michael H, Schwarzer K, Potapov A, Choi C, Kel-Margoulis O, Wingender E: TRANSPATH: an information resource for storing and visualizing signalling pathways and their pathological aberrations. Nucleic Acids Res. 2006, 34: D546-551. 10.1093/nar/gkj107

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  44. Matys V, Kel-Margoulis OV, Fricke E, Liebich I, Land S, Barre-Dirrie A, Reuter I, Chekmenev D, Krull M, Hornischer K, Voss N, Stegmaier P, Lewicki-Potapov B, Saxel H, Kel AE, Wingender E: TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes. Nucleic Acids Res. 2006, 34: D108-110. 10.1093/nar/gkj143

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Schomburg I, Chang A, Ebeling C, Gremse M, Heldt C, Huhn G, Schomburg D: BRENDA, the enzyme database: updates and major new developments. Nucleic Acids Res. 2004, 32: D431-433. 10.1093/nar/gkh081

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  46. Laplante JP, Borckmans P, Dewel G, Gimenez M, Micheau JC: Relaxation kinetics near the hysteresis limit of a bistable chemical system: The chlorite-iodide reaction in a CSTR. J Phys Chem. 1987, 91: 3401-3405. 10.1021/j100296a060.

    Article  CAS  Google Scholar 

  47. Dutt AK: Non-equilibrium thermodynamics of a model bistable chemical system. Chem Phys Lett. 2000, 322: 73-77. 10.1016/S0009-2614(00)00283-9.

    Article  CAS  Google Scholar 

  48. Strizhak PE, Demjanchyk I, Fecher F, Schneider FW, Muenster AF: Stochastic Resonance in a bistable chemical system: The oxidation of ascorbic acid by oxygen catalyzed by Copper(II) Ions. Angew Chem Int Ed. 2000, 39: 4573-4576. 10.1002/1521-3773(20001215)39:24<4573::AID-ANIE4573>3.0.CO;2-0.

    Article  CAS  Google Scholar 

  49. Goldblatt C, Lenton TM, Watson AJ: Bistability of atmospheric oxygen and the Great Oxidation. Nature. 2006, 443: 683-686. 10.1038/nature05169

    Article  CAS  PubMed  Google Scholar 

  50. Goldbeter A: Computational approaches to cellular rhythms. Nature. 2002, 420: 238-245. 10.1038/nature01259

    Article  CAS  PubMed  Google Scholar 

  51. Tsai TY-C, Choi YS, Ma W, Pomerening JR, Tang C, Ferrel JE: Robust, tunable biological oscillations from interlinked positive and negative feedback loops. Science. 2008, 321: 126-129. 10.1126/science.1156951

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  52. Goldbeter A: Biochemical oscillations and cellular rhythms: The molecular bases of periodic and chaotic behaviour. 2008, Cambridge Univ Press, Cambridge

    Google Scholar 

  53. Kruse K, Howard M, Margolin W: An experimentalist's guide to computational modelling of the Min system. Mol Microbiol. 2007, 63: 1279-1284. 10.1111/j.1365-2958.2007.05607.x

    Article  CAS  PubMed  Google Scholar 

  54. Wilhelm T, Heinrich R: Mathematical analysis of the smallest chemical reaction system with Hopf bifurcation. J Math Chem. 1996, 19: 111-130. 10.1007/BF01165179.

    Article  CAS  Google Scholar 

  55. Wilhelm T, Schuster S, Heinrich R: Kinetic and thermodynamic analyses of the reversible version of the smallest chemical reaction system with Hopf bifurcation. Nonlinear World. 1997, 4: 295-321.

    Google Scholar 

  56. Schneider KR, Wilhelm T: Model reduction by extended quasi-steady-state approximation. J Math Biol. 2000, 40: 443-450. 10.1007/s002850000026

    Article  CAS  PubMed  Google Scholar 

  57. Wilhelm T: Chemical systems consisting only of elementary steps - a paradigma for nonlinear behaviour. J Math Chem. 2000, 27: 71-88. 10.1023/A:1019131221994.

    Article  CAS  Google Scholar 

  58. Arnold VI: Ordinary differential equations. 1980, MIT Press, Cambridge

    Google Scholar 

  59. Murray JD: Mathematical Biology. 1993, Springer, Berlin

    Book  Google Scholar 

Download references


I thank Klaus R. Schneider for referring to the lemma of the index sum. The work was supported by a BBSRC Core Strategic Grant for IFR.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Thomas Wilhelm.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

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 (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Wilhelm, T. The smallest chemical reaction system with bistability. BMC Syst Biol 3, 90 (2009).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: