Metabolic flux analysis in Escherichia coli by integrating isotopic dynamic and isotopic stationary 13C labeling data

The novel concept of isotopic dynamic 13C metabolic flux analysis (ID-13C MFA) enables integrated analysis of isotopomer data from isotopic transient and/or isotopic stationary phase of a 13C labeling experiment, short-time experiments, and an extended range of applications of 13C MFA. In the presented work, an experimental and computational framework consisting of short-time 13C labeling, an integrated rapid sampling procedure, a LC-MS analytical method, numerical integration of the system of isotopomer differential equations, and estimation of metabolic fluxes was developed and applied to determine intracellular fluxes in glycolysis, pentose phosphate pathway (PPP), and citric acid cycle (TCA) in Escherichia coli grown in aerobic, glucose-limited chemostat culture at a dilution rate of D = 0.10 h(-1). Intracellular steady state concentrations were quantified for 12 metabolic intermediates. A total of 90 LC-MS mass isotopomers were quantified at sampling times t = 0, 91, 226, 346, 589 s and at isotopic stationary conditions. Isotopic stationarity was reached within 10 min in glycolytic and PPP metabolites. Consistent flux solutions were obtained by ID-13C MFA using isotopic dynamic and isotopic stationary 13C labeling data and by isotopic stationary 13C MFA (IS-13C MFA) using solely isotopic stationary data. It is demonstrated that integration of dynamic 13C labeling data increases the sensitivity of flux estimation, particularly at the glucose-6-phosphate branch point. The identified split ratio between glycolysis and PPP was 55%:44%. These results were confirmed by IS-13C MFA additionally using labeling data in proteinogenic amino acids (GC-MS) obtained after 5 h from sampled biomass.     

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.     

Rapid simulation and analysis of isotopomer distributions using constraints based on enzyme mechanisms: an example from HT29 cancer cells

MOTIVATION: Addition of labeled substrates and the measurement of the subsequent distribution of the labels in isotopomers in reaction networks provide a unique method for assessing metabolic fluxes in whole cells. However, owing to insufficiency of information, attempts to quantify the fluxes often yield multiple possible sets of solutions that are consistent with a given experimental pattern of isotopomers. In the study of the pentose phosphate pathways, the need to consider isotope exchange reactions of transketolase (TK) and transaldolase (TA) (which in past analyses have often been ignored) magnifies this problem; but accounting for the interrelation between the fluxes known from biochemical studies and kinetic modeling solves it. The mathematical relationships between kinetic and equilibrium constants restrict the domain of estimated fluxes to the ones compatible not only with a given set of experimental data, but also with other biochemical information. METHOD: We present software that integrates kinetic modeling with isotopomer distribution analysis. It solves the ordinary differential equations for total concentrations (accounting for the kinetic mechanisms) as well as for all isotopomers in glycolysis and the pentose phosphate pathway (PPP). In the PPP the fluxes created in the TK and TA reactions are expressed through unitary rate constants. The algorithms that account for all the kinetic and equilbrium constant constraints are integrated with the previously developed algorithms, which have been further optimized. The most time-consuming calculations were programmed directly in assembly language; this gave an order of magnitude decrease in the computation time, thus allowing analysis of more complex systems. The software was developed as C-code linked to a program written in Mathematica (Wolfram Research, Champaign, IL), and also as a C++ program independent from Mathematica. RESULTS: Implementing constraints imposed by kinetic and equilibrium constants in the isotopomer distribution analysis in the data from the cancer cells eliminated estimates of fluxes that were inconsistent with the kinetic mechanisms of TK and TA. Fluxes measured experimentally in cells can be used to estimate better the kinetics of TK and TA as they operate in situ. Thus, our approach of integrating various methods for in situ flux analysis opens up the possibility of designing new types of experiments to probe metabolic interrelationships, including the incorporation of additional biochemical information. AVAILABILITY: Software is available freely at: http://www.bq.ub.es/bioqint/selivanov.htm CONTACT: martacascante@ub.edu     

