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.

An Extended Super-Twisting Algorithm for simultaneous estimation of reaction rates and input disturbance in bioprocess

Examining a Pretty Simple and Low Cost Method for Modeling of Biogas Production from Biodegradable Solids

Robust asymptotic super twisting sliding mode observer for non-linear uncertain biochemical systems

Production Optimization for Concentration and Volume-Limited Fed-Batch Reactors in Biochemical Processes

Simulation results for on-line optimization of a batch bioreactor using nonlinear filtering and optimal control

Monitoring Ethanol Fermentation in Real Time by a Robust State Observer for Uncertainties

Global Stabilizing Control of a Continuous Ethanol Fermentation Process Starting from Batch Mode Production

State estimation and nonlinear tracking control simulation approach. Application to a bioethanol production system

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

On the Robustness of a Modified Super-Twisting Algorithm With Prescribed-Time Convergence

Modeling and estimation of production rate for the production phase of non-growth-associated high cell density processes

Fractional Lotka-Volterra model with time-delay and delayed controller for a bioreactor

Exponential super-twisting control for nonlinear systems with unknown polynomial perturbations

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

Disturbance-observer-based finite time sliding mode controller with unmatched uncertainties utilizing improved cubature Kalman filter

Adjoint Sensitivity Analysis on Multi-Scale Bioprocess Stochastic Reaction Network

Extremum seeking via continuation techniques for optimizing biogas production in the chemostat

Identification in an anaerobic batch system: global sensitivity analysis, multi-start strategy and optimization criterion selection

An identification algorithm of switched Box-Jenkins systems in the presence of bounded disturbances: An approach for approximating complex biological wastewater treatment models

Neural Network Model-Based On-Line Re-Optimisation Control of Fed-Batch Processes Using a Modified Iterative Dynamic Programming Algorithm

An Adaptive Generalized Super-Twisting Algorithm via Event-Triggered Control