Derivation of a quantitative minimal model from a detailed elementary-step mechanism supported by mathematical coupling analysis.

Accurate experimental data increasingly allow the development of detailed elementary-step mechanisms for complex chemical and biochemical reaction systems. Model reduction techniques are widely applied to obtain representations in lower-dimensional phase space which are more suitable for mathematical analysis, efficient numerical simulation, and model-based control tasks. Here, we exploit a recently implemented numerical algorithm for error-controlled computation of the minimum dimension required for a still accurate reduced mechanism based on automatic time scale decomposition and relaxation of fast modes. We determine species contributions to the active (slow) dynamical modes of the reaction system and exploit this information in combination with quasi-steady-state and partial-equilibrium approximations for explicit model reduction of a novel detailed chemical mechanism for the Ru-catalyzed light-sensitive Belousov-Zhabotinsky reaction. The existence of a minimum dimension of seven is demonstrated to be mandatory for the reduced model to show good quantitative consistency with the full model in numerical simulations. We derive such a maximally reduced seven-variable model from the detailed elementary-step mechanism and demonstrate that it reproduces quantitatively accurately the dynamical features of the full model within a given accuracy tolerance.

[1]  Kenneth Showalter,et al.  Reaction Mechanism for Light Sensitivity of the Ru(bpy)32+-Catalyzed Belousov−Zhabotinsky Reaction , 1997 .

[2]  Alexander N. Gorban,et al.  Reduced description in the reaction kinetics , 2000 .

[3]  Gene H. Golub,et al.  Matrix computations , 1983 .

[4]  S. Lam,et al.  The CSP method for simplifying kinetics , 1994 .

[5]  M. Diehl,et al.  Real-time optimization and nonlinear model predictive control of processes governed by differential-algebraic equations , 2000 .

[6]  Michael J. Davis,et al.  Geometrical Simplification of Complex Kinetic Systems , 2001 .

[7]  Simon J. Fraser,et al.  Geometrical picture of reaction in enzyme kinetics , 1989 .

[8]  Marc R. Roussel,et al.  Accurate steady-state approximations : implications for kinetics experiments and mechanism , 1991 .

[9]  L. Petzold,et al.  Model reduction for chemical kinetics: an optimization approach , 1999 .

[10]  Alison S. Tomlin,et al.  The use of algebraic sets in the approximation of inertial manifolds and lumping in chemical kinetic systems , 1995 .

[11]  A. N. Gorban,et al.  Constructive methods of invariant manifolds for kinetic problems , 2003 .

[12]  A. Zhabotinsky,et al.  Concentration Wave Propagation in Two-dimensional Liquid-phase Self-oscillating System , 1970, Nature.

[13]  M. Mavrovouniotis,et al.  Simplification of Mathematical Models of Chemical Reaction Systems. , 1998, Chemical reviews.

[14]  T. Amemiya,et al.  Experimental and Model Studies of Oscillations, Photoinduced Transitions, and Steady States in the Ru(bpy)32+-Catalyzed Belousov−Zhabotinsky Reaction under Different Solute Compositions , 2002 .

[15]  L. Kuhnert,et al.  A new optical photochemical memory device in a light-sensitive chemical active medium , 1986, Nature.

[16]  Simon J. Fraser,et al.  The steady state and equilibrium approximations: A geometrical picture , 1988 .

[17]  R. H. Snow,et al.  A Chemical Kinetics Computer Program for Homogeneous and Free-Radical Systems of Reactions , 1966 .

[18]  Tamás Turányi,et al.  On the error of the quasi-steady-state approximation , 1993 .

[19]  J. L. Hudson,et al.  Chemical complexity: Spontaneous and engineered structures , 2003 .

[20]  Daniel Walgraef,et al.  Spatio-temporal pattern formation , 1996 .

[21]  J. Demas,et al.  An oscillating chemical reaction with a luminescent indicator , 1973 .

[22]  M. Jinguji,et al.  Primary process of illumination effect on the tris(bipyridine)ruthenium(2+)-catalyzed Belousov-Zhabotinskii reaction , 1992 .

[23]  Ursula Kummer,et al.  A New Time-Dependent Complexity Reduction Method for Biochemical Systems , 2005, Trans. Comp. Sys. Biology.

[24]  T. Turányi.  Sensitivity analysis of complex kinetic systems. Tools and applications , 1990 .

[25]  Grégoire Nicolis,et al.  Self-Organization in nonequilibrium systems , 1977 .

[26]  K. Chou,et al.  Graphic rules in steady and non-steady state enzyme kinetics. , 1989, The Journal of biological chemistry.

[27]  Michael Menzinger,et al.  Nonlinear Dynamics of the BZ Reaction: A Simple Experiment that Illustrates Limit Cycles, Chaos, Bifurcations, and Noise , 1996 .

[28]  E. L. King,et al.  A Schematic Method of Deriving the Rate Laws for Enzyme-Catalyzed Reactions , 1956 .

[29]  Herschel Rabitz,et al.  New approaches to determination of constrained lumping schemes for a reaction system in the whole composition space , 1991 .

[30]  Ulrich Maas,et al.  Simplifying chemical kinetics: Intrinsic low-dimensional manifolds in composition space , 1992 .

[31]  P. Deuflhard,et al.  Extrapolation Integrators for Quasilinear Implicit ODEs , 1987 .

[32]  István Szalai,et al.  Oscillations in the Bromomalonic Acid/Bromate System Catalyzed by [Ru(bipy)3]2+ , 1995 .

[33]  H. Kitano,et al.  Computational systems biology , 2002, Nature.

[34]  János Tóth,et al.  A general analysis of exact nonlinear lumping in chemical kinetics , 1994 .

[35]  L. Kuhnert,et al.  Analysis of the modified complete Oregonator accounting for oxygen sensitivity and photosensitivity of Belousov-Zhabotinskii systems , 1990 .

[36]  D Lebiedz,et al.  Computing minimal entropy production trajectories: an approach to model reduction in chemical kinetics. , 2004, The Journal of chemical physics.

[37]  Marc R. Roussel,et al.  Geometry of the steady-state approximation: Perturbation and accelerated convergence methods , 1990 .

[38]  D J Park,et al.  The hierarchical structure of metabolic networks and the construction of efficient metabolic simulators. , 1974, Journal of theoretical biology.

[39]  Dimitris A. Goussis,et al.  Asymptotic Solution of Stiff PDEs with the CSP Method: The Reaction Diffusion Equation , 1998, SIAM J. Sci. Comput..

[40]  M. Kramer,et al.  Sensitivity Analysis in Chemical Kinetics , 1983 .

[41]  U. Maas.  Efficient calculation of intrinsic low-dimensional manifolds for the simplification of chemical kinetics , 1998 .

[42]  Kenneth Showalter,et al.  Design and Control of Wave Propagation Patterns in Excitable Media , 2002, Science.