Numerical bias estimation for mass spectrometric mass isotopomer analysis

Mass spectrometric (MS) isotopomer analysis has become a standard tool for investigating biological systems using stable isotopes. In particular, metabolic flux analysis uses mass isotopomers of metabolic products typically formed from ¹³C-labeled substrates to quantitate intracellular pathway fluxes. In the current work, we describe a model-driven method of numerical bias estimation regarding MS isotopomer analysis. Correct bias estimation is crucial for measuring statistical qualities of measurements and obtaining reliable fluxes. The model we developed for bias estimation corrects a priori unknown systematic errors unique for each individual mass isotopomer peak. For validation, we carried out both computational simulations and experimental measurements. From stochastic simulations, it was observed that carbon mass isotopomer distributions and measurement noise can be determined much more precisely only if signals are corrected for possible systematic errors. By removing the estimated background signals, the residuals resulting from experimental measurement and model expectation became consistent with normality, experimental variability was reduced, and data consistency was improved. The method is useful for obtaining systematic error-free data from ¹³C tracer experiments and can also be extended to other stable isotopes. As a result, the reliability of metabolic fluxes that are typically computed from mass isotopomer measurements is increased.     

Quantitative analysis of metabolic fluxes in Escherichia coli, using two-dimensional NMR spectroscopy and complete isotopomer models.

Complete isotopomer models that simulate distribution of label in 13C tracer experiments are applied to the quantification of metabolic fluxes in the primary carbon metabolism of E. coli under aerobic and anaerobic conditions. The concept of isotopomer mapping matrices (IMMs) is used to simplify the formulation of isotopomer mass balances by expressing all isotopomer mass balances of a metabolite pool in a single matrix equation. A numerically stable method to calculate the steady-state isotopomer distribution in metabolic networks in introduced. Net values of intracellular fluxes and the degree of reversibility of enzymatic steps are estimated by minimization of the deviations between experimental and simulated measurements. The metabolic model applied includes the Embden-Meyerhof-Parnas and the pentose phosphate pathway, the tricarboxylic acid cycle, anaplerotic reaction sequences and pathways involved in amino acid synthesis. The study clearly demonstrates the value of complete isotopomer models for maximizing the information obtainable from 13C tracer experiments. The approach applied here offers a completely general and comprehensive analysis of carbon tracer experiments where any set of experimental data on the labeling state and extracellular fluxes can be used for the quantification of metabolic fluxes in complex metabolic networks.     

Determination of metabolic fluxes in a non-steady-state system

Estimation of fluxes through metabolic networks from redistribution patterns of (13)C has become a well developed technique in recent years. However, the approach is currently limited to systems at metabolic steady-state; dynamic changes in metabolic fluxes cannot be assessed. This is a major impediment to understanding the behaviour of metabolic networks, because steady-state is not always experimentally achievable and a great deal of information about the control hierarchy of the network can be derived from the analysis of flux dynamics. To address this issue, we have developed a method for estimating non-steady-state fluxes based on the mass-balance of mass isotopomers. This approach allows multiple mass-balance equations to be written for the change in labelling of a given metabolite pool and thereby permits over-determination of fluxes. We demonstrate how linear regression methods can be used to estimate non-steady-state fluxes from these mass balance equations. The approach can be used to calculate fluxes from both mass isotopomer and positional isotopomer labelling information and thus has general applicability to data generated from common spectrometry- or NMR-based analytical platforms. The approach is applied to a GC-MS time-series dataset of (13)C-labelling of metabolites in a heterotrophic Arabidopsis cell suspension culture. Threonine biosynthesis is used to demonstrate that non-steady-state fluxes can be successfully estimated from such data while organic acid metabolism is used to highlight some common issues that can complicate flux estimation. These include multiple pools of the same metabolite that label at different rates and carbon skeleton rearrangements.     

An optimized algorithm for flux estimation from isotopomer distribution in glucose metabolites

