Metabolic flux elucidation for large-scale models using 13C labeled isotopes

Tracking the Metabolic Fate of Exogenous Arachidonic Acid in Ferroptosis Using Dual-Isotope Labeling Lipidomics

Lipid metabolism is implicated in a variety of diseases, including cancer, cell death, and inflammation, but lipidomics has proven to be challenging due to the vast structural diversity over a narrow range of mass and polarity of lipids. Isotope labeling is often used in metabolomics studies to follow the metabolism of exogenously added labeled compounds because they can be differentiated from endogenous compounds by the mass shift associated with the label. The application of isotope labeling to lipidomics has also been explored as a method to track the metabolism of lipids in various disease states. However, it can be difficult to differentiate a single isotopically labeled lipid from the rest of the lipidome due to the variety of endogenous lipids present over the same mass range. Here we report the development of a dual-isotope deuterium labeling method to track the metabolic fate of exogenous polyunsaturated fatty acids, e.g., arachidonic acid, in the context of ferroptosis using hydrophilic interaction-ion mobility-mass spectrometry (HILIC-IM-MS). Ferroptosis is a type of cell death that is dependent on lipid peroxidation. The use of two isotope labels rather than one enables the identification of labeled species by a signature doublet peak in the resulting mass spectra. A Python-based software, D-Tracer, was developed to efficiently extract metabolites with dual-isotope labels. The labeled species were then identified with LiPydomics based on their retention times, collision cross section, and m/z values. Changes in exogenous AA incorporation in the absence and presence of a ferroptosis inducer were elucidated.

Interpreting metabolic complexity via isotope-assisted metabolic flux analysis

Isotope-assisted metabolic flux analysis (iMFA) is a powerful method to mathematically determine the metabolic fluxome from experimental isotope labeling data and a metabolic network model. While iMFA was originally developed for industrial biotechnological applications, it is increasingly used to analyze eukaryotic cell metabolism in physiological and pathological states. In this review, we explain how iMFA estimates the intracellular fluxome, including data and network model (inputs), the optimization-based data fitting (process), and the flux map (output). We then describe how iMFA enables analysis of metabolic complexities and discovery of metabolic pathways. Our goal is to expand the use of iMFA in metabolism research, which is essential to maximizing the impact of metabolic experiments and continuing to advance iMFA and biocomputational techniques.

CeCaFLUX: the first web server for standardized and visual instationary 13C metabolic flux analysis

SUMMARY

The number of instationary 13C-metabolic flux (INST-MFA) studies grows every year, making it more important than ever to ensure the clarity, standardization and reproducibility of each study. We proposed CeCaFLUX, the first user-friendly web server that derives metabolic flux distribution from instationary 13C-labeled data. Flux optimization and statistical analysis are achieved through an evolutionary optimization in a parallel manner. It can visualize the flux optimizing process in real time and the ultimate flux outcome. It will also function as a database to enhance the consistency and to facilitate sharing of flux studies.

AVAILABILITY AND IMPLEMENTATION

CeCaFLUX is freely available at https://www.cecaflux.net, the source code can be downloaded at https://github.com/zhzhd82/CeCaFLUX.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Discovering Essential Multiple Gene Effects Through Large Scale Optimization: An Application to Human Cancer Metabolism

Computational modelling of metabolic processes has proven to be a useful approach to formulate our knowledge and improve our understanding of core biochemical systems that are crucial to maintaining cellular functions. Towards understanding the broader role of metabolism on cellular decision-making in health and disease conditions, it is important to integrate the study of metabolism with other core regulatory systems and omics within the cell, including gene expression patterns. After quantitatively integrating gene expression profiles with a genome-scale reconstruction of human metabolism, we propose a set of combinatorial methods to reverse engineer gene expression profiles and to find pairs and higher-order combinations of genetic modifications that simultaneously optimize multi-objective cellular goals. This enables us to suggest classes of transcriptomic profiles that are most suitable to achieve given metabolic phenotypes. We demonstrate how our techniques are able to compute beneficial, neutral or "toxic" combinations of gene expression levels. We test our methods on nine tissue-specific cancer models, comparing our outcomes with the corresponding normal cells, identifying genes as targets for potential therapies. Our methods open the way to a broad class of applications that require an understanding of the interplay among genotype, metabolism, and cellular behaviour, at scale.

Comprehensive Isotopic Targeted Mass Spectrometry: Reliable Metabolic Flux Analysis with Broad Coverage

Metabolic flux analysis (MFA) is highly relevant to understanding metabolic mechanisms of various biological processes. While the pace of methodology development in MFA has been rapid, a major challenge the field continues to witness is limited metabolite coverage, often restricted to a small to moderate number of well-known compounds. In addition, isotopic peaks from an enriched metabolite tend to have low abundances, which makes liquid chromatography tandem mass spectrometry (LC-MS/MS) highly useful in MFA due to its high sensitivity and specificity. Previously we have built large-scale LC-MS/MS approaches that can be routinely used for measurement of up to ∼1,900 metabolite/feature levels [Gu et al. Anal. Chem. 2015, 87, 12355-12362. Shi et al. Anal. Chem. 2019, 91, 13737-13745.]. In this study, we aim to expand our previous studies focused on metabolite level measurements to flux analysis and establish a novel comprehensive isotopic targeted mass spectrometry (CIT-MS) method for reliable MFA analysis with broad coverage. As a proof-of-principle, we have applied CIT-MS to compare the steady-state enrichment of metabolites between Myc(oncogene)-On and Myc-Off Tet21N human neuroblastoma cells cultured with U-13C6-glucose medium. CIT-MS is operationalized using multiple reaction monitoring (MRM) mode and is able to perform MFA of 310 identified metabolites (142 reliably detected, 46 kinetically profiled) selected from >35 metabolic pathways of strong biological significance. Further, we developed a novel concept of relative flux, which eliminates the requirement of absolute quantitation in traditional MFA and thus enables comparative MFA under the pseudosteady state. As a result, CIT-MS was shown to possess the advantages of broad coverage, easy implementation, fast throughput, and more importantly, high fidelity and accuracy in MFA. In principle, CIT-MS can be easily adapted to track the flux of other labeled tracers (such as 15N-tracers) in any metabolite detectable by LC-MS/MS and in various biological models (such as mice). Therefore, CIT-MS has great potential to bring new insights to both basic and clinical metabolism research.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For ¹³C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

Genome-Scale 13C Fluxomics Modeling for Metabolic Engineering of Saccharomyces cerevisiae

