Parallel labeling experiments and metabolic flux analysis: Past, present and future methodologies

The Carbon Flow Shifts from Primary to Secondary Metabolism during Xylem Vessel Cell Differentiation in Arabidopsis thaliana

Xylem vessel cell differentiation is characterized by the deposition of a secondary cell wall (SCW) containing cellulose, hemicellulose and lignin. VASCULAR-RELATED NAC-DOMAIN7 (VND7), a plant-specific NAC (NAM, ATAF1/2, and CUC2) transcription factor, is a master regulator of xylem vessel cell differentiation in Arabidopsis (Arabidopsis thaliana). Previous metabolome analysis using the VND7-inducible system in tobacco BY-2 cells successfully revealed significant quantitative changes in primary metabolites during xylem vessel cell differentiation. However, the flow of primary metabolites is not yet well understood. Here, we performed a metabolomic analysis of VND7-inducible Arabidopsis T87 suspension cells. Capillary electrophoresis-time-of-flight mass spectrometry quantified 57 metabolites, and subsequent data analysis highlighted active changes in the levels of UDP-glucose and phenylalanine, which are building blocks of cellulose and lignin, respectively. In a metabolic flow analysis using stable carbon 13 (13C) isotope, the 13C-labeling ratio specifically increased in 3-phosphoglycerate after 12 h of VND7 induction, followed by an increase in shikimate after 24 h of induction, while the inflow of 13C into lactate from pyruvate was significantly inhibited, indicating an active shift of carbon flow from glycolysis to the shikimate pathway during xylem vessel cell differentiation. In support of this notion, most glycolytic genes involved in the downstream of glyceraldehyde 3-phosphate were downregulated following the induction of xylem vessel cell differentiation, whereas genes for the shikimate pathway and phenylalanine biosynthesis were upregulated. These findings provide evidence for the active shift of carbon flow from primary metabolic pathways to the SCW polymer biosynthetic pathway at specific points during xylem vessel cell differentiation.

Utilization of lignocellulosic hydrolysates for photomixotrophic chemical production in Synechococcus elongatus PCC 7942

To meet the need for environmentally friendly commodity chemicals, feedstocks for biological chemical production must be diversified. Lignocellulosic biomass are an carbon source with the potential for effective use in a large scale and cost-effective production systems. Although the use of lignocellulosic biomass lysates for heterotrophic chemical production has been advancing, there are challenges to overcome. Here we aim to investigate the obligate photoautotroph cyanobacterium Synechococcus elongatus PCC 7942 as a chassis organism for lignocellulosic chemical production. When modified to import monosaccharides, this cyanobacterium is an excellent candidate for lysates-based chemical production as it grows well at high lysate concentrations and can fix CO2 to enhance carbon efficiency. This study is an important step forward in enabling the simultaneous use of two sugars as well as lignocellulosic lysate. Incremental genetic modifications enable catabolism of both sugars concurrently without experiencing carbon catabolite repression. Production of 2,3-butanediol is demonstrated to characterize chemical production from the sugars in lignocellulosic hydrolysates. The engineered strain achieves a titer of 13.5 g L-1 of 2,3-butanediol over 12 days under shake-flask conditions. This study can be used as a foundation for industrial scale production of commodity chemicals from a combination of sunlight, CO2, and lignocellulosic sugars.

A pathogen-specific isotope tracing approach reveals metabolic activities and fluxes of intracellular Salmonella

Pathogenic bacteria proliferating inside mammalian host cells need to rapidly adapt to the intracellular environment. How they achieve this and scavenge essential nutrients from the host has been an open question due to the difficulties in distinguishing between bacterial and host metabolites in situ. Here, we capitalized on the inability of mammalian cells to metabolize mannitol to develop a stable isotopic labeling approach to track Salmonella enterica metabolites during intracellular proliferation in host macrophage and epithelial cells. By measuring label incorporation into Salmonella metabolites with liquid chromatography-mass spectrometry (LC-MS), and combining it with metabolic modeling, we identify relevant carbon sources used by Salmonella, uncover routes of their metabolization, and quantify relative reaction rates in central carbon metabolism. Our results underline the importance of the Entner-Doudoroff pathway (EDP) and the phosphoenolpyruvate carboxylase for intracellularly proliferating Salmonella. More broadly, our metabolic labeling strategy opens novel avenues for understanding the metabolism of pathogens inside host cells.

Stochastic simulation algorithm for isotope-based dynamic flux analysis

Carbon isotope labeling method is a standard metabolic engineering tool for flux quantification in living cells. To cope with the high dimensionality of isotope labeling systems, diverse algorithms have been developed to reduce the number of variables or operations in metabolic flux analysis (MFA), but lacks generalizability to non-stationary metabolic conditions. In this study, we present a stochastic simulation algorithm (SSA) derived from the chemical master equation of the isotope labeling system. This algorithm allows to compute the time evolution of isotopomer concentrations in non-stationary conditions, with the valuable property that computational time does not scale with the number of isotopomers. The efficiency and limitations of the algorithm is benchmarked for the forward and inverse problems of 13C-DMFA in the pentose phosphate pathways, and is compared with EMU-based methods for NMFA and MFA including the central carbon metabolism. Overall, SSA constitutes an alternative class to deterministic approaches for metabolic flux analysis that is well adapted to comprehensive dataset including parallel labeling experiments, and whose limitations associated to the sampling size can be overcome by using Monte Carlo sampling approaches.

Dietary intake of deuterium oxide decreases cochlear metabolism and oxidative stress levels in a mouse model of age-related hearing loss

Age-related hearing loss (ARHL) is the most prevalent sensory disorder in the elderly. Currently, no treatment can effectively prevent or reverse ARHL. Aging auditory organs are often accompanied by exacerbated oxidative stress and metabolic deterioration. Here, we report the effect of deuterated oxygen (D2O), also known as "heavy water", mouse models of ARHL. Supplementing the normal mouse diet with 10% D2O from 4 to 9 weeks of age lowered hearing thresholds at selected frequencies in treated mice compared to untreated control group. Oxidative stress levels were significantly reduced and in the cochlear duct of treated vs. untreated mice. Through metabolic flux analysis, we found that D2O mainly slowed down catabolic reactions, and may delay metabolic deterioration related to aging to a certain extent. Experiments confirmed that the Nrf2/HO-1/glutathione axis was down-regulated in treated mice. Thus, D2O supplementation can hinder ARHL progression in mouse models by slowing the pace of metabolism and reducing endogenous oxidative stress production in the cochlea. These findings open new avenues for protecting the cochlea from oxidative stress and regulating metabolism to prevent ARHL.