MOTIVATION: Analysis of the conversion of (13)C glucose within the metabolic network allows the evaluation of the biochemical fluxes in interconnecting metabolic pathways. Such analyses require solving hundreds of equations with respect to individual isotopomer concentrations, and this assumes applying special software even for constructing the equations. The algorithm, proposed by others could be improved. METHOD: A C-code linked to the program written in Mathematica (Wolfram Research Inc.), constructs and solves differential equations for all isotopomer concentrations, using the general enzyme characteristics (K(m), equilibrium constant, etc.). This code uses innovative algorithm of determination for the isotopomers-products, thus essentially decreasing the computation time. Feasible metabolic fluxes are provided by the parameters of enzyme kinetics found from the data fitting. RESULTS: The software effectively evaluates metabolic fluxes based on the measured isotopomer distribution, as was illustrated by the analysis of glycolysis and pentose phosphate cycle. The mechanism of transketolase and transaldolase catalysis was shown to induce a specific kind of isotopomer re-distribution, which, despite the significance of its effect, usually is not taken into account. AVAILABILITY: The software could be freely downloaded from the site: http://bq.ub.es/bioqint/label_distribution/.     

GC-MS analysis of amino acids rapidly provides rich information for isotopomer balancing

Gas chromatography-mass spectrometry (GC-MS) is a rapid method that provides rich information on isotopomer distributions for metabolic flux analysis. First, we established a fast and reliable experimental protocol for GC-MS analysis of amino acids from total biomass hydrolyzates, and common experimental pitfalls are discussed. Second, a suitable interface for the use of mass isotopomer data is presented. For this purpose, a general, matrix-based correction procedure that accounts for naturally occurring isotopes is introduced. Simulated and experimentally determined mass distributions of unlabeled amino acids showed a deviation of less than 3% for 89% of the analyzed fragments. Third, to investigate general properties of GC-MS-based isotopomer balancing, altered flux ratios through glycolysis and pentose phosphate pathway and/or exchange fluxes were simulated. Different fluxes were found to exert specific and significant influence on the mass isotopomer distributions, thus indicating that GC-MS data contain valuable information for metabolic flux analysis. Fourth, comparison of different methods revealed that GC-MS analysis provides the largest number of independent constraints on amino acid isotopomer abundance, followed by correlation spectroscopy and fractional enrichment analysis by nuclear magnetic resonance.     

Identification of optimal measurement sets for complete flux elucidation in metabolic flux analysis experiments

Metabolic flux analysis (MFA) methods use external flux and isotopic measurements to quantify the magnitude of metabolic flows in metabolic networks. A key question in this analysis is choosing a set of measurements that is capable of yielding a unique flux distribution (identifiability). In this article, we introduce an optimization-based framework that uses incidence structure analysis to determine the smallest (or most cost-effective) set of measurements leading to complete flux elucidation. This approach relies on an integer linear programming formulation OptMeas that allows for the measurement of external fluxes and the complete (or partial) enumeration of the isotope forms of metabolites without requiring any of these to be chosen in advance. We subsequently query and refine the measurement sets suggested by OptMeas for identifiability and optimality. OptMeas is first tested on small to medium-size demonstration examples. It is subsequently applied to a large-scale E. coli isotopomer mapping model with more than 17,000 isotopomers. A number of additional measurements are identified leading to maximum flux elucidation in an amorphadiene producing E. coli strain.     

Evaluation of regression models in metabolic physiology: predicting fluxes from isotopic data without knowledge of the pathway

