Analysis of metabolic flux in <Emphasis Type="Italic">Escherichia coli</Emphasis> expressing human-like collagen in fed-batch culture

Metabolic flux distributions of recombinant Escherichia coli BL21 expressing human-like collagen were determined by means of a stoichiometric network and metabolic balancing. At the batch growth stage, the fluxes of the pentose phosphate pathway were higher than the fluxes of the fed-batch growth phase and the production stage. After the temperature was increased, there was a substantially elevated energy demand for synthesizing human-like collagen and heat-shock proteins, which resulted in changes in metabolic fluxes. The activities of the Embden-Meyerhof-Parnas pathway and the tricarboxylic acid cycle were significantly enhanced, leading to a reduction in the fluxes of the pentose phosphate pathway and other anabolic pathways. The temperature upshift also caused an increase in NADPH production by isocitrate dehydrogenase in the tricarboxylic acid cycle. The metabolic model predicted the involvement of a transhydrogenase that generates additional NADH from NADPH, thereby increasing ATP regeneration in the respiratory chain. These data indicated that the maintenance energy for cellular activity increased with the increase in biomass in fed-batch culture, and that cell growth and synthesis of human-like collagen could clearly represent the changes in metabolic fluxes. At the production stage, more NADPH was used to synthesize human-like collagen than for maintaining cellular activity, cell growth, and cell propagation. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Modeling and optimization of a multi-product biosynthesis factory for multiple objectives

Genetic algorithms and optimization in general, enable us to probe deeper into the metabolic pathway recipe for multi-product biosynthesis. An augmented model for optimizing serine and tryptophan flux ratios simultaneously in Escherichia coli, was developed by linking the dynamic tryptophan operon model and aromatic amino acid-tryptophan biosynthesis pathways to the central carbon metabolism model. Six new kinetic parameters of the augmented model were estimated with considerations of available experimental data and other published works. Major differences between calculated and reference concentrations and fluxes were explained. Sensitivities and underlying competition among fluxes for carbon sources were consistent with intuitive expectations based on metabolic network and previous results. Biosynthesis rates of serine and tryptophan were simultaneously maximized using the augmented model via concurrent gene knockout and manipulation. The optimization results were obtained using the elitist non-dominant sorting genetic algorithm (NSGA-II) supported by pattern recognition heuristics. A range of Pareto-optimal enzyme activities regulating the amino acids biosynthesis was successfully obtained and elucidated wherever possible vis-à-vis fermentation work based on recombinant DNA technology. The predicted potential improvements in various metabolic pathway recipes using the multi-objective optimization strategy were highlighted and discussed in detail.     

Metabolic flux analysis for serine alkaline protease fermentation by Bacillus licheniformis in a defined medium: effects of the oxygen transfer rate.

The metabolic fluxes through the central carbon pathways in the bioprocess for serine alkaline protease (SAP) production by Bacillus licheniformis were calculated by the metabolic flux-based stoichiometric model based on the proposed metabolic network that contains 102 metabolites and 133 reaction fluxes using the time profiles of citrate, dry cell, organic acids, amino acids, and SAP as the constraints. The model was solved by minimizing the SAP accumulation rate in the cell. The effects of the oxygen-transfer rate (OTR) on the metabolic fluxes were investigated in a defined medium where citrate was used as the sole carbon source. The central pathways were active for the growth and the SAP synthesis in all the periods of the bioprocess at low (LOT), medium (MOT), and high (HOT) oxygen-transfer conditions. The flux partitioning in the TCA cycle at alpha-ketoglutarate towards glutamate group and at oxalacetate (OA) toward aspartic acid group amino acids were dependent on the OTR. The flux of the anaplerotic reaction that connects the TCA cycle either from malate or OA to the gluconeogenesis pathway via the main branch point pyruvate (Pyr) was also influenced by the OTR. With the decrease in the OTR, the intracellular flux values after glycerate 3-phosphate (PG3) in the gluconeogenesis pathway and the specific growth rate decreased. The total ATP-generation rate increased with the increase in OTR. The pathway towards the aspartic acid family amino acids which is important for sporulation that precedes the SAP synthesis were all active throughout the bioprocess. Metabolic flux analysis results at LOT, MOT, and HOT conditions encourage the design of an oxygen-transfer strategy in the bioreactor; moreover, asparagine synthetase or aspartate kinase could be the potential metabolic engineering sites due to the low value of the flux from the branch point aspartate toward asparagine.     