Statin therapy inhibits fatty acid synthase via dynamic protein modifications

Statins are a class of drug widely prescribed for the prevention of cardiovascular disease, with pleiotropic cellular effects. Statins inhibit HMG-CoA reductase (HMGCR), which converts the metabolite HMG-CoA into mevalonate. Recent discoveries have shown HMG-CoA is a reactive metabolite that can non-enzymatically modify proteins and impact their activity. Therefore, we predicted that inhibition of HMGCR by statins might increase HMG-CoA levels and protein modifications. Upon statin treatment, we observe a strong increase in HMG-CoA levels and modification of only a single protein. Mass spectrometry identifies this protein as fatty acid synthase (FAS), which is modified on active site residues and, importantly, on non-lysine side-chains. The dynamic modifications occur only on a sub-pool of FAS that is located near HMGCR and alters cellular signaling around the ER and Golgi. These results uncover communication between cholesterol and lipid biosynthesis by the substrate of one pathway inhibiting another in a rapid and reversible manner.

Isotopically nonstationary 13C metabolic flux analysis in resting and activated human platelets

Platelet metabolism is linked to platelet hyper- and hypoactivity in numerous human diseases. Developing a detailed understanding of the link between metabolic shifts and platelet activation state is integral to improving human health. Here, we show the first application of isotopically nonstationary ¹³C metabolic flux analysis to quantitatively measure carbon fluxes in both resting and thrombin activated platelets. Metabolic flux analysis results show that resting platelets primarily metabolize glucose to lactate via glycolysis, while acetate is oxidized to fuel the tricarboxylic acid cycle. Upon activation with thrombin, a potent platelet agonist, platelets increase their uptake of glucose 3-fold. This results in an absolute increase in flux throughout central metabolism, but when compared to resting platelets they redistribute carbon dramatically. Activated platelets decrease relative flux to the oxidative pentose phosphate pathway and TCA cycle from glucose and increase relative flux to lactate. These results provide the first report of reaction-level carbon fluxes in platelets and allow us to distinguish metabolic fluxes with much higher resolution than previous studies.

How to Tackle Underdeterminacy in Metabolic Flux Analysis? A Tutorial and Critical Review

Metabolic flux analysis is often (not to say almost always) faced with system underdeterminacy. Indeed, the linear algebraic system formed by the steady-state mass balance equations around the intracellular metabolites and the equality constraints related to the measurements of extracellular fluxes do not define a unique solution for the distribution of intracellular fluxes, but instead a set of solutions belonging to a convex polytope. Various methods have been proposed to tackle this underdeterminacy, including flux pathway analysis, flux balance analysis, flux variability analysis and sampling. These approaches are reviewed in this article and a toy example supports the discussion with illustrative numerical results.

Utilizing tandem mass spectrometry for metabolic flux analysis

Metabolic flux analysis (MFA) aims at revealing the metabolic reaction rates in a complex biochemical network. To do so, MFA uses the input of stable isotope labeling patterns of the intracellular metabolites. Elementary metabolic unit (EMU) is the computational framework to simulate the metabolite labeling patterns in a network, which was originally designed for simulating mass isotopomer distributions (MIDs) at the MS1 level. Recently, the EMU framework is expanded to simulate tandem mass spectrometry data. Tandem mass spectrometry has emerged as a new experimental approach to provide information on the positional isotope labeling of metabolites and therefore greatly improves the precision of MFA. In this review, we will discuss the new EMU framework that can accommodate the tandem mass isotopomer distributions (TMIDs) data. We will also analyze the improvement on the MFA precision by using TMID. Our analysis shows that combining the MIDs of the parent and daughter ions and the TMID for the MFA is more powerful than using TMID alone.

Considerations for using isolated cell systems to understand cardiac metabolism and biology

Changes in myocardial metabolic activity are fundamentally linked to cardiac health and remodeling. Primary cardiomyocytes, induced pluripotent stem cell-derived cardiomyocytes, and transformed cardiomyocyte cell lines are common models used to understand how (patho)physiological conditions or stimuli contribute to changes in cardiac metabolism. These cell models are helpful also for defining metabolic mechanisms of cardiac dysfunction and remodeling. Although technical advances have improved our capacity to measure cardiomyocyte metabolism, there is often heterogeneity in metabolic assay protocols and cell models, which could hinder data interpretation and discernment of the mechanisms of cardiac (patho)physiology. In this review, we discuss considerations for integrating cardiomyocyte cell models with techniques that have become relatively common in the field, such as respirometry and extracellular flux analysis. Furthermore, we provide overviews of metabolic assays that complement XF analyses and that provide information on not only catabolic pathway activity, but biosynthetic pathway activity and redox status as well. Cultivating a more widespread understanding of the advantages and limitations of metabolic measurements in cardiomyocyte cell models will continue to be essential for the development of coherent metabolic mechanisms of cardiac health and pathophysiology.

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

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

Metabolic Modelling as a Framework for Metabolomics Data Integration and Analysis

Metabolic networks are regulated to ensure the dynamic adaptation of biochemical reaction fluxes to maintain cell homeostasis and optimal metabolic fitness in response to endogenous and exogenous perturbations. To this end, metabolism is tightly controlled by dynamic and intricate regulatory mechanisms involving allostery, enzyme abundance and post-translational modifications. The study of the molecular entities involved in these complex mechanisms has been boosted by the advent of high-throughput technologies. The so-called omics enable the quantification of the different molecular entities at different system layers, connecting the genotype with the phenotype. Therefore, the study of the overall behavior of a metabolic network and the omics data integration and analysis must be approached from a holistic perspective. Due to the close relationship between metabolism and cellular phenotype, metabolic modelling has emerged as a valuable tool to decipher the underlying mechanisms governing cell phenotype. Constraint-based modelling and kinetic modelling are among the most widely used methods to study cell metabolism at different scales, ranging from cells to tissues and organisms. These approaches enable integrating metabolomic data, among others, to enhance model predictive capabilities. In this review, we describe the current state of the art in metabolic modelling and discuss future perspectives and current challenges in the field.

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.