This study explores the ability of regression models, with no knowledge of the underlying physiology, to estimate physiological parameters relevant for metabolism and endocrinology. Four regression models were compared: multiple linear regression (MLR), principal component regression (PCR), partial least-squares regression (PLS) and regression using artificial neural networks (ANN). The pathway of mammalian gluconeogenesis was analyzed using [U−¹³C]glucose as tracer. A set of data was simulated by randomly selecting physiologically appropriate metabolic fluxes for the 9 steps of this pathway as independent variables. The isotope labeling patterns of key intermediates in the pathway were then calculated for each set of fluxes, yielding 29 dependent variables. Two thousand sets were created, allowing independent training and test data. Regression models were asked to predict the nine fluxes, given only the 29 isotopomers. For large training sets (>50) the artificial neural network model was superior, capturing 95% of the variability in the gluconeogenic flux, whereas the three linear models captured only 75%. This reflects the ability of neural networks to capture the inherent non-linearities of the metabolic system. The effect of error in the variables and the addition of random variables to the data set was considered. Model sensitivities were used to find the isotopomers that most influenced the predicted flux values. These studies provide the first test of multivariate regression models for the analysis of isotopomer flux data. They provide insight for metabolomics and the future of isotopic tracers in metabolic research where the underlying physiology is complex or unknown.We acknowledge the support of NIH Grant DK58533 and the DuPont-MIT Alliance.     

Systematic quantification of complex metabolic flux networks using stable isotopes and mass spectrometry.

Metabolic fluxes provide a detailed metric of the cellular metabolic phenotype. Fluxes are estimated indirectly from available measurements and various methods have been developed for this purpose. Of particular interest are methods making use of stable isotopic tracers as they enable the estimation of fluxes at a high resolution. In this paper, we present data validating the use of mass spectrometry (MS) for the quantification of complex metabolic flux networks. In the context of the lysine biosynthesis flux network of Corynebacterium glutamicum (ATCC 21799) under glucose limitation in continuous culture, operating at 0.1 x h(-1) after the introduction of 50% [1-13C]glucose, we deploy a bioreaction network analysis methodology for flux determination from mass isotopomer measurements of biomass hydrolysates, while thoroughly addressing the issues of measurement accuracy, flux observability and data reconciliation. The analysis enabled the resolution of the involved anaplerotic activity of the microorganism using only one labeled substrate, the determination of the range of most of the exchange fluxes and the validation of the flux estimates through satisfaction of redundancies. Specifically, we determined that phosphoenolpyruvate carboxykinase and synthase do not carry flux at these experimental conditions and identified a high futile cycle between oxaloacetate and pyruvate, indicating a highly active in vivo oxaloacetate decarboxylase. Both results validated previous in vitro activity measurements. The flux estimates obtained passed the chi2 statistical test. This is a very important result considering that prior flux analyses of extensive metabolic networks from isotopic measurements have failed criteria of statistical consistency.     

An improved method for statistical analysis of metabolic flux analysis using isotopomer mapping matrices with analytical expressions

Analytical expressions were derived for calculating the sensitivities of isotopomer distribution vectors, the weighted output matrix with respect to the fluxes, and the covariance matrix for the metabolic flux analysis based on isotopomers mapping matrices (IMM). These expressions allow us to implement efficient statistical analysis, avoiding the time-consuming Monte Carlo techniques for estimating the confidence interval of the fluxes. The analytical expressions are also useful in implementing a faster design of experiment, which requires repetitive computation of the covariance matrix that is not straightforward to make in practice with the numerical techniques based on the conventional IMM. The proposed method was applied for analyzing the central carbon metabolism of the mixotrophically cultivated Synechocystis sp. PCC6803, and the confidence intervals of all its fluxes were computed based on the isotopomer distribution measured using NMR and GC-MS. It was found that the best feasible mixture for labeling experiment is 70% unlabeled, 10% [U-13C] and 20% [1,2-13C2] labeled glucose to obtain the most reliable metabolic fluxes.     

The relationships between the metabolic fluxes and 13C-labeled isotopomer distribution for the flux analysis of the main metabolic pathways

