Reduction of multiscale stochastic biochemical reaction networks using exact moment derivation

Biochemical reaction networks (BRNs) in a cell frequently consist of reactions with disparate timescales. The stochastic simulations of such multiscale BRNs are prohibitively slow due to high computational cost for the simulations of fast reactions. One way to resolve this problem uses the fact that fast species regulated by fast reactions quickly equilibrate to their stationary distribution while slow species are unlikely to be changed. Thus, on a slow timescale, fast species can be replaced by their quasi-steady state (QSS): their stationary conditional expectation values for given slow species. As the QSS are determined solely by the state of slow species, such replacement leads to a reduced model, where fast species are eliminated. However, it is challenging to derive the QSS in the presence of nonlinear reactions. While various approximation schemes for the QSS have been developed, they often lead to considerable errors. Here, we propose two classes of multiscale BRNs which can be reduced by deriving an exact QSS rather than approximations. Specifically, if fast species constitute either a feedforward network or a complex balanced network, the reduced model based on the exact QSS can be derived. Such BRNs are frequently observed in a cell as the feedforward network is one of fundamental motifs of gene or protein regulatory networks. Furthermore, complex balanced networks also include various types of fast reversible bindings such as bindings between transcriptional factors and gene regulatory sites. The reduced models based on exact QSS, which can be calculated by the computational packages provided in this work, accurately approximate the slow scale dynamics of the original full model with much lower computational cost.

[1]  Matthew West,et al.  Binning for Efficient Stochastic Multiscale Particle Simulations , 2013, Multiscale Model. Simul..

[2]  Sookkyung Lim,et al.  Mathematical modeling and validation of glucose compensation of the neurospora circadian clock. , 2015, Biophysical journal.

[3]  Niraj Kumar,et al.  Exact distributions for stochastic gene expression models with bursting and feedback. , 2014, Physical review letters.

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

[5]  U. Alon,et al.  Search for computational modules in the C. elegans brain , 2004, BMC Biology.

[6]  Deena R. Schmidt,et al.  Steady-state fluctuations of a genetic feedback loop: an exact solution. , 2012, The Journal of chemical physics.

[7]  Jaques Reifman,et al.  Mathematical Modeling of the Heat-Shock Response in HeLa Cells , 2015, Biophysical journal.

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

[9]  D. Gillespie The chemical Langevin equation , 2000 .

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

[11]  Jae Kyoung Kim,et al.  Reduction for Stochastic Biochemical Reaction Networks with Multiscale Conservations , 2017, Multiscale Model. Simul..

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

[13]  Eduardo Sontag,et al.  Exact Moment Dynamics for Feedforward Nonlinear Chemical Reaction Networks , 2015, IEEE Life Sciences Letters.

[14]  D. Gillespie,et al.  Deterministic limit of stochastic chemical kinetics. , 2009, The journal of physical chemistry. B.

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

[16]  Doron Zeilberger,et al.  A symbolic computation approach to a problem involving multivariate Poisson distributions , 2010, Adv. Appl. Math..

[17]  Ioannis G Kevrekidis,et al.  A constrained approach to multiscale stochastic simulation of chemically reacting systems. , 2011, The Journal of chemical physics.

[18]  Krešimir Josić,et al.  Molecular mechanisms that regulate the coupled period of the mammalian circadian clock. , 2014, Biophysical journal.

[19]  Peter G Wolynes,et al.  Influence of decoys on the noise and dynamics of gene expression. , 2012, Physical review. E, Statistical, nonlinear, and soft matter physics.

[20]  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.

[21]  S. Meyn,et al.  Stability of Markovian processes III: Foster–Lyapunov criteria for continuous-time processes , 1993, Advances in Applied Probability.

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

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

[24]  Richard Bertram,et al.  Geometric Singular Perturbation Analysis of Bursting Oscillations in Pituitary Cells , 2015 .

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

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

[27]  Xiufen Zou,et al.  Negative feedback contributes to the stochastic expression of the interferon-β gene in virus-triggered type I interferon signaling pathways. , 2015, Mathematical biosciences.

[28]  E. O’Shea,et al.  Global analysis of protein expression in yeast , 2003, Nature.

[29]  Gilles Gnacadja,et al.  Monotonicity of interleukin-1 receptor-ligand binding with respect to antagonist in the presence of decoy receptor. , 2007, Journal of theoretical biology.

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

[31]  Linda R. Petzold,et al.  Stochastic modelling of gene regulatory networks , 2005 .

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

[33]  Per Lötstedt,et al.  Dimensional Reduction of the Fokker-Planck Equation for Stochastic Chemical Reactions , 2006, Multiscale Model. Simul..

