Bidirectional reaction steps in metabolic networks: II. Flux estimation and statistical analysis

Metabolic carbon labelling experiments enable a large amount of extracellular fluxes and intracellular carbon isotope enrichments to be measured. Since the relation between the measured quantities and the unknown intracellular metabolic fluxes is given by bilinear balance equations, flux determination from this data set requires the numerical solution of a nonlinear inverse problem. To this end, a general algorithm for flux estimation from metabolic carbon labelling experiments based on the least squares approach is developed in this contribution and complemented by appropriate tools for statistical analysis. The linearization technique usually applied for the computation of nonlinear confidence regions is shown to be inappropriate in the case of large exchange fluxes. For this reason a sophisticated compactification transformation technique for nonlinear statistical analysis is developed. Statistical analysis is then performed by computing appropriate statistical quality measures like output sensitivities, parameter sensitivities and the parameter covariance matrix. This allows one to determine the order of magnitude of exchange fluxes in most practical situations. An application study with a large data set from lysine-producing Corynebacterium glutamicum demonstrates the power and limitations of the carbon-labelling technique. It is shown that all intracellular fluxes in central metabolism can be quantitated without assumptions on intracellular energy yields. At the same time several exchange fluxes are determined which is invaluable information for metabolic engineering. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55: 118-135, 1997.     

Modeling isotopomer distributions in biochemical networks using isotopomer mapping matrices

Within the last decades NMR spectroscopy has undergone tremendous development and has become a powerful analytical tool for the investigation of intracellular flux distributions in biochemical networks using (13)C-labeled substrates. Not only are the experiments much easier to conduct than experiments employing radioactive tracer elements, but NMR spectroscopy also provides additional information on the labeling pattern of the metabolites. Whereas the maximum amount of information obtainable with (14)C-labeled substrates is the fractional enrichment in the individual carbon atom positions, NMR spectroscopy can also provide information on the degree of labeling at neighboring carbon atom positions by analyzing multiplet patterns in NMR spectra or using 2-dimensional NMR spectra. It is possible to quantify the mole fractions of molecules that show a specific labeling pattern, i.e., information of the isotopomer distribution in metabolite pools can be obtained. The isotopomer distribution is the maximum amount of information that in theory can be obtained from (13)C-tracer studies. The wealth of information contained in NMR spectra frequently leads to overdetermined algebraic systems. Consequently, fluxes must be estimated by nonlinear least squares analysis, in which experimental labeling data is compared with simulated steady state isotopomer distributions. Hence, mathematical models are required to compute the steady state isotopomer distribution as a function of a given set of steady state fluxes. Because 2(n) possible labeling patterns exist in a molecule of n carbon atoms, and each pattern corresponds to a separate state in the isotopomer model, these models are inherently complex. Model complexity, so far, has restricted usage of isotopomer information to relatively small metabolic networks. A general methodology for the formulation of isotopomer models is described. The model complexity of isotopomer models is reduced to that of classical metabolic models by expressing the 2(n) isotopomer mass balances of a metabolite pool in a single matrix equation. Using this approach an isotopomer model has been implemented that describes label distribution in primary carbon metabolism, i.e., in a metabolic network including the Embden-Meyerhof-Parnas and pentose phosphate pathway, the tricarboxylic acid cycle, and selected anaplerotic reaction sequences. The model calculates the steady state label distribution in all metabolite pools as a function of the steady state fluxes and is applied to demonstrate the effect of selected anaplerotic fluxes on the labeling pattern of the pathway intermediates. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 55:831-840, 1997.     

13C tracer experiments and metabolite balancing for metabolic flux analysis: comparing two approaches

Conventional metabolic flux analysis uses the information gained from determination of measurable fluxes and a steady-state assumption for intracellular metabolites to calculate the metabolic fluxes in a given metabolic network. The determination of intracellular fluxes depends heavily on the correctness of the assumed stoichiometry including the presence of all reactions with a noticeable impact on the model metabolite balances. Determination of fluxes in complex metabolic networks often requires the inclusion of NADH and NADPH balances, which are subject to controversial debate. Transhydrogenation reactions that transfer reduction equivalents from NADH to NADPH or vice versa can usually not be included in the stoichiometric model, because they result in singularities in the stoichiometric matrix. However, it is the NADPH balance that, to a large extent, determines the calculated flux through the pentose phosphate pathway. Hence, wrong assumptions on the presence or activity of transhydrogenation reactions will result in wrong estimations of the intracellular flux distribution. Using 13C tracer experiments and NMR analysis, flux analysis can be performed on the basis of only well established stoichiometric equations and measurements of the labeling state of intracellular metabolites. Neither NADH/NADPH balancing nor assumptions on energy yields need to be included to determine the intracellular fluxes. Because metabolite balancing methods and the use of 13C labeling measurements are two different approaches to the determination of intracellular fluxes, both methods can be used to verify each other or to discuss the origin and significance of deviations in the results. Flux analysis based entirely on metabolite balancing and flux analysis, including labeling information, have been performed independently for a wild-type strain of Aspergillus oryzae producing alpha-amylase. Two different nitrogen sources, NH4+ and NO3-, have been used to investigate the influence of the NADPH requirements on the intracellular flux distribution. The two different approaches to the calculation of fluxes are compared and deviations in the results are discussed. Copyright 1998 John Wiley & Sons, Inc.     

