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 Optimization Based Study of Dimerization Kinetics

Dynamic Optimization of Batch Free Radical Polymerization with Conditional Modeling Formulation Through the Adaptive Smoothing Strategy

Adaptable Parallel Acceleration Strategy For Dynamic Monte Carlo Simulations Of Polymerization With Microscopic Resolution

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

Frequency Adaptive Metadynamics for the Calculation of Rare-Event Kinetics.

Optimisation of Simulations of Stochastic Processes by Removal of Opposing Reactions

Parallel in Time Simulation of Multiscale Stochastic Chemical Kinetics

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

Lattice Models for Magnetic Fluids: Correlations Between Order Parameters

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

Extraction of kinetics from equilibrium distributions of states using the Metropolis Monte Carlo method

Computer Aided Recipe Design: Optimization of Polydisperse Chemical Mixtures Using Molecular Descriptors

Learning effective amino acid interactions through iterative stochastic techniques

Determining the DNA stability parameters for the breathing dynamics of heterogeneous DNA by stochastic optimization

Advanced isoconversional kinetic analysis: insight in mechanisms and simulations. Successes and future

Simulation and non-linear optimization of kinetic models for solid-state processes

Optimization-based Simultaneous Modelling of Stoichiometries and Kinetics in Complex Organic Reaction System

Simulating Protein Conformations through Global Optimization

Simultaneous Optimization of Biomolecular Energy Functions on Features from Small Molecules and Macromolecules.

Variational Selection of Features for Molecular Kinetics