Citrate production from n-alkane by Candida lipolytica in reference to carbon fluxes in vivo

Summary Citrate production by Candida lipolytica from n-alkane as sole carbon source was studied with a NH ₄ ⁺ -limited chemostat culture. This organism produced citrate and isocitrate simultaneously, the ratio of these acids remaining nearly unchanged within the range of dilution rates used. However, the specific production rates of these acids were maximal around a particular dilution rate. For each steady-state culture, specific activities of some key enzymes in tricarboxylate and glyoxylate cycles were measured in vitro. In order to correlate the enzyme activites with the change in specific production rate of each acid, the intracellular carbon fluxes were assessed by solving algebraic equations derived from carbon balance with respect to a simplified metabolic map of n-alkane metabolism. By comparing the fluxes assessed in vivo with the enzyme activities measured in vitro, the metabolic pathway(s) most likely to control the production of acids and cell material were investigated.     

Shewanella oneidensis MR-1 fluxome under various oxygen conditions

The central metabolic fluxes of Shewanella oneidensis MR-1 were examined under carbon-limited (aerobic) and oxygen-limited (microaerobic) chemostat conditions, using 13C-labeled lactate as the sole carbon source. The carbon labeling patterns of key amino acids in biomass were probed using both gas chromatography-mass spectrometry (GC-MS) and 13C nuclear magnetic resonance (NMR). Based on the genome annotation, a metabolic pathway model was constructed to quantify the central metabolic flux distributions. The model showed that the tricarboxylic acid (TCA) cycle is the major carbon metabolism route under both conditions. The Entner-Doudoroff and pentose phosphate pathways were utilized primarily for biomass synthesis (with a flux below 5% of the lactate uptake rate). The anaplerotic reactions (pyruvate to malate and oxaloacetate to phosphoenolpyruvate) and the glyoxylate shunt were active. Under carbon-limited conditions, a substantial amount (9% of the lactate uptake rate) of carbon entered the highly reversible serine metabolic pathway. Under microaerobic conditions, fluxes through the TCA cycle decreased and acetate production increased compared to what was found for carbon-limited conditions, and the flux from glyoxylate to glycine (serine-glyoxylate aminotransferase) became measurable. Although the flux distributions under aerobic, microaerobic, and shake flask culture conditions were different, the relative flux ratios for some central metabolic reactions did not differ significantly (in particular, between the shake flask and aerobic-chemostat groups). Hence, the central metabolism of S. oneidensis appears to be robust to environmental changes. Our study also demonstrates the merit of coupling GC-MS with 13C NMR for metabolic flux analysis to reduce the use of 13C-labeled substrates and to obtain more-accurate flux values.     

Optimization-driven identification of genetic perturbations accelerates the convergence of model parameters in ensemble modeling of metabolic networks

The ensemble modeling (EM) approach has shown promise in capturing kinetic and regulatory effects in the modeling of metabolic networks. Efficacy of the EM procedure relies on the identification of model parameterizations that adequately describe all observed metabolic phenotypes upon perturbation. In this study, we propose an optimization-based algorithm for the systematic identification of genetic/enzyme perturbations to maximally reduce the number of models retained in the ensemble after each round of model screening. The key premise here is to design perturbations that will maximally scatter the predicted steady-state fluxes over the ensemble parameterizations. We demonstrate the applicability of this procedure for an Escherichia coli metabolic model of central metabolism by successively identifying single, double, and triple enzyme perturbations that cause the maximum degree of flux separation between models in the ensemble. Results revealed that optimal perturbations are not always located close to reaction(s) whose fluxes are measured, especially when multiple perturbations are considered. In addition, there appears to be a maximum number of simultaneous perturbations beyond which no appreciable increase in the divergence of flux predictions is achieved. Overall, this study provides a systematic way of optimally designing genetic perturbations for populating the ensemble of models with relevant model parameterizations.     