Metabolic flux analysis of CHO cell metabolism in the late non‐growth phase

Chinese hamster ovary (CHO) cell cultures are commonly used for production of recombinant human therapeutic proteins. Often the goal of such a process is to separate the growth phase of the cells, from the non‐growth phase where ideally the cells are diverting resources to produce the protein of interest. Characterizing the way that the cells use nutrients in terms of metabolic fluxes as a function of culture conditions can provide a deeper understanding of the cell biology offering guidance for process improvements. To evaluate the fluxes, metabolic flux analysis of the CHO cell culture in the non‐growth phase was performed by a combination of steady‐state isotopomer balancing and stoichiometric modeling. Analysis of the glycolytic pathway and pentose phosphate pathway (PPP) indicated that almost all of the consumed glucose is diverted towards PPP with a high NADPH production; with even recycle from PPP to G6P in some cases. Almost all of the pyruvate produced from glycolysis entered the TCA cycle with little or no lactate production. Comparison of the non‐growth phase against previously reported fluxes from growth phase cultures indicated marked differences in the fluxes, in terms of the split between glycolysis and PPP, and also around the pyruvate node. Possible reasons for the high NADPH production are also discussed. Evaluation of the fluxes indicated that the medium strength, carbon dioxide level, and temperature with dissolved oxygen have statistically significant impacts on different nodes of the flux network. Biotechnol. Bioeng. 2011; 108:82–92. © 2010 Wiley Periodicals, Inc.     

Response of the central metabolism of Corynebacterium glutamicum to different flux burdens

To evaluate the importance of reactions within the central metabolism under different flux burdens the fluxes within the pentose phosphate pathway (PPP), as well as the other reactions of the central metabolism, were intensively analyzed and quantitated. For this purpose, Corynebacterium glutamicum was grown with [1-(13)C]glucose to metabolic and isotopic steady state and the fractional enrichments in precursor metabolites (e.g., pentose 5-phosphate) were quantified. Matrix calculus was used to express these data together with metabolite mass data. The detailed analysis of the dependence of (13)C enrichments on exchange fluxes enabled the transketolase-catalyzed exchange rate (2 pentose 5-phosphate <--> sedoheptulose 7-phosphate + glyceraldehyde 3-phosphate) to be quantified as 74.3% (molar metabolite flux) at a net flux of 10.3% and the exchange rate (pentose 5-phosphate + erythrose 4-phosphate <--> fructose 6-phosphate + glyceraldehyde 3-phosphate) to be quantified as 5.6% at a net flux of 8.1%. The flux entering the tricarboxylic acid cycle was 93.3%. The same comprehensive flux analysis as performed for the nonexcreting condition was done with the identical strain that had been forced to excrete L-glutamate. Because we had already quantified the fluxes for L-lysine excretion with an isogenic strain, three directly comparable flux situations are thus available. Consequently, this comparison permits a direct cause-and-effect relationship to be specified. In response to the different flux burdens of the cell, the PPP flux decreased from a maximum of 67% to 26%, with the glycolytic flux increasing accordingly. The carbon flux through isocitrate dehydrogenase increased from 20% to 36%. The bidirectional carbon flux between pyruvate and oxaloacetate decreased from 36% to 9%. Since the cause of the three different flux states was the allelic exchange in the final L-lysine assembling pathway or the glutamate export activity, respectively, the flexible response is the effect. This shows conclusively the enormous flexibility within the central metabolism of C. glutamicum to supply precursors upon their withdrawal for the synthesis of amino acids. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 56: 168-180, 1997.     

Possible pitfalls of flux calculations based on (13)C-labeling

