Comparison of different moment-closure approximations for stochastic chemical kinetics

Abstract

In recent years moment-closure approximations (MA) of the chemical master equation have become a popular method for the study of stochastic effects in chemical reaction systems. Several different MA methods have been proposed and applied in the literature, but it remains unclear how they perform with respect to each other. In this paper we study the normal, Poisson, log-normal and central-moment-neglect MAs by applying them to understand the stochastic properties of chemical systems whose deterministic rate equations show the properties of bistability, ultrasensitivity and oscillatory behaviour. Our results suggest that the normal MA is favourable over the other studied MAs. In particular we found that (i) the size of the region of parameter space where a closure gives physically meaningful results, e.g. positive mean and variance, is considerably larger for the normal closure than for the other three closures; (ii) the accuracy of the predictions of the four closures (relative to simulations using the stochastic simulation algorithm) is comparable in those regions of parameter space where all closures give physically meaningful results; (iii) the Poisson and log-normal MAs are not uniquely defined for systems involving conservation laws in molecule numbers. We also describe the new software package MOCA which enables the automated numerical analysis of various MA methods in a graphical user interface and which was used to perform the comparative analysis presented in this paper. MOCA allows the user to develop novel closure methods and can treat polynomial, non-polynomial, as well as time-dependent propensity functions, thus being applicable to virtually any chemical reaction system.

36 pages, 14 figures

References in corpus (1)

• Validity conditions for moment closure approximations in stochastic chemical kinetics

Cited by in corpus (21)

• Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
• Master equations and the theory of stochastic path integrals
• Generalized Method of Moments for Estimating Parameters of Stochastic Reaction Networks
• Parameter estimation for biochemical reaction networks using Wasserstein distances
• Rigorous bounds on the stationary distributions of the chemical master equation via mathematical programming
• A variational approach to moment-closure approximations for the kinetics of biomolecular reaction networks
• Efficient Low-Order Approximation of First-Passage Time Distributions
• Statistical physics approaches to subnetwork dynamics in biochemical systems
• Flows, scaling, and the control of moment hierarchies for stochastic chemical reaction networks
• Solving Moment Hierarchies for Chemical Reaction Networks
• Expectation propagation for continuous time stochastic processes
• Uniformisation techniques for stochastic simulation of chemical reaction networks
• Optimal control of protein copy number
• Autoregressive Point-Processes as Latent State-Space Models: a Moment-Closure Approach to Fluctuations and Autocorrelations
• Regulation of stem cell dynamics through volume exclusion
• Interval Analysis of Worst-case Stationary Moments for Stochastic Chemical Reactions with Uncertain Parameters
• Gaussian Closure Scheme in the Quasi-Linkage Equilibrium Regime of Evolving Genome Populations
• Insights into the relation between noise and biological complexity
• Improved estimations of stochastic chemical kinetics by finite state expansion
• An algebraic method to calculate parameter regions for constrained steady-state distribution in stochastic reaction networks
• Logistic dynamics of small populations with demographic stochasticity