Unbiased tau-leap methods for stochastic simulation of chemically reacting systems.

The tau-leap method first developed by Gillespie [D. T. Gillespie, J. Chem. Phys. 115, 1716 (2001)] can significantly speed up stochastic simulation of certain chemically reacting systems with acceptable losses in accuracy. Recently, several improved tau-leap methods, including the binomial, multinomial, and modified tau-leap methods, have been developed. However, in all these tau-leap methods, the mean of the number of times, K(m), that the mth reaction channel fires during a leap is not equal to the true mean. Therefore, all existing tau-leap methods produce biased simulation results, which limit the simulation accuracy and speed. In this paper, we analyze the mean of K(m) based on the chemical master equation. Using this analytical result, we develop unbiased Poisson and binomial tau-leap methods. Moreover, we analyze the variance of K(m), and then develop an unbiased Poisson/Gaussian/binomial tau-leap method to correct the errors in both the mean and variance of K(m). Simulation results demonstrate that our unbiased tau-leap method can significantly improve simulation accuracy without sacrificing speed.

[1]  D. Gillespie,et al.  Avoiding negative populations in explicit Poisson tau-leaping. , 2005, The Journal of chemical physics.

[2]  Voratas Kachitvichyanukul,et al.  Binomial random variate generation , 1988, CACM.

[3]  B. Kholodenko,et al.  Quantification of Short Term Signaling by the Epidermal Growth Factor Receptor* , 1999, The Journal of Biological Chemistry.

[4]  L. Devroye.  Generating the maximum of independent identically distributed random variables , 1980 .

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

[6]  K. Burrage,et al.  Binomial leap methods for simulating stochastic chemical kinetics. , 2004, The Journal of chemical physics.

[7]  William H. Press,et al.  Numerical recipes in C , 2002 .

[8]  Haluk Resat,et al.  Multinomial tau-leaping method for stochastic kinetic simulations. , 2007, The Journal of chemical physics.

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

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

[11]  H. Wiley,et al.  An integrated model of epidermal growth factor receptor trafficking and signal transduction. , 2003, Biophysical journal.

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

[13]  Xiaodong Cai,et al.  K-leap method for accelerating stochastic simulation of coupled chemical reactions. , 2007, The Journal of chemical physics.

[14]  Linda R Petzold,et al.  Efficient step size selection for the tau-leaping simulation method. , 2006, The Journal of chemical physics.

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

[16]  D. Vlachos,et al.  Binomial distribution based tau-leap accelerated stochastic simulation. , 2005, The Journal of chemical physics.

[17]  Anne Auger,et al.  R-leaping: accelerating the stochastic simulation algorithm by reaction leaps. , 2006, The Journal of chemical physics.

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

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