Metabolic engineers have enthusiastically adopted the (13)C-labeling technique as a powerful tool for elucidating fluxes in metabolic networks. This tracer technique makes it possible to determine fluxes that are unobservable using only metabolite balances and allows the elimination of doubtful cofactor balances that are indispensable in flux analysis based on metabolite balancing alone. The (13)C-labeling technique, however, relies on a number of assumptions that are not free from uncertainties. Two possible errors in the models that are needed to determine the metabolic fluxes from labeling data are omitted reactions and ignored occurrence of channeling. By means of two representative examples it is shown that these modeling errors may lead to serious errors in the calculated flux distributions despite the use of labeling data. A complicating fact is that the model errors are not always easily detected as poor models may still yield good fits of experimental data. Results of (13)C-labeling experiments should therefore be interpreted with appropriate caution.     

Metabolic isotopomer labeling systems. Part II: structural flux identifiability analysis

Metabolic flux analysis using carbon labeling experiments (CLEs) is an important tool in metabolic engineering where the intracellular fluxes have to be computed from the measured extracellular fluxes and the partially measured distribution of 13C labeling within the intracellular metabolite pools. The relation between unknown fluxes and measurements is described by an isotopomer labeling system (ILS) (see Part I [Math. Biosci. 169 (2001) 173]). Part II deals with the structural flux identifiability of measured ILSs in the steady state. The central question is whether the measured data contains sufficient information to determine the unknown intracellular fluxes. This question has to be decided a priori, i.e. before the CLE is carried out. In structural identifiability analysis the measurements are assumed to be noise-free. A general theory of structural flux identifiability for measured ILSs is presented and several algorithms are developed to solve the identifiability problem. In the particular case of maximal measurement information, a symbolical algorithm is presented that decides the identifiability question by means of linear methods. Several upper bounds of the number of identifiable fluxes are derived, and the influence of the chosen inputs is evaluated. By introducing integer arithmetic this algorithm can even be applied to large networks. For the general case of arbitrary measurement information, identifiability is decided by a local criterion. A new algorithm based on integer arithmetic enables an a priori local identifiability analysis to be performed for networks of arbitrary size. All algorithms have been implemented and flux identifiability is investigated for the network of the central metabolic pathways of a microorganism. Moreover, several small examples are worked out to illustrate the influence of input metabolite labeling and the paradox of information loss due to network simplification.     

A reliable simulator for dynamic flux balance analysis

Dynamic flux balance analysis (DFBA) provides a platform for detailed design, control and optimization of biochemical process technologies. It is a promising modeling framework that combines genome‐scale metabolic network analysis with dynamic simulation of the extracellular environment. Dynamic flux balance analysis assumes that the intracellular species concentrations are in equilibrium with the extracellular environment. The resulting underdetermined stoichiometric model is solved under the assumption of a biochemical objective such as growth rate maximization. The model of the metabolism is coupled with the dynamic mass balance equations of the extracellular environment via expressions for the rates of substrate uptake and product excretion, which imposes additional constraints on the linear program (LP) defined by growth rate maximization of the metabolism. The linear program is embedded into the dynamic model of the bioreactor, and together with the additional constraints this provides an accurate model of the substrate consumption, product secretion, and biomass production during operation. A DFBA model consists of a system of ordinary differential equations for which the evaluation of the right‐hand side requires not only function evaluations, but also the solution of one or more linear programs. The numerical tool presented here accurately and efficiently simulates large‐scale dynamic flux balance models. The main advantages that this approach has over existing implementation are that the integration scheme has a variable step size, that the linear program only has to be solved when qualitative changes in the optimal flux distribution of the metabolic network occur, and that it can reliably simulate behavior near the boundary of the domain where the model is defined. This is illustrated through large‐scale examples taken from the literature. Biotechnol. Bioeng. 2013; 110: 792–802. © 2012 Wiley Periodicals, Inc.     

Improvements in metabolic flux analysis using carbon bond labeling experiments: bondomer balancing and Boolean function mapping

The biosynthetically directed fractional (13)C labeling method for metabolic flux evaluation relies on performing a 2-D [(13)C, (1)H] NMR experiment on extracts from organisms cultured on a uniformly labeled carbon substrate. This article focuses on improvements in the interpretation of data obtained from such an experiment by employing the concept of bondomers. Bondomers take into account the natural abundance of (13)C; therefore many bondomers in a real network are zero, and can be precluded a priori--thus resulting in fewer balances. Using this method, we obtained a set of linear equations which can be solved to obtain analytical formulas for NMR-measurable quantities in terms of fluxes in glycolysis and the pentose phosphate pathways. For a specific case of this network with four degrees of freedom, a priori identifiability of the fluxes was shown possible for any set of fluxes. For a more general case with five degrees of freedom, the fluxes were shown identifiable for a representative set of fluxes. Minimal sets of measurements which best identify the fluxes are listed. Furthermore, we have delineated Boolean function mapping, a new method to iteratively simulate bondomer abundances or efficiently convert carbon skeleton rearrangement information to mapping matrices. The efficiency of this method is expected to be valuable while analyzing metabolic networks which are not completely known (such as in plant metabolism) or while implementing iterative bondomer balancing methods.     

