Optimization and model reduction in the high dimensional parameter space of a budding yeast cell cycle model

Abstract

BackgroundParameter estimation from experimental data is critical for mathematical modeling of protein regulatory networks. For realistic networks with dozens of species and reactions, parameter estimation is an especially challenging task. In this study, we present an approach for parameter estimation that is effective in fitting a model of the budding yeast cell cycle (comprising 26 nonlinear ordinary differential equations containing 126 rate constants) to the experimentally observed phenotypes (viable or inviable) of 119 genetic strains carrying mutations of cell cycle genes.ResultsStarting from an initial guess of the parameter values, which correctly captures the phenotypes of only 72 genetic strains, our parameter estimation algorithm quickly improves the success rate of the model to 105–111 of the 119 strains. This success rate is comparable to the best values achieved by a skilled modeler manually choosing parameters over many weeks. The algorithm combines two search and optimization strategies. First, we use Latin hypercube sampling to explore a region surrounding the initial guess. From these samples, we choose ∼20 different sets of parameter values that correctly capture wild type viability. These sets form the starting generation of differential evolution that selects new parameter values that perform better in terms of their success rate in capturing phenotypes. In addition to producing highly successful combinations of parameter values, we analyze the results to determine the parameters that are most critical for matching experimental outcomes and the most competitive strains whose correct outcome with a given parameter vector forces numerous other strains to have incorrect outcomes. These “most critical parameters” and “most competitive strains” provide biological insights into the model. Conversely, the “least critical parameters” and “least competitive strains” suggest ways to reduce the computational complexity of the optimization.ConclusionsOur approach proves to be a useful tool to help systems biologists fit complex dynamical models to large experimental datasets. In the process of fitting the model to the data, the tool identifies suggestive correlations among aspects of the model and the data.

Quantitative determinants of aerobic glycolysis identify flux through the enzyme GAPDH as a limiting step

“…Such a sampling was chosen to be large enough to cover all feasible solutions of aerobic glycolysis. The sampling was carried out using a Latin Hypercube sampling method (Oguz et al, 2013). For each chosen set of random parameter vectors the model was simulated between 2000 and 5000 times and its output was calculated.…”

A Mathematical Framework for Understanding Four-Dimensional Heterogeneous Differentiation of $$\hbox {CD4}^{+}$$ CD4 + T Cells

“…Starting with this population, we next implemented Differential Evolution (DE). The two-stage optimization approach based on LHS and DE has been presented previously by Oguz et al(Oguz et al 2013). Details of LHS and DE are provided in the Supplementary Text.…”

Conditional robustness analysis for fragility discovery and target identification in biochemical networks and in cancer systems biology

“…Once the pathway structure is drawn, the corresponding equations are relatively easy to write down using widely accepted kinetic laws, such as the law of mass action, the Michaelis-Menten law or the Hill functions. The generated model will depend on several parameters, and the corresponding identification and model selection problems are relevant issues (see [ 18 , 19 ] and the references therein).…”