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 review of Girsanov Reweighting and of Square Root Approximation for building molecular Markov State Models

Girsanov reweighting for simulations of underdamped Langevin dynamics. Theory

Implementation of Girsanov Reweighting in OpenMM and Deeptime

Tensor-SqRA: Modeling the Transition Rates of Interacting Molecular Systems in terms of Potential Energies

Markov-type state models to describe non-Markovian dynamics

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

Numerical implementation of some reweighted path integral methods

Accurate estimates of dynamical statistics using memory

RiteWeight: Randomized Iterative Trajectory Reweighting for Steady-State Distributions Without Discretization Error

Optimized Markov State Models for Metastable Systems

Arbitrarily accurate, nonparametric coarse graining with Markov renewal processes and the Mori-Zwanzig formulation

Markov State Models with Weighted Ensemble Simulation: How to Eliminate the Trajectory Merging Bias

Reweighting in nonequilibrium simulations

Reweighted Jarzynski sampling: Acceleration of rare events and free energy calculation with a bias potential learned from nonequilibrium work

Towards a Benchmark for Markov State Models: The Folding of HP35

Improved Reweighting of Accelerated Molecular Dynamics Simulations for Free Energy Calculation

Multiscale molecular kinetics by coupling Markov state models and reaction-diffusion dynamics

Estimation and uncertainty of reversible Markov models