Bioreaction network topology and metabolic flux ratio analysis by biosynthetic fractional 13C labeling and two-dimensional NMR spectroscopy

Biosynthetically directed fractional 13C labeling of the proteinogenic amino acids is achieved by feeding a mixture of uniformly 13C-labeled and unlabeled carbon source compounds into a bioreaction network. Analysis of the resulting labeling pattern enables both a comprehensive characterization of the network topology and the determination of metabolic flux ratios. Attractive features with regard to routine applications are (i) an inherently small demand for 13C-labeled source compounds and (ii) the high sensitivity of two-dimensional [13C,1H]-correlation nuclear magnetic resonance spectroscopy for analysis of 13C-labeling patterns. A user-friendly program, FCAL, is available to allow rapid data analysis. This novel approach, which recently also has been employed in conjunction with metabolic flux balancing to obtain reliable estimates of in vivo fluxes, enables efficient support of metabolic engineering and biotechnology process design.     

High-throughput metabolic flux analysis based on gas chromatography-mass spectrometry derived 13C constraints

13C-constrained flux balancing analysis based on gas chromatography-mass spectrometry data is presented here as a simple and robust method for the estimation of intracellular carbon fluxes. In this approach, the underdetermined system of metabolite balances deduced from stoichiometric relations and measured extracellular rates is complemented with 13C constraints from metabolic flux ratio analysis. Fluxes in central carbon metabolism of exponentially growing Escherichia coli were estimated by 13C-constrained flux balancing from three different 13C-labeled glucose experiments. The best resolution of the network was achieved using 13C constraints derived from [U-13C]glucose and [1-13C]glucose experiments. The corresponding flux estimate was in excellent agreement with a solution that was independently obtained with a comprehensive isotopomer model. This new methodology was also demonstrated to faithfully capture the intracellular flux distribution in E. coli shake flasks and 1-ml deep-well microtiter plates. Due to its simplicity, speed, and robustness, 13C-constrained metabolic flux balancing is promising for routine and high-throughput analysis on a miniaturized scale.     

Dynamic flux responses in riboflavin overproducing Bacillus subtilis to increasing glucose limitation in fed‐batch culture

How do intracellular fluxes respond to dynamically increasing glucose limitation when the physiology changes from strong overflow metabolism near to exclusively maintenance metabolism? Here we investigate this question in a typical industrial, glucose‐limited fed‐batch cultivation with a riboflavin overproducing Bacillus subtilis strain. To resolve dynamic flux changes, a novel approach to ¹³C flux analysis was developed that is based on recording ¹³C labeling patterns in free intracellular amino acids. Fluxes are then estimated with stationary flux ratio and iterative isotopomer balancing methods, for which a decomposition of the process into quasi‐steady states and estimation of isotopic steady state ¹³C labeling patterns was necessary. By this approach, we achieve a temporal resolution of 30–60 min that allows us to resolve the slow metabolic transients that typically occur in such cultivations. In the late process phase we found, most prominently, almost exclusive respiratory metabolism, significantly increased pentose phosphate pathway contribution and a strongly decreased futile cycle through the PEP carboxykinase. As a consequence, higher catabolic NADPH formation occurred than was necessary to satisfy the anabolic demands, suggesting a transhydrogenase‐like mechanism to close the balance of reducing equivalents. Biotechnol. Bioeng. 2010. 105: 795–804. © 2009 Wiley Periodicals, Inc.     

A transient isotopic labeling methodology for 13C metabolic flux analysis of photoautotrophic microorganisms

Metabolic flux analysis is increasingly recognized as an integral component of systems biology. However, techniques for experimental measurement of system-wide metabolic fluxes in purely photoautotrophic systems (growing on CO(2) as the sole carbon source) have not yet been developed due to the unique problems posed by such systems. In this paper, we demonstrate that an approach that balances positional isotopic distributions transiently is the only route to obtaining system-wide metabolic flux maps for purely autotrophic metabolism. The outlined transient (13)C-MFA methodology enables measurement of fluxes at a metabolic steady-state, while following changes in (13)C-labeling patterns of metabolic intermediates as a function of time, in response to a step-change in (13)C-label input. We use mathematical modeling of the transient isotopic labeling patterns of central intermediates to assess various experimental requirements for photoautotrophic MFA. This includes the need for intracellular metabolite concentration measurements and isotopic labeling measurements as a function of time. We also discuss photobioreactor design and operation in order to measure fluxes under precise environmental conditions. The transient MFA technique can be used to measure and compare fluxes under different conditions of light intensity, nitrogen sources or compare strains with various mutations or gene deletions and additions.     