Analysis of flux estimates based on (13)C-labelling experiments.

Modelling of the fluxes in central metabolism can be performed by combining labelling experiments with metabolite balancing. Using this approach, multiple samples from a cultivation of Saccharomyces cerevisiae in metabolic and isotopic steady state were analysed, and the metabolic fluxes in central metabolism were estimated. In the various samples, the estimates of the central metabolic pathways, the tricarboxylic acid cycle, the oxidative pentose phosphate pathway and the anaplerotic pathway, showed an unprecedented reproducibility. The high reproducibility was obtained with fractional labellings of individual carbon atoms as the calculational base, illustrating that the more complex modelling using isotopomers is not necessarily superior with respect to reproducibility of the flux estimates. Based on these results some general difficulties in flux estimation are discussed.     

Metabolic carbon fluxes and biosynthesis of polyhydroxyalkanoates in Ralstonia eutropha on short chain fatty acids

Short chain fatty acids such as acetic, propionic, and butyric acids can be synthesized into polyhydroxyalkanoates (PHAs) by Ralstonia eutropha. Metabolic carbon fluxes of the acids in living cells have significant effect on the yield, composition, and thermomechanical properties of PHA bioplastics. Based on the general knowledge of central metabolism pathways and the unusual metabolic pathways in R. eutropha, a metabolic network of 41 bioreactions is constructed to analyze the carbon fluxes on utilization of the short chain fatty acids. In fed-batch cultures with constant feeding of acid media, carbon metabolism and distribution in R. eutropha were measured involving CO2, PHA biopolymers, and residual cell mass. As the cells underwent unsteady state metabolism and PHA biosynthesis under nitrogen-limited conditions, accumulative carbon balance was applied for pseudo-steady-state analysis of the metabolic carbon fluxes. Cofactor NADP/NADPH balanced between PHA synthesis and the C3/C4 pathway provided an independent constraint for solution of the underdetermined metabolic network. A major portion of propionyl-CoA was directed to pyruvate via the 2-methylcitrate cycle and further decarboxylated to acetyl-CoA. Only a small amount of propionate carbon (<15% carbon) was directly condensed with acetyl-CoA for 3-hydroxyvalerate. The ratio of glyoxylate shunt to TCA cycle varies from 0 to 0.25, depending on the intracellular acetyl-CoA level and acetic acid in the medium. Malate is the node of the C3/C4 pathway and TCA cycle and its decarboxylation to dehydrogenation ranges from 0.33 to 1.28 in response to the demands on NADPH and oxaloacetate for short chain fatty acids utilization.     

Metabolic and morphological differences between rapidly proliferating cancerous and normal breast epithelial cells

The metabolic and morphological characteristics of two human epithelial breast cell populations--MCF7 cells, a cancerous cell line, and 48R human mammary epithelial cells (48R HMECs), a noncancerous, finite lifespan cell strain--were compared at identical growth rates. Both cell types were induced to grow rapidly in nutrient-rich media containing 13C-labeled glucose, and the isotopic enrichment of cellular metabolites was quantified to calculate metabolic fluxes in key pathways. Despite their similar growth rates, the cells exhibited distinctly different metabolic and morphological profiles. MCF7 cells have an 80% smaller exposed surface area and contain 26% less protein per cell than the 48R cells. Surprisingly, rapidly proliferating 48R cells exhibited a 225% higher per-cell glucose consumption rate, a 250% higher per-cell lactate production rate, and a nearly identical per-cell glutamine consumption rate relative to the cancer cell line. However, when fluxes were considered on the basis of exposed area, the cancer cells were observed to have higher glucose, lactate, and glutamine fluxes, demonstrating superior transport capabilities per unit area of cell membrane. MCF7 cells also consumed amino acids at rates much higher than are generally required for protein synthesis, whereas 48R cells generally did not. Pentose phosphate pathway activity was higher in MCF7 cells, and the flux of glutamine to glutamate was less reversible. Energy efficiency was significantly higher in MCF7 cells, as a result of a combination of their smaller size and greater reliance on the TCA cycle than the 48R cells. These observations support evolutionary models of cancer cell metabolism and suggest targets for metabolic drugs in metastatic breast cancers.     