Synthetic biology is a rapidly developing field that pursues the application of engineering principles and development approaches to biological engineering. Synthetic biology is poised to change the way biology is practiced, and has important practical applications: for example, building genetically engineered organisms to produce biofuels, medicines, and other chemicals. Traditionally, synthetic biology has focused on manipulating a few genes (e.g., in a single pathway or genetic circuit), but its combination with systems biology holds the promise of creating new cellular architectures and constructing complex biological systems from the ground up. Enabling this merge of synthetic and systems biology will require greater predictive capability for modeling the behavior of cellular systems, and more comprehensive data sets for building and calibrating these models. The so-called "-omics" data sets can now be generated via high throughput techniques in the form of genomic, proteomic, transcriptomic, and metabolomic information on the engineered biological system. Of particular interest with respect to the engineering of microbes capable of producing biofuels and other chemicals economically and at scale are metabolomic datasets, and their insights into intracellular metabolic fluxes. Metabolic fluxes provide a rapid and easy to understand picture of how carbon and energy flow throughout the cell. Here, we present a detailed guide to performing metabolic flux analysis and modeling using the open source JBEI Quantitative Metabolic Modeling (jQMM) library. This library allows the user to transform metabolomics data in the form of isotope labeling data from a ¹³C labeling experiment into a determination of cellular fluxes that can be used to develop genetic engineering strategies for metabolic engineering.The jQMM library presents a complete toolbox for performing a range of different tasks of interest in metabolic engineering. Various different types of flux analysis and modeling can be performed such as flux balance analysis, ¹³C metabolic flux analysis, and two-scale ¹³C metabolic flux analysis (2S-¹³C MFA). 2S-¹³C MFA is a novel method that determines genome-scale fluxes without the need of every single carbon transition in the metabolic network. In addition to several other capabilities, the jQMM library can make model based predictions for how various genetic engineering strategies can be incorporated toward bioengineering goals: it can predict the effects of reaction knockouts on metabolism using both the MoMA and ROOM methodologies. In this chapter, we will illustrate the use of the jQMM library through a step-by-step demonstration of flux determination and knockout prediction in a complex eukaryotic model organism: Saccharomyces cerevisiae (S. cerevisiae). Included with this chapter is a digital Jupyter Notebook file that provides a computable appendix showing a self-contained example of jQMM usage, which can be changed to fit the user's specific needs. As an open source software project, users can modify and extend the code base to make improvements at will, allowing them to share their development work and contribute back to the jQMM modeling community.

Prediction of metabolic fluxes from gene expression data with Huber penalty convex optimization function

As one of the critical parameters of a metabolic pathway, the metabolic flux in a metabolic network serves as an essential role in physiology and pathology. Constraint-based metabolic models are the widely used frameworks for predicting metabolic fluxes in genome-scale metabolic networks. Integrating the transcriptomic data into the constraint-based metabolic models can effectively predict context-specific fluxes across different conditions. However, these methods always need user-defined thresholds to identify the expression levels of metabolic genes or restrain the rate of biomass production, and the predictive results are sensitive to the thresholds. In this work, we present the Huber penalty convex optimization function (HPCOF) combined with the flux minimization principle to predict metabolic fluxes. Our HPCOF method integrates gene expression profiles into the genome-scale metabolic models (GEMs) to reduce the sensitivity to outliers, and uses continuous expression data to avoid selection of arbitrary threshold parameters. In the case studies of Saccharomyces cerevisiae (S. cerevisiae) and Escherichia coli (E. coli) strains under different conditions, the results show that our HPCOF method has a better fit to the experimentally measured values, and has a higher Pearson correlation coefficient, a smaller P-value and a lower sum of squared error than other methods. The HPCOF code can be freely downloaded from for academic users.

Two-Scale 13C Metabolic Flux Analysis for Metabolic Engineering

Accelerating the Design-Build-Test-Learn (DBTL) cycle in synthetic biology is critical to achieving rapid and facile bioengineering of organisms for the production of, e.g., biofuels and other chemicals. The Learn phase involves using data obtained from the Test phase to inform the next Design phase. As part of the Learn phase, mathematical models of metabolic fluxes give a mechanistic level of comprehension to cellular metabolism, isolating the principle drivers of metabolic behavior from the peripheral ones, and directing future experimental designs and engineering methodologies. Furthermore, the measurement of intracellular metabolic fluxes is specifically noteworthy as providing a rapid and easy-to-understand picture of how carbon and energy flow throughout the cell. Here, we present a detailed guide to performing metabolic flux analysis in the Learn phase of the DBTL cycle, where we show how one can take the isotope labeling data from a ¹³C labeling experiment and immediately turn it into a determination of cellular fluxes that points in the direction of genetic engineering strategies that will advance the metabolic engineering process.For our modeling purposes we use the Joint BioEnergy Institute (JBEI) Quantitative Metabolic Modeling (jQMM) library, which provides an open-source, python-based framework for modeling internal metabolic fluxes and making actionable predictions on how to modify cellular metabolism for specific bioengineering goals. It presents a complete toolbox for performing different types of flux analysis such as Flux Balance Analysis, ¹³C Metabolic Flux Analysis, and it introduces the capability to use ¹³C labeling experimental data to constrain comprehensive genome-scale models through a technique called two-scale ¹³C Metabolic Flux Analysis (2S-¹³C MFA) [1]. In addition to several other capabilities, the jQMM is also able to predict the effects of knockouts using the MoMA and ROOM methodologies. The use of the jQMM library is illustrated through a step-by-step demonstration, which is also contained in a digital Jupyter Notebook format that enhances reproducibility and provides the capability to be adopted to the user's specific needs. As an open-source software project, users can modify and extend the code base and make improvements at will, providing a base for future modeling efforts.

Simultaneous tracing of carbon and nitrogen isotopes in human cells

Stable isotope tracing is a powerful method for interrogating metabolic enzyme activities across the metabolic network of living cells. However, most studies of mammalian cells have used (13)C-labeled tracers only and focused on reactions in central carbon metabolism. Cellular metabolism, however, involves other biologically important elements, including nitrogen, hydrogen, oxygen, phosphate and sulfur. Tracing stable isotopes of such elements may help shed light on poorly understood metabolic pathways. Here, we demonstrate the use of high-resolution mass spectrometry to simultaneously trace carbon and nitrogen metabolism in human cells cultured with (13)C- and (15)N-labeled glucose and glutamine. To facilitate interpretation of the complex isotopomer data generated, we extend current methods for metabolic flux analysis to handle multivariate mass isotopomer distributions (MMIDs). We find that observed MMIDs are broadly consistent with known biochemical pathways. Whereas measured (13)C MIDs were informative for central carbon metabolism, (15)N isotopes provided evidence for nitrogen-carrying reactions in amino acid and nucleotide metabolism. This computational and experimental methodology expands the scope of metabolic flux analysis beyond carbon metabolism, and may prove important to understanding metabolic phenotypes in health and disease.

An analytic approximation of the feasible space of metabolic networks