Simultaneous tracers and a unified model of positional and mass isotopomers for quantification of metabolic flux in liver

Computational models based on the metabolism of stable isotope tracers can yield valuable insight into the metabolic basis of disease. The complexity of these models is limited by the number of tracers and the ability to characterize tracer labeling in downstream metabolites. NMR spectroscopy is ideal for multiple tracer experiments since it precisely detects the position of tracer nuclei in molecules, but it lacks sensitivity for detecting low-concentration metabolites. GC-MS detects stable isotope mass enrichment in low-concentration metabolites, but lacks nuclei and positional specificity. We performed liver perfusions and in vivo infusions of ²H and ¹³C tracers, yielding complex glucose isotopomers that were assigned by NMR and fit to a newly developed metabolic model. Fluxes regressed from ²H and ¹³C NMR positional isotopomer enrichments served to validate GC-MS-based flux estimates obtained from the same experimental samples. NMR-derived fluxes were largely recapitulated by modeling the mass isotopomer distributions of six glucose fragment ions measured by GC-MS. Modest differences related to limited fragmentation coverage of glucose C1-C3 were identified, but fluxes such as gluconeogenesis, glycogenolysis, cataplerosis and TCA cycle flux were tightly correlated between the methods. Most importantly, modeling of GC-MS data could assign fluxes in primary mouse hepatocytes, an experiment that is impractical by ²H or ¹³C NMR.

Bioenergetics and translational metabolism: implications for genetics, physiology and precision medicine

It is now becoming clear that human metabolism is extremely plastic and varies substantially between healthy individuals. Understanding the biochemistry that underlies this physiology will enable personalized clinical interventions related to metabolism. Mitochondrial quality control and the detailed mechanisms of mitochondrial energy generation are central to understanding susceptibility to pathologies associated with aging including cancer, cardiac and neurodegenerative diseases. A precision medicine approach is also needed to evaluate the impact of exercise or caloric restriction on health. In this review, we discuss how technical advances in assessing mitochondrial genetics, cellular bioenergetics and metabolomics offer new insights into developing metabolism-based clinical tests and metabolotherapies. We discuss informatics approaches, which can define the bioenergetic-metabolite interactome and how this can help define healthy energetics. We propose that a personalized medicine approach that integrates metabolism and bioenergetics with physiologic parameters is central for understanding the pathophysiology of diseases with a metabolic etiology. New approaches that measure energetics and metabolomics from cells isolated from human blood or tissues can be of diagnostic and prognostic value to precision medicine. This is particularly significant with the development of new metabolotherapies, such as mitochondrial transplantation, which could help treat complex metabolic diseases.

Metabolic flux ratio analysis by parallel 13C labeling of isoprenoid biosynthesis in Rhodobacter sphaeroides

Metabolic engineering for increased isoprenoid production often benefits from the simultaneous expression of the two naturally available isoprenoid metabolic routes, namely the 2-methyl-D-erythritol 4-phosphate (MEP) pathway and the mevalonate (MVA) pathway. Quantification of the contribution of these pathways to the overall isoprenoid production can help to obtain a better understanding of the metabolism within a microbial cell factory. Such type of investigation can benefit from ¹³C metabolic flux ratio studies. Here, we designed a method based on parallel labeling experiments (PLEs), using [1-¹³C]- and [4-¹³C]glucose as tracers to quantify the metabolic flux ratios in the glycolytic and isoprenoid pathways. By just analyzing a reporter isoprenoid molecule and employing only four equations, we could describe the metabolism involved from substrate catabolism to product formation. These equations infer ¹³C atom incorporation into the universal isoprenoid building blocks, isopentenyl-pyrophosphate (IPP) and dimethylallyl-pyrophosphate (DMAPP). Therefore, this renders the method applicable to the study of any of isoprenoid of interest. As proof of principle, we applied it to study amorpha-4,11-diene biosynthesis in the bacterium Rhodobacter sphaeroides. We confirmed that in this species the Entner-Doudoroff pathway is the major pathway for glucose catabolism, while the Embden-Meyerhof-Parnas pathway contributes to a lesser extent. Additionally, we demonstrated that co-expression of the MEP and MVA pathways caused a mutual enhancement of their metabolic flux capacity. Surprisingly, we also observed that the isoprenoid flux ratio remains constant under exponential growth conditions, independently from the expression level of the MVA pathway. Apart from proposing and applying a tool for studying isoprenoid biosynthesis within a microbial cell factory, our work reveals important insights from the co-expression of MEP and MVA pathways, including the existence of a yet unclear interaction between them.

Integration of flux measurements and pharmacological controls to optimize stable isotope-resolved metabolomics workflows and interpretation

Stable isotope-resolved metabolomics (SIRM) provides information regarding the relative activity of numerous metabolic pathways and the contribution of nutrients to specific metabolite pools; however, SIRM experiments can be difficult to execute, and data interpretation is challenging. Furthermore, standardization of analytical procedures and workflows remain significant obstacles for widespread reproducibility. Here, we demonstrate the workflow of a typical SIRM experiment and suggest experimental controls and measures of cross-validation that improve data interpretation. Inhibitors of glycolysis and oxidative phosphorylation as well as mitochondrial uncouplers serve as pharmacological controls, which help define metabolic flux configurations that occur under well-controlled metabolic states. We demonstrate how such controls and time course labeling experiments improve confidence in metabolite assignments as well as delineate metabolic pathway relationships. Moreover, we demonstrate how radiolabeled tracers and extracellular flux analyses integrate with SIRM to improve data interpretation. Collectively, these results show how integration of flux methodologies and use of pharmacological controls increase confidence in SIRM data and provide new biological insights.

