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.

A SOLUTION TO DISCONTINUOUS PARAMETER INDENTIFICATION OF PARTIAL DIFFERENTIAL EQUATIONS INVERSE PROBLEMS

Improved gene expression programming algorithm for solving inverse problems in partial differential equations

Parameter identification in PDEs by the solution of monotone inclusion problems

A Nonlinear Multigrid Method for the Parameter Identification Problem of Partial Differential Equations with Constraints

Multi-parameters Identification Problem for a Degenerate Parabolic Equation.

A New Method For Parameter Identifiability Of Dae Systems

An Improved Population-Based Incremental Learning Method For Inverse Problems

Inverse problems for evolutionary quasi-variational hemivariational inequalities with application to mixed boundary value problems

EXISTENCE AND NONEXISTENCE OF SOLUTIONS FOR A CLASS OF P-Laplacian-type FRACTIONAL FOUR-POINT BOUNDARY-VALUE PROBLEMS WITH A PARAMETER

A hybrid regularization method for identifying the source term and the initial value simultaneously for fractional pseudo-parabolic equation with involution

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

Inverse problem of reconstructing source term for a class of non‐divergence parabolic equations

Evolutionary Quasi-variational Hemivariational Inequalities: Existence and Parameter Identification

Parameter Identification for Partial Differential Equations with Spatiotemporal Varying Coefficients

A Restrained Optimal Perturbation Method for Solving a Kind of Inverse Problem with Variable Coefficient.

Parameter Identification of Discrete-time Linear Time-invariant Systems Using State and Input Data

A PARAMETERIZATION FIXED POINT METHOD FOR DIRECT SOLUTION OF OPTIMAL CONTROL PROBLEMS

Numerical Method for Parameter Inference of Systems of Nonlinear Ordinary Differential Equations with Partial Observations.

Fast numerical solvers for parameter identification problems in mathematical biology

Changes in human eighth nerve action potential as a function of stimulation rate.

Local Parameter Identifiability: Case of Discrete Infinite-Dimensional Parameter