Genome-scale model for Clostridium acetobutylicum: Part I. Metabolic network resolution and analysis

A genome-scale metabolic network reconstruction for Clostridium acetobutylicum (ATCC 824) was carried out using a new semi-automated reverse engineering algorithm. The network consists of 422 intracellular metabolites involved in 552 reactions and includes 80 membrane transport reactions. The metabolic network illustrates the reliance of clostridia on the urea cycle, intracellular L-glutamate solute pools, and the acetylornithine transaminase for amino acid biosynthesis from the 2-oxoglutarate precursor. The semi-automated reverse engineering algorithm identified discrepancies in reaction network databases that are major obstacles for fully automated network-building algorithms. The proposed semi-automated approach allowed for the conservation of unique clostridial metabolic pathways, such as an incomplete TCA cycle. A thermodynamic analysis was used to determine the physiological conditions under which proposed pathways (e.g., reverse partial TCA cycle and reverse arginine biosynthesis pathway) are feasible. The reconstructed metabolic network was used to create a genome-scale model that correctly characterized the butyrate kinase knock-out and the asolventogenic M5 pSOL1 megaplasmid degenerate strains. Systematic gene knock-out simulations were performed to identify a set of genes encoding clostridial enzymes essential for growth in silico.     

Development of a kinetic metabolic model: application to Catharanthus roseus hairy root

A kinetic metabolic model describing Catharanthus roseus hairy root growth and nutrition was developed. The metabolic network includes glycolysis, pentose-phosphate pathway, TCA cycle and the catabolic reactions leading to cell building blocks such as amino acids, organic acids, organic phosphates, lipids and structural hexoses. The central primary metabolic network was taken at pseudo-steady state and metabolic flux analysis technique allowed reducing from 31 metabolic fluxes to 20 independent pathways. Hairy root specific growth rate was described as a function of intracellular concentration in cell building blocks. Intracellular transport and accumulation kinetics for major nutrients were included. The model uses intracellular nutrients as well as energy shuttles to describe metabolic regulation. Model calibration was performed using experimental data obtained from batch and medium exchange liquid cultures of C. roseus hairy root using a minimal medium in Petri dish. The model is efficient in estimating the growth rate.     

Phosphoketolase Pathway for Xylose Catabolism in Clostridium acetobutylicum Revealed by 13C Metabolic Flux Analysis

Solvent-producing clostridia are capable of utilizing pentose sugars, including xylose and arabinose; however, little is known about how pentose sugars are catabolized through the metabolic pathways in clostridia. In this study, we identified the xylose catabolic pathways and quantified their fluxes in Clostridium acetobutylicum based on [1-(13)C]xylose labeling experiments. The phosphoketolase pathway was found to be active, which contributed up to 40% of the xylose catabolic flux in C. acetobutylicum. The split ratio of the phosphoketolase pathway to the pentose phosphate pathway was markedly increased when the xylose concentration in the culture medium was increased from 10 to 20 g liter(-1). To our knowledge, this is the first time that the in vivo activity of the phosphoketolase pathway in clostridia has been revealed. A phosphoketolase from C. acetobutylicum was purified and characterized, and its activity with xylulose-5-P was verified. The phosphoketolase was overexpressed in C. acetobutylicum, which resulted in slightly increased xylose consumption rates during the exponential growth phase and a high level of acetate accumulation.     

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.     

Genetically constrained metabolic flux analysis

Significant progress has been made in using existing metabolic databases to estimate metabolic fluxes. Traditional metabolic flux analysis generally starts with a predetermined metabolic network. This approach has been employed successfully to analyze the behaviors of recombinant strains by manually adding or removing the corresponding pathway(s) in the metabolic map. The current work focuses on the development of a new framework that utilizes genomic and metabolic databases, including available genetic/regulatory network structures and gene chip expression data, to constrain metabolic flux analysis. The genetic network consisting of the sensing/regulatory circuits will activate or deactivate a specific set of genes in response to external stimulus. The activation and/or repression of this set of genes will result in different gene expression levels that will in turn change the structure of the metabolic map. Hence, the metabolic map will automatically "adapt" to the external stimulus as captured by the genetic network. This adaptation selects a subnetwork from the pool of feasible reactions and so performs what we term "environmentally driven dimensional reduction." The Escherichia coli oxygen and redox sensing/regulatory system, which controls the metabolic patterns connected to glycolysis and the TCA cycle, was used as a model system to illustrate the proposed approach.     

