Dataset Information

Parameter optimization in S-system models.

ABSTRACT:

Background

The inverse problem of identifying the topology of biological networks from their time series responses is a cornerstone challenge in systems biology. We tackle this challenge here through the parameterization of S-system models. It was previously shown that parameter identification can be performed as an optimization based on the decoupling of the differential S-system equations, which results in a set of algebraic equations.

Results

A novel parameterization solution is proposed for the identification of S-system models from time series when no information about the network topology is known. The method is based on eigenvector optimization of a matrix formed from multiple regression equations of the linearized decoupled S-system. Furthermore, the algorithm is extended to the optimization of network topologies with constraints on metabolites and fluxes. These constraints rejoin the system in cases where it had been fragmented by decoupling. We demonstrate with synthetic time series why the algorithm can be expected to converge in most cases.

Conclusion

A procedure was developed that facilitates automated reverse engineering tasks for biological networks using S-systems. The proposed method of eigenvector optimization constitutes an advancement over S-system parameter identification from time series using a recent method called Alternating Regression. The proposed method overcomes convergence issues encountered in alternate regression by identifying nonlinear constraints that restrict the search space to computationally feasible solutions. Because the parameter identification is still performed for each metabolite separately, the modularity and linear time characteristics of the alternating regression method are preserved. Simulation studies illustrate how the proposed algorithm identifies the correct network topology out of a collection of models which all fit the dynamical time series essentially equally well.

SUBMITTER: Vilela M

PROVIDER: S-EPMC2333970 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:Ordinary differential equation (ODE) models are widely used to study biochemical reactions in cellular networks since they effectively describe the temporal evolution of these networks using mass action kinetics. The parameters of these models are rarely known a priori and must instead be estimated by calibration using experimental data. Optimization-based calibration of ODE models on is often challenging, even for low-dimensional problems. Multiple hypotheses have been advanced to explain why biochemical model calibration is challenging, including non-identifiability of model parameters, but there are few comprehensive studies that test these hypotheses, likely because tools for performing such studies are also lacking. Nonetheless, reliable model calibration is essential for uncertainty analysis, model comparison, and biological interpretation. We implemented an established trust-region method as a modular Python framework (fides) to enable systematic comparison of different approaches to ODE model calibration involving a variety of Hessian approximation schemes. We evaluated fides on a recently developed corpus of biologically realistic benchmark problems for which real experimental data are available. Unexpectedly, we observed high variability in optimizer performance among different implementations of the same mathematical instructions (algorithms). Analysis of possible sources of poor optimizer performance identified limitations in the widely used Gauss-Newton, BFGS and SR1 Hessian approximation schemes. We addressed these drawbacks with a novel hybrid Hessian approximation scheme that enhances optimizer performance and outperforms existing hybrid approaches. When applied to the corpus of test models, we found that fides was on average more reliable and efficient than existing methods using a variety of criteria. We expect fides to be broadly useful for ODE constrained optimization problems in biochemical models and to be a foundation for future methods development.

Project description:As countries attach importance to environmental protection, clean energy has become a hot topic. Among them, solar energy, as one of the efficient and easily accessible clean energy sources, has received widespread attention. An essential component in converting solar energy into electricity are solar cells. However, a major optimization difficulty remains in precisely and effectively calculating the parameters of photovoltaic (PV) models. In this regard, this study introduces an improved rime optimization algorithm (RIME), namely ERINMRIME, which integrates the Nelder-Mead simplex (NMs) with the environment random interaction (ERI) strategy. In the later phases of ERINMRIME, the ERI strategy serves as a complementary mechanism for augmenting the solution space exploration ability of the agent. By facilitating external interactions, this method improves the algorithm's efficacy in conducting a global search by keeping it from becoming stuck in local optima. Moreover, by incorporating NMs, ERINMRIME enhances its ability to do local searches, leading to improved space exploration. To evaluate ERINMRIME's optimization performance on PV models, this study conducted experiments on four different models: the single diode model (SDM), the double diode model (DDM), the three-diode model (TDM), and the photovoltaic (PV) module model. The experimental results show that ERINMRIME reduces root mean square error for SDM, DDM, TDM, and PV module models by 46.23%, 59.32%, 61.49%, and 23.95%, respectively, compared with the original RIME. Furthermore, this study compared ERINMRIME with nine improved classical algorithms. The results show that ERINMRIME is a remarkable competitor. Ultimately, this study evaluated the performance of ERINMRIME across three distinct commercial PV models, while considering varying irradiation and temperature conditions. The performance of ERINMRIME is superior to existing similar algorithms in different irradiation and temperature conditions. Therefore, ERINMRIME is an algorithm with great potential in identifying and recognizing unknown parameters of PV models.