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.

Parameter estimation in modulated, unbranched reaction chains within biochemical systems

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Complexity analysis and parameter estimation of dynamic metabolic systems

Bayesian Parameter Estimation for Biochemical Reaction Networks Using Region-Based Adaptive Parallel Tempering

Estimation of parameters in rational reaction rates of molecular biological systems via weighted least squares

Parameter Estimation for Dynamical Systems with Discrete Events and Logical Operations

Data2Dynamics: a Modeling Environment Tailored to Parameter Estimation in Dynamical Systems.

Bayesian inference of biochemical kinetic parameters using the linear noise approximation

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

Modeling, Simulating, and Parameter Fitting of Biochemical Kinetic Experiments

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

Global parameter identification of stochastic reaction networks from single trajectories

Parameter estimation in large-scale systems biology models: a parallel and self-adaptive cooperative strategy

Robust Parameter Estimation for Biological Systems: A Study on the Dynamics of Microbial Communities

A parallelizable sampling method for parameter inference of large biochemical reaction models

Inferring Catalysis in Biological Systems

High-dimensional Bayesian Parameter Estimation: Case Study for a Model of JAK2/STAT5 Signaling.

Data-driven Bayesian estimation of Monod kinetics

Separable Nonlinear Least-Squares Parameter Estimation for Complex Dynamic Systems

Efficient probabilistic inference in biochemical networks

Parameter inference for stochastic biochemical models from perturbation experiments parallelised at the single cell level