The ethylmalonyl-CoA pathway is used in place of the glyoxylate cycle by Methylobacterium extorquens AM1 during growth on acetate

Acetyl-CoA assimilation was extensively studied in organisms harboring the glyoxylate cycle. In this study, we analyzed the metabolism of the facultative methylotroph Methylobacterium extorquens AM1, which lacks isocitrate lyase, the key enzyme in the glyoxylate cycle, during growth on acetate. MS/MS-based proteomic analysis revealed that the protein repertoire of M. extorquens AM1 grown on acetate is similar to that of cells grown on methanol and includes enzymes of the ethylmalonyl-CoA (EMC) pathway that were recently shown to operate during growth on methanol. Dynamic 13C labeling experiments indicate the presence of distinct entry points for acetate: the EMC pathway and the TCA cycle. 13C steady-state metabolic flux analysis showed that oxidation of acetyl-CoA occurs predominantly via the TCA cycle and that assimilation occurs via the EMC pathway. Furthermore, acetyl-CoA condenses with the EMC pathway product glyoxylate, resulting in malate formation. The latter, also formed by the TCA cycle, is converted to phosphoglycerate by a reaction sequence that is reversed with respect to the serine cycle. Thus, the results obtained in this study reveal the utilization of common pathways during the growth of M. extorquens AM1 on C1 and C2 compounds, but with a major redirection of flux within the central metabolism. Furthermore, our results indicate that the metabolic flux distribution is highly complex in this model methylotroph during growth on acetate and is fundamentally different from organisms using the glyoxylate cycle.     

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

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

Intracellular pyruvate levels positively correlate with cytokine production capacity in tolerant monocytes from patients with pneumonia

BackgroundCommunity-acquired pneumonia (CAP) is responsible for a high morbidity and mortality worldwide. Monocytes are essential for pathogen recognition and the initiation of an innate immune response. Immune cells induce intracellular glycolysis upon activation to support several functions.ObjectiveTo obtain insight in the metabolic profile of blood monocytes during CAP, with a focus on glycolysis and branching metabolic pathways, and to determine a possible association between intracellular metabolite levels and monocyte function.MethodsMonocytes were isolated from blood of patients with CAP within 24 h of hospital admission and from control subjects matched for age, sex and chronic comorbidities. Changes in glycolysis, oxidative phosphorylation (OXPHOS), tricarboxylic acid (TCA) cycle and the pentose phosphate pathway were investigated through RNA sequencing and metabolomics measurements. Monocytes were stimulated ex vivo with lipopolysaccharide (LPS) to determine their capacity to produce tumor necrosis factor (TNF), interleukin (IL)-1β and IL-10.Results50 patients with CAP and 25 non-infectious control subjects were studied. When compared with control monocytes, monocytes from patients showed upregulation of many genes involved in glycolysis, including PKM, the gene encoding pyruvate kinase, the rate limiting enzyme for pyruvate production. Gene set enrichment analysis of OXPHOS, the TCA cycle and the pentose phosphate pathway did not reveal differences between monocytes from patients and controls. Patients' monocytes had elevated intracellular levels of pyruvate and the TCA cycle intermediate α-ketoglutarate. Monocytes from patients were less capable of producing cytokines upon LPS stimulation. Intracellular pyruvate (but not α-ketoglutarate) concentrations positively correlated with IL-1β and IL-10 levels released by patients' (but not control) monocytes upon exposure to LPS.ConclusionThese results suggest that elevated intracellular pyruvate levels may partially maintain cytokine production capacity of hyporesponsive monocytes from patients with CAP.     

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 ¹³C 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.     

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.     

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.     

Role of Central Metabolism in the Osmoadaptation of the Halophilic Bacterium Chromohalobacter salexigens

