Metabolic flux and metabolic network analysis of Penicillium chrysogenum using 2D [13C, 1H] COSY NMR measurements and cumulative bondomer simulation

At present two alternative methods are available for analyzing the fluxes in a metabolic network: (1) combining measurements of net conversion rates with a set of metabolite balances including the cofactor balances, or (2) leaving out the cofactor balances and fitting the resulting free fluxes to measured (13)C-labeling data. In this study these two approaches are applied to the fluxes in the glycolysis and pentose phosphate pathway of Penicillium chrysogenum growing on either ammonia or nitrate as the nitrogen source, which is expected to give different pentose phosphate pathway fluxes. The presented flux analyses are based on extensive sets of 2D [(13)C, (1)H] COSY data. A new concept is applied for simulation of this type of (13)C-labeling data: cumulative bondomer modeling. The outcomes of the (13)C-labeling based flux analysis substantially differ from those of the pure metabolite balancing approach. The fluxes that are determined using (13)C-labeling data are shown to be highly dependent on the chosen metabolic network. Extending the traditional nonoxidative pentose phosphate pathway with additional transketolase and transaldolase reactions, extending the glycolysis with a fructose 6-phosphate aldolase/dihydroxyacetone kinase reaction sequence or adding a phosphoenolpyruvate carboxykinase reaction to the model considerably improves the fit of the measured and the simulated NMR data. The results obtained using the extended version of the nonoxidative pentose phosphate pathway model show that the transketolase and transaldolase reactions need not be assumed reversible to get a good fit of the (13)C-labeling data. Strict statistical testing of the outcomes of (13)C-labeling based flux analysis using realistic measurement errors is demonstrated to be of prime importance for verifying the assumed metabolic model.     

A priori analysis of metabolic flux identifiability from (13)C-labeling data

The (13)C-labeling technique was introduced in the field of metabolic engineering as a tool for determining fluxes that could not be found using the 'classical' method of flux balancing. An a priori flux identifiability analysis is required in order to determine whether a (13)C-labeling experiment allows the identification of all the fluxes. In this article, we propose a method for identifiability analysis that is based on the recently introduced 'cumomer' concept. The method improves upon previous identifiability methods in that it provides a way of systematically reducing the metabolic network on the basis of structural elements that constitute a network and to use the implicit function theorem to analytically determine whether the fluxes in the reduced network are theoretically identifiable for various types of real measurement data. Application of the method to a realistic flux identification problem shows both the potential of the method in yielding new, interesting conclusions regarding the identifiability and its practical limitations that are caused by the fact that symbolic calculations grow fast with the dimension of the studied system.     

Metabolic flux analysis in Escherichia coli by integrating isotopic dynamic and isotopic stationary 13C labeling data

The novel concept of isotopic dynamic 13C metabolic flux analysis (ID-13C MFA) enables integrated analysis of isotopomer data from isotopic transient and/or isotopic stationary phase of a 13C labeling experiment, short-time experiments, and an extended range of applications of 13C MFA. In the presented work, an experimental and computational framework consisting of short-time 13C labeling, an integrated rapid sampling procedure, a LC-MS analytical method, numerical integration of the system of isotopomer differential equations, and estimation of metabolic fluxes was developed and applied to determine intracellular fluxes in glycolysis, pentose phosphate pathway (PPP), and citric acid cycle (TCA) in Escherichia coli grown in aerobic, glucose-limited chemostat culture at a dilution rate of D = 0.10 h(-1). Intracellular steady state concentrations were quantified for 12 metabolic intermediates. A total of 90 LC-MS mass isotopomers were quantified at sampling times t = 0, 91, 226, 346, 589 s and at isotopic stationary conditions. Isotopic stationarity was reached within 10 min in glycolytic and PPP metabolites. Consistent flux solutions were obtained by ID-13C MFA using isotopic dynamic and isotopic stationary 13C labeling data and by isotopic stationary 13C MFA (IS-13C MFA) using solely isotopic stationary data. It is demonstrated that integration of dynamic 13C labeling data increases the sensitivity of flux estimation, particularly at the glucose-6-phosphate branch point. The identified split ratio between glycolysis and PPP was 55%:44%. These results were confirmed by IS-13C MFA additionally using labeling data in proteinogenic amino acids (GC-MS) obtained after 5 h from sampled biomass.     

13C代谢流量分析发展30年

13C代谢流量分析(13C metabolic flux analysis,13C-MFA),是通过标记实验分析蛋白氨基酸或胞内代谢物同位素标记异构体的分布情况,从而准确定量胞内反应速率.该技术在系统理解细胞代谢特性、指导代谢工程改造和揭示病理生理学等方面起着重要作用,引起研究者的广泛重视.文中重点综述了代谢流分析30...     