Integrating systemic and molecular levels to infer key drivers sustaining metabolic adaptations

Metabolic adaptations to complex perturbations, like the response to pharmacological treatments in multifactorial diseases such as cancer, can be described through measurements of part of the fluxes and concentrations at the systemic level and individual transporter and enzyme activities at the molecular level. In the framework of Metabolic Control Analysis (MCA), ensembles of linear constraints can be built integrating these measurements at both systemic and molecular levels, which are expressed as relative differences or changes produced in the metabolic adaptation. Here, combining MCA with Linear Programming, an efficient computational strategy is developed to infer additional non-measured changes at the molecular level that are required to satisfy these constraints. An application of this strategy is illustrated by using a set of fluxes, concentrations, and differentially expressed genes that characterize the response to cyclin-dependent kinases 4 and 6 inhibition in colon cancer cells. Decreases and increases in transporter and enzyme individual activities required to reprogram the measured changes in fluxes and concentrations are compared with down-regulated and up-regulated metabolic genes to unveil those that are key molecular drivers of the metabolic response.     

Metabolomics Analysis of Metabolic Effects of Nicotinamide Phosphoribosyltransferase (NAMPT) Inhibition on Human Cancer Cells

Nicotinamide phosphoribosyltransferase (NAMPT) plays an important role in cellular bioenergetics. It is responsible for converting nicotinamide to nicotinamide adenine dinucleotide, an essential molecule in cellular metabolism. NAMPT has been extensively studied over the past decade due to its role as a key regulator of nicotinamide adenine dinucleotide–consuming enzymes. NAMPT is also known as a potential target for therapeutic intervention due to its involvement in disease. In the current study, we used a global mass spectrometry–based metabolomic approach to investigate the effects of FK866, a small molecule inhibitor of NAMPT currently in clinical trials, on metabolic perturbations in human cancer cells. We treated A2780 (ovarian cancer) and HCT-116 (colorectal cancer) cell lines with FK866 in the presence and absence of nicotinic acid. Significant changes were observed in the amino acids metabolism and the purine and pyrimidine metabolism. We also observed metabolic alterations in glycolysis, the citric acid cycle (TCA), and the pentose phosphate pathway. To expand the range of the detected polar metabolites and improve data confidence, we applied a global metabolomics profiling platform by using both non-targeted and targeted hydrophilic (HILIC)-LC-MS and GC-MS analysis. We used Ingenuity Knowledge Base to facilitate the projection of metabolomics data onto metabolic pathways. Several metabolic pathways showed differential responses to FK866 based on several matches to the list of annotated metabolites. This study suggests that global metabolomics can be a useful tool in pharmacological studies of the mechanism of action of drugs at a cellular level.     

Batch culture characterization and metabolic flux analysis of succinate-producing Escherichia coli strains

