Selection of tracers for 13C-metabolic flux analysis using elementary metabolite units (EMU) basis vector methodology

Metabolic flux analysis (MFA) is a powerful technique for elucidating in vivo fluxes in microbial and mammalian systems. A key step in (13)C-MFA is the selection of an appropriate isotopic tracer to observe fluxes in a proposed network model. Despite the importance of MFA in metabolic engineering and beyond, current approaches for tracer experiment design are still largely based on trial-and-error. The lack of a rational methodology for selecting isotopic tracers prevents MFA from achieving its full potential. Here, we introduce a new technique for tracer experiment design based on the concept of elementary metabolite unit (EMU) basis vectors. We demonstrate that any metabolite in a network model can be expressed as a linear combination of so-called EMU basis vectors, where the corresponding coefficients indicate the fractional contribution of the EMU basis vector to the product metabolite. The strength of this approach is the decoupling of substrate labeling, i.e. the EMU basis vectors, from the dependence on free fluxes, i.e. the coefficients. In this work, we demonstrate that flux observability inherently depends on the number of independent EMU basis vectors and the sensitivities of coefficients with respect to free fluxes. Specifically, the number of independent EMU basis vectors places hard limits on how many free fluxes can be determined in a model. This constraint is used as a guide for selecting feasible substrate labeling. In three example models, we demonstrate that by maximizing the number of independent EMU basis vectors the observability of a system is improved. Inspection of sensitivities of coefficients with respect to free fluxes provides additional constraints for proper selection of tracers. The present contribution provides a fresh perspective on an important topic in metabolic engineering, and gives practical guidelines and design principles for a priori selection of isotopic tracers for (13)C-MFA studies.     

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.     

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

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

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 ¹³C-metabolic flux analysis (¹³C-MFA). Here, cells were grown in parallel cultures with [1-¹³C]glucose and [U-¹³C]glucose as tracers and ¹³C-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 ¹³C-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 ¹³C-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 ¹³C-flux measurements in C. acetobutylicum.     

Elementary metabolite units (EMU): a novel framework for modeling isotopic distributions

Metabolic flux analysis (MFA) has emerged as a tool of great significance for metabolic engineering and mammalian physiology. An important limitation of MFA, as carried out via stable isotope labeling and GC/MS and nuclear magnetic resonance (NMR) measurements, is the large number of isotopomer or cumomer equations that need to be solved, especially when multiple isotopic tracers are used for the labeling of the system. This restriction reduces the ability of MFA to fully utilize the power of multiple isotopic tracers in elucidating the physiology of realistic situations comprising complex bioreaction networks. Here, we present a novel framework for the modeling of isotopic labeling systems that significantly reduces the number of system variables without any loss of information. The elementary metabolite unit (EMU) framework is based on a highly efficient decomposition method that identifies the minimum amount of information needed to simulate isotopic labeling within a reaction network using the knowledge of atomic transitions occurring in the network reactions. The functional units generated by the decomposition algorithm, called EMUs, form the new basis for generating system equations that describe the relationship between fluxes and stable isotope measurements. Isotopomer abundances simulated using the EMU framework are identical to those obtained using the isotopomer and cumomer methods, however, require significantly less computation time. For a typical (13)C-labeling system the total number of equations that needs to be solved is reduced by one order-of-magnitude (100s EMUs vs. 1000s isotopomers). As such, the EMU framework is most efficient for the analysis of labeling by multiple isotopic tracers. For example, analysis of the gluconeogenesis pathway with (2)H, (13)C, and (18)O tracers requires only 354 EMUs, compared to more than two million isotopomers.     

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

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

13C metabolic flux analysis of microbial and mammalian systems is enhanced with GC-MS measurements of glycogen and RNA labeling