A Pareto approach to resolve the conflict between information gain and experimental costs: Multiple-criteria design of carbon labeling experiments

Science revolves around the best way of conducting an experiment to obtain insightful results. Experiments with maximal information content can be found by computational experimental design (ED) strategies that identify optimal conditions under which to perform the experiment. Several criteria have been proposed to measure the information content, each emphasizing different aspects of the design goal, i.e., reduction of uncertainty. Where experiments are complex or expensive, second sight is at the budget governing the achievable amount of information. In this context, the design objectives cost and information gain are often incommensurable, though dependent. By casting the ED task into a multiple-criteria optimization problem, a set of trade-off designs is derived that approximates the Pareto-frontier which is instrumental for exploring preferable designs. In this work, we present a computational methodology for multiple-criteria ED of information-rich experiments that accounts for virtually any set of design criteria. The methodology is implemented for the case of ¹³C metabolic flux analysis (MFA), which is arguably the most expensive type among the ‘omics’ technologies, featuring dozens of design parameters (tracer composition, analytical platform, measurement selection etc.). Supported by an innovative visualization scheme, we demonstrate with two realistic showcases that the use of multiple criteria reveals deep insights into the conflicting interplay between information carriers and cost factors that are not amendable to single-objective ED. For instance, tandem mass spectrometry turns out as best-in-class with respect to information gain, while it delivers this information quality cheaper than the other, routinely applied analytical technologies. Therewith, our Pareto approach to ED offers the investigator great flexibilities in the conception phase of a study to balance costs and benefits.

Designing experiments is obligatory in the biosciences to valorize their scientific outcome. When the experiments are expensive, unfortunately, in practice often the costs emerge to be showstoppers. In this situation the question arises: How to get the most out of the experiment for your invest in terms of time and money? We approach this question by formulating the design task as a multiple-criteria optimization problem. Its solution produces a set of Pareto-optimal design proposals that feature the trade-off between information gain, as measured by different metrics, and the costs. Then, exploration of the design proposals allows us to make the best decision on information-economic experiments under given circumstances. Implemented in the field of isotope-based metabolic flux analysis, practical application of the Pareto approach provides detailed insight into the tight interplay of plenty of information carriers and cost factors. Supported by an innovative tailored visual representation scheme, the investigator is enabled to explore the options before conducting the experiment. With a practical showcase at hand, our computational study highlights the benefits of incorporating multiple information criteria apart from the costs, balancing the shortcomings of conventional single-objective experimental design strategies.

SWATH Tandem Mass Spectrometry Workflow for Quantification of Mass Isotopologue Distribution of Intracellular Metabolites and Fragments Labeled with Isotopic 13C Carbon

Accurate quantification of mass isotopologue distribution (MID) of metabolites is a prerequisite for ¹³C-metabolic flux analysis. Currently used mass spectrometric (MS) techniques based on multiple reaction monitoring (MRM) place limitations on the number of MIDs that can be analyzed in a single run. Moreover, the deconvolution step results in amplification of error. Here, we demonstrate that SWATH MS/MS, a data independent acquisition (DIA) technique allows quantification of a large number of precursor and product MIDs in a single run. SWATH sequentially fragments all precursor ions in stacked mass isolation windows. Co-fragmentation of all precursor isotopologues in a single SWATH window yields higher sensitivity enabling quantification of MIDs of fragments with low abundance and lower systematic and random errors. We quantify the MIDs of 53 precursor and product ions corresponding to 19 intracellular metabolites from a dynamic ¹³C-labeling of a model cyanobacterium, Synechococcus sp. PCC 7002. The use of product MIDs resulted in an improved precision of many measured fluxes compared to when only precursor MIDs were used for flux analysis. The approach is truly untargeted and allows additional metabolites to be quantified from the same data.

A guide to 13C metabolic flux analysis for the cancer biologist

Cancer metabolism is significantly altered from normal cellular metabolism allowing cancer cells to adapt to changing microenvironments and maintain high rates of proliferation. In the past decade, stable-isotope tracing and network analysis have become powerful tools for uncovering metabolic pathways that are differentially activated in cancer cells. In particular, ¹³C metabolic flux analysis (¹³C-MFA) has emerged as the primary technique for quantifying intracellular fluxes in cancer cells. In this review, we provide a practical guide for investigators interested in getting started with ¹³C-MFA. We describe best practices in ¹³C-MFA, highlight potential pitfalls and alternative approaches, and conclude with new developments that can further enhance our understanding of cancer metabolism.

Tracing tagged molecules can help researchers understand the altered metabolism of cancer cells. The abilities of cancer cells to multiply rapidly and invade new tissues are supported by metabolic alterations, which can be investigated by feeding tagged molecules to cells and tracing how they are metabolized. These techniques, such as ¹³C metabolic flux analysis (¹³C-MFA), have been perceived as difficult to use, but recent advances are making them more accessible. Maciek Antoniewicz, University of Delaware, Newark, USA, has published a practical guide for researchers wanting to use 13C-MFA. The review includes best practices, pitfalls, alternative approaches, and new developments, especially new user-friendly software that allows researchers without extensive training in mathematics, statistics, or coding to perform 13C-MFA. Broadening access to tools for investigating altered metabolic pathways may spur development of new cancer therapies targeting these pathways.

From correlation to causation: analysis of metabolomics data using systems biology approaches

INTRODUCTION

Metabolomics is a well-established tool in systems biology, especially in the top-down approach. Metabolomics experiments often results in discovery studies that provide intriguing biological hypotheses but rarely offer mechanistic explanation of such findings. In this light, the interpretation of metabolomics data can be boosted by deploying systems biology approaches.

OBJECTIVES

This review aims to provide an overview of systems biology approaches that are relevant to metabolomics and to discuss some successful applications of these methods.

METHODS

We review the most recent applications of systems biology tools in the field of metabolomics, such as network inference and analysis, metabolic modelling and pathways analysis.

RESULTS

We offer an ample overview of systems biology tools that can be applied to address metabolomics problems. The characteristics and application results of these tools are discussed also in a comparative manner.