It has been known that ¹³C-labeling technique is quite useful in estimating the metabolic fluxes. Although the program-based flux analysis is powerful, it is not easy to be confident with the result obtained without experiences and exhaustive trial and errors based on statistical analysis for the confidence intervals in practice. It is, therefore, quite important to grasp the relationship between the fluxes and the ¹³C-labeled isotopomer distribution to get deeper insight into the metabolic flux analysis. In the present research, it was shown explicitly how the isotopomer distribution changes with respect to the fluxes in relation to the labeling patterns of the substrate, where either labeled glucose, acetate, or pyruvate was used as a carbon source. Some of the analytical expressions were derived based on the matrix representation, and they were utilized for analysis. It was shown that the isotopomer pattern does not necessarily change uniformly with respect to fluxes, but changes in a complicated way in particular for the case of using pyruvate as a carbon source where some isotopomers do not necessarily change monotonically. It was shown to be quite important to grasp how the isotopomer pattern changes with respect to fluxes and the labeling pattern of the substrate for flux determination and the experimental design. It was also shown that the mixture of [1-¹³C] acetate and [2-¹³C] acetate should not be used from the information index point of view. Some of the experimental data were evaluated from the present approach. It was also shown that the isotopomer distribution is less sensitive to the bidirectional fluxes in the reversible pathway.     

13C-Labeled metabolic flux analysis of a fed-batch culture of elutriated Saccharomyces cerevisiae

This study addresses the question of whether observable changes in fluxes in the primary carbon metabolism of Saccharomyces cerevisiae occur between the different phases of the cell division cycle. To detect such changes by metabolic flux analysis, a 13C-labeling experiment was performed with a fed-batch culture inoculated with a partially synchronized cell population obtained through centrifugal elutriation. Such a culture exhibits dynamic changes in the fractions of cells in different cell cycle phases over time. The mass isotopomer distributions of free intracellular metabolites in central carbon metabolism were measured by liquid chromatography-mass spectrometry. For four time points during the culture, these distributions were used to obtain the best estimates for the metabolic fluxes. The obtained flux fits suggested that the optimally fitted split ratio for the pentose phosphate pathway changed by almost a factor of 2 up and down around a value of 0.27 during the experiment. Statistical analysis revealed that some of the fitted flux distributions for different time points were significantly different from each other, indicating that cell cycle-dependent variations in cytosolic metabolic fluxes indeed occurred.     

Modeling and optimization of a multi-product biosynthesis factory for multiple objectives

Genetic algorithms and optimization in general, enable us to probe deeper into the metabolic pathway recipe for multi-product biosynthesis. An augmented model for optimizing serine and tryptophan flux ratios simultaneously in Escherichia coli, was developed by linking the dynamic tryptophan operon model and aromatic amino acid-tryptophan biosynthesis pathways to the central carbon metabolism model. Six new kinetic parameters of the augmented model were estimated with considerations of available experimental data and other published works. Major differences between calculated and reference concentrations and fluxes were explained. Sensitivities and underlying competition among fluxes for carbon sources were consistent with intuitive expectations based on metabolic network and previous results. Biosynthesis rates of serine and tryptophan were simultaneously maximized using the augmented model via concurrent gene knockout and manipulation. The optimization results were obtained using the elitist non-dominant sorting genetic algorithm (NSGA-II) supported by pattern recognition heuristics. A range of Pareto-optimal enzyme activities regulating the amino acids biosynthesis was successfully obtained and elucidated wherever possible vis-à-vis fermentation work based on recombinant DNA technology. The predicted potential improvements in various metabolic pathway recipes using the multi-objective optimization strategy were highlighted and discussed in detail.     

Current status of 13C-metabolic flux analysis and future perspectives

