Bistability and oscillations in chemical reaction networks.

*(English)*Zbl 1311.92088Summary: Bifurcation theory is one of the most widely used approaches for analysis of dynamical behaviour of chemical and biochemical reaction networks. Some of the interesting qualitative behaviour that are analyzed are oscillations and bistability (a situation where a system has at least two coexisting stable equilibria). Both phenomena have been identified as central features of many biological and biochemical systems. This paper, using the theory of stoichiometric network analysis (SNA) and notions from algebraic geometry, presents sufficient conditions for a reaction network to display bifurcations associated with these phenomena. The advantage of these conditions is that they impose fewer algebraic conditions on model parameters than conditions associated with standard bifurcation theorems. To derive the new conditions, a coordinate transformation will be made that will guarantee the existence of branches of positive equilibria in the system. This is particularly useful in mathematical biology, where only positive variable values are considered to be meaningful. The first part of the paper will be an extended introduction to SNA and algebraic geometry-related methods which are used in the coordinate transformation and set up of the theorems. In the second part of the paper we will focus on the derivation of bifurcation conditions using SNA and algebraic geometry. Conditions will be derived for three bifurcations: the saddle-node bifurcation, a simple branching point, both linked to bistability, and a simple Hopf bifurcation. The latter is linked to oscillatory behaviour. The conditions derived are sufficient and they extend earlier results from stoichiometric network analysis as can be found in [B. D. Aguda and B. L. Clarke, “Bistability in chemical reaction networks: theory and application to the peroxidase-oxidase reaction”, J. Chem. Phys. 87, No. 6, 3461–3470 (1987); B. L. Clarke and W. Jiang, “Method for deriving Hopf and saddle-node bifurcation hypersurfaces and application to a model of the Belousov-Zhabotinskii reaction”, J. Chem. Phys. 99, No. 6, 4464–4476 (1993); K. Gatermann et al., J. Symb. Comput. 40, No. 6, 1361–1382 (2005; Zbl 1120.13033)]. In these papers, some necessary conditions for two of these bifurcations were given. A set of examples will illustrate that algebraic conditions arising from given sufficient bifurcation conditions are not more difficult to interpret nor harder to calculate than those arising from necessary bifurcation conditions. Hence an increasing amount of information is gained at no extra computational cost. The theory can also be used in a second step for a systematic bifurcation analysis of larger reaction networks.

##### MSC:

92C45 | Kinetics in biochemical problems (pharmacokinetics, enzyme kinetics, etc.) |

92E20 | Classical flows, reactions, etc. in chemistry |

34C23 | Bifurcation theory for ordinary differential equations |

37G10 | Bifurcations of singular points in dynamical systems |

37N25 | Dynamical systems in biology |

PDF
BibTeX
Cite

\textit{M. Domijan} and \textit{M. Kirkilionis}, J. Math. Biol. 59, No. 4, 467--501 (2009; Zbl 1311.92088)

Full Text:
DOI

##### References:

[1] | Aguda BD, Clarke BL (1987) Bistability in chemical reaction networks: theory and application to the peroxidase–oxidase reaction. J Chem Phys 87: 3461–3470 |

[2] | Arnold EA (2003) Modular algorithms for computing Gröbner bases. J Symb Comput 35: 403–419 · Zbl 1046.13018 |

[3] | Bayer D, Stillman M (1987) A theorem on refining division orders by the reverse lexicographic order. Duke J Math 55: 321–328 · Zbl 0638.06003 |

[4] | Bosma W, Cannon J, Playoust C (1997) The Magma algebra system. I. The user language. J Symb Comput 24(3–4): 235–265 · Zbl 0898.68039 |

[5] | Bronnikova TV, Fed’kina VR, Schaffer WM, Olsen LF (1995) Period-doubling bifurcations in a detailed model of the peroxidase–oxidase reaction. J Phys Chem 99: 9309–9312 |

[6] | Capani A, Niesi G, Robbiano L (2000) CoCoA: A System for Doing Computations in Commutative Algebra, 4.0 edn, available via anonymous ftp from http://cocoa.dima.unige.it/ |

[7] | Clarke BL (1975) Theorems on chemical network stability. J Chem Phys 62: 773–775 |

[8] | Clarke BL (1975) Stability of topologically similar chemical networks. J Chem Phys 62: 3726–3738 |

[9] | Clarke BL (1980) Stability of complex reaction networks In: Prigogine I, Rice S (eds) Adv Chem Phys Wiley, New York, vol 43, pp 1–216 |

[10] | Clarke BL, Jiang W (1993) Method for deriving Hopf and saddle-node bifurcation hypersurfaces and application to a model of the Belousov–Zhabotinskii reaction. J Chem Phys 99: 4464–4476 |

[11] | Cox D, Little J, O’Shea D (2007) Ideals, varieties, and algorithms, an introduction to computational algebraic geometry and commutative algebra, 3rd edn. Undergraduate Texts in Mathematics. Springer, New York · Zbl 1118.13001 |

[12] | Craciun G, Feinberg M (2005) Multiple equilibria in complex chemical reaction networks. I. The injectivity property. SIAM J Appl Math 65: 1526–1546 · Zbl 1094.80005 |

[13] | Craciun G, Tang Y, Feinberg M (2006) Understanding bistability in complex enzyme-driven reaction networks. PNAS 30(103): 8697–8702 · Zbl 1254.93116 |

[14] | Domijan M, Kirkilionis M (2008) Graph theory and qualitative analysis of reaction networks. Warwick Preprint 13/2007. Accepted: Networks and Heterogeneous Media · Zbl 1141.92021 |

