Abstract:Mechanistic dynamic models have become an essential tool for understanding biomolecular networks and other biological systems. Biochemical stochasticity can be extremely important in some situations, e.g., at the single-cell level where there is a low copy number of the species involved. In these scenarios, deterministic models are not suitable to characterize the dynamics, so stochastic dynamic models should be considered. Here, we address the challenging problem of parameter estimation in stochastic dynamic models. Despite recent advances, this area is considerably less mature than its deterministic counterpart. We present a novel strategy based on two components: (i) global optimization via a hybrid stochastic-deterministic approach, and (ii) stochastic simulation techniques tailored to the sparsity of the available experimental data. Regarding the latter, for cases of dense population data we make use of a novel approach using a Partial Integro-Differential Equation (PIDE) model solved using a semilagrangian method. In order to further speed up the simulations, we also present efficient parallel implementations for multi-core CPUs and also for graphical processing units (GPUs). Importantly, whereas SDE and Fokker Planck approximations of the Chemical Master Equation (CME) apply when the reactant populations are sufficiently large, the PIDE approximation to the CME is valid for very low copy numbers, and therefore they enable us to tackle parameter estimation for systems with large intrinsic molecular noise, (highly stochastic regimes far from the thermodynamic limit). We test our strategy with four challenging problems: a Lotka-Volterra system, a polarization system in S. cerevisiae, a genetic toggle switch, and a genetic circadian oscillator. Our method could successfully solve these problems in very reasonable computation times (often a few minutes for the first two problems) using standard low-cost hardware, showing very significant speedups with respect to recent alternative methods. The code used to obtain the results reported here is available at https://doi.org/10.5281/zenodo.5195408.

Stochastic simulation of biochemical systems with randomly fluctuating rate constants

Stochastic Simulation of Chemical Chaos

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

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

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

A general method for numerically simulating the stochastic time evolution of coupled chemical reactions

Dynamic disorder in simple enzymatic reactions induces stochastic amplification of substrate

Stochastic Simulation of Biomolecular Networks in Dynamic Environments

Optimisation of Simulations of Stochastic Processes by Removal of Opposing Reactions

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

Modelling non-Markovian dynamics in biochemical reactions

Stochastic quasi-steady state approximations for asymptotic solutions of the chemical master equation

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

Simulation and inference algorithms for stochastic biochemical reaction networks: from basic concepts to state-of-the-art

Stochastic theory of large-scale enzyme-reaction networks: Finite copy number corrections to rate equation models

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

Learning heterogenous reaction rates from stochastic simulations

Enzyme kinetic reactions as interacting particle systems: Stochastic averaging and parameter inference

Stochastic Differential Equations for Practical Simulation of Gene Circuits

Stochasticity or Noise in Biochemical Reactions

Global Optimization Approach for Parameter Estimation in Stochastic Dynamic Models of Biosystems