CONCLUSIONS

Systems biology-enhanced analysis of metabolomics data can provide insights into the molecular mechanisms originating the observed metabolic profiles and enhance the scientific impact of metabolomics studies.

Applications of Metabolomics in Cancer Studies

Since the start of metabolomics as a field of research, the number of studies related to cancer has grown to such an extent that cancer metabolomics now represents its own discipline. In this chapter, the applications of metabolomics in cancer studies are explored. Different approaches and analytical platforms can be employed for the analysis of samples depending on the goal of the study and the aspects of the cancer metabolome being investigated. Analyses have concerned a range of cancers including lung, colorectal, bladder, breast, gastric, oesophageal and thyroid, amongst others. Developments in these strategies and methodologies that have been applied are discussed, in addition to exemplifying the use of cancer metabolomics in the discovery of biomarkers and in the assessment of therapy (both pharmaceutical and nutraceutical). Finally, the application of cancer metabolomics in personalised medicine is presented.

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.

Synergizing 13C Metabolic Flux Analysis and Metabolic Engineering for Biochemical Production

Metabolic engineering of industrial microorganisms to produce chemicals, fuels, and drugs has attracted increasing interest as it provides an environment-friendly and renewable route that does not depend on depleting petroleum sources. However, the microbial metabolism is so complex that metabolic engineering efforts often have difficulty in achieving a satisfactory yield, titer, or productivity of the target chemical. To overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been developed to investigate rigorously the cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, ¹³C-MFA has been widely used in academic labs and the biotechnology industry to pinpoint the key issues related to microbial-based chemical production and to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this chapter we introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied to synergize with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production.

Comprehensive metabolic modeling of multiple 13C-isotopomer data sets to study metabolism in perfused working hearts

In many forms of cardiomyopathy, alterations in energy substrate metabolism play a key role in disease pathogenesis. Stable isotope tracing in rodent heart perfusion systems can be used to determine cardiac metabolic fluxes, namely those relative fluxes that contribute to pyruvate, the acetyl-CoA pool, and pyruvate anaplerosis, which are critical to cardiac homeostasis. Methods have previously been developed to interrogate these relative fluxes using isotopomer enrichments of measured metabolites and algebraic equations to determine a predefined metabolic flux model. However, this approach is exquisitely sensitive to measurement error, thus precluding accurate relative flux parameter determination. In this study, we applied a novel mathematical approach to determine relative cardiac metabolic fluxes using ¹³C-metabolic flux analysis (¹³C-MFA) aided by multiple tracer experiments and integrated data analysis. Using ¹³C-MFA, we validated a metabolic network model to explain myocardial energy substrate metabolism. Four different ¹³C-labeled substrates were queried (i.e., glucose, lactate, pyruvate, and oleate) based on a previously published study. We integrated the analysis of the complete set of isotopomer data gathered from these mouse heart perfusion experiments into a single comprehensive network model that delineates substrate contributions to both pyruvate and acetyl-CoA pools at a greater resolution than that offered by traditional methods using algebraic equations. To our knowledge, this is the first rigorous application of ¹³C-MFA to interrogate data from multiple tracer experiments in the perfused heart. We anticipate that this approach can be used widely to study energy substrate metabolism in this and other similar biological systems.

Optimal tracers for parallel labeling experiments and 13C metabolic flux analysis: A new precision and synergy scoring system

¹³C-Metabolic flux analysis (¹³C-MFA) is a widely used approach in metabolic engineering for quantifying intracellular metabolic fluxes. The precision of fluxes determined by ¹³C-MFA depends largely on the choice of isotopic tracers and the specific set of labeling measurements. A recent advance in the field is the use of parallel labeling experiments for improved flux precision and accuracy. However, as of today, no systemic methods exist for identifying optimal tracers for parallel labeling experiments. In this contribution, we have addressed this problem by introducing a new scoring system and evaluating thousands of different isotopic tracer schemes. Based on this extensive analysis we have identified optimal tracers for ¹³C-MFA. The best single tracers were doubly ¹³C-labeled glucose tracers, including [1,6-¹³C]glucose, [5,6-¹³C]glucose and [1,2-¹³C]glucose, which consistently produced the highest flux precision independent of the metabolic flux map (here, 100 random flux maps were evaluated). Moreover, we demonstrate that pure glucose tracers perform better overall than mixtures of glucose tracers. For parallel labeling experiments the optimal isotopic tracers were [1,6-¹³C]glucose and [1,2-¹³C]glucose. Combined analysis of [1,6-¹³C]glucose and [1,2-¹³C]glucose labeling data improved the flux precision score by nearly 20-fold compared to widely use tracer mixture 80% [1-¹³C]glucose +20% [U-¹³C]glucose.

Characterization of physiological responses to 22 gene knockouts in Escherichia coli central carbon metabolism

Understanding the impact of gene knockouts on cellular physiology, and metabolism in particular, is centrally important to quantitative systems biology and metabolic engineering. Here, we present a comprehensive physiological characterization of wild-type Escherichia coli and 22 knockouts of enzymes in the upper part of central carbon metabolism, including the PTS system, glycolysis, pentose phosphate pathway and Entner-Doudoroff pathway. Our results reveal significant metabolic changes that are affected by specific gene knockouts. Analysis of collective trends and correlations in the data using principal component analysis (PCA) provide new, and sometimes surprising, insights into E. coli physiology. Additionally, by comparing the data-to-model predictions from constraint-based approaches such as FBA, MOMA and RELATCH we demonstrate the important role of less well-understood kinetic and regulatory effects in central carbon metabolism.

Evidence for transketolase-like TKTL1 flux in CHO cells based on parallel labeling experiments and (13)C-metabolic flux analysis

