Implicit Simulation Methods for Stochastic Chemical Kinetics

In biochemical systems some of the chemical species are present with only small numbers of molecules. In this situation discrete and stochastic simulation approaches are more relevant than continuous and deterministic ones. The fundamental Gillespie's stochastic simulation algorithm (SSA) accounts for every reaction event, which occurs with a probability determined by the configuration of the system. This approach requires a considerable computational effort for models with many reaction channels and chemical species. In order to improve efficiency, tau-leaping methods represent multiple firings of each reaction during a simulation step by Poisson random variables. For stiff systems the mean of this variable is treated implicitly in order to ensure numerical stability. This paper develops fully implicit tau-leaping-like algorithms that treat implicitly both the mean and the variance of the Poisson variables. The construction is based on adapting weakly convergent discretizations of stochastic differential equations to stochastic chemical kinetic systems. Theoretical analyses of accuracy and stability of the new methods are performed on a standard test problem. Numerical results demonstrate the performance of the proposed tau-leaping methods.

[1]  Matthias Dehmer,et al.  Information Theory and Statistical Learning , 2010 .

[2]  T. Faniran Numerical Solution of Stochastic Differential Equations , 2015 .

[3]  D. Sherrington Stochastic Processes in Physics and Chemistry , 1983 .

[4]  W. Ebeling Stochastic Processes in Physics and Chemistry , 1995 .

[5]  Dan ie l T. Gil lespie A rigorous derivation of the chemical master equation , 1992 .

[6]  Linda R. Petzold,et al.  Improved leap-size selection for accelerated stochastic simulation , 2003 .

[7]  Bin Min,et al.  A weak second order tau-leaping method for chemical kinetic systems. , 2011, The Journal of chemical physics.

[8]  Tianhai Tian,et al.  Implicit Taylor methods for stiff stochastic differential equations , 2001 .

[9]  Hikmet Budak BIOCOMP'06: the 2006 International Conference on Bioinformatics & Computational Biology , 2006 .

[10]  A. Sandu,et al.  Implicit Second Order Weak Taylor Tau-Leaping Methods for the Stochastic Simulation of Chemical Kinetics , 2011, ICCS.

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

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

[13]  E. Renshaw,et al.  STOCHASTIC DIFFERENTIAL EQUATIONS , 1974 .

[14]  John J. Tyson,et al.  Stochastic exit from mitosis in budding yeast , 2011, Cell cycle.

[15]  Tatiana T Marquez-Lago,et al.  Binomial tau-leap spatial stochastic simulation algorithm for applications in chemical kinetics. , 2007, The Journal of chemical physics.

[16]  Tiejun Li,et al.  Analysis of Explicit Tau-Leaping Schemes for Simulating Chemically Reacting Systems , 2007, Multiscale Model. Simul..

[17]  Clifford A. Shaffer,et al.  Stochastic cell cycle modeling for budding yeast , 2009, SpringSim '09.

[18]  Adrian Sandu,et al.  Fully implicit tau-leaping methods for the stochastic simulation of chemical kinetics , 2011, SpringSim.

[19]  Clifford A. Shaffer,et al.  Cell Cycle Modeling for Budding Yeast with Stochastic Simulation Algorithms , 2008 .

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

[21]  Adrian Sandu,et al.  Parallel stochastic simulations of budding yeast cell cycle: load balancing strategies and theoretical analysis , 2010, BCB '10.

[22]  Solomon Kullback,et al.  Information Theory and Statistics , 1960 .

[23]  S. Shreve,et al.  Stochastic differential equations , 1955, Mathematical Proceedings of the Cambridge Philosophical Society.

[24]  N. Kampen,et al.  Stochastic processes in physics and chemistry , 1981 .

[25]  J. Elf,et al.  Spontaneous separation of bi-stable biochemical systems into spatial domains of opposite phases. , 2004, Systems biology.

[26]  A. Arkin,et al.  Stochastic mechanisms in gene expression. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[27]  Muruhan Rathinam,et al.  Consistency and Stability of Tau-Leaping Schemes for Chemical Reaction Systems , 2005, Multiscale Model. Simul..

[28]  Layne T. Watson,et al.  Stochastic Simulation Algorithms for Chemical Reactions , 2008, BIOCOMP.

[29]  Sheldon M. Ross Introduction to Probability Models. , 1995 .

[30]  Sheldon M. Ross,et al.  Introduction to Probability Models, Eighth Edition , 1972 .

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

[32]  Hong Li,et al.  Efficient formulation of the stochastic simulation algorithm for chemically reacting systems. , 2004, The Journal of chemical physics.

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