Metabolic flux analysis at ultra short time scale: isotopically non-stationary 13C labeling experiments

A novel approach to (13)C metabolic flux analysis (MFA) is presented using cytosolic metabolite pool sizes and their (13)C labeling data from an isotopically non-stationary (13)C labeling experiment (INST-CLE). The procedure is demonstrated with an E. coli wild type strain grown at fed batch conditions. The intra cellular labeling dynamics are excited by a sudden step increase of the (13)C portion in the substrate feed. Due to unchanged saturation of the substrate uptake system, the metabolic fluxes remain constant during the following sampling time period of only 16s, in which 20 samples are taken by an automated rapid sampling device immediately stopping metabolism by methanol quenching. Subsequent cell disruptive sample preparation and LC-MS/MS enabled simultaneous determination of pool sizes and mass isotopomers of intra cellular metabolites requiring detection limits in the nM range. Based on this data the new computational flux analysis tool 13CFLUX/INST is used to determine the intra cellular fluxes based on a complex carbon labeling network model. The measured data is in good agreement with the model predictions, thus proving the applicability of the new isotopically non-stationary (13)C metabolic flux analysis (INST-(13)C-MFA) concept. Moreover, it is shown that significant new information with respect to flux identifiability, non-measurable pool sizes, data consistency, or large storage pools can be taken from the novel kind of experimental data. This offers new insight into the biological operation of the metabolic network in vivo.     

Efficient Modeling of MS/MS Data for Metabolic Flux Analysis

Metabolic flux analysis (MFA) is a widely used method for quantifying intracellular metabolic fluxes. It works by feeding cells with isotopic labeled nutrients, measuring metabolite isotopic labeling, and computationally interpreting the measured labeling data to estimate flux. Tandem mass-spectrometry (MS/MS) has been shown to be useful for MFA, providing positional isotopic labeling data. Specifically, MS/MS enables the measurement of a metabolite tandem mass-isotopomer distribution, representing the abundance in which certain parent and product fragments of a metabolite have different number of labeled atoms. However, a major limitation in using MFA with MS/MS data is the lack of a computationally efficient method for simulating such isotopic labeling data. Here, we describe the tandemer approach for efficiently computing metabolite tandem mass-isotopomer distributions in a metabolic network, given an estimation of metabolic fluxes. This approach can be used by MFA to find optimal metabolic fluxes, whose induced metabolite labeling patterns match tandem mass-isotopomer distributions measured by MS/MS. The tandemer approach is applied to simulate MS/MS data in a small-scale metabolic network model of mammalian methionine metabolism and in a large-scale metabolic network model of E. coli. It is shown to significantly improve the running time by between two to three orders of magnitude compared to the state-of-the-art, cumomers approach. We expect the tandemer approach to promote broader usage of MS/MS technology in metabolic flux analysis. Implementation is freely available at www.cs.technion.ac.il/~tomersh/methods.html     

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.     

A Method to Constrain Genome-Scale Models with 13C Labeling Data

Current limitations in quantitatively predicting biological behavior hinder our efforts to engineer biological systems to produce biofuels and other desired chemicals. Here, we present a new method for calculating metabolic fluxes, key targets in metabolic engineering, that incorporates data from ¹³C labeling experiments and genome-scale models. The data from ¹³C labeling experiments provide strong flux constraints that eliminate the need to assume an evolutionary optimization principle such as the growth rate optimization assumption used in Flux Balance Analysis (FBA). This effective constraining is achieved by making the simple but biologically relevant assumption that flux flows from core to peripheral metabolism and does not flow back. The new method is significantly more robust than FBA with respect to errors in genome-scale model reconstruction. Furthermore, it can provide a comprehensive picture of metabolite balancing and predictions for unmeasured extracellular fluxes as constrained by ¹³C labeling data. A comparison shows that the results of this new method are similar to those found through ¹³C Metabolic Flux Analysis (¹³C MFA) for central carbon metabolism but, additionally, it provides flux estimates for peripheral metabolism. The extra validation gained by matching 48 relative labeling measurements is used to identify where and why several existing COnstraint Based Reconstruction and Analysis (COBRA) flux prediction algorithms fail. We demonstrate how to use this knowledge to refine these methods and improve their predictive capabilities. This method provides a reliable base upon which to improve the design of biological systems.