Fast stochastic simulation of biochemical reaction systems by alternative formulations of the chemical Langevin equation

Bence Mélykúti, Kevin Burrage, Konstantinos C. Zygalakis

Research output: Contribution to journalArticlepeer-review

35 Scopus citations

Abstract

The 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.
Original languageEnglish (US)
Pages (from-to)164109
JournalThe Journal of Chemical Physics
Volume132
Issue number16
DOIs
StatePublished - Apr 26 2010
Externally publishedYes

Fingerprint

Dive into the research topics of 'Fast stochastic simulation of biochemical reaction systems by alternative formulations of the chemical Langevin equation'. Together they form a unique fingerprint.

Cite this