Validity conditions for moment closure approximations in stochastic chemical kinetics

Abstract

Approximations based on moment-closure (MA) are commonly used to obtain estimates of the mean molecule numbers and of the variance of fluctuations in the number of molecules of chemical systems. The advantage of this approach is that it can be far less computationally expensive than exact stochastic simulations of the chemical master equation. Here we numerically study the conditions under which the MA equations yield results reflecting the true stochastic dynamics of the system. We show that for bistable and oscillatory chemical systems with deterministic initial conditions, the solution of the MA equations can be interpreted as a valid approximation to the true moments of the CME, only when the steady-state mean molecule numbers obtained from the chemical master equation fall within a certain finite range. The same validity criterion for monostable systems implies that the steady-state mean molecule numbers obtained from the chemical master equation must be above a certain threshold. For mean molecule numbers outside of this range of validity, the MA equations lead to either qualitatively wrong oscillatory dynamics or to unphysical predictions such as negative variances in the molecule numbers or multiple steady-state moments of the stationary distribution as the initial conditions are varied. Our results clarify the range of validity of the MA approach and show that pitfalls in the interpretation of the results can only be overcome through the systematic comparison of the solutions of the MA equations of a certain order with those of higher orders.

arXiv comment27 pages, 13 figures; accepted at Journal of Chemical Physics

References in corpus (3)

• Analytical distributions for stochastic gene expression
• Steady-state fluctuations of a genetic feedback loop: an exact solution
• A general moment expansion method for stochastic kinetic models

Cited by in corpus (16)

• Approximation and inference methods for stochastic biochemical kinetics - a tutorial review
• Master equations and the theory of stochastic path integrals
• Comparison of different moment-closure approximations for stochastic chemical kinetics
• Dissipation in noisy chemical networks: The role of deficiency
• Stochastic reaction networks in dynamic compartment populations
• Parameter estimation for biochemical reaction networks using Wasserstein distances
• Molecular finite-size effects in stochastic models of equilibrium chemical systems
• 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
• Expectation propagation for continuous time stochastic processes
• Autoregressive Point-Processes as Latent State-Space Models: a Moment-Closure Approach to Fluctuations and Autocorrelations
• Insights into the relation between noise and biological complexity
• Splitting Nodes and Linking Channels: A Method for Assembling Biocircuits from Stochastic Elementary Units
• 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