Parallel Labeling Experiments Research Articles

Abstract 6612: Therapy-resistant endocycling cancer cells alter mitochondrial structure and metabolism to survive therapeutic stress

Abstract Therapy resistance is responsible for 90% of cancer-related deaths. We identified a novel mechanism of therapy resistance in which cancer cells enter an endocycling cell state, undergoing repeated S and G phases without dividing. As a result, resistant endocycling cells are 40 times larger than mitotic cancer cell controls. Interestingly, we observed that these non-proliferative endocycling cells can eventually give rise to proliferative progeny weeks after chemotherapy is removed, modeling cancer recurrence. We hypothesize that this endocycling cancer cell state is a critical driver of therapy resistance and cancer lethality. We evaluated oxidative stress in endocycling cells that survived the days and weeks following chemotherapy treatment using DCF-DA staining. We found increased levels of reactive oxygen species (ROS) in endocycling prostate cancer cells compared to untreated mitotic cells when normalized to cell size. Given that mitochondria are major producers and regulators of ROS, we hypothesize that cells surviving in this resistant state alter mitochondrial structure and function in response to increased oxidative stress. We performed immunofluorescence and confocal microscopy on mitochondria in untreated cells and cells following chemotherapy treatment to observe their morphology. Surviving endocycling cells had increased mitochondrial fragmentation when compared to the untreated group, a result that was further supported by reduced expression of phospho-DRP1 Ser637 via western blot. To assess changes in glucose metabolism, we then utilized isotope tracing-based metabolomics with [U-13C]glucose on untreated and surviving endocycling cells following chemotherapy. Endocycling cells had decreased labeling of 13C in TCA cycle metabolites, suggesting a decrease in mitochondrial function when compared to untreated cells. Ongoing work involves further studying metabolic reprogramming as cells enter this resistant state with parallel labeling experiments using [1,2-13C]glucose and [U-13C]glutamine to perform mathematical modeling-based 13C-metabolic flux analysis. In addition, we engineered a cell line to track mitochondria in real-time to quantify mitochondrial dynamics underlying the morphology observed in untreated and surviving cells. Our goal is to characterize how this therapy-resistant cell state alters its mitochondrial structure and function to survive therapeutic stress. Identifying vulnerabilities of this phenotype will enable new approaches for eliminating drug resistance in cancer. Citation Format: Melvin Li, Sarah R. Amend, Kenneth J. Pienta. Therapy-resistant endocycling cancer cells alter mitochondrial structure and metabolism to survive therapeutic stress [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 6612.

13C-Metabolic flux analysis of 3T3-L1 adipocytes illuminates its core metabolism under hypoxia

Hypoxia has been identified as a major factor in the pathogenesis of adipose tissue inflammation, which is a hallmark of obesity and obesity-linked type 2 diabetes mellitus. In this study, we have investigated the impact of hypoxia (1% oxygen) on the physiology and metabolism of 3T3-L1 adipocytes, a widely used cell culture model of adipose. Specifically, we applied parallel labeling experiments, isotopomer spectral analysis, and 13C-metabolic flux analysis to quantify the impact of hypoxia on adipogenesis, de novo lipogenesis and metabolic flux reprogramming in adipocytes. We found that 3T3-L1 cells can successfully differentiate into lipid-accumulating adipocytes under hypoxia, although the production of lipids was reduced by about 40%. Quantitative flux analysis demonstrated that short-term (1 day) and long-term (7 days) exposure to hypoxia resulted in similar reprogramming of cellular metabolism. Overall, we found that hypoxia: 1) reduced redox and energy generation by more than 2-fold and altered the patterns of metabolic pathway contributions to production and consumption of energy and redox cofactors; 2) redirected glucose metabolism from pentose phosphate pathway and citric acid cycle to lactate production; 3) rewired glutamine metabolism, from net glutamine production to net glutamine catabolism; 4) suppressed branched chain amino acid consumption; and 5) reduced biosynthesis of odd-chain fatty acids and mono-unsaturated fatty acids, while synthesis of saturated even-chain fatty acids was not affected. Together, these results highlight the profound impact of extracellular microenvironment on adipocyte metabolic activity and function.

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.

Coordinated reprogramming of metabolism and cell function in adipocytes from proliferation to differentiation

