Validation-based model selection for 13C metabolic flux analysis with uncertain measurement errors

Accurate measurements of metabolic fluxes in living cells are central to metabolism research and metabolic engineering. The gold standard method is model-based metabolic flux analysis (MFA), where fluxes are estimated indirectly from mass isotopomer data with the use of a mathematical model of the metabolic network. A critical step in MFA is model selection: choosing what compartments, metabolites, and reactions to include in the metabolic network model. Model selection is often done informally during the modelling process, based on the same data that is used for model fitting (estimation data). This can lead to either overly complex models (overfitting) or too simple ones (underfitting), in both cases resulting in poor flux estimates. Here, we propose a method for model selection based on independent validation data. We demonstrate in simulation studies that this method consistently chooses the correct model in a way that is independent on errors in measurement uncertainty. This independence is beneficial, since estimating the true magnitude of these errors can be difficult. In contrast, commonly used model selection methods based on the χ²-test choose different model structures depending on the believed measurement uncertainty; this can lead to errors in flux estimates, especially when the magnitude of the error is substantially off. We present a new approach for quantification of prediction uncertainty of mass isotopomer distributions in other labelling experiments, to check for problems with too much or too little novelty in the validation data. Finally, in an isotope tracing study on human mammary epithelial cells, the validation-based model selection method identified pyruvate carboxylase as a key model component. Our results argue that validation-based model selection should be an integral part of MFA model development.     

Improving the accuracy of flux balance analysis through the implementation of carbon availability constraints for intracellular reactions

Constraint-based modeling methods, such as Flux Balance Analysis (FBA), have been extensively used to decipher complex, information rich -omics datasets to elicit system-wide behavioral patterns of cellular metabolism. FBA has been successfully used to gain insight in a wide range of applications, such as range of substrate utilization, product yields and to design metabolic engineering strategies to improve bioprocess performance. A well-known challenge associated with large genome-scale metabolic networks is that they result in underdetermined problem formulations. Consequently, rather than unique solutions, FBA and related methods examine ranges of reaction flux values that are consistent with the studied physiological conditions. The wider the reported flux ranges, the higher the uncertainty in the determination of basic reaction properties, limiting interpretability of and confidence in the results. Herein, we propose a new, computationally efficient approach that refines flux range predictions by constraining reaction fluxes on the basis of the elemental balance of carbon. We compared carbon constraint FBA (ccFBA) against experimentally-measured intracellular fluxes using the latest CHO GEM (iCHO1766) and were able to substantially improve the accuracy of predicted flux values compared with FBA. ccFBA can be used as a stand-alone method but is also compatible with and complimentary to other constraint-based approaches.     

INCA 2.0: A tool for integrated,dynamic modeling of NMR- and MS-based isotopomer measurements and rigorous metabolic flux analysis

Metabolic flux analysis (MFA) combines experimental measurements and computational modeling to determine biochemical reaction rates in live biological systems. Advancements in analytical instrumentation, such as nuclear magnetic resonance (NMR) spectroscopy and mass spectrometry (MS), have facilitated chemical separation and quantification of isotopically enriched metabolites. However, no software packages have been previously described that can integrate isotopomer measurements from both MS and NMR analytical platforms and have the flexibility to estimate metabolic fluxes from either isotopic steady-state or dynamic labeling experiments. By applying physiologically relevant cardiac and hepatic metabolic models to assess NMR isotopomer measurements, we herein test and validate new modeling capabilities of our enhanced flux analysis software tool, INCA 2.0. We demonstrate that INCA 2.0 can simulate and regress steady-state ¹³C NMR datasets from perfused hearts with an accuracy comparable to other established flux assessment tools. Furthermore, by simulating the infusion of three different ¹³C acetate tracers, we show that MFA based on dynamic ¹³C NMR measurements can more precisely resolve cardiac fluxes compared to isotopically steady-state flux analysis. Finally, we show that estimation of hepatic fluxes using combined ¹³C NMR and MS datasets improves the precision of estimated fluxes by up to 50%. Overall, our results illustrate how the recently added NMR data modeling capabilities of INCA 2.0 can enable entirely new experimental designs that lead to improved flux resolution and can be applied to a wide range of biological systems and measurement time courses.     

