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 Scaling Approach to Stochastic Chemical Reaction Networks

Analysis of Stochastic Chemical Reaction Networks with a Hierarchy of Timescales

Stochastic Chemical Reaction Networks with Discontinuous Limits and AIMD processes

Multiscale Analysis of Reaction Networks

A Stein's Method Approach to the Linear Noise Approximation for Stationary Distributions of Chemical Reaction Networks

Classification and threshold dynamics of stochastic reaction networks

A Scalable Computational Framework for Establishing Long-Term Behavior of Stochastic Reaction Networks

Reduction for stochastic biochemical reaction networks with multiscale conservations

Asymptotic analysis for stationary distributions of scaled reaction networks

Asymptotology of Chemical Reaction Networks

An Algebraic Approach to Product-form Stationary Distributions for Some Reaction Networks

Product-Form Poisson-Like Distributions and Complex Balanced Reaction Systems

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

Complex balanced distributions for chemical reaction networks

Universal Slow Dynamics of Chemical Reaction Networks

Transitions and Thermodynamics on Species Graphs of Chemical Reaction Networks

Fast reactions with non-interacting species in stochastic reaction networks

Stochastic analysis of chemical reactions in multi-component interacting systems at criticality

Stochastic reaction networks within interacting compartments

A network dynamics approach to chemical reaction networks