The pentose phosphate pathway (PPP) is a fundamental component of cellular metabolism. It provides precursors for the biosynthesis of nucleotides and contributes to the production of reducing power in the form of NADPH. It has been hypothesized that mammalian cells may contain a hidden reaction in PPP catalyzed by transketolase-like protein 1 (TKTL1) that is closely related to the classical transketolase enzyme; however, until now there has been no direct experimental evidence for this reaction. In this work, we have applied state-of-the-art techniques in (13)C metabolic flux analysis ((13)C-MFA) based on parallel labeling experiments and integrated flux fitting to estimate the TKTL1 flux in CHO cells. We identified a set of three parallel labeling experiments with [1-(13)C]glucose+[4,5,6-(13)C]glucose, [2-(13)C]glucose+[4,5,6-(13)C]glucose, and [3-(13)C]glucose+[4,5,6-(13)C]glucose and developed a new method to measure (13)C-labeling of fructose 6-phosphate by GC-MS that allows intuitive interpretation of mass isotopomer distributions to determine key fluxes in the model, including glycolysis, oxidative PPP, non-oxidative PPP, and the TKTL1 flux. Using these tracers we detected a significant TKTL1 flux in CHO cells at the stationary phase. The flux results suggest that the main function of oxidative PPP in CHO cells at the stationary phase is to fuel the TKTL1 reaction. Overall, this study demonstrates for the first time that carbon atoms can be lost in the PPP, by means other than the oxidative PPP, and that this loss of carbon atoms is consistent with the hypothesized TKTL1 reaction in mammalian cells.

Co-utilization of glucose and xylose by evolved Thermus thermophilus LC113 strain elucidated by (13)C metabolic flux analysis and whole genome sequencing

We evolved Thermus thermophilus to efficiently co-utilize glucose and xylose, the two most abundant sugars in lignocellulosic biomass, at high temperatures without carbon catabolite repression. To generate the strain, T. thermophilus HB8 was first evolved on glucose to improve its growth characteristics, followed by evolution on xylose. The resulting strain, T. thermophilus LC113, was characterized in growth studies, by whole genome sequencing, and (13)C-metabolic flux analysis ((13)C-MFA) with [1,6-(13)C]glucose, [5-(13)C]xylose, and [1,6-(13)C]glucose+[5-(13)C]xylose as isotopic tracers. Compared to the starting strain, the evolved strain had an increased growth rate (~2-fold), increased biomass yield, increased tolerance to high temperatures up to 90°C, and gained the ability to grow on xylose in minimal medium. At the optimal growth temperature of 81°C, the maximum growth rate on glucose and xylose was 0.44 and 0.46h(-1), respectively. In medium containing glucose and xylose the strain efficiently co-utilized the two sugars. (13)C-MFA results provided insights into the metabolism of T. thermophilus LC113 that allows efficient co-utilization of glucose and xylose. Specifically, (13)C-MFA revealed that metabolic fluxes in the upper part of metabolism adjust flexibly to sugar availability, while fluxes in the lower part of metabolism remain relatively constant. Whole genome sequence analysis revealed two large structural changes that can help explain the physiology of the evolved strain: a duplication of a chromosome region that contains many sugar transporters, and a 5x multiplication of a region on the pVV8 plasmid that contains xylose isomerase and xylulokinase genes, the first two enzymes of xylose catabolism. Taken together, (13)C-MFA and genome sequence analysis provided complementary insights into the physiology of the evolved strain.

Comparison of gas chromatography-mass spectrometry and gas chromatography-combustion-isotope ratio mass spectrometry analysis for in vivo estimates of metabolic fluxes

Gas chromatography-mass spectrometry (GC-MS) was compared with gas chromatography-combustion-isotope ratio mass spectrometry (GC-C-IRMS) for measurements of cholesterol (13)C enrichment after infusion of labeled precursor ([(13)C1,2]acetate). Paired results were significantly correlated, although GC-MS was less accurate than GC-C-IRMS for higher enrichments. Nevertheless, only GC-MS was able to provide information on isotopologue distribution, bringing new insights to lipid metabolism. Therefore, we assessed the isotopologue distribution of cholesterol in humans and dogs known to present contrasted cholesterol metabolic pathways. The labeled tracer incorporation was different in both species, highlighting the subsidiarity of GC-MS and GC-C-IRMS to analyze in vivo stable isotope studies.

Catabolism of Branched Chain Amino Acids Contributes Significantly to Synthesis of Odd-Chain and Even-Chain Fatty Acids in 3T3-L1 Adipocytes

The branched chain amino acids (BCAA) valine, leucine and isoleucine have been implicated in a number of diseases including obesity, insulin resistance, and type 2 diabetes mellitus, although the mechanisms are still poorly understood. Adipose tissue plays an important role in BCAA homeostasis by actively metabolizing circulating BCAA. In this work, we have investigated the link between BCAA catabolism and fatty acid synthesis in 3T3-L1 adipocytes using parallel ¹³C-labeling experiments, mass spectrometry and model-based isotopomer data analysis. Specifically, we performed parallel labeling experiments with four fully ¹³C-labeled tracers, [U-¹³C]valine, [U-¹³C]leucine, [U-¹³C]isoleucine and [U-¹³C]glutamine. We measured mass isotopomer distributions of fatty acids and intracellular metabolites by GC-MS and analyzed the data using the isotopomer spectral analysis (ISA) framework. We demonstrate that 3T3-L1 adipocytes accumulate significant amounts of even chain length (C14:0, C16:0 and C18:0) and odd chain length (C15:0 and C17:0) fatty acids under standard cell culture conditions. Using a novel GC-MS method, we demonstrate that propionyl-CoA acts as the primer on fatty acid synthase for the production of odd chain fatty acids. BCAA contributed significantly to the production of all fatty acids. Leucine and isoleucine contributed at least 25% to lipogenic acetyl-CoA pool, and valine and isoleucine contributed 100% to lipogenic propionyl-CoA pool. Our results further suggest that low activity of methylmalonyl-CoA mutase and mass action kinetics of propionyl-CoA on fatty acid synthase result in high rates of odd chain fatty acid synthesis in 3T3-L1 cells. Overall, this work provides important new insights into the connection between BCAA catabolism and fatty acid synthesis in adipocytes and underscores the high capacity of adipocytes for metabolizing BCAA.

13C-Metabolic Flux Analysis: An Accurate Approach to Demystify Microbial Metabolism for Biochemical Production