¹³C metabolic flux analysis (¹³C-MFA) is a widely used tool for quantitative analysis of microbial and mammalian metabolism. Until now, ¹³C-MFA was based mainly on measurements of isotopic labeling of amino acids derived from hydrolyzed biomass proteins and isotopic labeling of extracted intracellular metabolites. Here, we demonstrate that isotopic labeling of glycogen and RNA, measured with gas chromatography-mass spectrometry (GC-MS), provides valuable additional information for ¹³C-MFA. Specifically, we demonstrate that isotopic labeling of glucose moiety of glycogen and ribose moiety of RNA greatly enhances resolution of metabolic fluxes in the upper part of metabolism; importantly, these measurements allow precise quantification of net and exchange fluxes in the pentose phosphate pathway. To demonstrate the practical importance of these measurements for ¹³C-MFA, we have used Escherichia coli as a model microbial system and CHO cells as a model mammalian system. Additionally, we have applied this approach to determine metabolic fluxes of glucose and xylose co-utilization in the E. coli ΔptsG mutant. The convenience of measuring glycogen and RNA, which are stable and abundant in microbial and mammalian cells, offers the following key advantages: reduced sample size, no quenching required, no extractions required, and GC-MS can be used instead of more costly LC-MS/MS techniques. Overall, the presented approach for ¹³C-MFA will have widespread applicability in metabolic engineering and biomedical research.     

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.     

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 ¹³C-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 ¹³C metabolic flux analysis (¹³C-MFA). Specifically, we quantified precise metabolic fluxes for each condition by performing parallel labeling experiments and analyzing the data through integrated ¹³C-MFA using the optimal tracers [1,2-¹³C]glucose, [1,6-¹³C]glucose, [1,2-¹³C]xylose and [5-¹³C]xylose. We also quantified changes in biomass composition and confirmed turnover of macromolecules by applying [U-¹³C]glucose and [U-¹³C]xylose tracers. We demonstrated that under anaerobic growth conditions there is significant turnover of lipids and that a significant portion of CO₂ 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/FADH₂, 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.     

Re-examination of metabolic fluxes in <Emphasis Type="Italic">Escherichia coli</Emphasis> during anaerobic fermentation of glucose using <Superscript>13</Superscript>C labeling experiments and 2-dimensional nuclear magnetic resonance (NMR) spectroscopy

Improved design of metabolic flux estimation using mixed label ¹³C labeling experiments and identifiability analysis motivated re-examination of metabolic fluxes during anaerobic fermentation in the Escherichia coli. Comprehensive metabolic flux maps were determined by using a mixture of differently labeled glucose and compared to conventional flux maps obtained using extracellular measurements and comprehensive metabolic flux maps obtained using only U-¹³C glucose as the substrate. As expected, conventional flux analysis performs poorly in comparison to ¹³C-MFA, especially in the Embden-Meyerhof-Parnas (EMP) and pentose phosphate (PP) pathways. Identifiability analysis indicated and experiments confirmed that a mixture of 10% U-^l3C glucose, 25% 1-¹³C glucose, and 65% naturally labeled glucose significantly improved the statistical quality of all calculated fluxes in the PP pathway, the EMP pathway, the anaplerotic reactions, and the tricarboxylic acid cycle. Modifying the network topology for the presence and absence of the Entner-Doudoroff pathway and the glyoxylate shunt did not affect the value or quality of estimated fluxes significantly. Extracellular measurement of formate production was necessary for the accurate estimation of the fluxes around the formate node.     

Metabolic flux analysis of CHO cells at growth and non-growth phases using isotopic tracers and mass spectrometry

