The complex chemical Langevin equation

Abstract

The chemical Langevin equation (CLE) is a popular simulation method to probe the stochastic dynamics of chemical systems. The CLE's main disadvantage is its break down in finite time due to the problem of evaluating square roots of negative quantities whenever the molecule numbers become sufficiently small. We show that this issue is not a numerical integration problem, rather in many systems it is intrinsic to all representations of the CLE. Various methods of correcting the CLE have been proposed which avoid its break down. We show that these methods introduce undesirable artefacts in the CLE's predictions. In particular, for unimolecular systems, these correction methods lead to CLE predictions for the mean concentrations and variance of fluctuations which disagree with those of the chemical master equation. We show that, by extending the domain of the CLE to complex space, break down is eliminated, and the CLE's accuracy for unimolecular systems is restored. Although the molecule numbers are generally complex, we show that the "complex CLE" predicts real-valued quantities for the mean concentrations, the moments of intrinsic noise, power spectra and first passage times, hence admitting a physical interpretation. It is also shown to provide a more accurate approximation of the chemical master equation of simple biochemical circuits involving bimolecular reactions than the various corrected forms of the real-valued CLE, the linear-noise approximation and a commonly used two moment-closure approximation.

47 pages, 8 figures

References in corpus (2)

• Steady-state fluctuations of a genetic feedback loop: an exact solution
• How accurate are the non-linear chemical Fokker-Planck and chemical Langevin equations?

Cited by in corpus (13)

• Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
• Noise-induced multistability in chemical systems: Discrete vs Continuum modeling
• Hybrid framework for the simulation of stochastic chemical kinetics
• Approximate probability distributions of the master equation
• Model reduction for stochastic chemical systems with abundant species
• From interacting agents to density-based modeling with stochastic PDEs
• A practical guide to pseudo-marginal methods for computational inference in systems biology
• Dynamical mean-field theory: from ecosystems to reaction networks
• Design principles for biochemical oscillations with limited energy resources
• The chemical Langevin equation: a path integral view of Gillespie's derivation
• Cooling down stochastic differential equations: almost sure convergence
• Expectation propagation for continuous time stochastic processes
• Insights into the relation between noise and biological complexity