2015 Volume 5 Issue 3
Article Contents

Tae-Hyuk Ahn, Xiaoying Han, Adrian Sandu. IMPLICIT SIMULATION METHODS FOR STOCHASTIC CHEMICAL KINETICS[J]. Journal of Applied Analysis & Computation, 2015, 5(3): 420-452. doi: 10.11948/2015034
Citation: Tae-Hyuk Ahn, Xiaoying Han, Adrian Sandu. IMPLICIT SIMULATION METHODS FOR STOCHASTIC CHEMICAL KINETICS[J]. Journal of Applied Analysis & Computation, 2015, 5(3): 420-452. doi: 10.11948/2015034

IMPLICIT SIMULATION METHODS FOR STOCHASTIC CHEMICAL KINETICS

  • Fund Project:
  • 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.
    MSC: 65C05;65C30
  • 加载中
  • [1] T.-H. Ahn, Y. Cao and L.T. Watson, Stochastic simulation algorithms for chemical reactions, BIOCOMP, 2008, 431-436.

    Google Scholar

    [2] H. R. Arabnia, M. Q. Yang and J.Y. Yang, International Conference on Bioinformatics & Computational Biology, BIOCOMP, CSREA Press, Las Vegas Nevada, 2008, 14-17.

    Google Scholar

    [3] T.-H. Ahn, P. Wang, L.T. Watson, Y. Cao, C.A. Shaffer and W.T. Baumann, Stochastic cell cycle modeling for budding yeast, Proceedings of the 2009 Spring Simulation Multiconference, 2009, 113:1-113:6.

    Google Scholar

    [4] T.-H. Ahn, L.T. Watson, Y. Cao, C.A. Shaffer and W.T. Baumann, Cell cycle modeling for budding yeast with stochastic simulation algorithms, Computer Modeling in Engineering and Sciences, 51(2009), 27-52.

    Google Scholar

    [5] T.-H. Ahn and A. Sandu, Parallel stochastic simulations of budding yeast cell cycle:load balancing strategies and theoretical analysis, Proceedings of the First ACM International Conference on Bioinformatics and Computational Biology, ser. BCB'10, 2010, 237-246.

    Google Scholar

    [6] D.A. Ball, T.-H. Ahn, P. Wang, K.C. Chen, Y. Cao, J.J. Tyson, J. Peccoud and W.T. Baumann, Stochastic exit from mitosis in budding yeast:model predictions and experimental observations, Cell Cycle, 10(2011), 1-11.

    Google Scholar

    [7] T. Ahn and A. Sandu, Fully implicit tau-leaping methods for the stochastic simulation of chemical kinetics, Proceedings of the 19th High Performance Computing Symposia, 2011, 118-125.

    Google Scholar

    [8] T.-H. Ahn and A. Sandu, Implicit second order weak Taylor Tau-Leaping methods for the stochastic simulation of chemical kinetics, Procedia Computer Science, 2011, 2297-2306.

    Google Scholar

    [9] Y. Cao, H. Li and L. Petzold, Efficient formulation of the stochastic simulation algorithm for chemically reacting systems, Journal of Chemical Physics, 121(2004), 4059-4067.

    Google Scholar

    [10] Y. Cao, L. Petzold, M. Rathinam, D. Gillespie, The numerical stability of leaping methods for stochastic simulation of chemically reacting systems, Journal of Chemical Physics, 121(2004), 12169-12178.

    Google Scholar

    [11] Y. Cao and L. Petzold, Trapezoidal Tau-leaping formula for the stochastic simulation of biochemical systems, Proceedings of Foundations of Systems Biology in Engineering (FOSBE 2005), 2005, 149-52.

    Google Scholar

    [12] J. Elf and M. Ehrenberg, Spontaneous separation of bi-stable biochemical systems into spatial domains of opposite phases, Systems Biology, IEE Proceedings, 1(2004), 230-236.

    Google Scholar

    [13] F.E. Streib and M. Dehmer, Information Theory and Statistical Learning, New York, NY, Springer, 2008.

    Google Scholar

    [14] M.A. Gibson and J. Bruck, Efficient exact stochastic simulation of chemical systems with many species and many channels, Journal of Physical Chemistry, 104(2000), 1876-1889.

    Google Scholar

    [15] I.I. Gikhman and A.V. Skorokhod, Stochastic Differential Equations, New York, NY, Springer, 1972.

    Google Scholar

    [16] D.T. Gillespie, Exact stochastic simulation of coupled chemical reactions, Journal of Physical Chemistry, 81(1977), 2340-2361.

    Google Scholar

    [17] D.T. Gillespie, A rigorous derivation of the chemical master equation, Physica A, 188(1992), 404-425.

    Google Scholar

    [18] D.T. Gillespie, Approximate accelerated stochastic simulation of chemically reacting systems, Journal of Chemical Physics, 115(2001), 1716-1733.

    Google Scholar

    [19] D.T. Gillespie and L.R. Petzold, Improved leap-size selection for accelerated stochastic simulation, Journal of Chemical Physics, 119(2003), 8229-8234.

    Google Scholar

    [20] N.G.van Kampen, Stochastic Processes in Physics and Chemistry, North Holland, Netherlands, Third Edition, 2007.

    Google Scholar

    [21] P.E. Kloeden and E. Platen, Numerical Solution of Stochastic Differential Equations, New York, NY, Springer, 1999.

    Google Scholar

    [22] H.H. McAdams and A. Arkin, Stochastic mechanisms in gene expression, Proc. Natl. Acad. Sci., 94(1997), 814-819.

    Google Scholar

    [23] M. Rathinam, L.R. Petzold, Y. Cao and D.T. Gillespie, Stiffness in stochastic chemically reacting systems:the implicit Tau-leaping method, Journal of Chemical Physics, 119(2003), 12784-12794.

    Google Scholar

    [24] R. Muruhan, L.R. Petzold, Y. Cao and D.T. Gillespie, Consistency and stability of Tau Leaping Schemes for Chemical Reaction Systems, SIAM Journal of Multiscale Modeling and Simulation, 4(2005), 867-895.

    Google Scholar

    [25] S.M. Ross, Introduction to Probability Models, Ninth Edition, Academic Press, Orlando, FL, USA, Inc., 2006.

    Google Scholar

    [26] T. Tian and K. Burrage, Implicit Taylor methods for stiff stochastic differential equations, Applied Numerical Mathematics, 38(2001), 167-185.

    Google Scholar

    [27] T. Li, Analysis of explicit Tau-Leaping schemes for simulating chemically reacting systems, Multiscale Modeling and Simulation, 6(2007), 417-436.

    Google Scholar

    [28] Y. Hu, T. Li and B. Min, A weak second order tau-leaping method for chemical kinetic systems, Journal of Chemical Physics, 135(2011), 024113.

    Google Scholar

    [29] Tatiana T. Marquez-Lago and Kevin Burrage Binomial tau-leap spatial stochastic simulation algorithm for applications in chemical kinetics, Journal of Chemical Physics, 127(2007), 104101.

    Google Scholar

Article Metrics

Article views(1937) PDF downloads(614) Cited by(0)

Access History

Other Articles By Authors

Catalog

    /

    DownLoad:  Full-Size Img  PowerPoint