Metabolic engineering of various industrial microorganisms to produce chemicals, fuels, and drugs has raised interest since it is environmentally friendly, sustainable, and independent of nonrenewable resources. However, microbial metabolism is so complex that only a few metabolic engineering efforts have been able to achieve a satisfactory yield, titer or productivity of the target chemicals for industrial commercialization. In order to overcome this challenge, ¹³C Metabolic Flux Analysis (¹³C-MFA) has been continuously developed and widely applied to rigorously investigate cell metabolism and quantify the carbon flux distribution in central metabolic pathways. In the past decade, many ¹³C-MFA studies have been performed in academic labs and biotechnology industries to pinpoint key issues related to microbe-based chemical production. Insightful information about the metabolic rewiring has been provided to guide the development of the appropriate metabolic engineering strategies for improving the biochemical production. In this review, we will introduce the basics of ¹³C-MFA and illustrate how ¹³C-MFA has been applied via integration with metabolic engineering to identify and tackle the rate-limiting steps in biochemical production for various host microorganisms

Complete genome sequence, metabolic model construction and phenotypic characterization of Geobacillus LC300, an extremely thermophilic, fast growing, xylose-utilizing bacterium

We have isolated a new extremely thermophilic fast-growing Geobacillus strain that can efficiently utilize xylose, glucose, mannose and galactose for cell growth. When grown aerobically at 72 °C, Geobacillus LC300 has a growth rate of 2.15 h(-1) on glucose and 1.52 h(-1) on xylose (doubling time less than 30 min). The corresponding specific glucose and xylose utilization rates are 5.55 g/g/h and 5.24 g/g/h, respectively. As such, Geobacillus LC300 grows 3-times faster than E. coli on glucose and xylose, and has a specific xylose utilization rate that is 3-times higher than the best metabolically engineered organism to date. To gain more insight into the metabolism of Geobacillus LC300 its genome was sequenced using PacBio's RS II single-molecule real-time (SMRT) sequencing platform and annotated using the RAST server. Based on the genome annotation and the measured biomass composition a core metabolic network model was constructed. To further demonstrate the biotechnological potential of this organism, Geobacillus LC300 was grown to high cell-densities in a fed-batch culture, where cells maintained a high xylose utilization rate under low dissolved oxygen concentrations. All of these characteristics make Geobacillus LC300 an attractive host for future metabolic engineering and biotechnology applications.

Parallel labeling experiments for pathway elucidation and (13)C metabolic flux analysis

Metabolic pathway models provide the foundation for quantitative studies of cellular physiology through the measurement of intracellular metabolic fluxes. For model organisms metabolic models are well established, with many manually curated genome-scale model reconstructions, gene knockout studies and stable-isotope tracing studies. However, for non-model organisms a similar level of knowledge is often lacking. Compartmentation of cellular metabolism in eukaryotic systems also presents significant challenges for quantitative (13)C-metabolic flux analysis ((13)C-MFA). Recently, innovative (13)C-MFA approaches have been developed based on parallel labeling experiments, the use of multiple isotopic tracers and integrated data analysis, that allow more rigorous validation of pathway models and improved quantification of metabolic fluxes. Applications of these approaches open new research directions in metabolic engineering, biotechnology and medicine.

A roadmap for interpreting (13)C metabolite labeling patterns from cells

Measuring intracellular metabolism has increasingly led to important insights in biomedical research. (13)C tracer analysis, although less information-rich than quantitative (13)C flux analysis that requires computational data integration, has been established as a time-efficient method to unravel relative pathway activities, qualitative changes in pathway contributions, and nutrient contributions. Here, we review selected key issues in interpreting (13)C metabolite labeling patterns, with the goal of drawing accurate conclusions from steady state and dynamic stable isotopic tracer experiments.

(13)C-metabolic flux analysis of co-cultures: A novel approach

In this work, we present a novel approach for performing (13)C metabolic flux analysis ((13)C-MFA) of co-culture systems. We demonstrate for the first time that it is possible to determine metabolic flux distributions in multiple species simultaneously without the need for physical separation of cells or proteins, or overexpression of species-specific products. Instead, metabolic fluxes for each species in a co-culture are estimated directly from isotopic labeling of total biomass obtained using conventional mass spectrometry approaches such as GC-MS. In addition to determining metabolic fluxes, this approach estimates the relative population size of each species in a mixed culture and inter-species metabolite exchange. As such, it enables detailed studies of microbial communities including species dynamics and interactions between community members. The methodology is experimentally validated here using a co-culture of two E. coli knockout strains. Taken together, this work greatly extends the scope of (13)C-MFA to a large number of multi-cellular systems that are of significant importance in biotechnology and medicine.

(13)C metabolic flux analysis of the extremely thermophilic, fast growing, xylose-utilizing Geobacillus strain LC300

Thermophiles are increasingly used as versatile hosts in the biotechnology industry. One of the key advantages of thermophiles is the potential to achieve high rates of feedstock conversion at elevated temperatures. The recently isolated Geobacillus strain LC300 grows extremely fast on xylose, with a doubling time of less than 30 min. In the accompanying paper, the genome of Geobacillus LC300 was sequenced and annotated. In this work, we have experimentally validated the metabolic network model using parallel (13)C-labeling experiments and applied (13)C-metabolic flux analysis to quantify precise metabolic fluxes. Specifically, the complete set of singly labeled xylose tracers, [1-(13)C], [2-(13)C], [3-(13)C], [4-(13)C], and [5-(13)C]xylose, was used for the first time. Isotopic labeling of biomass amino acids was measured by gas chromatography mass spectrometry (GC-MS). Isotopic labeling of carbon dioxide in the off-gas was also measured by an on-line mass spectrometer. The (13)C-labeling data was then rigorously integrated for flux elucidation using the COMPLETE-MFA approach. The results provided important new insights into the metabolism of Geobacillus LC300, its efficient xylose utilization pathways, and the balance between carbon, redox and energy fluxes. The pentose phosphate pathway, glycolysis and TCA cycle were found to be highly active in Geobacillus LC300. The oxidative pentose phosphate pathway was also active and contributed significantly to NADPH production. No transhydrogenase activity was detected. Results from this work provide a solid foundation for future studies of this strain and its metabolic engineering and biotechnological applications.

