Production-Passage-Time Approximation: A New Approximation Method to Accelerate the Simulation Process of Enzymatic Reactions

Given the substantial computational requirements of stochastic simulation, approximation is essential for efficient analysis of any realistic biochemical system. This paper introduces a new approximation method to reduce the computational cost of stochastic simulations of an enzymatic reaction scheme which in biochemical systems often includes rapidly changing fast reactions with enzyme and enzyme-substrate complex molecules present in very small counts. Our new method removes the substrate dissociation reaction by approximating the passage time of the formation of each enzyme-substrate complex molecule which is destined to a production reaction. This approach skips the firings of unimportant yet expensive reaction events, resulting in a substantial acceleration in the stochastic simulations of enzymatic reactions. Additionally, since all the parameters used in our new approach can be derived by the Michaelis-Menten parameters which can actually be measured from experimental data, applications of this approximation can be practical even without having full knowledge of the underlying enzymatic reaction. Furthermore, since our approach does not require a customized simulation procedure for enzymatic reactions, it allows biochemical systems that include such reactions to still take advantage of standard stochastic simulation tools. Here, we apply this new method to various enzymatic reaction systems, resulting in a speedup of orders of magnitude in temporal behavior analysis without any significant loss in accuracy.

[1]  A Finney,et al.  Systems biology markup language: Level 2 and beyond. , 2003, Biochemical Society transactions.

[2]  Sidney Yip,et al.  Handbook of Materials Modeling , 2005 .

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

[4]  S Achimescu,et al.  Signal propagation in nonlinear stochastic gene regulatory networks. , 2005, Systems biology.

[5]  D. Gillespie Markov Processes: An Introduction for Physical Scientists , 1991 .

[6]  D. Gillespie The chemical Langevin equation , 2000 .

[7]  N. Friedman,et al.  Stochastic protein expression in individual cells at the single molecule level , 2006, Nature.

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

[9]  A. Arkin,et al.  Stochastic kinetic analysis of developmental pathway bifurcation in phage lambda-infected Escherichia coli cells. , 1998, Genetics.

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

[11]  Jeffrey W. Smith,et al.  Stochastic Gene Expression in a Single Cell , .

[12]  A. Arkin,et al.  Stochastic amplification and signaling in enzymatic futile cycles through noise-induced bistability with oscillations. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[13]  C. Rao,et al.  Control, exploitation and tolerance of intracellular noise , 2002, Nature.

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

[15]  J. Derisi,et al.  Single-cell proteomic analysis of S. cerevisiae reveals the architecture of biological noise , 2006, Nature.

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

[17]  Axel Groß,et al.  Quantum Theory of Reactive Scattering and Adsorption at Surfaces , 2005 .

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

[19]  Michael S. Samoilov,et al.  Automated Abstraction Methodology for Genetic Regulatory Networks , 2006, Trans. Comp. Sys. Biology.

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

[21]  D. Gillespie,et al.  Accelerated stochastic simulation of the stiff enzyme-substrate reaction. , 2005, The Journal of chemical physics.

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

[23]  C. Rao,et al.  Stochastic chemical kinetics and the quasi-steady-state assumption: Application to the Gillespie algorithm , 2003 .

[24]  A. van Oudenaarden,et al.  Noise Propagation in Gene Networks , 2005, Science.

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

[26]  Michael Hucka,et al.  Systems Biology Markup Language (SBML) Level 2: Structures and Facilities for Model Definitions , 2007, WWW 2007.