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.

Fast Evaluation of Fluctuations in Biochemical Networks With the Linear Noise Approximation

Extending the linear-noise approximation to biochemical systems influenced by intrinsic noise and slow lognormally distributed extrinsic noise

Inference for reaction networks using the Linear Noise Approximation

LNA++: Linear Noise Approximation with First and Second Order Sensitivities

Stochastic enzyme kinetics and the quasi-steady-state reductions: Application of the slow scale linear noise approximation à la Fenichel

Bayesian inference of biochemical kinetic parameters using the linear noise approximation

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

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

Stochastic lumping analysis for linear kinetics and its application to the fluctuation relations between hierarchical kinetic networks

Noise Effects in Nonlinear Biochemical Signaling

Propagation of Fluctuations in Biochemical Systems, I: Linear SSC Networks

Coarse-graining stochastic biochemical networks: quasi-stationary approximation and fast simulations using a stochastic path integral technique

Structured Projection-Based Model Reduction with Application to Stochastic Biochemical Networks

Noise and critical phenomena in biochemical signaling cycles at small molecule numbers

CERENA: ChEmical REaction Network Analyzer—A Toolbox for the Simulation and Analysis of Stochastic Chemical Kinetics

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

Stochastic simulation of biochemical systems with randomly fluctuating rate constants

Distribution of GABA immunoreactivity in the central and peripheral nervous system of amphioxus (Branchiostoma lanceolatum pallas)

Mass Action Dynamics of Coupled Reactions using Fluctuation Theory

Investigating the two-moment characterisation of subcellular biochemical networks

Steady-state fluctuations of a genetic feedback loop with fluctuating rate parameters using the unified colored noise approximation