Interpretation of metabolic flux maps by limitation potentials and constrained limitation sensitivities

Two new concepts, "Limitation Potential" and "Constraint Limitation Sensitivity" are introduced that use definitions derived from metabolic flux analysis (MFA) and metabolic network analysis (MNA). They are applied to interpret a measured flux distribution in the context of all possible flux distributions and thus combine MFA with MNA. The proposed measures are used to quantify and compare the influence of intracellular fluxes on the production yield. The methods are purely based on the stoichiometry of the network and constraints that are given from irreversible fluxes. In contrast to metabolic control analysis (MCA), within this approach no information about the kinetic mechanisms are needed. A limitation potential (LP) is defined as the reduction of the reachable (theoretical) maximum by a measured flux. Measured fluxes that strongly narrow the reachable maximum are assumed to be limiting as the network has no ability to counterbalance the restriction due to the observed flux. In a second step, the sensitivity of the reduced maximum is regarded. This measure provides information about the necessitated changes to reach higher yields. The methods are applied to interpret the capabilities of a network based on measured fluxes for a L-phenylalanine producer. The strain was examined by a series of experiments and three flux maps of the production phase are analyzed. It can be shown that the reachable yield is drastically reduced by the measured efflux into the TCA cycle, while the oxidative pentose-phosphate pathway only plays a secondary role on the reachable maximum.     

Incorporating metabolic flux ratios into constraint-based flux analysis by using artificial metabolites and converging ratio determinants

One of the well-established approaches for the quantitative characterization of large-scale underdetermined metabolic network is constraint-based flux analysis, which quantifies intracellular metabolic fluxes to characterize the metabolic status. The system is typically underdetermined, and thus usually is solved by linear programming with the measured external fluxes as constraints. Thus, the intracellular flux distribution calculated may not represent the true values. (13)C-constrained flux analysis allows more accurate determination of internal fluxes, but is currently limited to relatively small metabolic networks due to the requirement of complicated mathematical formulation and limited parameters available. Here, we report a strategy of employing such partial information obtained from the (13)C-labeling experiments as additional constraints during the constraint-based flux analysis. A new methodology employing artificial metabolites and converging ratio determinants (CRDs) was developed for improving constraint-based flux analysis. The CRDs were determined based on the metabolic flux ratios obtained from (13)C-labeling experiments, and were incorporated into the mass balance equations for the artificial metabolites. These new mass balance equations were used as additional constraints during the constraint-based flux analysis with genome-scale E. coli metabolic model, which allowed more accurate determination of intracellular metabolic fluxes.     

A scalable algorithm to explore the Gibbs energy landscape of genome-scale metabolic networks

The integration of various types of genomic data into predictive models of biological networks is one of the main challenges currently faced by computational biology. Constraint-based models in particular play a key role in the attempt to obtain a quantitative understanding of cellular metabolism at genome scale. In essence, their goal is to frame the metabolic capabilities of an organism based on minimal assumptions that describe the steady states of the underlying reaction network via suitable stoichiometric constraints, specifically mass balance and energy balance (i.e. thermodynamic feasibility). The implementation of these requirements to generate viable configurations of reaction fluxes and/or to test given flux profiles for thermodynamic feasibility can however prove to be computationally intensive. We propose here a fast and scalable stoichiometry-based method to explore the Gibbs energy landscape of a biochemical network at steady state. The method is applied to the problem of reconstructing the Gibbs energy landscape underlying metabolic activity in the human red blood cell, and to that of identifying and removing thermodynamically infeasible reaction cycles in the Escherichia coli metabolic network (iAF1260). In the former case, we produce consistent predictions for chemical potentials (or log-concentrations) of intracellular metabolites; in the latter, we identify a restricted set of loops (23 in total) in the periplasmic and cytoplasmic core as the origin of thermodynamic infeasibility in a large sample (10(6)) of flux configurations generated randomly and compatibly with the prior information available on reaction reversibility.     

