Abstract:Motivation: Stochastic molecular processes are a leading cause of cell‐to‐cell variability. Their dynamics are often described by continuous‐time discrete‐state Markov chains and simulated using stochastic simulation algorithms. As these stochastic simulations are computationally demanding, ordinary differential equation models for the dynamics of the statistical moments have been developed. The number of state variables of these approximating models, however, grows at least quadratically with the number of biochemical species. This limits their application to small‐ and medium‐sized processes. Results: In this article, we present a scalable moment‐closure approximation (sMA) for the simulation of statistical moments of large‐scale stochastic processes. The sMA exploits the structure of the biochemical reaction network to reduce the covariance matrix. We prove that sMA yields approximating models whose number of state variables depends predominantly on local properties, i.e. the average node degree of the reaction network, instead of the overall network size. The resulting complexity reduction is assessed by studying a range of medium‐ and large‐scale biochemical reaction networks. To evaluate the approximation accuracy and the improvement in computational efficiency, we study models for JAK2/STAT5 signalling and NF&kgr;B signalling. Our method is applicable to generic biochemical reaction networks and we provide an implementation, including an SBML interface, which renders the sMA easily accessible. Availability and implementation: The sMA is implemented in the open‐source MATLAB toolbox CERENA and is available from https://github.com/CERENADevelopers/CERENA. Contact: jan.hasenauer@helmholtz‐muenchen.de or atefeh.kazeroonian@tum.de Supplementary information: Supplementary data are available at Bioinformatics online.

A local polynomial moment approximation for compartmentalised biochemical systems

Modeling of Stochastic Biological Processes with Non-Polynomial Propensities Using Non-Central Conditional Moment Equation

Method of Conditional Moments (MCM) for the Chemical Master Equation

Inference for Stochastic Chemical Kinetics Using Moment Equations and System Size Expansion

A Scalable Moment-Closure Approximation for Large-Scale Biochemical Reaction Networks

An approach to the residence time distribution for stochastic multi-compartment models

Probabilistic analysis of a class of compartmental models formulated by random differential equations

Bounds on stationary moments in stochastic chemical kinetics

Moment-based analysis of biochemical networks in a heterogeneous population of communicating cells

Learning Moment Closure in Reaction-Diffusion Systems with Spatial Dynamic Boltzmann Distributions

Moment-based parameter inference with error guarantees for stochastic reaction networks

Tighter bounds on transient moments of stochastic chemical systems

Investigating the two-moment characterisation of subcellular biochemical networks

Deterministic Bounding Systems for Stochastic Compartmental Spreading Processes

Moment Propagation of Polynomial Systems Through Carleman Linearization for Probabilistic Safety Analysis

A New Moment Method for Solving the Coagulation Equation for Particles in Brownian Motion

Effective one-dimension reduction of multi-compartment complex systems dynamics

Approximation and inference methods for stochastic biochemical kinetics - a tutorial review

Exact and Approximate Moment Derivation for Probabilistic Loops With Non-Polynomial Assignments

Inferring the dynamics of quasi-reaction systems via nonlinear local mean-field approximations

Modelling non-Markovian dynamics in biochemical reactions