Chinese hamster ovary (CHO) cells are the main platform for production of biotherapeutics in the biopharmaceutical industry. However, relatively little is known about the metabolism of CHO cells in cell culture. In this work, metabolism of CHO cells was studied at the growth phase and early stationary phase using isotopic tracers and mass spectrometry. CHO cells were grown in fed-batch culture over a period of six days. On days 2 and 4, [1,2-¹³C] glucose was introduced and the labeling of intracellular metabolites was measured by gas chromatography-mass spectrometry (GC–MS) at 6, 12 and 24 h following the introduction of tracer. Intracellular metabolic fluxes were quantified from measured extracellular rates and ¹³C-labeling dynamics of intracellular metabolites using non-stationary ¹³C-metabolic flux analysis (¹³C-MFA). The flux results revealed significant rewiring of intracellular metabolic fluxes in the transition from growth to non-growth, including changes in energy metabolism, redox metabolism, oxidative pentose phosphate pathway and anaplerosis. At the exponential phase, CHO cell metabolism was characterized by a high flux of glycolysis from glucose to lactate, anaplerosis from pyruvate to oxaloacetate and from glutamate to α-ketoglutarate, and cataplerosis though malic enzyme. At the stationary phase, the flux map was characterized by a reduced flux of glycolysis, net lactate uptake, oxidative pentose phosphate pathway flux, and reduced rate of anaplerosis. The fluxes of pyruvate dehydrogenase and TCA cycle were similar at the exponential and stationary phase. The results presented here provide a solid foundation for future studies of CHO cell metabolism for applications such as cell line development and medium optimization for high-titer production of recombinant proteins.     

Automated Bayesian model development for frequency detection in biological time series

Background

The ability to perform quantitative studies using isotope tracers and metabolic flux analysis (MFA) is critical for detecting pathway bottlenecks and elucidating network regulation in biological systems, especially those that have been engineered to alter their native metabolic capacities. Mathematically, MFA models are traditionally formulated using separate state variables for reaction fluxes and isotopomer abundances. Analysis of isotope labeling experiments using this set of variables results in a non-convex optimization problem that suffers from both implementation complexity and convergence problems.

Results

This article addresses the mathematical and computational formulation of ¹³C MFA models using a new set of variables referred to as fluxomers. These composite variables combine both fluxes and isotopomer abundances, which results in a simply-posed formulation and an improved error model that is insensitive to isotopomer measurement normalization. A powerful fluxomer iterative algorithm (FIA) is developed and applied to solve the MFA optimization problem. For moderate-sized networks, the algorithm is shown to outperform the commonly used 13CFLUX cumomer-based algorithm and the more recently introduced OpenFLUX software that relies upon an elementary metabolite unit (EMU) network decomposition, both in terms of convergence time and output variability.

Conclusions

Substantial improvements in convergence time and statistical quality of results can be achieved by applying fluxomer variables and the FIA algorithm to compute best-fit solutions to MFA models. We expect that the fluxomer formulation will provide a more suitable basis for future algorithms that analyze very large scale networks and design optimal isotope labeling experiments.     

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