A kinetic model of Escherichia coli core metabolism satisfying multiple sets of mutant flux data

In contrast to stoichiometric-based models, the development of large-scale kinetic models of metabolism has been hindered by the challenge of identifying kinetic parameter values and kinetic rate laws applicable to a wide range of environmental and/or genetic perturbations. The recently introduced ensemble modeling (EM) procedure provides a promising remedy to address these challenges by decomposing metabolic reactions into elementary reaction steps and incorporating all phenotypic observations, upon perturbation, in its model parameterization scheme. Here, we present a kinetic model of Escherichia coli core metabolism that satisfies the fluxomic data for wild-type and seven mutant strains by making use of the EM concepts. This model encompasses 138 reactions, 93 metabolites and 60 substrate-level regulatory interactions accounting for glycolysis/gluconeogenesis, pentose phosphate pathway, TCA cycle, major pyruvate metabolism, anaplerotic reactions and a number of reactions in other parts of the metabolism. Parameterization is performed using a formal optimization approach that minimizes the discrepancies between model predictions and flux measurements. The predicted fluxes by the model are within the uncertainty range of experimental flux data for 78% of the reactions (with measured fluxes) for both the wild-type and seven mutant strains. The remaining flux predictions are mostly within three standard deviations of reported ranges. Converting the EM-based parameters into a Michaelis–Menten equivalent formalism revealed that 35% of K_m and 77% of k_cat parameters are within uncertainty range of the literature-reported values. The predicted metabolite concentrations by the model are also within uncertainty ranges of metabolomic data for 68% of the metabolites. A leave-one-out cross-validation test to evaluate the flux prediction performance of the model showed that metabolic fluxes for the mutants located in the proximity of mutations used for training the model can be predicted more accurately. The constructed model and the parameterization procedure presented in this study pave the way for the construction of larger-scale kinetic models with more narrowly distributed parameter values as new metabolomic/fluxomic data sets are becoming available for E. coli and other organisms.     

Isotopomer analysis of myocardial substrate metabolism: a systems biology approach

The increasing accessibility of mass isotopomer data via GC-MS and NMR technology has necessitated the use of a systematic and reliable method to take advantage of such data for flux analysis. Here we applied a nonlinear, optimization-based method to study substrate metabolism in cardiomyocytes using (13)C data from perfused mouse hearts. The myocardial metabolic network used in this study accounts for 257 reactions and 240 metabolites, which are further compartmentalized into extracellular space, cytosol, and mitochondrial matrix. Analysis of the perfused mouse heart showed that the steady-state ATP production rate was 16.6 +/- 2.3 micromol/min . gww, with 30% of the ATP coming from glycolysis. Of the four substrates available in the perfusate (glucose, pyruvate, lactate, and oleate), exogenous glucose forms the majority of cytosolic pyruvate. Pyruvate decaboxylation is significantly higher than carboxylation, suggesting that anaplerosis is low in the perfused heart. Exchange fluxes were predicted to be high for reversible enzymes in the citric acid cycle (CAC), but low in the glycolytic pathway. Pseudoketogenesis amounted to approximately 50% of the net ketone body uptake. Sensitivity analysis showed that the estimated flux distributions were relatively insensitive to experimental errors. The application of isotopomer data drastically improved the estimation of reaction fluxes compared to results computed with respect to reaction stoichiometry alone. Further study of 12 commonly used (13)C glucose mixtures showed that the mixtures of 20% [U-(13)C(6)] glucose, 80% [3 (13)C] glucose and 20% [U-(13)C(6)] glucose, 80% [4 (13)C] were best for resolving fluxes in the current network.     

Application of hybrid cybernetic model in simulating myeloma cell culture co-consuming glucose and glutamine with mixed consumption patterns