Adipose tissue plays a major role in regulating lipid and energy homeostasis by storing excess nutrients, releasing energetic substrates through lipolysis, and regulating metabolism of other tissues and organs through endocrine and paracrine signaling. Adipocytes within fat tissues store excess nutrients through increased cell number (hyperplasia), increased cell size (hypertrophy), or both. The differentiation of pre-adipocytes into mature lipid-accumulating adipocytes requires a complex interaction of metabolic pathways that is still incompletely understood. Here, we applied parallel labeling experiments and 13C-metabolic flux analysis to quantify precise metabolic fluxes in proliferating and differentiated 3T3-L1 cells, a widely used model to study adipogenesis. We found that morphological and biomass composition changes in adipocytes were accompanied by significant shifts in metabolic fluxes, encompassing all major metabolic pathways. In contrast to proliferating cells, differentiated adipocytes 1) increased glucose uptake and redirected glucose utilization from lactate production to lipogenesis and energy generation; 2) increased pathway fluxes through glycolysis, oxidative pentose phosphate pathway and citric acid cycle; 3) reduced lactate secretion, resulting in increased ATP generation via oxidative phosphorylation; 4) rewired glutamine metabolism, from glutaminolysis to de novo glutamine synthesis; 5) increased cytosolic NADPH production, driven mostly by increased cytosolic malic enzyme flux; 6) increased production of monounsaturated C16:1; and 7) activated a mitochondrial pyruvate cycle through simultaneous activity of pyruvate carboxylase, malate dehydrogenase and malic enzyme. Taken together, these results quantitatively highlight the complex interplay between pathway fluxes and cell function in adipocytes, and suggest a functional role for metabolic reprogramming in adipose differentiation and lipogenesis.

Non-conventional pathways enable pennycress (Thlaspi arvense L.) embryos to achieve high efficiency of oil biosynthesis.

Pennycress (Thlaspi arvense L.) accumulates oil up to 35% of the total seed biomass, and its overall fatty acid composition is suitable for aviation fuel. However, for this plant to become economically viable, its oil production needs to be improved. In vivo culture conditions that resemble the development of pennycress embryos in planta were developed based on the composition of the liquid endosperm. Then, substrate uptake rates and biomass accumulation were measured from cultured pennycress embryos, revealing a biosynthetic efficiency of 93%, which is one of the highest in comparison with other oilseeds to date. Additionally, the ratio of carbon in oil to CO2 indicated that non-conventional pathways are likely to be responsible for such a high carbon conversion efficiency. To identify the reactions enabling this phenomenon, parallel labeling experiments with 13C-labeled substrates were conducted in pennycress embryos. The main findings of these labeling experiments include: (i) the occurrence of the oxidative reactions of the pentose phosphate pathway in the cytosol; (ii) the reversibility of isocitrate dehydrogenase; (iii) the operation of the plastidic NADP-dependent malic enzyme; and (iv) the refixation of CO2 by Rubisco. These reactions are key providers of carbon and reductant for fatty acid synthesis and elongation.

13C-Metabolic Flux Analysis in Developing Flax (Linum usitatissinum L.) Embryos to Understand Storage Lipid Biosynthesis.

Flax (Linum usitatissinum L.) oil is an important source of α-linolenic (C18:3 ω-3). This polyunsaturated fatty acid is well known for its nutritional role in human and animal diets. Understanding storage lipid biosynthesis in developing flax embryos can lead to an increase in seed yield via marker-assisted selection. While a tremendous amount of work has been done on different plant species to highlight their metabolism during embryo development, a comprehensive analysis of metabolic flux in flax is still lacking. In this context, we have utilized in vitro cultured developing embryos of flax and determined net fluxes by performing three complementary parallel labeling experiments with 13C-labeled glucose and glutamine. Metabolic fluxes were estimated by computer-aided modeling of the central metabolic network including 11 cofactors of 118 reactions of the central metabolism and 12 pseudo-fluxes. A focus on lipid storage biosynthesis and the associated pathways was done in comparison with rapeseed, arabidopsis, maize and sunflower embryos. In our hands, glucose was determined to be the main source of carbon in flax embryos, leading to the conversion of phosphoenolpyruvate to pyruvate. The oxidative pentose phosphate pathway (OPPP) was identified as the producer of NADPH for fatty acid biosynthesis. Overall, the use of 13C-metabolic flux analysis provided new insights into the flax embryo metabolic processes involved in storage lipid biosynthesis. The elucidation of the metabolic network of this important crop plant reinforces the relevance of the application of this technique to the analysis of complex plant metabolic systems.

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 13C metabolic flux ratio studies. Here, we designed a method based on parallel labeling experiments (PLEs), using [1-13C]- and [4-13C]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 13C 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.

High-resolution 13C metabolic flux analysis.