Current ¹³C-metabolic flux analysis methods were reviewed as well as the weakness of the conventional metabolic flux analysis without ¹³C-labeled experiments. Although it has been recognized that ¹³C-labeling technique is powerful in estimating the metabolic fluxes, and the program-based flux analysis is necessary, one may not be confident with the result obtained without experiences and exhaustive trial and errors in practice due to its black box nature. In the present article, we call attention to the importance of investigating the relationships between fluxes and isotopomer or mass isotopomer distributions to understand the mechanism of generating specific isotopomers. Then, the experimental design for the preferred mixture of the specific ¹³C-labeled substrate was discussed. The effect of the reversibility in the bidirectional flux on the isotopomer distribution was also mentioned, and it was shown why the reliability of the bidirectional fluxes becomes lower. Moreover, by noting that recent development of measurement techniques enables us to measure the isotopomer patterns of intracellular metabolites instead of proteinogenic amino acids, it is mentioned that this enables us to estimate the flux changes during time-variant batch culture. Some future perspectives are discussed in relation to the integration of different levels of information in the cell.     

Metabolic flux analysis in a nonstationary system: fed-batch fermentation of a high yielding strain of E. coli producing 1,3-propanediol

Metabolic fluxes estimated from stable-isotope studies provide a key to understanding cell physiology and regulation of metabolism. A limitation of the classical method for metabolic flux analysis (MFA) is the requirement for isotopic steady state. To extend the scope of flux determination from stationary to nonstationary systems, we present a novel modeling strategy that combines key ideas from isotopomer spectral analysis (ISA) and stationary MFA. Isotopic transients of the precursor pool and the sampled products are described by two parameters, D and G parameters, respectively, which are incorporated into the flux model. The G value is the fraction of labeled product in the sample, and the D value is the fractional contribution of the feed for the production of labeled products. We illustrate the novel modeling strategy with a nonstationary system that closely resembles industrial production conditions, i.e. fed-batch fermentation of Escherichia coli that produces 1,3-propanediol (PDO). Metabolic fluxes and the D and G parameters were estimated by fitting labeling distributions of biomass amino acids measured by GC/MS to a model of E. coli metabolism. We obtained highly consistent fits from the data with 82 redundant measurements. Metabolic fluxes were estimated for 20 time points during course of the fermentation. As such we established, for the first time, detailed time profiles of in vivo fluxes. We found that intracellular fluxes changed significantly during the fed-batch. The intracellular flux associated with PDO pathway increased by 10%. Concurrently, we observed a decrease in the split ratio between glycolysis and pentose phosphate pathway from 70/30 to 50/50 as a function of time. The TCA cycle flux, on the other hand, remained constant throughout the fermentation. Furthermore, our flux results provided additional insight in support of the assumed genotype of the organism.     

K-FIT: An accelerated kinetic parameterization algorithm using steady-state fluxomic data

Kinetic models predict the metabolic flows by directly linking metabolite concentrations and enzyme levels to reaction fluxes. Robust parameterization of organism-level kinetic models that faithfully reproduce the effect of different genetic or environmental perturbations remains an open challenge due to the intractability of existing algorithms. This paper introduces Kinetics-based Fluxomics Integration Tool (K-FIT), a robust kinetic parameterization workflow that leverages a novel decomposition approach to identify steady-state fluxes in response to genetic perturbations followed by a gradient-based update of kinetic parameters until predictions simultaneously agree with the fluxomic data in all perturbed metabolic networks. The applicability of K-FIT to large-scale models is demonstrated by parameterizing an expanded kinetic model for E. coli (307 reactions and 258 metabolites) using fluxomic data from six mutants. The achieved thousand-fold speed-up afforded by K-FIT over meta-heuristic approaches is transformational enabling follow-up robustness of inference analyses and optimal design of experiments to inform metabolic engineering strategies.     

Metabolic flux analysis in biotechnology processes

Metabolic flux analysis (MFA) has become a fundamental tool of metabolic engineering to elucidate the metabolic state of the cell and has been applied to various biotechnological processes. In recent years, considerable technical advances have been made. Developments of analytical instruments allow us to determine ¹³C labeling distribution of intracellular metabolites with high accuracy and sensitivity. Moreover, kinetic information of intracellular label distribution during isotopic instationary enables us to calculate metabolic fluxes with shortened experimental time and decreased amount of labeled substrate. The ¹³C MFA may be one of the most promising approaches for the target estimation to improve strain performances and production processes.