Radioactive and stable isotopes have been applied for decades to elucidate metabolic pathways and quantify carbon flow in cellular systems using mass and isotope balancing approaches. Isotope-labeling experiments can be conducted as a single tracer experiment, or as parallel labeling experiments. In the latter case, several experiments are performed under identical conditions except for the choice of substrate labeling. In this review, we highlight robust approaches for probing metabolism and addressing metabolically related questions though parallel labeling experiments. In the first part, we provide a brief historical perspective on parallel labeling experiments, from the early metabolic studies when radioisotopes were predominant to present-day applications based on stable-isotopes. We also elaborate on important technical and theoretical advances that have facilitated the transition from radioisotopes to stable-isotopes. In the second part of the review, we focus on parallel labeling experiments for ¹³C-metabolic flux analysis (¹³C-MFA). Parallel experiments offer several advantages that include: tailoring experiments to resolve specific fluxes with high precision; reducing the length of labeling experiments by introducing multiple entry-points of isotopes; validating biochemical network models; and improving the performance of ¹³C-MFA in systems where the number of measurements is limited. We conclude by discussing some challenges facing the use of parallel labeling experiments for ¹³C-MFA and highlight the need to address issues related to biological variability, data integration, and rational tracer selection.     

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 ¹³C metabolic flux analysis (¹³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-¹³C]glucose+[4,5,6-¹³C]glucose, [2-¹³C]glucose+[4,5,6-¹³C]glucose, and [3-¹³C]glucose+[4,5,6-¹³C]glucose and developed a new method to measure ¹³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 13C 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 ¹³C-metabolic flux analysis (¹³C-MFA) with [1,6-¹³C]glucose, [5-¹³C]xylose, and [1,6-¹³C]glucose+[5-¹³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.46 h⁻¹, respectively. In medium containing glucose and xylose the strain efficiently co-utilized the two sugars. ¹³C-MFA results provided insights into the metabolism of T. thermophilus LC113 that allows efficient co-utilization of glucose and xylose. Specifically, ¹³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, ¹³C-MFA and genome sequence analysis provided complementary insights into the physiology of the evolved strain.     

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.     

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 ¹³C metabolic flux analysis (¹³C-MFA) to estimate intracellular fluxes using parallel labeling experiments with the optimal tracers [1,2−¹³C]glucose and [1,6−¹³C]glucose. During exponential growth in glucose minimal medium, V. natriegens had a growth rate of 1.70 1/h (doubling time of 24 min) and a glucose uptake rate of 3.90 g/g/h, which is more than two 2-fold faster than E. coli, although slower than the fast-growing thermophile Geobacillus LC300. ¹³C-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.     

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

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

Induced pluripotent stem cells can utilize lactate as a metabolic substrate to support proliferation

Human-induced pluripotent stem cells (iPSCs) hold the promise to improve cell-based therapies. Yet, to meet rising demands and become clinically impactful, sufficient high-quality iPSCs in quantity must be generated, a task that exceeds current capabilities. In this study, K3 iPSCs cultures were examined using parallel-labeling metabolic flux analysis (¹³C-MFA) to quantify intracellular fluxes at relevant bioprocessing stages: glucose concentrations representative of initial media concentrations and high lactate concentrations representative of fed-batch culture conditions, prior to and after bolus glucose feeds. The glucose and lactate concentrations are also representative of concentrations that might be encountered at different locations within 3D cell aggregates. Furthermore, a novel method was developed to allow the isotopic tracer [U-¹³C₃] lactate to be used in the ¹³C-MFA model. The results indicated that high extracellular lactate concentrations decreased glucose consumption and lactate production, while glucose concentrations alone did not affect rates of aerobic glycolysis. Moreover, for the high lactate cultures, lactate was used as a metabolic substrate to support oxidative mitochondrial metabolism. These results demonstrate that iPSCs have metabolic flexibility and possess the capacity to metabolize lactate to support exponential growth, and that high lactate concentrations alone do not adversely impact iPSC proliferation.     

Application of radiolabeled tracers to biocatalytic flux analysis.

Radiolabeled tracers can provide valuable information about the structure of and flux distributions in biocatalytic reaction networks. This method derives from prior studies of glucose metabolism in mammalian systems and is implemented by pulsing a culture with a radiolabeled metabolite that can be transported into the cells and subsequently measuring the radioactivity of all network metabolites following separation by liquid chromatography. Intracellular fluxes can be directly determined from the transient radioactivity count data by tracking the depletion of the radiolabeled metabolite and/or the accompanying accumulation of any products formed. This technique differs from previous methods in that it is applied within a systems approach to the problem of flux determination. It has been used for the investigation of the indene bioconversion network expressed in Rhodococcus sp. KY1. Flux estimates obtained by radioactive tracers were confirmed by macroscopic metabolite balancing and showed that indene oxidation in steady state chemostat cultures proceeds primarily through a monooxygenase activity forming (1S,2R)-indan oxide, with no dehydrogenation of trans-(1R,2R)-indandiol. These results confirmed the significance of indan oxide formation and identified the hydrolysis of indan oxide as a key step in maximizing the production of (2R)-indandiol, a chiral precursor of the HIV protease inhibitor, Crixivan.