Stochastic modeling of biochemical reactions taking place at the cellular level has become the subject of intense research in recent years. Molecular interactions in a single cell exhibit random fluctuations. These fluctuations may be significant when small populations of some reacting species are present and then a stochastic description of the cellular dynamics is required. Often, the biochemically reacting systems encountered in applications consist of many species interacting through many reaction channels. Also, the dynamics of such systems is typically non-linear and presents multiple time-scales. Consequently, the stochastic mathematical models of biochemical systems can be quite complex and their analysis challenging. In this paper, we present a method to reduce a stochastic continuous model of well-stirred biochemical systems, the Chemical Langevin Equation, while preserving the overall behavior of the system. Several tests of our method on models of practical interest gave excellent results.
H. Huang, M. Fairweather, J. F. Griffiths, A. S. Tomlin and R. B. Brad, “A Systematic Lumping Approach for the Reduction of Comprehensive Kinetic Models,” Proceedings of the Combustion Institute, Vol. 30, No.1, 2004, pp. 1309-1316. doi:10.1016/j.proci.2004.08.001
A. Kumar, P. D. Christofides and P. Daoutidis, “Singular Perturbation Modeling of Non-Linear Processes with Non-Explicit Time-Scale Separation,” Chemical Engineering Science, Vol. 53, No. 8, 1998, pp. 1491-1504.
C. V. Rao and A. P. Arkin, “Stochastic Chemical Kinetics and the Quasi-Steady-State Assumption: Application to the Gillespie Algorithm,” The Journal of Chemical Physics, Vol. 118, No. 11, 2003, pp. 4999-5010.
A. Samant and D. Vlachos, “Overcoming Stiffness in Stochastic Simulation Stemming from Partial Equilibrium: A Multiscale Monte-Carlo Algorithm,” The Journal of Chemical Physics, Vol. 123, No. 14, 2005, pp. 144114-144122. doi:10.1063/1.2046628
D. T. Gillespie, “A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions”, The Journal of Computational Physics, Vol. 22, No. 4, 1976, pp. 403-434.
M. Rathinam, P. W. Sheppard and M. Khammash, “Efficient Computation of Parameter Sensitivities of Discrete Stochastic Chemical Reaction Networks,” The Journal of Chemical Physics, Vol. 132, No. 3, 2010, pp. 034103-034116. doi:10.1063/1.3280166
P. W. Sheppard, M. Rathinam and M. Khammash, “A Pathwise Derivative Approach to the Computation of Parameter Sensitivities in Discrete Stochastic Chemical Systems,” The Journal of Chemical Physics, Vol. 136, No. 3, 2012, pp. 034115-034128. doi:10.1063/1.3677230
V. Sotiropoulos, M. N. Contou-Carrere, P. Daoutidis and Y. N. Kaznessis, “Model Reduction of Multiscale Chemical Langevin Equations: A Numerical Case Study,” IEEE/ACM Transactions on Computational Biology and Bioinformatics, Vol. 6, No. 3, 2009, pp. 470-482.
S. L. Cotter, K. C. Zygalakis, I. G. Kevrekidis and R. Erban, “A Constrained Approach to Multiscale Stochastic Simulation of Chemically Reacting Systems,” The Journal of Chemical Physics, Vol. 135, No. 9, 2011, pp. 094102-094114. doi:10.1063/1.3624333