Optimization-based metabolic control analysis

In this work, a novel optimization-based metabolic control analysis (OMCA) method is introduced for reducing data requirement for metabolic control analysis (MCA). It is postulated that using the optimal control approach, the fluxes in a metabolic network are correlated to metabolite concentrations and enzyme activities as a state-feedback control system that is optimal with respect to a homeostasis objective. It is then shown that the optimal feedback gains are directly related to the elasticity coefficients (ECs) of MCA. This approach requires determination of the relative "importance" of metabolites and fluxes for the system, which is possible with significantly reduced experimental data, as compared with typical MCA requirements. The OMCA approach is applied to a top-down control model of glycolysis in hepatocytes. It is statistically demonstrated that the OMCA model is capable of predicting the ECs observed experimentally with few exceptions. Further, an OMCA-based model reconciliation study shows that the modification of four assumed stoichiometric coefficients in the model can explain most of the discrepancies, with the exception of elasticities with respect to the NADH/NAD ratio.