This study presents an in-depth analysis of the anaerobic metabolic fluxes of various mutant strains of Escherichia coli overexpressing the Lactococcus lactis pyruvate carboxylase (PYC) for the production of succinate. Previously, a metabolic network design that includes an active glyoxylate pathway implemented in vivo increased succinate yield from glucose in an E. coli mutant to 1.6 mol/mol under fully anaerobic conditions. The design consists of a dual succinate synthesis route, which diverts required quantities of NADH through the traditional fermentative pathway and maximizes the carbon converted to succinate by balancing the carbon flux through the fermentative pathway and the glyoxylate pathway (which has a lower NADH requirement). Mutant strains previously constructed during the development of high-yield succinate-producing strains were selected for further characterization to understand their metabolic response as a result of several genetic manipulations and to determine the significance of the fermentative and the glyoxylate pathways in the production of succinate. Measured fluxes obtained under batch cultivation conditions were used to estimate intracellular fluxes and identify critical branch point flux split ratios. The comparison of changes in branch point flux split ratios to the glyoxylate pathway and the fermentative pathway at the oxaloacetate (OAA) node as a result of different mutations revealed the sensitivity of succinate yield to these manipulations. The most favorable split ratio to obtain the highest succinate yield was the fractional partition of OAA to glyoxylate of 0.32 and 0.68 to the fermentative pathway obtained in strains SBS550MG (pHL413) and SBS990MG (pHL413). The succinate yields achieved in these two strains were 1.6 and 1.7 mol/mol, respectively. In addition, an active glyoxylate pathway in an ldhA, adhE, ack-pta mutant strain is shown to be responsible for the high succinate yields achieved anaerobically. Furthermore, in vitro activity measurements of seven crucial enzymes involved in the pathways studied and intracellular measurements of key intermediate metabolite pools provided additional insights on the physiological perturbations caused by these mutations. The characterization of these recombinant mutant strains in terms of flux distribution pattern, in vitro enzyme activity and intracellular metabolite pools provides useful information for the rational modification of metabolic fluxes to improve succinate production.     

Parameter estimation and dynamic control analysis of central carbon metabolism in <Emphasis Type="Italic">Escherichia coli</Emphasis>

The central carbon metabolic pathway is the most important among metabolic pathways in all microorganisms since it produces energy and precursors for biosynthesis. In this study, a dynamic model for central carbon metabolism in Escherichia coli (E. coli) consisting of the phosphotransferase (PTS) system, glycolysis, pentose-phosphate pathway (PPP), and storage materials was obtained by ameliorating the model proposed by Chassagnole et al. (2002). In order to improve the performance of the model, principal parameters were estimated through the experimental measurements of intracellular concentrations of metabolites under transient conditions. Through dynamic metabolic control analysis (MCA), the tendencies of the metabolic fluxes at branch points were investigated, and the key parameters and enzyme kinetics that most dominantly affected the productivity of the desired metabolites were determined.     

Metabolic profiling by 13C-NMR spectroscopy: [1,2-13C2]glucose reveals a heterogeneous metabolism in human leukemia T cells