Assuming a steady-state condition within a cell, metabolic fluxes satisfy an underdetermined linear system of stoichiometric equations. Characterizing the space of fluxes that satisfy such equations along with given bounds (and possibly additional relevant constraints) is considered of utmost importance for the understanding of cellular metabolism. Extreme values for each individual flux can be computed with linear programming (as flux balance analysis), and their marginal distributions can be approximately computed with Monte Carlo sampling. Here we present an approximate analytic method for the latter task based on expectation propagation equations that does not involve sampling and can achieve much better predictions than other existing analytic methods. The method is iterative, and its computation time is dominated by one matrix inversion per iteration. With respect to sampling, we show through extensive simulation that it has some advantages including computation time, and the ability to efficiently fix empirically estimated distributions of fluxes.

Large-scale metabolic models of organisms from microbes to mammals can provide great insight into cellular function, but their analysis remains challenging. Here, the authors provide an approximate analytic method to estimate the feasible solution space for the flux vectors of metabolic networks, enabling more accurate analysis under a wide range of conditions of interest.

Estimation of flux ratios without uptake or release data: Application to serine and methionine metabolism

Model-based metabolic flux analysis (MFA) using isotope-labeled substrates has provided great insight into intracellular metabolic activities across a host of organisms. One challenge with applying MFA in mammalian systems, however, is the need for absolute quantification of nutrient uptake, biomass composition, and byproduct release fluxes. Such measurements are often not feasible in complex culture systems or in vivo. One way to address this issue is to estimate flux ratios, the fractional contribution of a flux to a metabolite pool, which are independent of absolute measurements and yet informative for cellular metabolism. Prior work has focused on "local" estimation of a handful of flux ratios for specific metabolites and reactions. Here, we perform systematic, model-based estimation of all flux ratios in a metabolic network using isotope labeling data, in the absence of uptake/release data. In a series of examples, we investigate what flux ratios can be well estimated with reasonably tight confidence intervals, and contrast this with confidence intervals on normalized fluxes. We find that flux ratios can provide useful information on the metabolic state, and is complementary to normalized fluxes: for certain metabolic reactions, only flux ratios can be well estimated, while for others normalized fluxes can be obtained. Simulation studies of a large human metabolic network model suggest that estimation of flux ratios is technically feasible for complex networks, but additional studies on data from actual isotopomer labeling experiments are needed to validate these results. Finally, we experimentally study serine and methionine metabolism in cancer cells using flux ratios. We find that, in these cells, the methionine cycle is truncated with little remethylation from homocysteine, and polyamine synthesis in the absence of methionine salvage leads to loss of 5-methylthioadenosine, suggesting a new mode of overflow metabolism in cancer cells. This work highlights the potential for flux ratio analysis in the absence of absolute quantification, which we anticipate will be important for both in vitro and in vivo studies of cancer metabolism.

Stable isotope studies reveal pathways for the incorporation of non-essential amino acids in Acyrthosiphon pisum (pea aphids)

Plant roots incorporate inorganic nitrogen into the amino acids glutamine, glutamic acid, asparagine and aspartic acid, which together serve as the primary metabolites of nitrogen transport to other tissues. Given the preponderance of these four amino acids, phloem sap is a nutritionally unbalanced diet for phloem-feeding insects. Therefore, aphids and other phloem feeders typically rely on microbial symbionts for the synthesis of essential amino acids. To investigate the metabolism of the four main transport amino acids by the pea aphid (Acyrthosiphon pisum), and its Buchnera aphidicola endosymbionts, aphids were fed defined diets with stable isotope-labeled glutamine, glutamic acid, asparagine or aspartic acid (U-(13)C, U-(15)N; U-(15)N; α-(15)N; or γ-(15)N). The metabolic fate of the dietary (15)N and (13)C was traced using gas chromatography-mass spectrometry (GC-MS). Nitrogen was the major contributor to the observed amino acid isotopomers with one additional unit mass (M+1). However, there was differential incorporation, with the amine nitrogen of asparagine being incorporated into other amino acids more efficiently than the amide nitrogen. Higher isotopomers (M+2, M+3 and M+4) indicated the incorporation of varying numbers of (13)C atoms into essential amino acids. GC-MS assays also showed that, even with an excess of dietary labeled glutamine, glutamic acid, asparagine or aspartic acid, the overall content of these amino acids in aphid bodies was mostly the product of catabolism of dietary amino acids and subsequent re-synthesis within the aphids. Thus, these predominant dietary amino acids are not passed directly to Buchnera endosymbionts for synthesis of essential amino acids, but are rather are produced de novo, most likely by endogenous aphid enzymes.

Estimating Metabolic Fluxes Using a Maximum Network Flexibility Paradigm

Motivation

Genome-scale metabolic networks can be modeled in a constraint-based fashion. Reaction stoichiometry combined with flux capacity constraints determine the space of allowable reaction rates. This space is often large and a central challenge in metabolic modeling is finding the biologically most relevant flux distributions. A widely used method is flux balance analysis (FBA), which optimizes a biologically relevant objective such as growth or ATP production. Although FBA has proven to be highly useful for predicting growth and byproduct secretion, it cannot predict the intracellular fluxes under all environmental conditions. Therefore, alternative strategies have been developed to select flux distributions that are in agreement with experimental “omics” data, or by incorporating experimental flux measurements. The latter, unfortunately can only be applied to a limited set of reactions and is currently not feasible at the genome-scale. On the other hand, it has been observed that micro-organisms favor a suboptimal growth rate, possibly in exchange for a more “flexible” metabolic network. Instead of dedicating the internal network state to an optimal growth rate in one condition, a suboptimal growth rate is used, that allows for an easier switch to other nutrient sources. A small decrease in growth rate is exchanged for a relatively large gain in metabolic capability to adapt to changing environmental conditions.

Results

Here, we propose Maximum Metabolic Flexibility (MMF) a computational method that utilizes this observation to find the most probable intracellular flux distributions. By mapping measured flux data from central metabolism to the genome-scale models of Escherichia coli and Saccharomyces cerevisiae we show that i) indeed, most of the measured fluxes agree with a high adaptability of the network, ii) this result can be used to further reduce the space of feasible solutions iii) this reduced space improves the quantitative predictions made by FBA and contains a significantly larger fraction of the measured fluxes compared to the flux space that was reduced by a uniform sampling approach and iv) MMF can be used to select reactions in the network that contribute most to the steady-state flux space. Constraining the selected reactions improves the quantitative predictions of FBA considerably more than adding an equal amount of flux constraints, selected using a more naïve approach. Our method can be applied to any cell type without requiring prior information.

Availability

MMF is freely available as a MATLAB plugin at: http://cs.ru.nl/~wmegchel/mmf.

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.