Experimental design principles for isotopically instationary 13C labeling experiments

13C metabolic flux analysis (MFA) is a well-established tool in Metabolic Engineering that found numerous applications in recent years. However, one strong limitation of the current method is the requirement of an-at least approximate-isotopic stationary state at sampling time. This requirement leads to a principle lower limit for the duration of a 13C labeling experiment. A new methodological development is based on repeated sampling during the instationary transient of the 13C labeling dynamics. The statistical and computational treatment of such instationary experiments is a completely new terrain. The computational effort is very high because large differential equations have to be solved and, moreover, the intracellular pool sizes play a significant role. For this reason, the present contribution works out principles and strategies for the experimental design of instationary experiments based on a simple example network. Hereby, the potential of isotopically instationary experiments is investigated in detail. Various statistical results on instationary flux identifiability are presented and possible pitfalls of experimental design are discussed. Finally, a framework for almost optimal experimental design of isotopically instationary experiments is proposed which provides a practical guideline for the analysis of large-scale networks.     

Metabolic flux analysis of CHO cells in perfusion culture by metabolite balancing and 2D [13C, 1H] COSY NMR spectroscopy

The physiological state of CHO cells in perfusion culture was quantified by determining fluxes through the bioreaction network using ¹³C glucose and 2D-NMR spectroscopy. CHO cells were cultivated in a 2.5 L perfusion bioreactor with glucose and glutamine as the primary carbon and energy sources. The reactor was inoculated at a cell density of 8×10⁶ cells/mL and operated at ～10×10⁶ cells/mL using unlabeled glucose for the first 13 days. The second phase lasted 12 days and the medium consisted of 10% [U-¹³C]glucose, 40% labeled [1-¹³C]glucose with the balance unlabeled. After the culture attained isotopic steady state, biomass samples from the last 3 days of cultivation were considered representative and used for flux estimation. They were hydrolyzed and analyzed by 2D [¹³C, ¹H] COSY measurements using the heteronuclear single quantum correlation sequence with gradients for artifacts suppression. Metabolic fluxes were determined using the 13C-Flux software package by minimizing the residuals between the experimental and the simulated NMR data. Normalized residuals exhibited a Gaussian distribution indicating good model fit to experimental data. The glucose consumption rate was 5-fold higher than that of glutamine with 41% of glucose channeled through the pentose phosphate pathway. The fluxes at the pyruvate branch point were almost equally distributed between lactate and the TCA cycle (55% and 45%, respectively). The anaplerotic conversion of pyruvate to oxaloacetate by pyruvate carboxylase accounted for 10% of the pyruvate flux with the remaining 90% entering the TCA cycle through acetyl-CoA. The conversion of malate to pyruvate catalyzed by the malic enzyme was 70% higher than that for the anaplerotic reaction catalyzed by pyruvate carboxylase. Most amino acid catabolic and biosynthetic fluxes were significantly lower than the glycolytic and TCA cycle fluxes. Metabolic flux data from NMR analysis validated a simplified model where metabolite balancing was used for flux estimation. In this reduced flux space, estimates from these two methods were in good agreement. This simplified model can routinely be used in bioprocess development experiments to estimate metabolic fluxes with much reduced analytical investment. The high resolution flux information from 2D-NMR spectroscopy coupled with the capability to validate a simplified metabolite balancing based model for routine use make ¹³C-isotopomer analysis an attractive bioprocess development tool for mammalian cell cultures.     

New tools for mass isotopomer data evaluation in (13)C flux analysis: mass isotope correction, data consistency checking, and precursor relationships

13C metabolic flux analysis (MFA) is based on carbon-labeling experiments where a specifically (13)C labeled substrate is fed. The labeled carbon atoms distribute over the metabolic network and the label enrichment of certain metabolic pools is measured by using different methods. Recently, MS methods have been dramatically improved-large and precise datasets are now available. MS data has to be preprocessed and corrected for natural stable mass isotopes. In this article we present (1). a new elegant method to correct MS measurement data for natural stable mass isotopes by infinite dimensional matrix calculus and (2). we statistically analyze and discuss a reconstruction of labeling pattern in metabolic precursors from biosynthesis molecules. Moreover, we establish a new method for consistency checking of MS spectra that can be applied for automatic error recognition in high-throughput flux analysis procedures. Preprocessing the measurement data changes their statistical properties which have to be considered in the subsequent parameter fitting process for (13)C MFA. We show that correcting for stable mass isotopes leads to rather small correlations. On the other hand, a direct reconstruction of a precursor labeling pattern from an aromatic amino acid measurement turns out to be critical. Reasonable results are only obtained if additional, independent information about the labeling of at least one precursor is available. A versatile MatLab tool for the rapid correction and consistency checking of MS spectra is presented. Practical examples for the described methods are also given.     