Metabolic profiling is defined as the simultaneous assessment of substrate fluxes within and among the different pathways of metabolite synthesis and energy production under various physiological conditions. The use of stable-isotope tracers and the analysis of the distribution of labeled carbons in various intermediates, by both mass spectrometry and NMR spectroscopy, allow the role of several metabolic processes in cell growth and death to be defined. In the present paper we describe the metabolic profiling of Jurkat cells by isotopomer analysis using (13)C-NMR spectroscopy and [1,2-(13)C(2)]glucose as the stable-isotope tracer. The isotopomer analysis of the lactate, alanine, glutamate, proline, serine, glycine, malate and ribose-5-phosphate moiety of nucleotides has allowed original integrated information regarding the pentose phosphate pathway, TCA cycle, and amino acid metabolism in proliferating human leukemia T cells to be obtained. In particular, the contribution of the glucose-6-phosphate dehydrogenase and transketolase activities to phosphoribosyl-pyrophosphate synthesis was evaluated directly by the determination of isotopomers of the [1'-(13)C], [4',5'-(13)C(2)]ribosyl moiety of nucleotides. Furthermore, the relative contribution of the glycolysis and pentose cycle to lactate production was estimated via analysis of lactate isotopomers. Interestingly, pyruvate carboxylase and pyruvate dehydrogenase flux ratios measured by glutamate isotopomers and the production of isotopomers of several metabolites showed that the metabolic processes described could not take place simultaneously in the same macrocompartments (cells). Results revealed a heterogeneous metabolism in an asynchronous cell population that may be interpreted on the basis of different metabolic phenotypes of subpopulations in relation to different cell cycle phases.     

Shewanella oneidensis MR-1 Fluxome under Various Oxygen Conditions

The central metabolic fluxes of Shewanella oneidensis MR-1 were examined under carbon-limited (aerobic) and oxygen-limited (microaerobic) chemostat conditions, using ¹³C-labeled lactate as the sole carbon source. The carbon labeling patterns of key amino acids in biomass were probed using both gas chromatography-mass spectrometry (GC-MS) and ¹³C nuclear magnetic resonance (NMR). Based on the genome annotation, a metabolic pathway model was constructed to quantify the central metabolic flux distributions. The model showed that the tricarboxylic acid (TCA) cycle is the major carbon metabolism route under both conditions. The Entner-Doudoroff and pentose phosphate pathways were utilized primarily for biomass synthesis (with a flux below 5% of the lactate uptake rate). The anaplerotic reactions (pyruvate to malate and oxaloacetate to phosphoenolpyruvate) and the glyoxylate shunt were active. Under carbon-limited conditions, a substantial amount (9% of the lactate uptake rate) of carbon entered the highly reversible serine metabolic pathway. Under microaerobic conditions, fluxes through the TCA cycle decreased and acetate production increased compared to what was found for carbon-limited conditions, and the flux from glyoxylate to glycine (serine-glyoxylate aminotransferase) became measurable. Although the flux distributions under aerobic, microaerobic, and shake flask culture conditions were different, the relative flux ratios for some central metabolic reactions did not differ significantly (in particular, between the shake flask and aerobic-chemostat groups). Hence, the central metabolism of S. oneidensis appears to be robust to environmental changes. Our study also demonstrates the merit of coupling GC-MS with ¹³C NMR for metabolic flux analysis to reduce the use of ¹³C-labeled substrates and to obtain more-accurate flux values.     

Metabolic flux analysis of Arthrobacter sp. CGMCC 3584 for cAMP production based on C tracer experiments and gas chromatography–mass spectrometry

Arthrobacter sp. CGMCC 3584 are able to produce cAMP from glucose by the purine synthesis pathway via de novo or salvage biosynthesis. In order to gain an improved understanding of its metabolism, ¹³C-labeling experiment and gas chromatography–mass spectrometry (GC–MS) analysis were employed to determine the metabolic network structure and estimate the intracellular fluxes. GC–MS analysis helps to reflect the activity of the intracellular pathways and reactions. The metabolic network mainly contains glycolytic and pentose phosphate pathways, the tricarboxylic acid cycle, and the inactive glyoxylate shunt. Hypoxanthine as a precursor of cAMP and sodium fluoride as an inhibitor of glycolysis were found to increase the cAMP production, as well as the flux through the PP pathway. The effects of adding hypoxanthine and sodium fluoride are discussed based on the enzyme assays and metabolic flux analysis. In conclusion, our results provide quantitative insights into how cells manipulate the metabolic network under different culture conditions and this may be of value in metabolic regulation for desirable production.     

Identifying causal genes and dysregulated pathways in complex diseases

In complex diseases, various combinations of genomic perturbations often lead to the same phenotype. On a molecular level, combinations of genomic perturbations are assumed to dys-regulate the same cellular pathways. Such a pathway-centric perspective is fundamental to understanding the mechanisms of complex diseases and the identification of potential drug targets. In order to provide an integrated perspective on complex disease mechanisms, we developed a novel computational method to simultaneously identify causal genes and dys-regulated pathways. First, we identified a representative set of genes that are differentially expressed in cancer compared to non-tumor control cases. Assuming that disease-associated gene expression changes are caused by genomic alterations, we determined potential paths from such genomic causes to target genes through a network of molecular interactions. Applying our method to sets of genomic alterations and gene expression profiles of 158 Glioblastoma multiforme (GBM) patients we uncovered candidate causal genes and causal paths that are potentially responsible for the altered expression of disease genes. We discovered a set of putative causal genes that potentially play a role in the disease. Combining an expression Quantitative Trait Loci (eQTL) analysis with pathway information, our approach allowed us not only to identify potential causal genes but also to find intermediate nodes and pathways mediating the information flow between causal and target genes. Our results indicate that different genomic perturbations indeed dys-regulate the same functional pathways, supporting a pathway-centric perspective of cancer. While copy number alterations and gene expression data of glioblastoma patients provided opportunities to test our approach, our method can be applied to any disease system where genetic variations play a fundamental causal role.