Hybrid discrete/continuum algorithms for stochastic reaction networks

Direct solutions of the Chemical Master Equation (CME) governing Stochastic Reaction Networks (SRNs) are generally prohibitively expensive due to excessive numbers of possible discrete states in such systems. To enhance computational efficiency we develop a hybrid approach where the evolution of states with low molecule counts is treated with the discrete CME model while that of states with large molecule counts is modeled by the continuum Fokker-Planck equation. The Fokker-Planck equation is discretized using a 2nd order finite volume approach with appropriate treatment of flux components. The numerical construction at the interface between the discrete and continuum regions implements the transfer of probability reaction by reaction according to the stoichiometry of the system. The performance of this novel hybrid approach is explored for a two-species circadian model with computational efficiency gains of about one order of magnitude.

[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]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[3]  Muruhan Rathinam,et al.  Stiffness in stochastic chemically reacting systems: The implicit tau-leaping method , 2003 .

[4]  Chi-Wang Shu.  Essentially non-oscillatory and weighted essentially non-oscillatory schemes for hyperbolic conservation laws , 1998 .

[5]  A. Harten.  High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[6]  R. LeVeque.  Finite Volume Methods for Hyperbolic Problems: Characteristics and Riemann Problems for Linear Hyperbolic Equations , 2002 .

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

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

[9]  Wang Chi-Shu,et al.  Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws , 1997 .

[10]  Per Lötstedt,et al.  Fokker–Planck approximation of the master equation in molecular biology , 2009 .

[11]  Roger B. Sidje,et al.  Inexact Uniformization Method for Computing Transient Distributions of Markov Chains , 2007, SIAM J. Sci. Comput..

[12]  Stefan Engblom.  A discrete spectral method for the chemical master equation , 2008 .

[13]  C. J. Burden,et al.  A solver for the stochastic master equation applied to gene regulatory networks , 2007 .

[14]  Adrian Sandu,et al.  A new look at the chemical master equation , 2013, Numerical Algorithms.

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

[16]  D. W. Scott,et al.  Multivariate Density Estimation, Theory, Practice and Visualization , 1992 .

[17]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[18]  Muruhan Rathinam,et al.  The numerical stability of leaping methods for stochastic simulation of chemically reacting systems. , 2004, The Journal of chemical physics.

[19]  S. Leibler,et al.  Mechanisms of noise-resistance in genetic oscillators , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Brian Munsky,et al.  The Finite State Projection Approach for the Analysis of Stochastic Noise in Gene Networks , 2008, IEEE Transactions on Automatic Control.

[21]  P. Lax.  Weak solutions of nonlinear hyperbolic equations and their numerical computation , 1954 .

[22]  Travis E. Oliphant,et al.  Python for Scientific Computing , 2007, Computing in Science & Engineering.

[23]  M. Khammash,et al.  The finite state projection algorithm for the solution of the chemical master equation. , 2006, The Journal of chemical physics.

[24]  Roger B. Sidje,et al.  Multiscale Modeling of Chemical Kinetics via the Master Equation , 2008, Multiscale Model. Simul..

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

[26]  Jingwei Zhang,et al.  Radial Basis Function Collocation for the Chemical Master Equation , 2010, BIOCOMP.

[27]  Sheng Wu,et al.  StochKit2: software for discrete stochastic simulation of biochemical systems with events , 2011, Bioinform..

[28]  C. W. Gardiner,et al.  Handbook of stochastic methods - for physics, chemistry and the natural sciences, Second Edition , 1986, Springer series in synergetics.

[29]  D. Gillespie.  A rigorous derivation of the chemical master equation , 1992 .

[30]  D. Gillespie.  A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions , 1976 .

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

[32]  Paul Sjöberg,et al.  PDE and Monte Carlo approaches to solving the master equation applied to gene regulation , 2007 .

[33]  Stefan Engblom.  Spectral approximation of solutions to the chemical master equation , 2009 .

[34]  Yang Cao,et al.  Multiscale stochastic simulation algorithm with stochastic partial equilibrium assumption for chemically reacting systems , 2005 .

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

[36]  B. Silverman.  Density estimation for statistics and data analysis , 1986 .

[37]  John D. Hunter,et al.  Matplotlib: A 2D Graphics Environment , 2007, Computing in Science & Engineering.

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

[39]  Jingwei Zhang,et al.  A modified uniformization method for the solution of the chemical master equation , 2010, Comput. Math. Appl..

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

[41]  Wael Hassan.  Simplified Wrapper and Interface Generator , 2000 .

[42]  Michael A. Gibson,et al.  Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels , 2000 .

[43]  Tobias Jahnke,et al.  Solving chemical master equations by adaptive wavelet compression , 2010, J. Comput. Phys..

[44]  Brian Munsky,et al.  A multiple time interval finite state projection algorithm for the solution to the chemical master equation , 2007, J. Comput. Phys..

[45]  Johan Paulsson,et al.  Near-critical phenomena in intracellular metabolite pools. , 2003, Biophysical journal.

[46]  Markus Hegland,et al.  An optimal Finite State Projection Method , 2010, ICCS.

[47]  D. Gillespie.  Approximate accelerated stochastic simulation of chemically reacting systems , 2001 .

[48]  E Weinan,et al.  Nested stochastic simulation algorithms for chemical kinetic systems with multiple time scales , 2007, J. Comput. Phys..

[49]  Desmond J. Higham,et al.  Chemical Master Equation and Langevin regimes for a gene transcription model , 2008, Theor. Comput. Sci..

[50]  Per Lötstedt,et al.  Adaptive, Conservative Solution of the Fokker-Planck Equation in Molecular Biology , 2004 .