Identification of metabolic fluxes in hepatic cells from transient 13C-labeling experiments: Part II. Flux estimation

This contribution addresses the identification of metabolic fluxes and metabolite concentrations in mammalian cells from transient (13)C-labeling experiments. Whilst part I describes experimental set-up and acquisition of required metabolite and (13)C-labeling data, part II focuses on setting up network models and the estimation of intracellular fluxes. Metabolic fluxes were determined in glycolysis, pentose-phosphate pathway (PPP), and citric acid cycle (TCA) in a hepatoma cell line grown in aerobic batch cultures. In glycolytic and PPP metabolite pools isotopic stationarity was observed within 30 min, whereas in the TCA cycle the labeling redistribution did not reach isotopic steady state even within 180 min. In silico labeling dynamics were in accordance with in vivo (13)C-labeling data. Split ratio between glycolysis and PPP was 57%:43%; intracellular glucose concentration was estimated at 101.6 nmol per 10(6) cells. In contrast to isotopic stationary (13)C-flux analysis, transient (13)C-flux analysis can also be applied to industrially relevant mammalian cell fed-batch and batch cultures.     

Innovations in generation and analysis of 2D [(13)C,(1)H] COSY NMR spectra for metabolic flux analysis purposes.

2D [(13)C,(1)H] COSY NMR is used by the metabolic engineering community for determining (13)C-(13)C connectivities in intracellular compounds that contain information regarding the steady-state fluxes in cellular metabolism. This paper proposes innovations in the generation and analysis of these specific NMR spectra. These include a computer tool that allows accurate determination of the relative peak areas and their complete covariance matrices even in very complex spectra. Additionally, a method is introduced for correcting the results for isotopic non-steady-state conditions. The proposed methods were applied to measured 2D [(13)C,(1)H] COSY NMR spectra. Peak intensities in a one-dimensional section of the spectrum are frequently not representative for relative peak volumes in the two-dimensional spectrum. It is shown that for some spectra a significant amount of additional information can be gained from long-range (13)C-(13)C scalar couplings in 2D [(13)C,(1)H] COSY NMR spectra. Finally, the NMR resolution enhancement by dissolving amino acid derivatives in a nonpolar solvent is demonstrated.     

Elucidating the role of copper in CHO cell energy metabolism using 13C metabolic flux analysis

¹³C‐metabolic flux analysis was used to understand copper deficiency‐related restructuring of energy metabolism, which leads to excessive lactate production in recombinant protein‐producing CHO cells. Stationary‐phase labeling experiments with U‐¹³C glucose were conducted on CHO cells grown under high and limiting copper in 3 L fed‐batch bioreactors. The resultant labeling patterns of soluble metabolites were measured by GC‐MS and used to estimate metabolic fluxes in the central carbon metabolism pathways using OpenFlux. Fluxes were evaluated 300 times from stoichiometrically feasible random guess values and their confidence intervals calculated by Monte Carlo simulations. Results from metabolic flux analysis exhibited significant carbon redistribution throughout the metabolic network in cells under Cu deficiency. Specifically, glycolytic fluxes increased (25%–79% relative to glucose uptake) whereas fluxes through the TCA and pentose phosphate pathway (PPP) were lower (15%–23% and 74%, respectively) compared with the Cu‐containing condition. Furthermore, under Cu deficiency, 33% of the flux entering TCA via the pyruvate node was redirected to lactate and malate production. Based on these results, we hypothesize that Cu deficiency disrupts the electron transport chain causing ATP deficiency, redox imbalance, and oxidative stress, which in turn drive copper‐deficient CHO cells to produce energy via aerobic glycolysis, which is associated with excessive lactate production, rather than the more efficient route of oxidative phosphorylation. © 2015 American Institute of Chemical Engineers Biotechnol. Prog., 31:1179–1186, 2015     

Photobioreactor design for isotopic non‐stationary 13C‐metabolic flux analysis (INST 13C‐MFA) under photoautotrophic conditions