While metabolic fluxes constitute the most direct window into a cell’s metabolism, their accurate measurement is non trivial. The gold standard for flux measurement involves providing a labeled feed where some of the carbon atoms have been substituted by isotopes with higher atomic mass (¹³C instead of ¹²C). The ensuing labeling found in intracellular metabolites is then used to computationally infer the metabolic fluxes that produced the observed pattern. However, this procedure is typically performed with small metabolic models encompassing only central carbon metabolism. The genomic revolution has afforded us easily available genomes and, with them, comprehensive genome-scale models of cellular metabolism. It would be desirable to use the ¹³C labeling experimental data to constrain genome-scale models: these data constrain fluxes very effectively and provide in the labeling data fit an obvious proof that the underlying model correctly explains measured quantities. Here, we introduce a rigorous, self-consistent method that uses the full amount of information contained in ¹³C labeling data to constrain fluxes for a genome-scale model where underlying assumptions are explicitly stated.

13C metabolic flux analysis at a genome-scale

Metabolic models used in 13C metabolic flux analysis generally include a limited number of reactions primarily from central metabolism. They typically omit degradation pathways, complete cofactor balances, and atom transition contributions for reactions outside central metabolism. This study addresses the impact on prediction fidelity of scaling-up mapping models to a genome-scale. The core mapping model employed in this study accounts for (75 reactions and 65 metabolites) primarily from central metabolism. The genome-scale metabolic mapping model (GSMM) (697 reaction and 595 metabolites) is constructed using as a basis the iAF1260 model upon eliminating reactions guaranteed not to carry flux based on growth and fermentation data for a minimal glucose growth medium. Labeling data for 17 amino acid fragments obtained from cells fed with glucose labeled at the second carbon was used to obtain fluxes and ranges. Metabolic fluxes and confidence intervals are estimated, for both core and genome-scale mapping models, by minimizing the sum of square of differences between predicted and experimentally measured labeling patterns using the EMU decomposition algorithm. Overall, we find that both topology and estimated values of the metabolic fluxes remain largely consistent between core and GSM model. Stepping up to a genome-scale mapping model leads to wider flux inference ranges for 20 key reactions present in the core model. The glycolysis flux range doubles due to the possibility of active gluconeogenesis, the TCA flux range expanded by 80% due to the availability of a bypass through arginine consistent with labeling data, and the transhydrogenase reaction flux was essentially unresolved due to the presence of as many as five routes for the inter-conversion of NADPH to NADH afforded by the genome-scale model. By globally accounting for ATP demands in the GSMM model the unused ATP decreased drastically with the lower bound matching the maintenance ATP requirement. A non-zero flux for the arginine degradation pathway was identified to meet biomass precursor demands as detailed in the iAF1260 model. Inferred ranges for 81% of the reactions in the genome-scale metabolic (GSM) model varied less than one-tenth of the basis glucose uptake rate (95% confidence test). This is because as many as 411 reactions in the GSM are growth coupled meaning that the single measurement of biomass formation rate locks the reaction flux values. This implies that accurate biomass formation rate and composition are critical for resolving metabolic fluxes away from central metabolism and suggests the importance of biomass composition (re)assessment under different genetic and environmental backgrounds. In addition, the loss of information associated with mapping fluxes from MFA on a core model to a GSM model is quantified.

Metabolic flux analysis of Escherichia coli MG1655 under octanoic acid (C8) stress

Systems metabolic engineering has made the renewable production of industrial chemicals a feasible alternative to modern operations. One major example of a renewable process is the production of carboxylic acids, such as octanoic acid (C8), from Escherichia coli, engineered to express thioesterase enzymes. C8, however, is toxic to E. coli above a certain concentration, which limits the final titer. (13)C metabolic flux analysis of E. coli was performed for both C8 stress and control conditions using NMR2Flux with isotopomer balancing. A mixture of labeled and unlabeled glucose was used as the sole carbon source for bacterial growth for (13)C flux analysis. By comparing the metabolic flux maps of the control condition and C8 stress condition, pathways that were altered under the stress condition were identified. C8 stress was found to reduce carbon flux in several pathways: the tricarboxylic acid (TCA) cycle, the CO2 production, and the pyruvate dehydrogenase pathway. Meanwhile, a few pathways became more active: the pyruvate oxidative pathway, and the extracellular acetate production. These results were statistically significant for three biological replicates between the control condition and C8 stress. As a working hypothesis, the following causes are proposed to be the main causes for growth inhibition and flux alteration for a cell under stress: membrane disruption, low activity of electron transport chain, and the activation of the pyruvate dehydrogenase regulator (PdhR).

OpenFLUX2: (13)C-MFA modeling software package adjusted for the comprehensive analysis of single and parallel labeling experiments

BACKGROUND

Steady-state (13)C-based metabolic flux analysis ((13)C-MFA) is the most powerful method available for the quantification of intracellular fluxes. These analyses include concertedly linked experimental and computational stages: (i) assuming the metabolic model and optimizing the experimental design; (ii) feeding the investigated organism using a chosen (13)C-labeled substrate (tracer); (iii) measuring the extracellular effluxes and detecting the (13)C-patterns of intracellular metabolites; and (iv) computing flux parameters that minimize the differences between observed and simulated measurements, followed by evaluating flux statistics. In its early stages, (13)C-MFA was performed on the basis of data obtained in a single labeling experiment (SLE) followed by exploiting the developed high-performance computational software. Recently, the advantages of parallel labeling experiments (PLEs), where several LEs are conducted under the conditions differing only by the tracer(s) choice, were demonstrated, particularly with regard to improving flux precision due to the synergy of complementary information. The availability of an open-source software adjusted for PLE-based (13)C-MFA is an important factor for PLE implementation.

RESULTS

The open-source software OpenFLUX, initially developed for the analysis of SLEs, was extended for the computation of PLE data. Using the OpenFLUX2, in silico simulation confirmed that flux precision is improved when (13)C-MFA is implemented by fitting PLE data to the common model compared with SLE-based analysis. Efficient flux resolution could be achieved in the PLE-mediated analysis when the choice of tracer was based on an experimental design computed to minimize the flux variances from different parts of the metabolic network. The analysis provided by OpenFLUX2 mainly includes (i) the optimization of the experimental design, (ii) the computation of the flux parameters from LEs data, (iii) goodness-of-fit testing of the model's adequacy, (iv) drawing conclusions concerning the identifiability of fluxes and construction of a contribution matrix reflecting the relative contribution of the measurement variances to the flux variances, and (v) precise determination of flux confidence intervals using a fine-tunable and convergence-controlled Monte Carlo-based method.

CONCLUSIONS

The developed open-source OpenFLUX2 provides a friendly software environment that facilitates beginners and existing OpenFLUX users to implement LEs for steady-state (13)C-MFA including experimental design, quantitative evaluation of flux parameters and statistics.