Precise quantification of metabolic pathway fluxes in biological systems is of major importance in guiding efforts in metabolic engineering, biotechnology, microbiology, human health, and cell culture. 13C metabolic flux analysis (13C-MFA) is the predominant technique used for determining intracellular fluxes. Here, we present a protocol for 13C-MFA that incorporates recent advances in parallel labeling experiments, isotopic labeling measurements, and statistical analysis, as well as best practices developed through decades of experience. Experimental design to ensure that fluxes are estimated with the highest precision is an integral part of the protocol. The protocol is based on growing microbes in two (or more) parallel cultures with 13C-labeled glucose tracers, followed by gas chromatography-mass spectrometry (GC-MS) measurements of isotopic labeling of protein-bound amino acids, glycogen-bound glucose, and RNA-bound ribose. Fluxes are then estimated using software for 13C-MFA, such as Metran, followed by comprehensive statistical analysis to determine the goodness of fit and calculate confidence intervals of fluxes. The presented protocol can be completed in 4 d and quantifies metabolic fluxes with a standard deviation of ≤2%, a substantial improvement over previous implementations. The presented protocol is exemplified using an Escherichia coli ΔtpiA case study with full supporting data, providing a hands-on opportunity to step through a complex troubleshooting scenario. Although applications to prokaryotic microbial systems are emphasized, this protocol can be easily adjusted for application to eukaryotic organisms.

Metabolism in dense microbial colonies: 13C metabolic flux analysis of E. coli grown on agar identifies two distinct cell populations with acetate cross-feeding

In this study, we have investigated for the first time the metabolism of E. coli grown on agar using 13C metabolic flux analysis (13C-MFA). To date, all 13C-MFA studies on microbes have been performed with cells grown in liquid culture. Here, we extend the scope of 13C-MFA to biological systems where cells are grown in dense microbial colonies. First, we identified new optimal 13C tracers to quantify fluxes in systems where the acetate yield cannot be easily measured. We determined that three parallel labeling experiments with the tracers [1,2–13C]glucose, [1,6–13C]glucose, and [4,5,6–13C]glucose permit precise estimation of not only intracellular fluxes, but also of the amount of acetate produced from glucose. Parallel labeling experiments were then performed with wild-type E. coli and E. coli ΔackA grown in liquid culture and on agar plates. Initial attempts to fit the labeling data from wild-type E. coli grown on agar did not produce a statistically acceptable fit. To resolve this issue, we employed the recently developed co-culture 13C-MFA approach, where two E. coli subpopulations were defined in the model that engaged in metabolite cross-feeding. The flux results identified two distinct E. coli cell populations, a dominant cell population (92% of cells) that metabolized glucose via conventional metabolic pathways and secreted a large amount of acetate (~40% of maximum theoretical yield), and a second smaller cell population (8% of cells) that consumed the secreted acetate without any glucose influx. These experimental results are in good agreement with recent theoretical simulations. Importantly, this study provides a solid foundation for future investigations of a wide range of problems involving microbial biofilms that are of great interest in biotechnology, ecology and medicine, where metabolite cross-feeding between cell populations is a core feature of the communities.

Metabolism of the fast-growing bacterium Vibrio natriegens elucidated by 13C metabolic flux analysis

Vibrio natriegens is a fast-growing, non-pathogenic bacterium that is being considered as the next-generation workhorse for the biotechnology industry. However, little is known about the metabolism of this organism which is limiting our ability to apply rational metabolic engineering strategies. To address this critical gap in current knowledge, here we have performed a comprehensive analysis of V. natriegens metabolism. We constructed a detailed model of V. natriegens core metabolism, measured the biomass composition, and performed high-resolution 13C metabolic flux analysis (13C-MFA) to estimate intracellular fluxes using parallel labeling experiments with the optimal tracers [1,2−13C]glucose and [1,6−13C]glucose. During exponential growth in glucose minimal medium, V. natriegens had a growth rate of 1.70 1/h (doubling time of 24min) and a glucose uptake rate of 3.90g/g/h, which is more than two 2-fold faster than E. coli, although slower than the fast-growing thermophile Geobacillus LC300. 13C-MFA revealed that the core metabolism of V. natriegens is similar to that of E. coli, with the main difference being a 33% lower normalized flux through the oxidative pentose phosphate pathway. Quantitative analysis of co-factor balances provided additional insights into the energy and redox metabolism of V. natriegens. Taken together, the results presented in this study provide valuable new information about the physiology of V. natriegens and establish a solid foundation for future metabolic engineering efforts with this promising microorganism.

