A multi-time-scale analysis of chemical reaction networks: II. Stochastic systems

We consider stochastic descriptions of chemical reaction networks in which there are both fast and slow reactions, and for which the time scales are widely separated. We develop a computational algorithm that produces the generator of the full chemical master equation for arbitrary systems, and show how to obtain a reduced equation that governs the evolution on the slow time scale. This is done by applying a state space decomposition to the full equation that leads to the reduced dynamics in terms of certain projections and the invariant distributions of the fast system. The rates or propensities of the reduced system are shown to be the rates of the slow reactions conditioned on the expectations of fast steps. We also show that the generator of the reduced system is a Markov generator, and we present an efficient stochastic simulation algorithm for the slow time scale dynamics. We illustrate the numerical accuracy of the approximation by simulating several examples. Graph-theoretic techniques are used throughout to describe the structure of the reaction network and the state-space transitions accessible under the dynamics.

[1]  H. Othmer.  A GRAPH-THEORETIC ANALYSIS OF CHEMICAL REACTION NETWORKS I. Invariants, Network Equivalence and Nonexistence of Various Types of Steady States ∗ , 1981 .

[2]  Eric Vanden-Eijnden,et al.  Nested stochastic simulation algorithm for chemical kinetic systems with disparate rates. , 2005, The Journal of chemical physics.

[3]  Hans G Othmer,et al.  A multi-time-scale analysis of chemical reaction networks: I. Deterministic systems , 2010, Journal of mathematical biology.

[4]  Vladimir A. Kazeev,et al.  Direct Solution of the Chemical Master Equation Using Quantized Tensor Trains , 2014, PLoS Comput. Biol..

[5]  Daniel T Gillespie,et al.  Stochastic simulation of chemical kinetics. , 2007, Annual review of physical chemistry.

[6]  Richard J. Boucherie,et al.  Product forms for queueing networks with state-dependent multiple job transitions , 1991, Advances in Applied Probability.

[7]  Hans G. Othmer,et al.  Nonuniqueness of equilibria in closed reacting systems , 1976 .

[8]  C. Rao,et al.  Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm , 2003 .

[9]  Wilhelm Huisinga,et al.  Hybrid Stochastic-Deterministic Solution of the Chemical Master Equation , 2012, Multiscale Model. Simul..

[10]  D. A. Mcquarrie.  Stochastic approach to chemical kinetics , 1967, Journal of Applied Probability.

[11]  Eric L Haseltine,et al.  The stochastic quasi-steady-state assumption: reducing the model but not the noise. , 2011, The Journal of chemical physics.

[12]  Can Huang,et al.  Strong Convergence and Speed up of Nested Stochastic Simulation Algorithm , 2014 .

[13]  Brian Munsky,et al.  Reduction and solution of the chemical master equation using time scale separation and finite state projection. , 2006, The Journal of chemical physics.

[14]  W. Huisinga,et al.  Solving the chemical master equation for monomolecular reaction systems analytically , 2006, Journal of mathematical biology.

[15]  J. Goutsias.  Quasiequilibrium approximation of fast reaction kinetics in stochastic biochemical systems. , 2005, The Journal of chemical physics.

[16]  Tosio Kato.  Perturbation theory for linear operators , 1966 .

[17]  S. Schnell,et al.  Reaction kinetics in intracellular environments with macromolecular crowding: simulations and rate laws. , 2004, Progress in biophysics and molecular biology.

[18]  David F. Anderson,et al.  Continuous Time Markov Chain Models for Chemical Reaction Networks , 2011 .

[19]  Linda R Petzold,et al.  The slow-scale stochastic simulation algorithm. , 2005, The Journal of chemical physics.

[20]  J. Goutsias,et al.  Markovian dynamics on complex reaction networks , 2012, 1205.5524.

[21]  David F. Anderson,et al.  Product-Form Stationary Distributions for Deficiency Zero Chemical Reaction Networks , 2008, Bulletin of mathematical biology.

[22]  João P. Hespanha,et al.  Equilibrium distributions of simple biochemical reaction systems for time-scale separation in stochastic reaction networks , 2014, Journal of The Royal Society Interface.

[23]  Eric L Haseltine,et al.  Two classes of quasi-steady-state model reductions for stochastic kinetics. , 2007, The Journal of chemical physics.

[24]  H. Othmer,et al.  The effects of cell density and metabolite flux on cellular dynamics , 1978, Journal of mathematical biology.

[25]  R. Bundschuh,et al.  Fluctuations and slow variables in genetic networks. , 2003, Biophysical Journal.

[26]  H. Othmer,et al.  A stochastic analysis of first-order reaction networks , 2005, Bulletin of mathematical biology.

[27]  Yiannis Kaznessis,et al.  Accurate hybrid stochastic simulation of a system of coupled chemical or biochemical reactions. , 2005, The Journal of chemical physics.

[28]  J. A. M. Janssen,et al.  The elimination of fast variables in complex chemical reactions. III. Mesoscopic level (irreducible case) , 1989 .

[29]  L. R. Petzold,et al.  Nested stochastic simulation algorithm for chemical kinetic systems with disparate rates ” [ JCP , 2007 .

[30]  Chang Hyeong Lee,et al.  A reduction method for multiple time scale stochastic reaction networks , 2009 .

[31]  R. Jackson,et al.  General mass action kinetics , 1972 .

[32]  Philipp Thomas,et al.  Communication: limitations of the stochastic quasi-steady-state approximation in open biochemical reaction networks. , 2011, The Journal of chemical physics.

[33]  T. B. Boffey,et al.  Applied Graph Theory , 1973 .

[34]  R. Aris.  Prolegomena to the rational analysis of systems of chemical reactions , 1965 .

[35]  H. El-Samad,et al.  A rigorous framework for multiscale simulation of stochastic cellular networks. , 2009, The Journal of chemical physics.

[36]  Peter Deuflhard,et al.  Adaptive Discrete Galerkin Methods Applied to the Chemical Master Equation , 2008, SIAM J. Sci. Comput..

[37]  Gianluca Setti,et al.  Design and Analysis of Biomolecular Circuits , 2011 .

[38]  Jifeng Hu,et al.  Stochastic Analysis of Reaction–Diffusion Processes , 2013, Bulletin of Mathematical Biology.

[39]  J. Rawlings,et al.  Approximate simulation of coupled fast and slow reactions for stochastic chemical kinetics , 2002 .

[40]  Krešimir Josić,et al.  The validity of quasi-steady-state approximations in discrete stochastic simulations. , 2014, Biophysical journal.

[41]  C. D. Meyer,et al.  Generalized inverses of linear transformations , 1979 .

[42]  Thomas Wilhelm,et al.  The smallest chemical reaction system with bistability , 2009, BMC Systems Biology.

[43]  Ramon Grima,et al.  Model reduction for stochastic chemical systems with abundant species. , 2015, The Journal of chemical physics.

[44]  Andreas Hellander,et al.  Hybrid method for the chemical master equation , 2007, J. Comput. Phys..

[45]  Simon L. Cotter,et al.  Constrained approximation of effective generators for multiscale stochastic reaction networks and application to conditioned path sampling , 2015, J. Comput. Phys..