The relationship between stochastic and deterministic quasi-steady state approximations

BackgroundThe quasi steady-state approximation (QSSA) is frequently used to reduce deterministic models of biochemical networks. The resulting equations provide a simplified description of the network in terms of non-elementary reaction functions (e.g. Hill functions). Such deterministic reductions are frequently a basis for heuristic stochastic models in which non-elementary reaction functions are used to define reaction propensities. Despite their popularity, it remains unclear when such stochastic reductions are valid. It is frequently assumed that the stochastic reduction can be trusted whenever its deterministic counterpart is accurate. However, a number of recent examples show that this is not necessarily the case.ResultsHere we explain the origin of these discrepancies, and demonstrate a clear relationship between the accuracy of the deterministic and the stochastic QSSA for examples widely used in biological systems. With an analysis of a two-state promoter model, and numerical simulations for a variety of other models, we find that the stochastic QSSA is accurate whenever its deterministic counterpart provides an accurate approximation over a range of initial conditions which cover the likely fluctuations from the quasi steady-state (QSS). We conjecture that this relationship provides a simple and computationally inexpensive way to test the accuracy of reduced stochastic models using deterministic simulations.ConclusionsThe stochastic QSSA is one of the most popular multi-scale stochastic simulation methods. While the use of QSSA, and the resulting non-elementary functions has been justified in the deterministic case, it is not clear when their stochastic counterparts are accurate. In this study, we show how the accuracy of the stochastic QSSA can be tested using their deterministic counterparts providing a concrete method to test when non-elementary rate functions can be used in stochastic simulations.

[1]  K. Burrage,et al.  Stochastic chemical kinetics and the total quasi-steady-state assumption: application to the stochastic simulation algorithm and chemical master equation. , 2008, The Journal of chemical physics.

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

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

[4]  K. Burrage,et al.  Stochastic models for regulatory networks of the genetic toggle switch. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[5]  Mark A. Stalzer,et al.  Efficient Formulations for Exact Stochastic Simulation of Chemical Systems , 2011, IEEE/ACM Transactions on Computational Biology and Bioinformatics.

[6]  Eric L Haseltine,et al.  On the origins of approximations for stochastic chemical kinetics. , 2005, The Journal of chemical physics.

[7]  Animesh Agarwal,et al.  On the precision of quasi steady state assumptions in stochastic dynamics. , 2012, The Journal of chemical physics.

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

[9]  Philipp Thomas,et al.  The slow-scale linear noise approximation: an accurate, reduced stochastic description of biochemical networks under timescale separation conditions , 2012, BMC Systems Biology.

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

[11]  Yiannis N Kaznessis,et al.  An equation-free probabilistic steady-state approximation: dynamic application to the stochastic simulation of biochemical reaction networks. , 2005, The Journal of chemical physics.

[12]  R. Erban,et al.  Stochastic modelling of reaction–diffusion processes: algorithms for bimolecular reactions , 2009, Physical biology.

[13]  Kevin R. Sanft,et al.  Legitimacy of the stochastic Michaelis-Menten approximation. , 2011, IET systems biology.

[14]  Michael A Schwemmer,et al.  Effects of moderate noise on a limit cycle oscillator: counterrotation and bistability. , 2014, Physical review letters.

[15]  Achim Kramer,et al.  Global parameter search reveals design principles of the mammalian circadian clock , 2008, BMC Systems Biology.

[16]  Cellular Physiology , 2022 .

[17]  Ovidiu Radulescu,et al.  Hybrid stochastic simplifications for multiscale gene networks , 2009, BMC Systems Biology.

[18]  T. Kepler,et al.  Stochasticity in transcriptional regulation: origins, consequences, and mathematical representations. , 2001, Biophysical journal.

[19]  B. Ingalls,et al.  Estimations of intrinsic and extrinsic noise in models of nonlinear genetic networks. , 2006, Chaos.

[20]  O Wolkenhauer,et al.  Family tree of Markov models in systems biology. , 2006, IET systems biology.

[21]  Andrzej M. Kierzek,et al.  STOCKS: STOChastic Kinetic Simulations of biochemical systems with Gillespie algorithm , 2002, Bioinform..

[22]  L. A. Segel,et al.  The Quasi-Steady-State Assumption: A Case Study in Perturbation , 1989, SIAM Rev..

[23]  T. Wager,et al.  Modeling and Validating Chronic Pharmacological Manipulation of Circadian Rhythms , 2013, CPT: pharmacometrics & systems pharmacology.

[24]  C. Daub,et al.  BMC Systems Biology , 2007 .

[25]  Yang Cao,et al.  Stochastic simulation of enzyme-catalyzed reactions with disparate timescales. , 2008, Biophysical journal.

[26]  Mudita Singhal,et al.  Simulation of Biochemical Networks using Copasi - A Complex Pathway Simulator , 2006, Proceedings of the 2006 Winter Simulation Conference.

[27]  José Halloy,et al.  How molecular should your molecular model be? On the level of molecular detail required to simulate biological networks in systems and synthetic biology. , 2011, Methods in enzymology.