Bacterial osmoadaptation involves the cytoplasmic accumulation of compatible solutes to counteract extracellular osmolarity. The halophilic and highly halotolerant bacterium Chromohalobacter salexigens is able to grow up to 3 m NaCl in a minimal medium due to the de novo synthesis of ectoines. This is an osmoregulated pathway that burdens central metabolic routes by quantitatively drawing off TCA cycle intermediaries. Consequently, metabolism in C. salexigens has adapted to support this biosynthetic route. Metabolism of C. salexigens is more efficient at high salinity than at low salinity, as reflected by lower glucose consumption, lower metabolite overflow, and higher biomass yield. At low salinity, by-products (mainly gluconate, pyruvate, and acetate) accumulate extracellularly. Using [1-¹³C]-, [2-¹³C]-, [6-¹³C]-, and [U-¹³C₆]glucose as carbon sources, we were able to determine the main central metabolic pathways involved in ectoines biosynthesis from glucose. C. salexigens uses the Entner-Doudoroff pathway rather than the standard glycolytic pathway for glucose catabolism, and anaplerotic activity is high to replenish the TCA cycle with the intermediaries withdrawn for ectoines biosynthesis. Metabolic flux ratios at low and high salinity were similar, revealing a certain metabolic rigidity, probably due to its specialization to support high biosynthetic fluxes and partially explaining why metabolic yields are so highly affected by salinity. This work represents an important contribution to the elucidation of specific metabolic adaptations in compatible solute-accumulating halophilic bacteria.     

Metabolic flux analysis using stoichiometric models for Aspergillus niger: comparison under glucoamylase-producing and non-producing conditions

Aspergillus niger AB1.13 cultures with glucoamylase production (with D-glucose as substrate) and without glucoamylase production (with D-xylose as substrate) were characterized by metabolic flux analysis. Two comprehensive metabolic models for d-glucose- as well as for D-xylose-consumption were used to quantify and compare the metabolic fluxes through the central pathways of carbon metabolism at different pH-values. The models consist of the most relevant metabolic pathways for A. niger including glycolysis, pentose-phosphate pathway, citrate cycle, energy metabolism and anaplerotic reactions comprising two intracellular compartments, the cytoplasm and mitochondrion. When D-xylose was used as the sole carbon source, the relative flux of the substrate through the oxidative pentose-phosphate pathway (PPP) via G6P-dehydrogenase was unaffected by the pH-value of the culture medium. About 30% of D-xylose consumed was routed through the oxidative PPP. In contrast, the flux of D-glucose (i.e., under glucoamylase-producing conditions) through the oxidative PPP was remarkably higher and, in addition was significantly affected by the pH-value of the culture medium (40% at pH 5.5, 56% at pH 3.7, respectively). Summarizing, the flux through the PPP under glucoamylase producing conditions was 30-90% higher than for non-producing conditions.     

Pathway confirmation and flux analysis of central metabolic pathways in Desulfovibrio vulgaris hildenborough using gas chromatography-mass spectrometry and Fourier transform-ion cyclotron resonance mass spectrometry

Flux distribution in central metabolic pathways of Desulfovibrio vulgaris Hildenborough was examined using 13C tracer experiments. Consistent with the current genome annotation and independent evidence from enzyme activity assays, the isotopomer results from both gas chromatography-mass spectrometry (GC-MS) and Fourier transform-ion cyclotron resonance mass spectrometry (FT-ICR MS) indicate the lack of an oxidatively functional tricarboxylic acid (TCA) cycle and an incomplete pentose phosphate pathway. Results from this study suggest that fluxes through both pathways are limited to biosynthesis. The data also indicate that >80% of the lactate was converted to acetate and that the reactions involved are the primary route of energy production [NAD(P)H and ATP production]. Independently of the TCA cycle, direct cleavage of acetyl coenzyme A to CO and 5,10-methyl tetrahydrofuran also leads to production of NADH and ATP. Although the genome annotation implicates a ferredoxin-dependent oxoglutarate synthase, isotopic evidence does not support flux through this reaction in either the oxidative or the reductive mode; therefore, the TCA cycle is incomplete. FT-ICR MS was used to locate the labeled carbon distribution in aspartate and glutamate and confirmed the presence of an atypical enzyme for citrate formation suggested in previous reports [the citrate synthesized by this enzyme is the isotopic antipode of the citrate synthesized by the (S)-citrate synthase]. These findings enable a better understanding of the relation between genome annotation and actual metabolic pathways in D. vulgaris and also demonstrate that FT-ICR MS is a powerful tool for isotopomer analysis, overcoming the problems with both GC-MS and nuclear magnetic resonance spectroscopy.     