[34]  Jae Kyoung Kim,et al.  The relationship between stochastic and deterministic quasi-steady state approximations , 2015, BMC Systems Biology.

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

[36]  M. Feinberg Chemical reaction network structure and the stability of complex isothermal reactors—I. The deficiency zero and deficiency one theorems , 1987 .

[37]  Elijah Roberts,et al.  Approximation and inference methods for stochastic biochemical kinetics—a tutorial review , 2017 .

[38]  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.

[39]  Ying Liu,et al.  Stochastic resonance in feedforward-loop neuronal network motifs in astrocyte field. , 2013, Journal of theoretical biology.

[40]  Domitilla Del Vecchio,et al.  Model order reduction for Linear Noise Approximation using time-scale separation , 2016, 2016 IEEE 55th Conference on Decision and Control (CDC).

[41]  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.

[42]  Jae Kyoung Kim Protein sequestration versus Hill-type repression in circadian clock models. , 2016, IET systems biology.

[43]  T. Hwa,et al.  Stochastic fluctuations in metabolic pathways , 2007, Proceedings of the National Academy of Sciences.

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

[45]  John J. Tyson,et al.  Model-driven experimental approach reveals the complex regulatory distribution of p53 by the circadian factor Period 2 , 2016, Proceedings of the National Academy of Sciences.

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

[47]  John J Tyson,et al.  Functional motifs in biochemical reaction networks. , 2010, Annual review of physical chemistry.

[48]  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.

[49]  M. Biggin Animal transcription networks as highly connected, quantitative continua. , 2011, Developmental cell.

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

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

[52]  D. Frishman,et al.  Protein abundance profiling of the Escherichia coli cytosol , 2008, BMC Genomics.

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

[54]  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.

[55]  D. Gillespie,et al.  Legitimacy of the stochastic Michaelis-Menten approximation. , 2011, IET systems biology.

[56]  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.

[57]  S. Mangan,et al.  Structure and function of the feed-forward loop network motif , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[58]  P J Staff A stochastic development of the reversible Michaelis-Menten mechanism. , 1970, Journal of theoretical biology.

[59]  Frank Kelly,et al.  Reversibility and Stochastic Networks , 1979 .

[60]  Eric Shea-Brown,et al.  The What and Where of Adding Channel Noise to the Hodgkin-Huxley Equations , 2011, PLoS Comput. Biol..

[61]  Guang Qiang Dong,et al.  Simplification of Stochastic Chemical Reaction Models with Fast and Slow Dynamics , 2007, Journal of biological physics.

[62]  N. Kampen,et al.  Elimination of fast variables , 1985 .

[63]  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.

[64]  J. Tóth,et al.  A full stochastic description of the Michaelis-Menten reaction for small systems. , 1977, Acta biochimica et biophysica; Academiae Scientiarum Hungaricae.

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

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

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

[68]  Hye-Won Kang A multiscale approximation in a heat shock response model of E. coli , 2011, BMC Systems Biology.

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

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

[71]  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..

[72]  Jean Mairesse,et al.  Deficiency Zero Petri Nets and Product Form , 2009, Petri Nets.

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

[74]  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.

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

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

[77]  Guilherme C P Innocentini,et al.  Comment on "Steady-state fluctuations of a genetic feedback loop: an exact solution" [J. Chem. Phys. 137, 035104 (2012)]. , 2013, The Journal of chemical physics.

[78]  Shinya Kuroda,et al.  Temporal coding of insulin action through multiplexing of the AKT pathway. , 2012, Molecular cell.

[79]  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.

[80]  Andreas Hellander,et al.  Accuracy of the Michaelis–Menten approximation when analysing effects of molecular noise , 2015, Journal of The Royal Society Interface.

[81]  Eric Vanden-Eijnden,et al.  NUMERICAL TECHNIQUES FOR MULTI-SCALE DYNAMICAL SYSTEMS WITH STOCHASTIC EFFECTS ⁄ , 2003 .

[82]  M. Feinberg The existence and uniqueness of steady states for a class of chemical reaction networks , 1995 .

[83]  Guido Sanguinetti,et al.  The complex chemical Langevin equation. , 2014, The Journal of chemical physics.

[84]  Min Zhou,et al.  A Period2 Phosphoswitch Regulates and Temperature Compensates Circadian Period. , 2015, Molecular cell.

[85]  Terence Hwa,et al.  Transcriptional regulation by the numbers: models. , 2004, Current opinion in genetics & development.

[86]  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.

[87]  T. Kurtz The Relationship between Stochastic and Deterministic Models for Chemical Reactions , 1972 .