[28]  Andrew J Black,et al.  Stochastic amplification in an epidemic model with seasonal forcing. , 2010, Journal of theoretical biology.

[29]  Linda R Petzold,et al.  Velocity Response Curves Support the Role of Continuous Entrainment in Circadian Clocks , 2010, Journal of biological rhythms.

[30]  Ovidiu Radulescu,et al.  Convergence of stochastic gene networks to hybrid piecewise deterministic processes , 2011, 1101.1431.

[31]  Mudita Singhal,et al.  COPASI - a COmplex PAthway SImulator , 2006, Bioinform..

[32]  John J. Tyson,et al.  Modeling Networks of Coupled Enzymatic Reactions Using the Total Quasi-Steady State Approximation , 2007, PLoS Comput. Biol..

[33]  Jae Kyoung Kim,et al.  On the Existence and Uniqueness of Biological Clock Models Matching Experimental Data , 2012, SIAM J. Appl. Math..

[34]  Sebastian Walcher,et al.  A constructive approach to quasi-steady state reductions , 2014, Journal of Mathematical Chemistry.

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

[36]  Robert E. O'Malley,et al.  Analyzing Multiscale Phenomena Using Singular Perturbation Methods , 1999 .

[37]  Bruce Tidor,et al.  Combined Model of Intrinsic and Extrinsic Variability for Computational Network Design with Application to Synthetic Biology , 2013, PLoS Comput. Biol..

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

[39]  M. Bennett,et al.  Transient dynamics of genetic regulatory networks. , 2007, Biophysical journal.

[40]  J Halloy,et al.  Deterministic Versus Stochastic Models for Circadian Rhythms , 2002, Journal of biological physics.

[41]  J. Collins,et al.  Combinatorial promoter design for engineering noisy gene expression , 2007, Proceedings of the National Academy of Sciences.

[42]  Philipp Thomas,et al.  Rigorous elimination of fast stochastic variables from the linear noise approximation using projection operators. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[43]  Daniel B. Forger,et al.  A mechanism for robust circadian timekeeping via stoichiometric balance , 2012, Molecular systems biology.

[44]  N. Berglund,et al.  Geometric singular perturbation theory for stochastic differential equations , 2002, math/0204008.

[45]  D. Gillespie.  Exact Stochastic Simulation of Coupled Chemical Reactions , 1977 .

[46]  A. van Oudenaarden,et al.  Noise Propagation in Gene Networks , 2005, Science.

[47]  J. Elf,et al.  Fast evaluation of fluctuations in biochemical networks with the linear noise approximation. , 2003, Genome research.

[48]  Vahid Shahrezaei,et al.  Analytical distributions for stochastic gene expression , 2008, Proceedings of the National Academy of Sciences.

[49]  Jae Kyoung Kim,et al.  Mechanisms That Enhance Sustainability of p53 Pulses , 2013, PloS one.

[50]  Eshel Ben-Jacob,et al.  Turning Oscillations Into Opportunities: Lessons from a Bacterial Decision Gate , 2013, Scientific Reports.

[51]  Melanie I. Stefan,et al.  Molecules for memory: modelling CaMKII , 2007, BMC Systems Biology.

[52]  M. L. Simpson,et al.  Frequency domain analysis of noise in autoregulated gene circuits , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[53]  A. Tzafriri,et al.  Michaelis-Menten kinetics at high enzyme concentrations , 2003, Bulletin of mathematical biology.

[54]  Krešimir Josić,et al.  Reduced models of networks of coupled enzymatic reactions. , 2011, Journal of theoretical biology.

[55]  J. Tyson,et al.  Computational Cell Biology , 2010 .

[56]  Hye-Won Kang,et al.  Central limit theorems and diffusion approximations for multiscale Markov chain models , 2012, 1208.3783.

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

[58]  Andrea Riba,et al.  A Combination of Transcriptional and MicroRNA Regulation Improves the Stability of the Relative Concentrations of Target Genes , 2013, PLoS Comput. Biol..

[59]  Samuel A. Isaacson,et al.  Incorporating Diffusion in Complex Geometries into Stochastic Chemical Kinetics Simulations , 2006, SIAM J. Sci. Comput..

[60]  T. Kurtz,et al.  Submitted to the Annals of Applied Probability ASYMPTOTIC ANALYSIS OF MULTISCALE APPROXIMATIONS TO REACTION NETWORKS , 2022 .

[61]  Michael B. Elowitz,et al.  Architecture-Dependent Noise Discriminates Functionally Analogous Differentiation Circuits , 2009, Cell.

[62]  M. Thattai,et al.  Intrinsic noise in gene regulatory networks , 2001, Proceedings of the National Academy of Sciences of the United States of America.

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

[64]  Hans G Othmer,et al.  A multi-time-scale analysis of chemical reaction networks: II. Stochastic systems , 2015, Journal of Mathematical Biology.

[65]  L Glass,et al.  Discontinuities in phase-resetting experiments. , 1984, The American journal of physiology.

[66]  M. Kaufman,et al.  From structure to dynamics: frequency tuning in the p53-Mdm2 network. II Differential and stochastic approaches. , 2010, Journal of theoretical biology.