The oxidative pentose phosphate pathway is the primary source of NADPH for lipid overproduction from glucose in Yarrowia lipolytica

Oleaginous microbes represent an attractive means of converting a diverse range of feedstocks into oils that can be transesterified to biodiesel. However, the mechanism of lipid overproduction in these organisms is incompletely understood, hindering the development of strategies for engineering superior biocatalysts for "single-cell oil" production. In particular, it is unclear which pathways are used to generate the large quantities of NADPH required for overproduction of the highly reduced fatty acid species. While early studies implicated malic enzyme as having a key role in production of lipogenic NADPH in oleaginous fungi, several recent reports have cast doubts as to whether malic enzyme may contribute to production of lipogenic NADPH in the model oleaginous yeast Yarrowia lipolytica. To address this problem we have used (13)C-Metabolic Flux Analysis to estimate the metabolic flux distributions during lipid accumulation in two Y. lipolytica strains; a control strain and a previously published engineered strain capable of producing lipids at roughly twice the yield. We observe a dramatic rearrangement of the metabolic flux distribution in the engineered strain which supports lipid overproduction. The NADPH-producing flux through the oxidative Pentose Phosphate Pathway is approximately doubled in the engineered strain in response to the roughly two-fold increase in fatty acid biosynthesis, while the flux through malic enzyme does not differ significantly between the two strains. Moreover, the estimated rate of NADPH production in the oxidative Pentose Phosphate Pathway is in good agreement with the estimated rate of NADPH consumption in fatty acid biosynthesis in both strains. These results suggest the oxidative Pentose Phosphate Pathway is the primary source of lipogenic NADPH in Y. lipolytica.

Methods and advances in metabolic flux analysis: a mini-review

Metabolic flux analysis (MFA) is one of the pillars of metabolic engineering. Over the past three decades, it has been widely used to quantify intracellular metabolic fluxes in both native (wild type) and engineered biological systems. Through MFA, changes in metabolic pathway fluxes are quantified that result from genetic and/or environmental interventions. This information, in turn, provides insights into the regulation of metabolic pathways and may suggest new targets for further metabolic engineering of the strains. In this mini-review, we discuss and classify the various methods of MFA that have been developed, which include stoichiometric MFA, (13)C metabolic flux analysis, isotopic non-stationary (13)C metabolic flux analysis, dynamic metabolic flux analysis, and (13)C dynamic metabolic flux analysis. For each method, we discuss key advantages and limitations and conclude by highlighting important recent advances in flux analysis approaches.

Integrated 13C-metabolic flux analysis of 14 parallel labeling experiments in Escherichia coli

The use of parallel labeling experiments for (13)C metabolic flux analysis ((13)C-MFA) has emerged in recent years as the new gold standard in fluxomics. The methodology has been termed COMPLETE-MFA, short for complementary parallel labeling experiments technique for metabolic flux analysis. In this contribution, we have tested the limits of COMPLETE-MFA by demonstrating integrated analysis of 14 parallel labeling experiments with Escherichia coli. An effort on such a massive scale has never been attempted before. In addition to several widely used isotopic tracers such as [1,2-(13)C]glucose and mixtures of [1-(13)C]glucose and [U-(13)C]glucose, four novel tracers were applied in this study: [2,3-(13)C]glucose, [4,5,6-(13)C]glucose, [2,3,4,5,6-(13)C]glucose and a mixture of [1-(13)C]glucose and [4,5,6-(13)C]glucose. This allowed us for the first time to compare the performance of a large number of isotopic tracers. Overall, there was no single best tracer for the entire E. coli metabolic network model. Tracers that produced well-resolved fluxes in the upper part of metabolism (glycolysis and pentose phosphate pathways) showed poor performance for fluxes in the lower part of metabolism (TCA cycle and anaplerotic reactions), and vice versa. The best tracer for upper metabolism was 80% [1-(13)C]glucose+20% [U-(13)C]glucose, while [4,5,6-(13)C]glucose and [5-(13)C]glucose both produced optimal flux resolution in the lower part of metabolism. COMPLETE-MFA improved both flux precision and flux observability, i.e. more independent fluxes were resolved with smaller confidence intervals, especially exchange fluxes. Overall, this study demonstrates that COMPLETE-MFA is a powerful approach for improving flux measurements and that this methodology should be considered in future studies that require very high flux resolution.

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.

(13)C metabolic flux analysis of recombinant expression hosts

Identifying host cell metabolic phenotypes that promote high recombinant protein titer is a major goal of the biotech industry. (13)C metabolic flux analysis (MFA) provides a rigorous approach to quantify these metabolic phenotypes by applying isotope tracers to map the flow of carbon through intracellular metabolic pathways. Recent advances in tracer theory and measurements are enabling more information to be extracted from (13)C labeling experiments. Sustained development of publicly available software tools and standardization of experimental workflows is simultaneously encouraging increased adoption of (13)C MFA within the biotech research community. A number of recent (13)C MFA studies have identified increased citric acid cycle and pentose phosphate pathway fluxes as consistent markers of high recombinant protein expression, both in mammalian and microbial hosts. Further work is needed to determine whether redirecting flux into these pathways can effectively enhance protein titers while maintaining acceptable glycan profiles.

Quantifying biomass composition by gas chromatography/mass spectrometry

We developed a set of methods for the quantification of four major components of microbial biomass using gas chromatography/mass spectrometry (GC/MS). Specifically, methods are described to quantify amino acids, RNA, fatty acids, and glycogen, which comprise an estimated 88% of the dry weight of Escherichia coli. Quantification is performed by isotope ratio analysis with fully (13)C-labeled biomass as internal standard, which is generated by growing E. coli on [U-(13)C]glucose. This convenient, reliable, and accurate single-platform (GC/MS) workflow for measuring biomass composition offers significant advantages over existing methods. We demonstrate the consistency, accuracy, precision, and utility of this procedure by applying it to three metabolically unique E. coli strains. The presented methods will have widespread applicability in systems microbiology and bioengineering.

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.