A dynamic model called hybrid cybernetic model (HCM) based on structured metabolic network is established for simulating mammalian cell metabolism featured with partially substitutable and partially complementary consumption patterns of two substrates, glucose and glutamine. Benefiting from the application of elementary mode analysis (EMA), the complicated metabolic network is decomposed into elementary modes (EMs) facilitating the employment of the hybrid cybernetic framework to investigate the external and internal flux distribution and the regulation mechanism among them. According to different substrate combination, two groups of EMs are obtained, i.e., EMs associated with glucose uptake and simultaneous uptake of glucose and glutamine. Uptake fluxes through various EMs are coupled together via cybernetic variables to maximize substrate uptake. External fluxes and internal fluxes could be calculated and estimated respectively, by the combination of the stoichiometrics of metabolic networks and fluxes through regulated EMs. The model performance is well validated via three sets of experimental data. Through parameter identification of limited number of experimental data, other external metabolites are precisely predicted. The obtained kinetic parameters of three experimental cultures have similar values, which indicates the robustness of the model. Furthermore, the prediction performance of the model is successfully validated based on identified parameters.     

Mass spectrometry for metabolic flux analysis

Mass spectrometry in combination with tracer experiments based on 13C substrates can serve as a powerful tool for the modeling and analysis of intracellular fluxes and the investigation of biochemical networks. The theoretical background for the application of mass spectrometry to metabolic flux analysis is discussed. Mass spectrometry methods are especially useful to determine mass distribution of metabolites. Additional information gained from fragmentation of metabolites, e.g., by electron impact ionization, allows further localization of labeling positions, up to complete resolution of isotopomer pools. To effectively handle mass distributions in simulation experiments, a matrix based general methodology is formulated. The natural isotope distribution of carbon, oxygen, hydrogen and nitrogen in the target metabolites is considered by introduction of correction matrices. It is shown by simulation results for the central carbon metabolism that neglecting natural isotope distributions causes significant errors in intracellular flux distributions. By varying relative fluxes into pentosephosphate pathway and pyruvate carboxylation reaction, marked changes in the mass distributions of metabolites result, which are illustrated for pyruvate, oxaloacetate, and alpha-ketoglutarate. In addition mass distributions of metabolites are significantly influenced over a broad range by the degree of reversibility of transaldolase and transketolase reactions in the pentosephosphate pathway. The mass distribution of metabolites is very sensitive towards intracellular flux patterns and can be measured with high accuracy by routine mass spectrometry methods. Copyright 1999 John Wiley & Sons, Inc.     

Development of a mathematical model for evaluating the dynamics of normal and apoptotic Chinese hamster ovary cells

A metabolic flux based methodology was developed for modeling the metabolism of a Chinese hamster ovary cell line. The elimination of insignificant fluxes resulted in a simplified metabolic network which was the basis for modeling the significant metabolites. Employing kinetic rate expressions for growing and non-growing subpopulations, a logistic model was developed for cell growth and dynamic models were formulated to describe culture composition and monoclonal antibody (MAb) secretion. The model was validated for a range of nutrient concentrations. Good agreement was obtained between model predictions and experimental data. The ultimate goal of this study is to establish a comprehensive dynamic model which may be used for model-based optimization of the cell culture for MAb production in both batch and fed-batch systems.     

A guide to metabolic flux analysis in metabolic engineering: Methods,tools and applications

The field of metabolic engineering is primarily concerned with improving the biological production of value-added chemicals, fuels and pharmaceuticals through the design, construction and optimization of metabolic pathways, redirection of intracellular fluxes, and refinement of cellular properties relevant for industrial bioprocess implementation. Metabolic network models and metabolic fluxes are central concepts in metabolic engineering, as was emphasized in the first paper published in this journal, “Metabolic fluxes and metabolic engineering” (Metabolic Engineering, 1: 1–11, 1999). In the past two decades, a wide range of computational, analytical and experimental approaches have been developed to interrogate the capabilities of biological systems through analysis of metabolic network models using techniques such as flux balance analysis (FBA), and quantify metabolic fluxes using constrained-based modeling approaches such as metabolic flux analysis (MFA) and more advanced experimental techniques based on the use of stable-isotope tracers, i.e. ¹³C-metabolic flux analysis (¹³C-MFA). In this review, we describe the basic principles of metabolic flux analysis, discuss current best practices in flux quantification, highlight potential pitfalls and alternative approaches in the application of these tools, and give a broad overview of pragmatic applications of flux analysis in metabolic engineering practice.     

