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.

Generative machine learning produces kinetic models that accurately characterize intracellular metabolic states

Toward Modeling Metabolic State from Single-Cell Transcriptomics.

Building kinetic models for metabolic engineering

Large-Scale Modeling Approach Reveals Functional Metabolic Shifts during Hepatic Differentiation

Efficient Parameter Estimation Enables the Prediction of Drug Response Using a Mechanistic Pan-Cancer Pathway Model

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Efficient parameterization of large-scale mechanistic models enables drug response prediction for cancer cell lines

Predicting cell morphological responses to perturbations using generative modeling

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

Automation on the generation of genome scale metabolic models

Perspective on automated predictive kinetics using estimates derived from large datasets

Enhancing the accuracy of genome-scale metabolic models with kinetic information

Hybrid physics-informed metabolic cybergenetics: process rates augmented with machine-learning surrogates informed by flux balance analysis

Deep Learning Basedkcatprediction Enables Improved Enzyme Constrained Model Reconstruction

COVRECON: Combining Genome-scale Metabolic Network Reconstruction and Data-driven Inverse Modeling to Reveal Changes in Metabolic Interaction Networks

Omics data for sampling thermodynamically feasible kinetic models

Context-Specific Metabolic Networks Are Consistent with Experiments

Prediction of intracellular metabolic states from extracellular metabolomic data

Scalable computation of intracellular metabolite concentrations

Modelling dynamic host-pathway interactions at the genome scale