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.

An all-in-one state-observer for protein refolding reactions using particle filters and delayed measurements

Online monitoring of protein refolding in inclusion body processing using intrinsic fluorescence

Advances in monitoring and control of refolding kinetics combining PAT and modeling

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

Harnessing Liquid Crystal Sensors for High-Throughput Real-Time Detection of Structural Changes in Lysozyme during Refolding Processes

Watching Single Unmodified Enzymes at Work

An Electrochemical Biosensor Architecture Based on Protein Folding Supports Direct Real‐Time Measurements in Whole Blood

Measurement of Protein-Protein Interaction Dynamics Using Microfluidics and Particle Diffusometry

Advances in Mixer Design and Detection Methods for Kinetics Studies of Macromolecular Folding and Binding on the Microsecond Time Scale

Filtered finite state projection method for the analysis and estimation of stochastic biochemical reaction networks

Real-time monitoring of quality attributes by in-line Fourier transform infrared spectroscopic sensors at ultrafiltration and diafiltration of bioprocess

Optimizing Bioprocessing Efficiency with OptFed: Dynamic Nonlinear Modeling Improves Product-to-Biomass Yield

Autonomous adaptive optimization of NMR experimental conditions for precise inference of minor conformational states of proteins based on chemical exchange saturation transfer

Exploration of Protein Unfolding by Modelling Calorimetry Data from Reheating

A Systematic Protein Refolding Screen Method using the DGR Approach Reveals that Time and Secondary TSA are Essential Variables

On-chip Measurements of Protein Unfolding from Direct Observations of Micron-Scale Diffusion.

Comparison of Expectation-Maximization based parameter estimation using Particle Filter, Unscented and Extended Kalman Filtering techniques

Protein phosphorylation reagent concentration curves via Raman Spectroscopy for real-time reaction monitoring in a microfluidic reactor

Real-time monitoring of gradient chromatography using dual Kalman-filters

An Evaluation of Force Field Accuracy for the Mini-Protein Chignolin using Markov State Models

Optofluidic Force Induction as a Process Analytical Technology