[15] | Eiswirth M, Freund A, Ross J (1991) Mechanistic classification of chemical oscillators and the role of species. Adv Chem Phys 80: 127–199 |

[16] | Eiswirth M, Bürger J, Strasser P, Ertl G (1996) Oscillating Langmuir–Hinshelwood mechanisms. J Phys Chem 100: 19118–19123 |

[17] | Field RJ, Körös E, Noyes RM (1972) Oscillations in chemical systems. 2. Thorough analysis of temporal oscillation in bromate–cerium–malonic acid system. J Am Chem Soc 94(25): 8649– 8664 |

[18] | Field RJ, Noyes RM (1974) Oscillations in chemical systems. IV. Limit cycle behavior in a model of a real chemical reaction. J Chem Phys 60: 1877–1884 |

[19] | FluxAnalyzer (2002) An interactive program in Matlab. Max-Planck Insitute Magdeburg. Available on request |

[20] | Gantmacher FR (1959) Application of the theory of matrices. Interscience Publishers Ltd., London · Zbl 0085.01001 |

[21] | Gatermann K, Eiswirth M, Sensse A (2005) Toric ideals and graph theory to analyze Hopf bifurcations in mass-action systems. J Symb Comput 40: 1361–1382 · Zbl 1120.13033 |

[22] | Goldbeter A, Dupont G (1990) Allosteric regulation, cooperativity, and biochemical oscillations. Biophys Chem 37: 341–353 |

[23] | Grayson DR, Stillman ME (1993) Macaulay 2, a software system for research in algebraic geometry. Available at http://www.math.uiuc.edu/Macaulay2/ |

[24] | Greuel G-M, Pfister G, Schönemann H (2007) Singular 3-0-4, A computer Algebra System for Polynomial Computations. Centre for Computer Algebra, University of Kaiserslautern. http://www.singular.uni-kl.de/index.html · Zbl 1344.13002 |

[25] | Guckenheimer J, Holmes JP (2002) Nonlinear oscillations, dynamical systems and bifurcations of vector fields. Applied Mathematics Sciences 42. Springer, Berlin |

[26] | Hashemi A, Pfister G, Schönemann H (2007) Singular library, modst.lib, A Singular library for computing Gröbner bases by modular methods |

[27] | Heinrich R, Schuster S (1996) The regulation of cellular processes. Chapman & Hall, London · Zbl 0895.92013 |

[28] | Hunt KLC, Hunt PM, Ross J (1990) Nonlinear dynamics and thermodynamics of chemical reactions far from equilibrium. Annu Rev Phys Chem 41: 409–439 |

[29] | Klonowski W (1983) Simplifying principles for chemical and enzyme reaction kinetics. Biophys Chem 18: 73–87 |

[30] | Krischer K, Eiswirth M, Ertl G (1992) Oscillatory CO oxidation on Pt(110): modeling of temporal self-organisation. J Chem Phys 96: 9161–9172 |

[31] | Kuznetsov Y (1998) Elements of applied bifurcation theory. Applied mathematical sciences, 2nd edn. Springer, New York, p 112 · Zbl 0914.58025 |

[32] | Liu W (1994) Criterion of Hopf bifurcations without using eigenvalues. J Math Anal Appl 182: 250–256 · Zbl 0794.34033 |

[33] | Mincheva M, Roussel MR (2007) Graph-theoretic methods for the analysis of chemical and biochemical networks. I. Multistability and oscillations in ordinary differential equation models. J Math Biol 55:61–86 · Zbl 1149.34031 |

[34] | Perelson AS, Wallwork D (1977) The arbitrary dynamic behavior of open chemical reaction systems. J Chem Phys 66: 4390–4394 |

[35] | Rockafeller RT (1970) Convex analysis. Princeton University Press, Princeton |

[36] | Selkov EE (1968) Self-oscillations in glycolysis. 1. A simple kinetic model. Eur J Biochem 4: 79–86 |

[37] | Sensse A (2002) Algebraic methods for the analysis of Hopf bifurcations in biochemical networks. Diploma thesis, Humboldt Universitat Berlin |

[38] | Sensse A (2005) Convex and toric geometry to analyze complex dynamics in chemical reaction systems. PhD thesis, Otto-von-Guericke-Universitat Magdeburg · Zbl 1120.13033 |

[39] | Sensse A, Eiswirth M (2005) Feedback loops for chaos in activator–inhibitor systems. J Chem Phys 122: 044516 · Zbl 1120.13033 |

[40] | Sensse A, Gatermann K, Eiswirth M (2005) Analytic solution for the electrocatalytic oxidation of formic acid. J Electroanal Chem 577: 35–46 |

[41] | Slepchenko BM, Terasaki M (2003) Cyclin aggregation and robustness of bio-switching. Mol Biol Cell 14: 4695–4706 |

[42] | Theis C (2000) Singular library, toric.lib, A Singular 3.0 library for computing standard basis of toric ideals |

[43] | Trott M (2006) The mathematica guidebook for symbolics. Springer, New York, pp 32–50 · Zbl 1095.68001 |

[44] | Tyson JJ, Chen K, Novak B (2001) Network dynamics and cell physiology. Nat Rev Mol Cell Biol 2: 908–916 |

[45] | Volpert A, Hudyaev S (1985) Analysis in classes of discontinuous functions and equations of mathematical physics (Chapter 12). Martinus Nijhoff, Dordrecht |

[46] | Walter W (1998) Ordinary differential equations. Springer, New York · Zbl 0991.34001 |

This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.