Estimating parameters for generalized mass action models with connectivity information

Abstract

BackgroundDetermining the parameters of a mathematical model from quantitative measurements is the main bottleneck of modelling biological systems. Parameter values can be estimated from steady-state data or from dynamic data. The nature of suitable data for these two types of estimation is rather different. For instance, estimations of parameter values in pathway models, such as kinetic orders, rate constants, flux control coefficients or elasticities, from steady-state data are generally based on experiments that measure how a biochemical system responds to small perturbations around the steady state. In contrast, parameter estimation from dynamic data requires time series measurements for all dependent variables. Almost no literature has so far discussed the combined use of both steady-state and transient data for estimating parameter values of biochemical systems.ResultsIn this study we introduce a constrained optimization method for estimating parameter values of biochemical pathway models using steady-state information and transient measurements. The constraints are derived from the flux connectivity relationships of the system at the steady state. Two case studies demonstrate the estimation results with and without flux connectivity constraints. The unconstrained optimal estimates from dynamic data may fit the experiments well, but they do not necessarily maintain the connectivity relationships. As a consequence, individual fluxes may be misrepresented, which may cause problems in later extrapolations. By contrast, the constrained estimation accounting for flux connectivity information reduces this misrepresentation and thereby yields improved model parameters.ConclusionThe method combines transient metabolic profiles and steady-state information and leads to the formulation of an inverse parameter estimation task as a constrained optimization problem. Parameter estimation and model selection are simultaneously carried out on the constrained optimization problem and yield realistic model parameters that are more likely to hold up in extrapolations with the model.

Parametric identifier of metabolic network associated to hydrogen production in Escherichia coli based on robust sliding-mode differentiation

“…In parameter characterization of the metabolic network model, one of the major problems is that distinctly different set of parameter values may fit to the experimental data (Ko et al 2009;Meskin et al 2013;Lillacci and Khammash 2010). The results obtained show that none of the errors obtained as a result of applying the algorithm using the estimated derivative is above 5.0 %.…”

Seeding-inspired chemotaxis genetic algorithm for the inference of biological systems

“…Because a genetic network is sparsely connected (Thieffry et al, 1998), researchers have introduced the concentration error, slope error and skeletalstructure penalty (Kikuchi et al, 2003) as pruning (penalty) terms. Wang et al used eÀconstrain to integrate these three penalties (Liu and Wang, 2008;Ko et al, 2009). Kimura et al sorted kinetic orders in a descending order and allowed the penalty to act only on specific genes that are regulated by many genes (Kimura et al, 2005).…”

The best models of metabolism

“…Uncounted reviews have been published describing these challenges, even within the field of biology alone, along with a wide variety of possible solutions. 148,[262][263][264][265][266][267][268][269][270][271][272][273][274][275][276][277][278][279][280] Unfortunately, powerful reverse engineering and parameter estimation algorithms from other fields are often not easily applicable as they require long, dense time series. 281,282 In spite of strong, concerted efforts, the scientific community is still awaiting effective solutions that work in most applications.…”