Fast stochastic simulation of biochemical reaction systems by alternative formulations of the chemical Langevin equation
KAUST Grant NumberKUK-C1-013-04
MetadataShow full item record
AbstractThe Chemical Langevin Equation (CLE), which is a stochastic differential equation driven by a multidimensional Wiener process, acts as a bridge between the discrete stochastic simulation algorithm and the deterministic reaction rate equation when simulating (bio)chemical kinetics. The CLE model is valid in the regime where molecular populations are abundant enough to assume their concentrations change continuously, but stochastic fluctuations still play a major role. The contribution of this work is that we observe and explore that the CLE is not a single equation, but a parametric family of equations, all of which give the same finite-dimensional distribution of the variables. On the theoretical side, we prove that as many Wiener processes are sufficient to formulate the CLE as there are independent variables in the equation, which is just the rank of the stoichiometric matrix. On the practical side, we show that in the case where there are m1 pairs of reversible reactions and m2 irreversible reactions there is another, simple formulation of the CLE with only m1 + m2 Wiener processes, whereas the standard approach uses 2 m1 + m2. We demonstrate that there are considerable computational savings when using this latter formulation. Such transformations of the CLE do not cause a loss of accuracy and are therefore distinct from model reduction techniques. We illustrate our findings by considering alternative formulations of the CLE for a human ether a-go-go related gene ion channel model and the Goldbeter-Koshland switch. © 2010 American Institute of Physics.
CitationMélykúti B, Burrage K, Zygalakis KC (2010) Fast stochastic simulation of biochemical reaction systems by alternative formulations of the chemical Langevin equation. J Chem Phys 132: 164109. Available: http://dx.doi.org/10.1063/1.3380661.
SponsorsThe authors wish to thank Thomas G. Kurtz for observing a mistake in an earlier version of this material which has been corrected subsequently. B.M.’s work was funded by the Engineering and Physical Sciences Research Council through a Doctoral Training Centre grant for the Life Sciences Interface Doctoral Training Centre, University of Oxford (Grant No. EP/E501605/1). K. C. Z. was supported by Award No. KUK-C1-013-04 made by the King Abdullah University of Science and Technology (KAUST).
JournalThe Journal of Chemical Physics
CollectionsPublications Acknowledging KAUST Support
- Developing Itô stochastic differential equation models for neuronal signal transduction pathways.
- Authors: Manninen T, Linne ML, Ruohonen K
- Issue date: 2006 Aug
- Accurate hybrid stochastic simulation of a system of coupled chemical or biochemical reactions.
- Authors: Salis H, Kaznessis Y
- Issue date: 2005 Feb 1
- Stochastic chemical kinetics and the total quasi-steady-state assumption: application to the stochastic simulation algorithm and chemical master equation.
- Authors: Macnamara S, Bersani AM, Burrage K, Sidje RB
- Issue date: 2008 Sep 7
- On the origins of approximations for stochastic chemical kinetics.
- Authors: Haseltine EL, Rawlings JB
- Issue date: 2005 Oct 22
- The complex chemical Langevin equation.
- Authors: Schnoerr D, Sanguinetti G, Grima R
- Issue date: 2014 Jul 14