13C metabolic flux analysis of three divergent extremely thermophilic bacteria: Geobacillus sp. LC300, Thermus thermophilus HB8, and Rhodothermus marinus DSM 4252

Thermophilic organisms are being increasingly investigated and applied in metabolic engineering and biotechnology. The distinct metabolic and physiological characteristics of thermophiles, including broad substrate range and high uptake rates, coupled with recent advances in genetic tool development, present unique opportunities for strain engineering. However, poor understanding of the cellular physiology and metabolism of thermophiles has limited the application of systems biology and metabolic engineering tools to these organisms. To address this concern, we applied high resolution 13C metabolic flux analysis to quantify fluxes for three divergent extremely thermophilic bacteria from separate phyla: Geobacillus sp. LC300, Thermus thermophilus HB8, and Rhodothermus marinus DSM 4252. We performed 18 parallel labeling experiments, using all singly labeled glucose tracers for each strain, reconstructed and validated metabolic network models, measured biomass composition, and quantified precise metabolic fluxes for each organism. In the process, we resolved many uncertainties regarding gaps in pathway reconstructions and elucidated how these organisms maintain redox balance and generate energy. Overall, we found that the metabolisms of the three thermophiles were highly distinct, suggesting that adaptation to growth at high temperatures did not favor any particular set of metabolic pathways. All three strains relied heavily on glycolysis and TCA cycle to generate key cellular precursors and cofactors. None of the investigated organisms utilized the Entner-Doudoroff pathway and only one strain had an active oxidative pentose phosphate pathway. Taken together, the results from this study provide a solid foundation for future model building and engineering efforts with these and related thermophiles.

Comprehensive analysis of glucose and xylose metabolism in Escherichia coli under aerobic and anaerobic conditions by 13C metabolic flux analysis

Glucose and xylose are the two most abundant sugars derived from the breakdown of lignocellulosic biomass. While aerobic glucose metabolism is relatively well understood in E. coli, until now there have been only a handful of studies focused on anaerobic glucose metabolism and no 13C-flux studies on xylose metabolism. In the absence of experimentally validated flux maps, constraint-based approaches such as MOMA and RELATCH cannot be used to guide new metabolic engineering designs. In this work, we have addressed this critical gap in current understanding by performing comprehensive characterizations of glucose and xylose metabolism under aerobic and anaerobic conditions, using recent state-of-the-art techniques in 13C metabolic flux analysis (13C-MFA). Specifically, we quantified precise metabolic fluxes for each condition by performing parallel labeling experiments and analyzing the data through integrated 13C-MFA using the optimal tracers [1,2-13C]glucose, [1,6-13C]glucose, [1,2-13C]xylose and [5-13C]xylose. We also quantified changes in biomass composition and confirmed turnover of macromolecules by applying [U-13C]glucose and [U-13C]xylose tracers. We demonstrated that under anaerobic growth conditions there is significant turnover of lipids and that a significant portion of CO2 originates from biomass turnover. Using knockout strains, we also demonstrated that β-oxidation is critical for anaerobic growth on xylose. Quantitative analysis of co-factor balances (NADH/FADH2, NADPH, and ATP) for different growth conditions provided new insights regarding the interplay of energy and redox metabolism and the impact on E. coli cell physiology.

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

13C-Metabolic flux analysis (13C-MFA) is a widely used approach in metabolic engineering for quantifying intracellular metabolic fluxes. The precision of fluxes determined by 13C-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 13C-MFA. The best single tracers were doubly 13C-labeled glucose tracers, including [1,6-13C]glucose, [5,6-13C]glucose and [1,2-13C]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-13C]glucose and [1,2-13C]glucose. Combined analysis of [1,6-13C]glucose and [1,2-13C]glucose labeling data improved the flux precision score by nearly 20-fold compared to widely use tracer mixture 80% [1-13C]glucose +20% [U-13C]glucose.

Evidence for transketolase-like TKTL1 flux in CHO cells based on parallel labeling experiments and 13C-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 13C metabolic flux analysis (13C-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-13C]glucose+[4,5,6-13C]glucose, [2-13C]glucose+[4,5,6-13C]glucose, and [3-13C]glucose+[4,5,6-13C]glucose and developed a new method to measure 13C-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.

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 13C-labeling experiments, mass spectrometry and model-based isotopomer data analysis. Specifically, we performed parallel labeling experiments with four fully 13C-labeled tracers, [U-13C]valine, [U-13C]leucine, [U-13C]isoleucine and [U-13C]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.

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