Adaptive metabolic behavior of photoautotrophic microorganisms toward genetic and environmental perturbations can be interpreted in a quantitative depiction of carbon flow through a biochemical reaction network using isotopic non‐stationary ¹³C‐metabolic flux analysis (INST ¹³C‐MFA). To evaluate ¹³C‐metabolic flux maps for Chlamydomonas reinhardtii, an original experimental framework was designed allowing rapid, reliable collection of high‐quality isotopomer data against time. It involved (i) a short‐time ¹³C labeling injection device based on mixing control in a torus‐shaped photobioreactor with plug‐flow hydrodynamics allowing a sudden step‐change in the ¹³C proportion in the substrate feed and (ii) a rapid sampling procedure using an automatic fast filtration method coupled to a manual rapid liquid nitrogen quenching step. ¹³C‐substrate labeling enrichment was controlled through the total dissolved inorganic carbon concentration in the pulsed solution. First results were obtained from steady‐state continuous culture measurements allowing the characterization of the kinetics of label incorporation into light‐limited growing cells cultivated in a photobioreactor operating at the maximal biomass productivity for an incident photon flux density of 200 µmol m^?2 s^?1. ¹³C label incorporation was measured for 21 intracellular metabolites using IC‐MS/MS in 58 samples collected across a labeling experiment duration of 7 min. The fastest labeling rate was observed for 2/3‐phosphoglycerate with an apparent isotopic stationary state reached after 300 s. The labeling rate was consistent with the optimized mixing time of about 4.9 s inside the reactor and the shortest reliable sampling period assessed at 5 s. Biotechnol. Bioeng. 2012; 109: 3030–3040. © 2012 Wiley Periodicals, Inc.     

Metabolic flux analysis using 13C peptide label measurements

¹³C metabolic flux analysis (MFA) has become the experimental method of choice to investigate the cellular metabolism of microbes, cell cultures and plant seeds. Conventional steady‐state MFA utilizes isotopic labeling measurements of amino acids obtained from protein hydrolysates. To retain spatial information in conventional steady‐state MFA, tissues or subcellular fractions must be dissected or biochemically purified. In contrast, peptides retain their identity in complex protein extracts, and may therefore be associated with a specific time of expression, tissue type and subcellular compartment. To enable ‘single‐sample’ spatially and temporally resolved steady‐state flux analysis, we investigated the suitability of peptide mass distributions (PMDs) as an alternative to amino acid label measurements. PMDs are the discrete convolution of the mass distributions of the constituent amino acids of a peptide. We investigated the requirements for the unique deconvolution of PMDs into amino acid mass distributions (AAMDs), the influence of peptide sequence length on parameter sensitivity, and how AAMD and flux estimates that are determined through deconvolution compare to estimates from a conventional GC–MS measurement‐based approach. Deconvolution of PMDs of the storage protein β–conglycinin of soybean (Glycine max) resulted in good AAMD and flux estimates if fluxes were directly fitted to PMDs. Unconstrained deconvolution resulted in inferior AAMD and flux estimates. PMD measurements do not include amino acid backbone fragments, which increase the information content in GC–MS‐derived analyses. Nonetheless, the resulting flux maps were of comparable quality due to the precision of Orbitrap quantification and the larger number of peptide measurements.     

13C-metabolic flux analysis of Clostridium ljungdahlii illuminates its core metabolism under mixotrophic culture conditions

Carbon dioxide-fixing acetogenic bacteria (acetogens) utilizing the Wood-Ljungdahl Pathway (WLP) play an important role in CO₂ fixation in the biosphere and in the development of biological processes – alone or in cocultures, under both autotrophic and mixotrophic conditions – for production of chemicals and fuels. To date, limited work has been reported in experimentally validating and quantifying reaction fluxes of their core metabolic pathways. Here, the core metabolic model of the acetogen Clostridium ljungdahlii was interrogated using ¹³C-metabolic flux analysis (¹³C-MFA), which required the development of a new defined culture medium. Autotrophic, heterotrophic, and mixotrophic growth in defined medium was possible by adding 1 mM methionine to replace yeast extract. Our ¹³C-MFA found an incomplete TCA cycle and inactive core pathways/reactions, notably those of the oxidative pentose phosphate pathway, Entner-Doudoroff pathway, and malate dehydrogenase. ¹³C-MFA during mixotrophic growth using the parallel tracers [1–¹³C]fructose, [1,2–¹³C]fructose, [1,2,3–¹³C]fructose, and [U–¹³C]asparagine found that externally supplied CO₂ contributed the majority of carbon consumed. All internally-produced CO₂ from the catabolism of asparagine and fructose was consumed by the WLP. While glycolysis of fructose was active, it was not a major contributor to overall production of ATP, NADH, and acetyl-CoA. Gluconeogenic reactions were active despite the availability of organic carbon. Asparagine was catabolized equally via conversion to threonine and subsequent cleavage to produce acetaldehyde and glycine, and via deamination to fumarate and then the anaplerotic conversion of malate to pyruvate. Both pathways for asparagine catabolism produced acetyl-CoA, either directly via pyruvate or indirectly via the WLP. Cofactor stoichiometry based on our data predicted an essentially zero flux through the ferredoxin-dependent transhydrogenase (Nfn) reaction. Instead, nearly all of NADPH generated from the hydrogenase reaction was consumed by the WLP. Reduced ferredoxin produced by the hydrogenase reaction and glycolysis was mostly used for ATP generation via the RNF/ATPase system, with the remainder consumed by the WLP. NADH produced by RNF/ATPase was entirely consumed via the WLP.