The principle of flux minimization and its application to estimate stationary fluxes in metabolic networks.

Cellular functions are ultimately linked to metabolic fluxes brought about by thousands of chemical reactions and transport processes. The synthesis of the underlying enzymes and membrane transporters causes the cell a certain 'effort' of energy and external resources. Considering that those cells should have had a selection advantage during natural evolution that enabled them to fulfil vital functions (such as growth, defence against toxic compounds, repair of DNA alterations, etc.) with minimal effort, one may postulate the principle of flux minimization, as follows: given the available external substrates and given a set of functionally important 'target' fluxes required to accomplish a specific pattern of cellular functions, the stationary metabolic fluxes have to become a minimum. To convert this principle into a mathematical method enabling the prediction of stationary metabolic fluxes, the total flux in the network is measured by a weighted linear combination of all individual fluxes whereby the thermodynamic equilibrium constants are used as weighting factors, i.e. the more the thermodynamic equilibrium lies on the right-hand side of the reaction, the larger the weighting factor for the backward reaction. A linear programming technique is applied to minimize the total flux at fixed values of the target fluxes and under the constraint of flux balance (= steady-state conditions) with respect to all metabolites. The theoretical concept is applied to two metabolic schemes: the energy and redox metabolism of erythrocytes, and the central metabolism of Methylobacterium extorquens AM1. The flux rates predicted by the flux-minimization method exhibit significant correlations with flux rates obtained by either kinetic modelling or direct experimental determination. Larger deviations occur for segments of the network composed of redundant branches where the flux-minimization method always attributes the total flux to the thermodynamically most favourable branch. Nevertheless, compared with existing methods of structural modelling, the principle of flux minimization appears to be a promising theoretical approach to assess stationary flux rates in metabolic systems in cases where a detailed kinetic model is not yet available.     

Parallel labeling experiments with [U-(13)C]glucose validate E. coli metabolic network model for (13)C metabolic flux analysis

(13)C-metabolic flux analysis (MFA) is a widely used method for measuring intracellular metabolic fluxes in living cells. (13)C MFA relies on several key assumptions: (1) the assumed metabolic network model is complete, in that it accounts for all significant enzymatic and transport reactions; (2) (13)C-labeling measurements are accurate and precise; and (3) enzymes and transporters do not discriminate between (12)C- and (13)C-labeled metabolites. In this study, we tested these inherent assumptions of (13)C MFA for wild-type E. coli by parallel labeling experiments with [U-(13)C]glucose as tracer. Cells were grown in six parallel cultures in custom-constructed mini-bioreactors, starting from the same inoculum, on medium containing different mixtures of natural glucose and fully labeled [U-(13)C]glucose, ranging from 0% to 100% [U-(13)C]glucose. Macroscopic growth characteristics of E. coli showed no observable kinetic isotope effect. The cells grew equally well on natural glucose, 100% [U-(13)C]glucose, and mixtures thereof. (13)C MFA was then used to determine intracellular metabolic fluxes for several metabolic network models: an initial network model from literature; and extended network models that accounted for potential dilution effects of isotopic labeling. The initial network model did not give statistically acceptable fits and produced inconsistent flux results for the parallel labeling experiments. In contrast, an extended network model that accounted for dilution of intracellular CO(2) by exchange with extracellular CO(2) produced statistically acceptable fits, and the estimated metabolic fluxes were consistent for the parallel cultures. This study illustrates the importance of model validation for (13)C MFA. We show that an incomplete network model can produce statistically unacceptable fits, as determined by a chi-square test for goodness-of-fit, and return biased metabolic fluxes. The validated metabolic network model for E. coli from this study can be used in future investigations for unbiased metabolic flux measurements.     