CeCaFDB: a curated database for the documentation, visualization and comparative analysis of central carbon metabolic flux distributions explored by 13C-fluxomics

The Central Carbon Metabolic Flux Database (CeCaFDB, available at http://www.cecafdb.org) is a manually curated, multipurpose and open-access database for the documentation, visualization and comparative analysis of the quantitative flux results of central carbon metabolism among microbes and animal cells. It encompasses records for more than 500 flux distributions among 36 organisms and includes information regarding the genotype, culture medium, growth conditions and other specific information gathered from hundreds of journal articles. In addition to its comprehensive literature-derived data, the CeCaFDB supports a common text search function among the data and interactive visualization of the curated flux distributions with compartmentation information based on the Cytoscape Web API, which facilitates data interpretation. The CeCaFDB offers four modules to calculate a similarity score or to perform an alignment between the flux distributions. One of the modules was built using an inter programming algorithm for flux distribution alignment that was specifically designed for this study. Based on these modules, the CeCaFDB also supports an extensive flux distribution comparison function among the curated data. The CeCaFDB is strenuously designed to address the broad demands of biochemists, metabolic engineers, systems biologists and members of the -omics community.

Central metabolic responses to the overproduction of fatty acids in Escherichia coli based on 13C-metabolic flux analysis

We engineered a fatty acid overproducing Escherichia coli strain through overexpressing tesA (“pull”) and fadR (“push”) and knocking out fadE (“block”). This “pull-push-block” strategy yielded 0.17 g of fatty acids (C12–C18) per gram of glucose (equivalent to 48% of the maximum theoretical yield) in batch cultures during the exponential growth phase under aerobic conditions. Metabolic fluxes were determined for the engineered E. coli and its control strain using tracer ([1,2-13C]glucose) experiments and 13C-metabolic flux analysis. Cofactor (NADPH) and energy (ATP) balances were also investigated for both strains based on estimated fluxes. Compared to the control strain, fatty acid overproduction led to significant metabolic responses in the central metabolism: (1) Acetic acid secretion flux decreased 10-fold; (2) Pentose phosphate pathway and Entner–Doudoroff pathway fluxes increased 1.5- and 2.0-fold, respectively; (3) Biomass synthesis flux was reduced 1.9-fold; (4) Anaplerotic phosphoenolpyruvate carboxylation flux decreased 1.7-fold; (5) Transhydrogenation flux converting NADH to NADPH increased by 1.7-fold. Real-time quantitative RT-PCR analysis revealed the engineered strain increased the transcription levels of pntA (encoding the membrane-bound transhydrogenase) by 2.1-fold and udhA (encoding the soluble transhydrogenase) by 1.4-fold, which is in agreement with the increased transhydrogenation flux. Cofactor and energy balances analyses showed that the fatty acid overproducing E. coli consumed significantly higher cellular maintenance energy than the control strain. We discussed the strategies to future strain development and process improvements for fatty acid production in E. coli.

A peptide-based method for 13C Metabolic Flux Analysis in microbial communities

The study of intracellular metabolic fluxes and inter-species metabolite exchange for microbial communities is of crucial importance to understand and predict their behaviour. The most authoritative method of measuring intracellular fluxes, ¹³C Metabolic Flux Analysis (¹³C MFA), uses the labeling pattern obtained from metabolites (typically amino acids) during ¹³C labeling experiments to derive intracellular fluxes. However, these metabolite labeling patterns cannot easily be obtained for each of the members of the community. Here we propose a new type of ¹³C MFA that infers fluxes based on peptide labeling, instead of amino acid labeling. The advantage of this method resides in the fact that the peptide sequence can be used to identify the microbial species it originates from and, simultaneously, the peptide labeling can be used to infer intracellular metabolic fluxes. Peptide identity and labeling patterns can be obtained in a high-throughput manner from modern proteomics techniques. We show that, using this method, it is theoretically possible to recover intracellular metabolic fluxes in the same way as through the standard amino acid based ¹³C MFA, and quantify the amount of information lost as a consequence of using peptides instead of amino acids. We show that by using a relatively small number of peptides we can counter this information loss. We computationally tested this method with a well-characterized simple microbial community consisting of two species.

Microbial communities underlie a variety of important biochemical processes ranging from underground cave formation to gold mining or the onset of obesity. Metabolic fluxes describe how carbon and energy flow through the microbial community and therefore provide insights that are rarely captured by other techniques, such as metatranscriptomics or metaproteomics. The most authoritative method to measure fluxes for pure cultures consists of feeding the cells a labeled carbon source and deriving the fluxes from the ensuing metabolite labeling pattern (typically amino acids). Since we cannot easily separate cells of metabolite for each species in a community, this approach is not generally applicable to microbial communities. Here we present a method to derive fluxes from the labeling of peptides, instead of amino acids. This approach has the advantage that peptides can be assigned to each species in a community in a high-throughput fashion through modern proteomic methods. We show that, by using this method, it is theoretically possible to recover the same amount of information as through the standard approach, if enough peptides are used. We computationally tested this method with a well-characterized simple microbial community consisting of two species.

The metabolic flux phenotype of heterotrophic Arabidopsis cells reveals a flexible balance between the cytosolic and plastidic contributions to carbohydrate oxidation in response to phosphate limitation

Understanding the mechanisms that allow plants to respond to variable and reduced availability of inorganic phosphate is of increasing agricultural importance because of the continuing depletion of the rock phosphate reserves that are used to combat inadequate phosphate levels in the soil. Changes in gene expression, protein levels, enzyme activities and metabolite levels all point to a reconfiguration of the central metabolic network in response to reduced availability of inorganic phosphate, but the metabolic significance of these changes can only be assessed in terms of the fluxes supported by the network. Steady-state metabolic flux analysis was used to define the metabolic phenotype of a heterotrophic Arabidopsis thaliana cell culture grown on a Murashige and Skoog medium containing 0, 1.25 or 5 mm inorganic phosphate. Fluxes through the central metabolic network were deduced from the redistribution of (13) C into metabolic intermediates and end products when cells were labelled with [1-(13) C], [2-(13) C], or [(13) C6 ]glucose, in combination with (14) C measurements of the rates of biomass accumulation. Analysis of the flux maps showed that reduced levels of phosphate in the growth medium stimulated flux through phosphoenolpyruvate carboxylase and malic enzyme, altered the balance between cytosolic and plastidic carbohydrate oxidation in favour of the plastid, and increased cell maintenance costs. We argue that plant cells respond to phosphate deprivation by reconfiguring the flux distribution through the pathways of carbohydrate oxidation to take advantage of better phosphate homeostasis in the plastid.

Systematic evaluation of methods for integration of transcriptomic data into constraint-based models of metabolism

Constraint-based models of metabolism are a widely used framework for predicting flux distributions in genome-scale biochemical networks. The number of published methods for integration of transcriptomic data into constraint-based models has been rapidly increasing. So far the predictive capability of these methods has not been critically evaluated and compared. This work presents a survey of recently published methods that use transcript levels to try to improve metabolic flux predictions either by generating flux distributions or by creating context-specific models. A subset of these methods is then systematically evaluated using published data from three different case studies in E. coli and S. cerevisiae. The flux predictions made by different methods using transcriptomic data are compared against experimentally determined extracellular and intracellular fluxes (from 13C-labeling data). The sensitivity of the results to method-specific parameters is also evaluated, as well as their robustness to noise in the data. The results show that none of the methods outperforms the others for all cases. Also, it is observed that for many conditions, the predictions obtained by simple flux balance analysis using growth maximization and parsimony criteria are as good or better than those obtained using methods that incorporate transcriptomic data. We further discuss the differences in the mathematical formulation of the methods, and their relation to the results we have obtained, as well as the connection to the underlying biological principles of metabolic regulation.

Elucidation of intrinsic biosynthesis yields using 13C-based metabolism analysis

This paper discusses the use of 13C-based metabolism analysis for the assessment of intrinsic product yields - the actual carbon contribution from a single carbon substrate to the final product via a specific biosynthesis route - in the following four cases. First, undefined nutrients (such as yeast extract) in fermentation may contribute significantly to product synthesis, which can be quantified through an isotopic dilution method. Second, product and biomass synthesis may be dependent on the co-metabolism of multiple-carbon sources. 13C labeling experiments can track the fate of each carbon substrate in the cell metabolism and identify which substrate plays a main role in product synthesis. Third, 13C labeling can validate and quantify the contribution of the engineered pathway (versus the native pathway) to the product synthesis. Fourth, the loss of catabolic energy due to cell maintenance (energy used for functions other than production of new cell components) and low P/O ratio (Phosphate/Oxygen Ratio) significantly reduces product yields. Therefore, 13C-metabolic flux analysis is needed to assess the influence of suboptimal energy metabolism on microbial productivity, and determine how ATP/NAD(P)H are partitioned among various cellular functions. Since product yield is a major determining factor in the commercialization of a microbial cell factory, we foresee that 13C-isotopic labeling experiments, even without performing extensive flux calculations, can play a valuable role in the development and verification of microbial cell factories.

optGpSampler: an improved tool for uniformly sampling the solution-space of genome-scale metabolic networks

Constraint-based models of metabolic networks are typically underdetermined, because they contain more reactions than metabolites. Therefore the solutions to this system do not consist of unique flux rates for each reaction, but rather a space of possible flux rates. By uniformly sampling this space, an estimated probability distribution for each reaction’s flux in the network can be obtained. However, sampling a high dimensional network is time-consuming. Furthermore, the constraints imposed on the network give rise to an irregularly shaped solution space. Therefore more tailored, efficient sampling methods are needed. We propose an efficient sampling algorithm (called optGpSampler), which implements the Artificial Centering Hit-and-Run algorithm in a different manner than the sampling algorithm implemented in the COBRA Toolbox for metabolic network analysis, here called gpSampler. Results of extensive experiments on different genome-scale metabolic networks show that optGpSampler is up to 40 times faster than gpSampler. Application of existing convergence diagnostics on small network reconstructions indicate that optGpSampler converges roughly ten times faster than gpSampler towards similar sampling distributions. For networks of higher dimension (i.e. containing more than 500 reactions), we observed significantly better convergence of optGpSampler and a large deviation between the samples generated by the two algorithms. Availability:optGpSampler for Matlab and Python is available for non-commercial use at: http://cs.ru.nl/~wmegchel/optGpSampler/.

Customization of ¹³C-MFA strategy according to cell culture system

(13)C-MFA is far from being a simple assay for quantifying metabolic activity. It requires considerable up-front experimental planning and familiarity with the cell culture system in question, as well as optimized analytics and adequate computation frameworks. The success of a (13)C-MFA experiment is ultimately rated by the ability to accurately quantify the flux of one or more reactions of interest. In this chapter, we describe the different (13)C-MFA strategies that have been developed for the various fermentation or cell culture systems, as well as the limitations of the respective strategies. The strategies are affected by many factors and the (13)C-MFA modeling and experimental strategy must be tailored to conditions. The prevailing philosophy in the computation process is that any metabolic processes that produce significant systematic bias in the labeling pattern of the metabolites being measured must be described in the model. It is equally important to plan a labeling strategy by analytical screening or by heuristics.

Refining carbon flux paths using atomic trace data

MOTIVATION

Pathway analysis tools are a powerful strategy to analyze 'omics' data in the field of systems biology. From a metabolic perspective, several pathway definitions can be found in the literature, each one appropriate for a particular study. Recently, a novel pathway concept termed carbon flux paths (CFPs) was introduced and benchmarked against existing approaches, showing a clear advantage for finding linear pathways from a given source to target metabolite. CFPs are simple paths in a metabolite-metabolite graph that satisfy typical constraints in stoichiometric models: mass balancing and thermodynamics (irreversibility). In addition, CFPs guarantee carbon exchange in each of their intermediate steps, but not between the source and the target metabolites and consequently false positive solutions may arise. These pathways often lack biological interest, particularly when studying biosynthetic or degradation routes of a metabolite. To overcome this issue, we amend the formulation in CFP, so as to account for atomic fate information. This approach is termed atomic CFP (aCFP).

RESULTS

By means of a side-by-side comparison in a medium scale metabolic network in Escherichia Coli, we show that aCFP provides more biologically relevant pathways than CFP, because canonical pathways are more easily recovered, which reflects the benefits of removing false positives. In addition, we demonstrate that aCFP can be successfully applied to genome-scale metabolic networks. As the quality of genome-scale atomic reconstruction is improved, methods such as the one presented here will undoubtedly be of value to interpret 'omics' data.

Computational evaluation of cellular metabolic costs successfully predicts genes whose expression is deleterious

Gene suppression and overexpression are both fundamental tools in linking genotype to phenotype in model organisms. Computational methods have proven invaluable in studying and predicting the deleterious effects of gene deletions, and yet parallel computational methods for overexpression are still lacking. Here, we present Expression-Dependent Gene Effects (EDGE), an in silico method that can predict the deleterious effects resulting from overexpression of either native or foreign metabolic genes. We first test and validate EDGE's predictive power in bacteria through a combination of small-scale growth experiments that we performed and analysis of extant large-scale datasets. Second, a broad cross-species analysis, ranging from microorganisms to multiple plant and human tissues, shows that genes that EDGE predicts to be deleterious when overexpressed are indeed typically down-regulated. This reflects a universal selection force keeping the expression of potentially deleterious genes in check. Third, EDGE-based analysis shows that cancer genetic reprogramming specifically suppresses genes whose overexpression impedes proliferation. The magnitude of this suppression is large enough to enable an almost perfect distinction between normal and cancerous tissues based solely on EDGE results. We expect EDGE to advance our understanding of human pathologies associated with up-regulation of particular transcripts and to facilitate the utilization of gene overexpression in metabolic engineering.

The ability of flux balance analysis to predict evolution of central metabolism scales with the initial distance to the optimum

The most powerful genome-scale framework to model metabolism, flux balance analysis (FBA), is an evolutionary optimality model. It hypothesizes selection upon a proposed optimality criterion in order to predict the set of internal fluxes that would maximize fitness. Here we present a direct test of the optimality assumption underlying FBA by comparing the central metabolic fluxes predicted by multiple criteria to changes measurable by a ¹³C-labeling method for experimentally-evolved strains. We considered datasets for three Escherichia coli evolution experiments that varied in their length, consistency of environment, and initial optimality. For ten populations that were evolved for 50,000 generations in glucose minimal medium, we observed modest changes in relative fluxes that led to small, but significant decreases in optimality and increased the distance to the predicted optimal flux distribution. In contrast, seven populations evolved on the poor substrate lactate for 900 generations collectively became more optimal and had flux distributions that moved toward predictions. For three pairs of central metabolic knockouts evolved on glucose for 600–800 generations, there was a balance between cases where optimality and flux patterns moved toward or away from FBA predictions. Despite this variation in predictability of changes in central metabolism, two generalities emerged. First, improved growth largely derived from evolved increases in the rate of substrate use. Second, FBA predictions bore out well for the two experiments initiated with ancestors with relatively sub-optimal yield, whereas those begun already quite optimal tended to move somewhat away from predictions. These findings suggest that the tradeoff between rate and yield is surprisingly modest. The observed positive correlation between rate and yield when adaptation initiated further from the optimum resulted in the ability of FBA to use stoichiometric constraints to predict the evolution of metabolism despite selection for rate.

The most common method of modeling genome-scale metabolism, flux balance analysis, involves using known stoichiometry to define feasible metabolic states and then choosing between these states by proposing that evolution has selected a metabolic flux that optimizes fitness. But does evolution optimize metabolism, and if so, what component of metabolism equates to fitness? We directly tested the underlying assumption of stoichiometric optimality by comparing predicted flux distributions with changes in fluxes that occurred following experimental evolution. Across three experiments ranging in length from a few hundred to fifty thousand generations, we found that substrate uptake – an input to the model – always increased, but supposed optimality criteria such as yield only increased sometimes. Despite this, there was a clear trend. Highly optimal ancestors evolved slightly lower yield in the course of increasing the overall rate, whereas more sub-optimal strains were able to increase both. These results suggest that flux balance analysis is capable of predicting either the initial metabolic behavior of strains or how they will evolve, but not both.

Error Analysis and Propagation in Metabolomics Data Analysis

Error analysis plays a fundamental role in describing the uncertainty in experimental results. It has several fundamental uses in metabolomics including experimental design, quality control of experiments, the selection of appropriate statistical methods, and the determination of uncertainty in results. Furthermore, the importance of error analysis has grown with the increasing number, complexity, and heterogeneity of measurements characteristic of 'omics research. The increase in data complexity is particularly problematic for metabolomics, which has more heterogeneity than other omics technologies due to the much wider range of molecular entities detected and measured. This review introduces the fundamental concepts of error analysis as they apply to a wide range of metabolomics experimental designs and it discusses current methodologies for determining the propagation of uncertainty in appropriate metabolomics data analysis. These methodologies include analytical derivation and approximation techniques, Monte Carlo error analysis, and error analysis in metabolic inverse problems. Current limitations of each methodology with respect to metabolomics data analysis are also discussed.

Systematic applications of metabolomics in metabolic engineering

The goals of metabolic engineering are well-served by the biological information provided by metabolomics: information on how the cell is currently using its biochemical resources is perhaps one of the best ways to inform strategies to engineer a cell to produce a target compound. Using the analysis of extracellular or intracellular levels of the target compound (or a few closely related molecules) to drive metabolic engineering is quite common. However, there is surprisingly little systematic use of metabolomics datasets, which simultaneously measure hundreds of metabolites rather than just a few, for that same purpose. Here, we review the most common systematic approaches to integrating metabolite data with metabolic engineering, with emphasis on existing efforts to use whole-metabolome datasets. We then review some of the most common approaches for computational modeling of cell-wide metabolism, including constraint-based models, and discuss current computational approaches that explicitly use metabolomics data. We conclude with discussion of the broader potential of computational approaches that systematically use metabolomics data to drive metabolic engineering.

Mathematical optimization applications in metabolic networks

Genome-scale metabolic models are increasingly becoming available for a variety of microorganisms. This has spurred the development of a wide array of computational tools, and in particular, mathematical optimization approaches, to assist in fundamental metabolic network analyses and redesign efforts. This review highlights a number of optimization-based frameworks developed towards addressing challenges in the analysis and engineering of metabolic networks. In particular, three major types of studies are covered here including exploring model predictions, correction and improvement of models of metabolism, and redesign of metabolic networks for the targeted overproduction of a desired compound. Overall, the methods reviewed in this paper highlight the diversity of queries, breadth of questions and complexity of redesigns that are amenable to mathematical optimization strategies.

Predicting outcomes of steady-state ¹³C isotope tracing experiments using Monte Carlo sampling

BACKGROUND

Carbon-13 (13C) analysis is a commonly used method for estimating reaction rates in biochemical networks. The choice of carbon labeling pattern is an important consideration when designing these experiments. We present a novel Monte Carlo algorithm for finding the optimal substrate input label for a particular experimental objective (flux or flux ratio). Unlike previous work, this method does not require assumption of the flux distribution beforehand.

RESULTS

Using a large E. coli isotopomer model, different commercially available substrate labeling patterns were tested computationally for their ability to determine reaction fluxes. The choice of optimal labeled substrate was found to be dependent upon the desired experimental objective. Many commercially available labels are predicted to be outperformed by complex labeling patterns. Based on Monte Carlo Sampling, the dimensionality of experimental data was found to be considerably less than anticipated, suggesting that effectiveness of 13C experiments for determining reaction fluxes across a large-scale metabolic network is less than previously believed.

CONCLUSIONS

While 13C analysis is a useful tool in systems biology, high redundancy in measurements limits the information that can be obtained from each experiment. It is however possible to compute potential limitations before an experiment is run and predict whether, and to what degree, the rate of each reaction can be resolved.

Toward biosynthetic design and implementation of Escherichia coli-derived paclitaxel and other heterologous polyisoprene compounds

Escherichia coli offers unparalleled engineering capacity in the context of heterologous natural product biosynthesis. However, as with other heterologous hosts, cellular metabolism must be designed or redesigned to support final compound formation. This task is at once complicated and aided by the fact that the cell does not natively produce an abundance of natural products. As a result, the metabolic engineer avoids complicated interactions with native pathways closely associated with the outcome of interest, but this convenience is tempered by the need to implement the required metabolism to allow functional biosynthesis. This review focuses on engineering E. coli for the purpose of polyisoprene formation, as it is related to isoprenoid compounds currently being pursued through a heterologous approach. In particular, the review features the compound paclitaxel and early efforts to design and overproduce intermediates through E. coli.

Metabolic pathway determination and flux analysis in nonmodel microorganisms through 13C-isotope labeling

C-isotope labeling is a commonly used technique for determining and quantifying pathways in microorganisms under various growth conditions. The experimental protocol consists of feeding the cell with a composition-defined substrate and measuring isotopic labeling patterns in the synthesized metabolites (often the amino acids). Not only can the labeling information be cross-referenced with genomic information to identify the novel pathways, but it can also be used to decipher absolute carbon fluxes through the metabolic network of interest. This technique can be widely used for functional characterization of nonmodel microbial species, and thus we provide a (13)C-pathway and flux analysis protocol. The five key procedures are: (1) growing cells using labeled substrates, (2) measuring extracellular metabolite and biomass component, (3) analyzing isotopic labeling patterns in amino acids and central metabolites using gas chromatography-mass spectrometry, (4) tracing (13)C carbon transitions in metabolites and discovering new pathways, and (5) estimating flux distributions based on isotopomer constraints. This protocol provides complementary information to the recently published protocol for (13)C-based metabolic flux analysis of the model species Escherichia coli (Nat Protoc 4:878-892, 2009).

Optimization of 13C isotopic tracers for metabolic flux analysis in mammalian cells

Mammalian cells consume and metabolize various substrates from their surroundings for energy generation and biomass synthesis. Glucose and glutamine, in particular, are the primary carbon sources for proliferating cancer cells. While this combination of substrates generates static labeling patterns for use in (13)C metabolic flux analysis (MFA), the inability of single tracers to effectively label all pathways poses an obstacle for comprehensive flux determination within a given experiment. To address this issue we applied a genetic algorithm to optimize mixtures of (13)C-labeled glucose and glutamine for use in MFA. We identified tracer combinations that minimized confidence intervals in an experimentally determined flux network describing central carbon metabolism in tumor cells. Additional simulations were used to determine the robustness of the [1,2-(13)C(2)]glucose/[U-(13)C(5)]glutamine tracer combination with respect to perturbations in the network. Finally, we experimentally validated the improved performance of this tracer set relative to glucose tracers alone in a cancer cell line. This versatile method allows researchers to determine the optimal tracer combination to use for a specific metabolic network, and our findings applied to cancer cells significantly enhance the ability of MFA experiments to precisely quantify fluxes in higher organisms.

Integration of in vivo and in silico metabolic fluxes for improvement of recombinant protein production

The filamentous fungus Aspergillus niger is an efficient host for the recombinant production of the glycosylated enzyme fructofuranosidase, a biocatalyst of commercial interest for the synthesis of pre-biotic sugars. In batch culture on a minimal glucose medium, the recombinant strain A. niger SKAn1015, expressing the fructofuranosidase encoding suc1 gene secreted 45U/mL of the target enzyme, whereas the parent wild type SKANip8 did not exhibit production. The production of the recombinant enzyme induced a significant change of in vivo fluxes in central carbon metabolism, as assessed by (13)C metabolic flux ratio analysis. Most notably, the flux redistribution enabled an elevated supply of NADPH via activation of the cytosolic pentose phosphate pathway (PPP) and mitochondrial malic enzyme, whereas the flux through energy generating TCA cycle was reduced. In addition, the overall possible flux space of fructofuranosidase producing A. niger was investigated in silico by elementary flux mode analysis. This provided theoretical flux distributions for multiple scenarios with differing production capacities. Subsequently, the measured flux changes linked to improved production performance were projected into the in silico flux space. This provided a quantitative evaluation of the achieved optimization and a priority ranked target list for further strain engineering. Interestingly, the metabolism was shifted largely towards the optimum flux pattern by sole expression of the recombinant enzyme, which seems an inherent attractive property of A. niger. Selected fluxes, however, changed contrary to the predicted optimum and thus revealed novel targets-including reactions linked to NADPH metabolism and gluconate formation.

Collisional fragmentation of central carbon metabolites in LC-MS/MS increases precision of ¹³C metabolic flux analysis

Experimental determination of fluxes by (13)C-tracers relies on detection of (13)C-patterns in metabolites or by-products. In the field of (13)C metabolic flux analysis, the most recent developments point toward recording labeling patterns by liquid chromatography (LC)-mass spectrometry (MS)/MS directly in intermediates in central carbon metabolism (CCM) to increase temporal resolution. Surprisingly, the flux studies published so far with LC-MS measurements were based on intact metabolic intermediates-thus neglected the potential benefits of using positional information to improve flux estimation. For the first time, we exploit collisional fragmentation to obtain more fine-grained (13)C-data on intermediates of CCM and investigate their impact in (13)C metabolic flux analysis. For the case study of Bacillus subtilis grown in mineral medium with (13)C-labeled glucose, we compare the flux estimates obtained by iterative isotopologue balancing of (13)C-data obtained either by LC-MS/MS for solely intact intermediates or LC-MS/MS for intact and fragmented intermediates of CCM. We show that with LC-MS/MS data, fragment information leads to more precise estimates of fluxes in pentose phosphate pathway, glycolysis, and to the tricarboxylic acid cycle. Additionally, we present an efficient analytical strategy to rapidly acquire large sets of (13)C-patterns by tandem MS, and an in-depth analysis of the collisional fragmentation of primary intermediates. In the future, this catalogue will enable comprehensive in silico calculability analyses to identify the most sensitive measurements and direct experimental design.

Experimental flux measurements on a network scale

Metabolic flux is a fundamental property of living organisms. In recent years, methods for measuring metabolic flux in plants on a network scale have evolved further. One major challenge in studying flux in plants is the complexity of the plant's metabolism. In particular, in the presence of parallel pathways in multiple cellular compartments, the core of plant central metabolism constitutes a complex network. Hence, a common problem with the reliability of the contemporary results of (13)C-Metabolic Flux Analysis in plants is the substantial reduction in complexity that must be included in the simulated networks; this omission partly is due to limitations in computational simulations. Here, I discuss recent emerging strategies that will better address these shortcomings.