Effect of sucA or sucC gene knockout on the metabolism in Escherichia coli based on gene expressions, enzyme activities, intracellular metabolite concentrations and metabolic fluxes by C-labeling experiments

The effect of sucA or sucC gene knockout on the metabolism in Escherichia coli was investigated for the aerobic cell growth in batch and continuous cultivations based on gene expressions, enzyme activities, intracellular metabolite concentrations and metabolic flux analysis. In the batch cultivation, the cell growth rate and the glucose uptake rate were lower for sucA mutant as compared with the parent strain, while it was not the case for sucC mutant. A significantly higher amount of acetate was produced, and it was not utilized in sucC mutant, while a little less acetate was produced in sucA mutant as compared with the parent strain. Unlike the parent strain and sucC mutant, sucA mutant excreted a little amount of l-glutamate. Enzyme activity results show that some of the glycolytic enzymes such as Tpi and Pgk were up-regulated, while Pfk, Fba and Pyk activities were down-regulated for sucA mutant as compared with the parent strain. For sucC mutant, the activities of Pfk, Fba, Tpi, GAPDH, Pgk and Pyk activities were down-regulated. As for the TCA cycle enzymes, the activities of CS and ICDH were down-regulated, while those of Icl, MS, Fum and MDH were up-regulated for sucA mutant. The activities of the oxidative pentose phosphate (PP) pathway enzymes such as G6PDH and 6PGDH and the gluconeogenic pathway enzyme such as Mez were up-regulated in sucA mutant. The Ack activity was down-regulated for sucA mutant, but not for sucC mutant. In continuous cultivation, the gene expression results indicate that the global regulatory genes such as fadR and iclR were slightly down-regulated in sucA mutant, which enhanced the expression of aceA gene and caused the up-regulation of the isocitrate lyase activity in sucA mutant, while fadR and iclR of sucC mutant changed little and no isocitrate lyase activation was observed for sucC mutant. Some other global regulatory genes such as arcA and fnr genes were down-regulated in both mutants, which caused some of the TCA cycle genes to be up-regulated. The effect of the sucA gene knockout on the metabolic flux distributions was investigated based on ¹H–¹³C NMR spectra and GC–MS signals obtained from ¹³C-labeling experiments. Flux analysis results indicate that the knockout of sucA gene caused the activation of PP pathway and the glyoxylate shunt. The fluxes through glycolysis and the TCA cycle were down-regulated in the sucA mutant. On the other hand, the fluxes through PP pathway and the anaplerotic reactions of Ppc-Pck and Mez increased.     

Evolution of metabolic pathways by chance assembly of enzyme proteins generated from sense and antisense strands of pre-existing genes.

In order to get an insight into the evolutionary aspect of metabolic pathways, especially of the ubiquitous glycolytic pathway, we have carried out an extensive search of sense-sense and sense-antisense similarities for enzyme proteins in the glycolytic pathway, the pentose phosphate cycle, alcohol and lactate fermentation pathways and the TCA cycle. This investigation of amino acid sequences reveals a curious pattern of similarity relations; no similarity can be found between the enzyme proteins in a section of the glycolytic pathway where the glyceraldehyde-3-phosphate or even glycerol-3-phosphate is converted into the pyruvate while many examples of sense-sense and sense-antisense similarities are found even between enzyme proteins in distant blocks, e.g. between the proteins in the TCA cycle and those in the pentose phosphate cycle, as well as between the functionally associated proteins in each of these blocks. Complementary to this characteristic pattern of amino acid sequence similarity, the search for similarities of nucleotide sequences also finds that the similarities of glycolytic enzyme genes, some sense-sense and others sense-antisense similarities, are concentrated on the nucleotide sequences of prokaryotic 16S or eukaryotic 18S ribosomal RNA gene with its flanks, although some of the copy sequences are also found in transfer RNA genes as well as in 23S or 26S ribosomal RNA gene. These results strongly suggest that the metabolic pathways have been developed by the chance assembly of enzyme proteins generated from the sense and antisense strands of pre-existing genes, e.g. the fermentation pathways and pentose phosphate cycle by the proteins from the genes of enzymes in the glycolytic pathway and the TCA cycle from all these successively increased genes, ascribing the origin of metabolic enzyme genes to the close relation between the glycolytic enzyme protein genes and the RNA gene cluster.