Abstract:Background Ordinary differential equation (ODE) models are widely used to describe (bio-)chemical and biological processes. To enhance the predictive power of these models, their unknown parameters are estimated from experimental data. These experimental data are mostly collected in perturbation experiments, in which the processes are pushed out of steady state by applying a stimulus. The information that the initial condition is a steady state of the unperturbed process provides valuable information, as it restricts the dynamics of the process and thereby the parameters. However, implementing steady-state constraints in the optimization often results in convergence problems. Results In this manuscript, we propose two new methods for solving optimization problems with steady-state constraints. The first method exploits ideas from optimization algorithms on manifolds and introduces a retraction operator, essentially reducing the dimension of the optimization problem. The second method is based on the continuous analogue of the optimization problem. This continuous analogue is an ODE whose equilibrium points are the optima of the constrained optimization problem. This equivalence enables the use of adaptive numerical methods for solving optimization problems with steady-state constraints. Both methods are tailored to the problem structure and exploit the local geometry of the steady-state manifold and its stability properties. A parameterization of the steady-state manifold is not required. The efficiency and reliability of the proposed methods is evaluated using one toy example and two applications. The first application example uses published data while the second uses a novel dataset for Raf/MEK/ERK signaling. The proposed methods demonstrated better convergence properties than state-of-the-art methods employed in systems and computational biology. Furthermore, the average computation time per converged start is significantly lower. In addition to the theoretical results, the analysis of the dataset for Raf/MEK/ERK signaling provides novel biological insights regarding the existence of feedback regulation. Conclusion Many optimization problems considered in systems and computational biology are subject to steady-state constraints. While most optimization methods have convergence problems if these steady-state constraints are highly nonlinear, the methods presented recover the convergence properties of optimizers which can exploit an analytical expression for the parameter-dependent steady state. This renders them an excellent alternative to methods which are currently employed in systems and computational biology.

Differential simulated annealing: a robust and efficient global optimization algorithm for parameter estimation of biological networks.

A New Strategy for Parameter Estimation of Dynamic Differential Equations Based on NSGA II.

A New DNA Genetic Algorithm and Its Application in Parameter Estimation

Scalable Parameter Estimation for Genome-Scale Biochemical Reaction Networks

Identification of Optimal Strategies for State Transition of Complex Biological Networks.

Tailored Parameter Optimization Methods for Ordinary Differential Equation Models with Steady-State Constraints

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

A Parameter Estimation Method for Biological Systems Modelled by ODE/DDE Models Using Spline Approximation and Differential Evolution Algorithm

Differential Evolution Improved with Self-Adaptive Control Parameters Based on Simulated Annealing.

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

The performance analysis of differential evolution algorithm and its application in parameter estimation of a dynamic model

Kriging-Based Parameter Estimation Algorithm for Metabolic Networks Combined with Single-Dimensional Optimization and Dynamic Coordinate Perturbation

Differential annealing for global optimization

Comparing AI versus Optimization Workflows for Simulation-Based Inference of Spatial-Stochastic Systems

An Improved Swarm Optimization for Parameter Estimation and Biological Model Selection

Hybrid Simulation Annealing and Evolution Strategy Algorithm in Non-Linear Parameter Estimation Application

Efficient parameter inference in networked dynamical systems via steady states: A surrogate objective function approach integrating mean-field and nonlinear least squares

Modeling of system of ordinary differential equations for dynamic systems: An experimental study based on different search techniques

Differential Evolution Algorithm Based On Simulated Annealing

Reducing local minima in fitness landscapes of parameter estimation by using piecewise evaluation and state estimation

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