The use of parallel labeling experiments for 13C metabolic flux analysis (13C-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-13C]glucose and mixtures of [1-13C]glucose and [U-13C]glucose, four novel tracers were applied in this study: [2,3-13C]glucose, [4,5,6-13C]glucose, [2,3,4,5,6-13C]glucose and a mixture of [1-13C]glucose and [4,5,6-13C]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-13C]glucose+20% [U-13C]glucose, while [4,5,6-13C]glucose and [5-13C]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.

BackgroundSteady-state 13C-based metabolic flux analysis (13C-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 13C-labeled substrate (tracer); (iii) measuring the extracellular effluxes and detecting the 13C-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, 13C-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 13C-MFA is an important factor for PLE implementation.ResultsThe 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 13C-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.ConclusionsThe developed open-source OpenFLUX2 provides a friendly software environment that facilitates beginners and existing OpenFLUX users to implement LEs for steady-state 13C-MFA including experimental design, quantitative evaluation of flux parameters and statistics.Electronic supplementary materialThe online version of this article (doi:10.1186/s12934-014-0152-x) contains supplementary material, which is available to authorized users.

Parallel labeling experiments validate Clostridium acetobutylicum metabolic network model for 13C metabolic flux analysis

In this work, we provide new insights into the metabolism of Clostridium acetobutylicum ATCC 824 obtained using a systematic approach for quantifying fluxes based on parallel labeling experiments and 13C-metabolic flux analysis (13C-MFA). Here, cells were grown in parallel cultures with [1-13C]glucose and [U-13C]glucose as tracers and 13C-MFA was used to quantify intracellular metabolic fluxes. Several metabolic network models were compared: an initial model based on current knowledge, and extended network models that included additional reactions that improved the fits of experimental data. While the initial network model did not produce a statistically acceptable fit of 13C-labeling data, an extended network model with five additional reactions was able to fit all data with 292 redundant measurements. The model was subsequently trimmed to produce a minimal network model of C. acetobutylicum for 13C-MFA, which could still reproduce all of the experimental data. The flux results provided valuable new insights into the metabolism of C. acetobutylicum. First, we found that TCA cycle was effectively incomplete, as there was no measurable flux between α-ketoglutarate and succinyl-CoA, succinate and fumarate, and malate and oxaloacetate. Second, an active pathway was identified from pyruvate to fumarate via aspartate. Third, we found that isoleucine was produced exclusively through the citramalate synthase pathway in C. acetobutylicum and that CAC3174 was likely responsible for citramalate synthase activity. These model predictions were confirmed in several follow-up tracer experiments. The validated metabolic network model established in this study can be used in future investigations for unbiased 13C-flux measurements in C. acetobutylicum.

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

BackgroundSteady-state 13C-based metabolic flux analysis (13C-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 13C-labeled substrate (tracer); (iii) measuring the extracellular effluxes and detecting the 13C-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, 13C-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 13C-MFA is an important factor for PLE implementation.ResultsThe 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 13C-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.ConclusionsThe developed open-source OpenFLUX2 provides a friendly software environment that facilitates beginners and existing OpenFLUX users to implement LEs for steady-state 13C-MFA including experimental design, quantitative evaluation of flux parameters and statistics.

COMPLETE-MFA: Complementary parallel labeling experiments technique for metabolic flux analysis

We have developed a novel approach for measuring highly accurate and precise metabolic fluxes in living cells, termed COMPLETE-MFA, short for complementary parallel labeling experiments technique for metabolic flux analysis. The COMPLETE-MFA method is based on combined analysis of multiple isotopic labeling experiments, where the synergy of using complementary tracers greatly improves the precision of estimated fluxes. In this work, we demonstrate the COMPLETE-MFA approach using all singly labeled glucose tracers, [1-13C], [2-13C], [3-13C], [4-13C], [5-13C], and [6-13C]glucose to determine precise metabolic fluxes for wild-type Escherichia coli. Cells were grown in six parallel cultures on defined medium with glucose as the only carbon source. Mass isotopomers of biomass amino acids were measured by gas chromatography–mass spectrometry (GC–MS). The data from all six experiments were then fitted simultaneously to a single flux model to determine accurate intracellular fluxes. We obtained a statistically acceptable fit with more than 300 redundant measurements. The estimated flux map is the most precise flux result obtained thus far for E. coli cells. To our knowledge, this is the first time that six isotopic labeling experiments have been successfully integrated for high-resolution 13C-flux analysis.