Reconstruction and flux analysis of coupling between metabolic pathways of astrocytes and neurons: application to cerebral hypoxia

Background

It is a daunting task to identify all the metabolic pathways of brain energy metabolism and develop a dynamic simulation environment that will cover a time scale ranging from seconds to hours. To simplify this task and make it more practicable, we undertook stoichiometric modeling of brain energy metabolism with the major aim of including the main interacting pathways in and between astrocytes and neurons.

Model

The constructed model includes central metabolism (glycolysis, pentose phosphate pathway, TCA cycle), lipid metabolism, reactive oxygen species (ROS) detoxification, amino acid metabolism (synthesis and catabolism), the well-known glutamate-glutamine cycle, other coupling reactions between astrocytes and neurons, and neurotransmitter metabolism. This is, to our knowledge, the most comprehensive attempt at stoichiometric modeling of brain metabolism to date in terms of its coverage of a wide range of metabolic pathways. We then attempted to model the basal physiological behaviour and hypoxic behaviour of the brain cells where astrocytes and neurons are tightly coupled.

Results

The reconstructed stoichiometric reaction model included 217 reactions (184 internal, 33 exchange) and 216 metabolites (183 internal, 33 external) distributed in and between astrocytes and neurons. Flux balance analysis (FBA) techniques were applied to the reconstructed model to elucidate the underlying cellular principles of neuron-astrocyte coupling. Simulation of resting conditions under the constraints of maximization of glutamate/glutamine/GABA cycle fluxes between the two cell types with subsequent minimization of Euclidean norm of fluxes resulted in a flux distribution in accordance with literature-based findings. As a further validation of our model, the effect of oxygen deprivation (hypoxia) on fluxes was simulated using an FBA-derivative approach, known as minimization of metabolic adjustment (MOMA). The results show the power of the constructed model to simulate disease behaviour on the flux level, and its potential to analyze cellular metabolic behaviour in silico.

Conclusion

The predictive power of the constructed model for the key flux distributions, especially central carbon metabolism and glutamate-glutamine cycle fluxes, and its application to hypoxia is promising. The resultant acceptable predictions strengthen the power of such stoichiometric models in the analysis of mammalian cell metabolism.     

Pathway analysis of Pichia pastoris to elucidate methanol metabolism and its regulation for production of recombinant proteins

This research rationally analyzes metabolic pathways of Pichia pastoris to study the metabolic flux responses of this yeast under methanol metabolism. A metabolic model of P. pastoris was constructed and analyzed by elementary mode analysis (EMA). EMA was used to comprehensively identify the cell's metabolic flux profiles and its underlying regulation mechanisms for the production of recombinant proteins from methanol. Change in phenotypes and flux profiles during methanol adaptation with varying feed mixture of glycerol and methanol was examined. EMA identified increasing and decreasing fluxes during the glycerol–methanol metabolic shift, which well agreed with experimental observations supporting the validity of the metabolic network model. Analysis of all the identified pathways also led to the determination of the metabolic capacities as well as the optimum metabolic pathways for recombinant protein synthesis during methanol induction. The network sensitivity analysis revealed that the production of proteins can be improved by manipulating the flux ratios at the pyruvate branch point. In addition, EMA suggested that protein synthesis is optimum under hypoxic culture conditions. The metabolic modeling and analysis presented in this study could potentially form a valuable knowledge base for future research on rational design and optimization of P. pastoris by determining target genes, pathways, and culture conditions for enhanced recombinant protein synthesis. The metabolic pathway analysis is also of considerable value for production of therapeutic proteins by P. pastoris in biopharmaceutical applications. © 2013 American Institute of Chemical Engineers Biotechnol. Prog., 30:28–37, 2014