Branch point control by the phosphorylation state of isocitrate dehydrogenase. A quantitative examination of fluxes during a regulatory transition

Merging Metabolic Modeling and Imaging for Screening Therapeutic Targets in Colorectal Cancer

Cancer-associated fibroblasts (CAFs) play a key role in metabolic reprogramming and are well-established contributors to drug resistance in colorectal cancer (CRC). To exploit this metabolic crosstalk, we integrated a systems biology approach that identified key metabolic targets in a data-driven method and validated them experimentally. This process involved high-throughput computational screening to investigate the effects of enzyme perturbations predicted by a computational model of CRC metabolism to understand system-wide effects efficiently. Our results highlighted hexokinase (HK) as one of the crucial targets, which subsequently became our focus for experimental validation using patient-derived tumor organoids (PDTOs). Through metabolic imaging and viability assays, we found that PDTOs cultured in CAF conditioned media exhibited increased sensitivity to HK inhibition. Our approach emphasizes the critical role of integrating computational and experimental techniques in exploring and exploiting CRC-CAF crosstalk.

Multiple strategies for metabolic engineering of Escherichia coli for efficient production of glycolate

Glycolate is a bulk chemical with wide applications in the textile, food processing, and pharmaceutical industries. Glycolate can be produced from glucose via the glycolysis and glyoxylate shunt pathways, followed by reduction to glycolate. However, two problems limit the productivity and yield of glycolate when using glucose as the sole carbon source. The first is a cofactor imbalance in the production of glycolate from glucose via the glycolysis pathway, since NADPH is required for glycolate production, while glycolysis generates NADH. To rectify this imbalance, the NADP⁺ -dependent glyceraldehyde 3-phosphate dehydrogenase GapC from Clostridium acetobutylicum was introduced to generate NADPH instead of NADH in the oxidation of glyceraldehyde 3-phosphate during glycolysis. The soluble transhydrogenase SthA was further eliminated to conserve NADPH by blocking its conversion into NADH. The second problem is an unfavorable carbon flux distribution between the tricarboxylic acid cycle and the glyoxylate shunt. To solve this problem, isocitrate dehydrogenase (ICDH) was eliminated to increase the carbon flux of glyoxylate and thereby improve the glycolate titer. After engineering through the integration of gapC, combined with the inactivation of ICDH, SthA, and by-product pathways, as well as the upregulation of the two key enzymes isocitrate lyase (encoding by aceA), and glyoxylate reductase (encoding by ycdW), the glycolate titer increased to 5.3 g/L with a yield of 1.89 mol/mol glucose. Moreover, an optimized fed-batch fermentation reached a titer of 41 g/L with a yield of 1.87 mol/mol glucose after 60 h.

Acetate as a potential feedstock for the production of value-added chemicals: Metabolism and applications

Acetate is regarded as a promising carbon feedstock in biological production owing to its possible derivation from C₁ gases such as CO, CO₂ and methane. To best use of acetate, comprehensive understanding of acetate metabolisms from genes and enzymes to pathways and regulations is needed. This review aims to provide an overview on the potential of acetate as carbon feedstock for industrial biotechnology. Biochemical, microbial and biotechnological aspects of acetate metabolism are described. Especially, the current state-of-the art in the production of value-added chemicals from acetate is summarized. Challenges and future perspectives are also provided.

Multi-omics analysis unravels a segregated metabolic flux network that tunes co-utilization of sugar and aromatic carbons in Pseudomonas putida

Pseudomonas species thrive in different nutritional environments and can catabolize divergent carbon substrates. These capabilities have important implications for the role of these species in natural and engineered carbon processing. However, the metabolic phenotypes enabling Pseudomonas to utilize mixed substrates remain poorly understood. Here, we employed a multi-omics approach involving stable isotope tracers, metabolomics, fluxomics, and proteomics in Pseudomonas putida KT2440 to investigate the constitutive metabolic network that achieves co-utilization of glucose and benzoate, respectively a monomer of carbohydrate polymers and a derivative of lignin monomers. Despite nearly equal consumption of both substrates, metabolite isotopologues revealed nonuniform assimilation throughout the metabolic network. Gluconeogenic flux of benzoate-derived carbons from the tricarboxylic acid cycle did not reach the upper Embden-Meyerhof-Parnas pathway nor the pentose-phosphate pathway. These latter two pathways were populated exclusively by glucose-derived carbons through a cyclic connection with the Entner-Doudoroff pathway. We integrated the ¹³C-metabolomics data with physiological parameters for quantitative flux analysis, demonstrating that the metabolic segregation of the substrate carbons optimally sustained biosynthetic flux demands and redox balance. Changes in protein abundance partially predicted the metabolic flux changes in cells grown on the glucose:benzoate mixture versus on glucose alone. Notably, flux magnitude and directionality were also maintained by metabolite levels and regulation of phosphorylation of key metabolic enzymes. These findings provide new insights into the metabolic architecture that affords adaptability of P. putida to divergent carbon substrates and highlight regulatory points at different metabolic nodes that may underlie the high nutritional flexibility of Pseudomonas species.

Machine learning applied to enzyme turnover numbers reveals protein structural correlates and improves metabolic models

Knowing the catalytic turnover numbers of enzymes is essential for understanding the growth rate, proteome composition, and physiology of organisms, but experimental data on enzyme turnover numbers is sparse and noisy. Here, we demonstrate that machine learning can successfully predict catalytic turnover numbers in Escherichia coli based on integrated data on enzyme biochemistry, protein structure, and network context. We identify a diverse set of features that are consistently predictive for both in vivo and in vitro enzyme turnover rates, revealing novel protein structural correlates of catalytic turnover. We use our predictions to parameterize two mechanistic genome-scale modelling frameworks for proteome-limited metabolism, leading to significantly higher accuracy in the prediction of quantitative proteome data than previous approaches. The presented machine learning models thus provide a valuable tool for understanding metabolism and the proteome at the genome scale, and elucidate structural, biochemical, and network properties that underlie enzyme kinetics.

Multiple Optimal Phenotypes Overcome Redox and Glycolytic Intermediate Metabolite Imbalances in Escherichia coli pgi Knockout Evolutions

A mechanistic understanding of how new phenotypes develop to overcome the loss of a gene product provides valuable insight on both the metabolic and regulatory functions of the lost gene. The pgi gene, whose product catalyzes the second step in glycolysis, was deleted in a growth-optimized Escherichia coli K-12 MG1655 strain. The initial knockout (KO) strain exhibited an 80% drop in growth rate that was largely recovered in eight replicate, but phenotypically distinct, cultures after undergoing adaptive laboratory evolution (ALE). Multi-omic data sets showed that the loss of pgi substantially shifted pathway usage, leading to a redox and sugar phosphate stress response. These stress responses were overcome by unique combinations of innovative mutations selected for by ALE. Thus, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.IMPORTANCE A mechanistic understanding of how microbes are able to overcome the loss of a gene through regulatory and metabolic changes is not well understood. Eight independent adaptive laboratory evolution (ALE) experiments with pgi knockout strains resulted in eight phenotypically distinct endpoints that were able to overcome the gene loss. Utilizing multi-omics analysis, the coordinated mechanisms from genome to metabolome that lead to multiple optimal phenotypes after the loss of a major gene product were revealed.

Gluconeogenic precursor availability regulates flux through the glyoxylate shunt in Pseudomonas aeruginosa

The glyoxylate shunt bypasses the oxidative decarboxylation steps of the tricarboxylic acid (TCA) cycle, thereby conserving carbon skeletons for gluconeogenesis and biomass production. In Escherichia coli, carbon flux is redirected through the first enzyme of the glyoxylate shunt, isocitrate lyase (ICL), following phosphorylation and inactivation of the TCA cycle enzyme, isocitrate dehydrogenase (ICD), by the kinase/phosphatase, AceK. In contrast, mycobacterial species lack AceK and employ a phosphorylation-insensitive isocitrate dehydrogenase (IDH), which is allosterically activated by the product of ICL activity, glyoxylate. However, Pseudomonas aeruginosa expresses IDH, ICD, ICL, and AceK, raising the question of how these enzymes are regulated to ensure proper flux distribution between the competing pathways. Here, we present the structure, kinetics, and regulation of ICL, IDH, and ICD from P. aeruginosa We found that flux partitioning is coordinated through reciprocal regulation of these enzymes, linking distribution of carbon flux to the availability of the key gluconeogenic precursors, oxaloacetate and pyruvate. Specifically, a greater abundance of these metabolites activated IDH and inhibited ICL, leading to increased TCA cycle flux. Regulation was also exerted through AceK-dependent phosphorylation of ICD; high levels of acetyl-CoA (which would be expected to accumulate when oxaloacetate is limiting) stimulated the kinase activity of AceK, whereas high levels of oxaloacetate stimulated its phosphatase activity. In summary, the TCA cycle-glyoxylate shunt branch point in P. aeruginosa has a complex enzymology that is profoundly different from those in other species characterized to date. Presumably, this reflects its predilection for consuming fatty acids, especially during infection scenarios.

Engineering of Escherichia coli for Krebs cycle-dependent production of malic acid

Background

Malate is a C4-dicarboxylic acid widely used as an acidulant in the food and beverage industry. Rational engineering has been performed in the past for the development of microbial strains capable of efficient production of this metabolite. However, as malate can be a precursor for specialty chemicals, such as 2,4-dihydroxybutyric acid, that require additional cofactors NADP(H) and ATP, we set out to reengineer Escherichia coli for Krebs cycle-dependent production of malic acid that can satisfy these requirements.

Results

We found that significant malate production required at least simultaneous deletion of all malic enzymes and dehydrogenases, and concomitant expression of a malate-insensitive PEP carboxylase. Metabolic flux analysis using ¹³C-labeled glucose indicated that malate-producing strains had a very high flux over the glyoxylate shunt with almost no flux passing through the isocitrate dehydrogenase reaction. The highest malate yield of 0.82 mol/mol was obtained with E. coli Δmdh Δmqo ΔmaeAB ΔiclR ΔarcA which expressed malate-insensitive PEP carboxylase Ppc^K620S and NADH-insensitive citrate synthase GltA^R164L. We also showed that inactivation of the dicarboxylic acid transporter DcuA strongly reduced malate production arguing for a pivotal role of this permease in malate export.

Conclusions

Since more NAD(P)H and ATP cofactors are generated in the Krebs cycle-dependent malate production when compared to pathways which depend on the function of anaplerotic PEP carboxylase or PEP carboxykinase enzymes, the engineered strain developed in this study can serve as a platform to increase biosynthesis of malate-derived metabolites such as 2,4-dihydroxybutyric acid.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-0959-y) contains supplementary material, which is available to authorized users.

Dynamics of benzoate metabolism in Pseudomonas putida KT2440

Soil microorganisms mineralize lignin-derived aromatic carbon sources using oxidative catabolic pathways, such as the β-ketoadipate pathway. Although this aromatic pathway is one of the best-studied pathways in biochemistry, the complete pathway, including its regulation by aromatic carbon sources, has not been integrated into the metabolic network. In particular, information about the in vivo operation (e.g., kinetics and flux capacity) of the pathway is lacking. In this contribution, we use kinetic modeling and thermodynamic analysis to evaluate the in vivo operation of this key aromatic multi-step pathway. The resulting ab initio deterministic model of benzoate degradation via the β-ketoadipate (ortho-cleavage) pathway in Pseudomonas putida KT2440 is presented. The kinetic model includes mechanistic rate expressions for the enzymes and transport processes. The design and experimental validation of the model are driven by data generated from short-term perturbation experiments in a benzoate-limited continuous culture. The results of rigorous modeling of the in vivo dynamics provide strong support for flux regulation by the benzoate transporter and the enzymes forming and cleaving catechol. Revisiting the β-ketoadipate pathway might be valuable for applications in different fields, such as biochemistry and metabolic engineering, that use lignin monomers as a carbon source.

•

We describe a kinetic model for the β-ketoadipate pathway.

•

Short term metabolic responses were tracked on metabolite level by rapid sampling.

•

The model captures steady state and dynamic conditions of the β-ketoadipate pathway.

•

Thermodynamic analysis revealed regulation points of the pathway.

•

The results are discussed in the context of metabolic network operation.

NADPH-generating systems in bacteria and archaea

Reduced nicotinamide adenine dinucleotide phosphate (NADPH) is an essential electron donor in all organisms. It provides the reducing power that drives numerous anabolic reactions, including those responsible for the biosynthesis of all major cell components and many products in biotechnology. The efficient synthesis of many of these products, however, is limited by the rate of NADPH regeneration. Hence, a thorough understanding of the reactions involved in the generation of NADPH is required to increase its turnover through rational strain improvement. Traditionally, the main engineering targets for increasing NADPH availability have included the dehydrogenase reactions of the oxidative pentose phosphate pathway and the isocitrate dehydrogenase step of the tricarboxylic acid (TCA) cycle. However, the importance of alternative NADPH-generating reactions has recently become evident. In the current review, the major canonical and non-canonical reactions involved in the production and regeneration of NADPH in prokaryotes are described, and their key enzymes are discussed. In addition, an overview of how different enzymes have been applied to increase NADPH availability and thereby enhance productivity is provided.

Making metabolism accessible and meaningful: is the definition of a central metabolic dogma within reach?

Intermediary metabolism, a dominant research area before the emergence of molecular biology, is attracting renewed interest for fundamental and applied reasons as documented here. Nonetheless, the field may appear to be a thicket precluding entry to all but the most determined. Here we present a metabolic overview that makes this important and fascinating area accessible to a broad range of the molecular biological and biotechnological communities that are being attracted to biological problems crying out for metabolic solutions. This is accomplished by identifying seven key concepts, a so-called metabolic central dogma, that provide a core understanding analogous to the "Central Dogma of Molecular Biology" which focused upon maintenance and flow of genetic information.

The functional structure of central carbon metabolism in Pseudomonas putida KT2440

What defines central carbon metabolism? The classic textbook scheme of central metabolism includes the Embden-Meyerhof-Parnas (EMP) pathway of glycolysis, the pentose phosphate pathway, and the citric acid cycle. The prevalence of this definition of central metabolism is, however, equivocal without experimental validation. We address this issue using a general experimental approach that combines the monitoring of transcriptional and metabolic flux changes between steady states on alternative carbon sources. This approach is investigated by using the model bacterium Pseudomonas putida with glucose, fructose, and benzoate as carbon sources. The catabolic reactions involved in the initial uptake and metabolism of these substrates are expected to show a correlated change in gene expressions and metabolic fluxes. However, there was no correlation for the reactions linking the 12 biomass precursor molecules, indicating a regulation mechanism other than mRNA synthesis for central metabolism. This result substantiates evidence for a (re)definition of central carbon metabolism including all reactions that are bound to tight regulation and transcriptional invariance. Contrary to expectations, the canonical Entner-Doudoroff and EMP pathways sensu stricto are not a part of central carbon metabolism in P. putida, as they are not regulated differently from the aromatic degradation pathway. The regulatory analyses presented here provide leads on a qualitative basis to address the use of alternative carbon sources by deregulation and overexpression at the transcriptional level, while rate improvements in central carbon metabolism require careful adjustment of metabolite concentrations, as regulation resides to a large extent in posttranslational and/or metabolic regulation.

A comparative analysis of industrial Escherichia coli K-12 and B strains in high-glucose batch cultivations on process-, transcriptome- and proteome level

Escherichia coli K-12 and B strains are among the most frequently used bacterial hosts for production of recombinant proteins on an industrial scale. To improve existing processes and to accelerate bioprocess development, we performed a detailed host analysis. We investigated the different behaviors of the E. coli production strains BL21, RV308, and HMS174 in response to high-glucose concentrations. Tightly controlled cultivations were conducted under defined environmental conditions for the in-depth analysis of physiological behavior. In addition to acquisition of standard process parameters, we also used DNA microarray analysis and differential gel electrophoresis (Ettan(TM) DIGE). Batch cultivations showed different yields of the distinct strains for cell dry mass and growth rate, which were highest for BL21. In addition, production of acetate, triggered by excess glucose supply, was much higher for the K-12 strains compared to the B strain. Analysis of transcriptome data showed significant alteration in 347 of 3882 genes common among all three hosts. These differentially expressed genes included, for example, those involved in transport, iron acquisition, and motility. The investigation of proteome patterns additionally revealed a high number of differentially expressed proteins among the investigated hosts. The subsequently selected 38 spots included proteins involved in transport and motility. The results of this comprehensive analysis delivered a full genomic picture of the three investigated strains. Differentially expressed groups for targeted host modification were identified like glucose transport or iron acquisition, enabling potential optimization of strains to improve yield and process quality. Dissimilar growth profiles of the strains confirm different genotypes. Furthermore, distinct transcriptome patterns support differential regulation at the genome level. The identified proteins showed high agreement with the transcriptome data and suggest similar regulation within a host at both levels for the identified groups. Such host attributes need to be considered in future process design and operation.

Flux analysis and metabolomics for systematic metabolic engineering of microorganisms

Rational engineering of metabolism is important for bio-production using microorganisms. Metabolic design based on in silico simulations and experimental validation of the metabolic state in the engineered strain helps in accomplishing systematic metabolic engineering. Flux balance analysis (FBA) is a method for the prediction of metabolic phenotype, and many applications have been developed using FBA to design metabolic networks. Elementary mode analysis (EMA) and ensemble modeling techniques are also useful tools for in silico strain design. The metabolome and flux distribution of the metabolic pathways enable us to evaluate the metabolic state and provide useful clues to improve target productivity. Here, we reviewed several computational applications for metabolic engineering by using genome-scale metabolic models of microorganisms. We also discussed the recent progress made in the field of metabolomics and (13)C-metabolic flux analysis techniques, and reviewed these applications pertaining to bio-production development. Because these in silico or experimental approaches have their respective advantages and disadvantages, the combined usage of these methods is complementary and effective for metabolic engineering.

Isocitrate dehydrogenase from Streptococcus mutans: biochemical properties and evaluation of a putative phosphorylation site at Ser102

Isocitrate deyhdrogenase (IDH) is a reversible enzyme in the tricarboxylic acid cycle that catalyzes the NAD(P)⁺-dependent oxidative decarboxylation of isocitrate to α-ketoglutarate (αKG) and the NAD(P)H/CO₂-dependent reductive carboxylation of αKG to isocitrate. The IDH gene from Streptococcus mutans was fused with the icd gene promoter from Escherichia coli to initiate its expression in the glutamate auxotrophic strain E. coli Δicd::kan^r of which the icd gene has been replaced by kanamycin resistance gene. The expression of S. mutans IDH (SmIDH) may restore the wild-type phenotype of the icd-defective strain on minimal medium without glutamate. The molecular weight of SmIDH was estimated to be 70 kDa by gel filtration chromatography, suggesting a homodimeric structure. SmIDH was divalent cation-dependent and Mn²⁺ was found to be the most effective cation. The optimal pH of SmIDH was 7.8 and the maximum activity was around 45°C. SmIDH was completely NAD⁺ dependent and its apparent K_m for NAD⁺ was 137 μM. In order to evaluate the role of the putative phosphorylation site at Ser102 in catalysis, two “stably phosphorylated” mutants were constructed by converting Ser102 into Glu102 or Asp102 in SmIDH to mimick a constitutively phosphorylated state. Meanwhile, the functional roles of another four amino acids (threonine, glycine, alanine and tyrosine) containing variant size of side chains were investigated. The replacement of Asp102 or Glu102 totally inactivated the enzyme, while the S102T, S102G, S102A and S102Y mutants decreased the affinity to isocitrate and only retained 16.0%, 2.8%, 3.3% and 1.1% of the original activity, respectively. These results reveal that Ser102 plays important role in substrate binding and is required for the enzyme function. Also, Ser102 in SmIDH is a potential phosphorylation site, indicating that the ancient NAD-dependent IDHs might be the underlying origin of “phosphorylation mechanism” used by their bacterial NADP-dependent homologs.

Resolving cell composition through simple measurements, genome-scale modeling, and a genetic algorithm

The biochemical composition of a cell is very complex and dynamic. It varies greatly among different organisms and environmental conditions. Inclusion of proper cell composition data is critical for accurate genome-scale metabolic flux modeling using flux balance analysis (FBA). However, determining cell composition experimentally is currently time-consuming and resource intensive. In this chapter, a method for predicting cell composition using a genome-scale model and "easy to measure" culture data (e.g., glucose uptake rate, and specific growth rate) is presented. The method makes use of a genetic algorithm for nonlinear optimization of a biomass equation (a mathematical description of cell composition). As a case study, the method was used to optimize a biomass equation for Escherichia coli MG1655 under multiple growth environments. The availability of experimentally determined (13)C flux data allowed a direct comparison with FBA predicted fluxes through the TCA cycle. Results showed dramatic improvement upon optimization of the biomass equation. In a second case study, biomass equation optimization was also applied to Clostridium acetobutylicum, an organism with less available biochemical cell composition data in the literature. The method produced a biomass equation highly similar to one determined experimentally for the closely related Gram-positive Bacillus subtilis.

Simple constrained-optimization view of acetate overflow in E. coli

The production of acetate by aerobically growing E. coli is examined. The problem is formulated in terms of a flow network that has as its objective maximal ATP synthesis. It is found that when loads are imposed and flux constraints exist either at the level of NADH turnover rate or the activity of a key Krebs cycle enzyme, switching to acetate overflow is predicted. Moreover, the result found for the latter constraint can be shown to be formally equivalent to a correlation experimentally determined for the specific rate of acetate production by E. coli K-12.

Ultrasensitive regulation of anapleurosis via allosteric activation of PEP carboxylase

Anapleurosis is the filling of the tricarboxylic acid cycle with four-carbon units. The common substrate for both anapleurosis and glucose phosphorylation in bacteria is the terminal glycolytic metabolite phosphoenolpyruvate (PEP). Here we show that Escherichia coli quickly and almost completely turns off PEP consumption upon glucose removal. The resulting buildup of PEP is used to quickly import glucose if it becomes available again. The switch-like termination of anapleurosis results from depletion of fructose-1,6-bisphosphate (FBP), an ultrasensitive allosteric activator of PEP carboxylase. E. coli expressing an FBP-insensitive point mutant of PEP carboxylase grow normally when glucose is steadily available. However, they fail to build up PEP upon glucose removal, grow poorly when glucose availability oscillates and suffer from futile cycling at the PEP node on gluconeogenic substrates. Thus, bacterial central carbon metabolism is intrinsically programmed with ultrasensitive allosteric regulation to enable rapid adaptation to changing environmental conditions.

Ligand binding and structural changes associated with allostery in yeast NAD(+)-specific isocitrate dehydrogenase

Yeast NAD(+)-specific isocitrate dehydrogenase (IDH) is an octameric enzyme composed of four each of regulatory IDH1 and catalytic IDH2 subunits that share 42% sequence identity. IDH2 contains catalytic isocitrate/Mg2+ and NAD+ binding sites whereas IDH1 contains homologous binding sites, respectively, for cooperative binding of isocitrate and for allosteric binding of AMP. Ligand binding is highly ordered in vitro, and IDH exhibits the unusual property of half-site binding for all ligands. The structures of IDH solved in the absence or presence of ligands have shown: (a) a heterodimer to be the basic structural/functional unit of the enzyme, (b) the organization of heterodimers to form tetramer and octamer structures, (c) structural differences that may underlie cooperative and allosteric regulatory mechanisms, and (d) the possibility for formation of a disulfide bond that could reduce catalytic activity. In vivo analyses of mutant enzymes have elucidated the physiological importance of catalytic activity and allosteric regulation of this tricarboxylic acid cycle enzyme. Other studies have established the importance of a disulfide bond in regulation of IDH activity in vivo, as well as contributions of this bond to the property of half-site ligand binding exhibited by the wild-type enzyme.

Basis for half-site ligand binding in yeast NAD(+)-specific isocitrate dehydrogenase

Yeast NAD(+)-specific isocitrate dehydrogenase is an allosterically regulated octameric enzyme composed of four heterodimers of a catalytic IDH2 subunit and a regulatory IDH1 subunit. Despite structural predictions that the enzyme would contain eight isocitrate binding sites, four NAD(+) binding sites, and four AMP binding sites, only half of the sites for each ligand can be measured in binding assays. On the basis of a potential interaction between side chains of Cys-150 residues in IDH2 subunits in each tetramer of the enzyme, ligand binding assays of wild-type (IDH1/IDH2) and IDH1/IDH2(C150S) octameric enzymes were conducted in the presence of dithiothreitol. These assays demonstrated the presence of eight isocitrate and four AMP binding sites for the wild-type enzyme in the presence of dithiothreitol and for the IDH1/IDH2(C150S) enzyme in the absence or presence of this reagent, suggesting that interactions between sulfhydryl side chains of IDH2 Cys-150 residues limit access to these sites. However, only two NAD(+) sites could be measured for either enzyme. A tetrameric form of IDH (an IDH1(G15D)/IDH2 mutant enzyme) demonstrated half-site binding for isocitrate (two sites) in the absence of dithiothreitol and full-site binding (four sites) in the presence of dithiothreitol. Only one NAD(+) site could be measured for the tetramer under both conditions. In the context of the structure of the enzyme, these results suggest that an observed asymmetry between heterotetramers in the holoenzyme contributes to interactions between IDH2 Cys-150 residues and to half-site binding of isocitrate, but that a form of negative cooperativity may limit access to apparently equivalent NAD(+) binding sites.

Construction and analyses of tetrameric forms of yeast NAD+-specific isocitrate dehydrogenase

Yeast NAD(+)-specific isocitrate dehydrogenase (IDH) is an octameric enzyme composed of four heterodimers of regulatory IDH1 and catalytic IDH2 subunits. The crystal structure suggested that the interactions between tetramers in the octamer are restricted to defined regions in IDH1 subunits from each tetramer. Using truncation and mutagenesis, we constructed three tetrameric forms of IDH. Truncation of five residues from the amino terminus of IDH1 did not alter the octameric form of the enzyme, but this truncation with an IDH1 G15D or IDH1 D168K residue substitution produced tetrameric enzymes as assessed by sedimentation velocity ultracentrifugation. The IDH1 G15D substitution in the absence of any truncation of IDH1 was subsequently found to be sufficient for production of a tetrameric enzyme. The tetrameric forms of IDH exhibited ∼50% reductions in V(max) and in cooperativity with respect to isocitrate relative to those of the wild-type enzyme, but they retained the property of allosteric activation by AMP. The truncated (-5)IDH1/IDH2 and tetrameric enzymes were much more sensitive than the wild-type enzyme to inhibition by the oxidant diamide and concomitant formation of a disulfide bond between IDH2 Cys-150 residues. Binding of ligands reduced the sensitivity of the wild-type enzyme to diamide but had no effect on inhibition of the truncated or tetrameric enzymes. These results suggest that the octameric structure of IDH has in part evolved for regulation of disulfide bond formation and activity by ensuring the proximity of the amino terminus of an IDH1 subunit of one tetramer to the IDH2 Cys-150 residues in the other tetramer.

13C-metabolic flux analysis for batch culture of Escherichia coli and its Pyk and Pgi gene knockout mutants based on mass isotopomer distribution of intracellular metabolites

Since most bio-production processes are conducted in a batch or fed-batch manner, the evaluation of metabolism with respect to time is highly desirable. Toward this aim, we applied (13)C-metabolic flux analysis to nonstationary conditions by measuring the mass isotopomer distribution of intracellular metabolites. We performed our analysis on batch cultures of wild-type Escherichia coli, as well as on Pyk and Pgi mutants, obtained the fluxes and metabolite concentrations as a function of time. Our results for the wild-type indicated that the TCA cycle flux tended to increase during growth on glucose. Following glucose exhaustion, cells controlled the branch ratio between the glyoxylate pathway and the TCA cycle, depending on the availability of acetate. In the Pyk mutant, the concentrations of glycolytic intermediates changed drastically over time due to the dumping and feedback inhibition caused by PEP accumulation. Nevertheless, the flux distribution and free amino acid concentrations changed little. The growth rate and the fluxes remained constant in the Pgi mutant and the glucose-6-phosphate dehydrogenase reaction was the rate-limiting step. The measured fluxes were compared with those predicted by flux balance analysis using maximization of biomass yield or ATP production. Our findings indicate that the objective function of biosynthesis became less important as time proceeds on glucose in the wild-type, while it remained highly important in the Pyk mutant. Furthermore, ATP production was the primary objective function in the Pgi mutant. This study demonstrates how cells adjust their metabolism in response to environmental changes and/or genetic perturbations in the batch cultivation.

Carbon flux rerouting during Mycobacterium tuberculosis growth arrest

A hallmark of the Mycobacterium tuberculosis life cycle is the pathogen's ability to switch between replicative and non-replicative states in response to host immunity. Transcriptional profiling by qPCR of ∼ 50 M. tuberculosis genes involved in central and lipid metabolism revealed a re-routing of carbon flow associated with bacterial growth arrest during mouse lung infection. Carbon rerouting was marked by a switch from metabolic pathways generating energy and biosynthetic precursors in growing bacilli to pathways for storage compound synthesis during growth arrest. Results of flux balance analysis using an in silico metabolic network were consistent with the transcript abundance data obtained in vivo. Similar transcriptional changes were seen in vitro when M. tuberculosis cultures were treated with bacteriostatic stressors under different nutritional conditions. Thus, altered expression of key metabolic genes reflects growth rate changes rather than changes in substrate availability. A model describing carbon flux rerouting was formulated that (i) provides a coherent interpretation of the adaptation of M. tuberculosis metabolism to immunity-induced stress and (ii) identifies features common to mycobacterial dormancy and stress responses of other organisms.

Biochemical production capabilities of Escherichia coli

Microbial metabolism provides at mechanism for the conversion of substrates into useful biochemicals. Utilization of microbes in industrial processes requires a modification of their natural metabolism in order to increase the efficiency of the desired conversion. Redirection of metabolic fluxes forms the basis of the newly defined field of metabolic engineering. In this study we use a flux balance based approach to study the biosynthesis of the 20 amino acids and 4 nucleotides as biochemical products. These amino acids and nucleotides are primary products of biosynthesis as well as important industrial products and precursors for the production of other biochemicals. The biosynthetic reactions of the bacterium Escherichia coli have been formulated into a metabolic network, and growth has been defined as a balanced drain on the metabolite pools corresponding to the cellular composition. Theoretical limits on the conversion of glucose, glycerol, and acetate substrates to biomass as well as the biochemical products have been computed. The substrate that results in the maximal carbon conversion to a particular product is identified. Criteria have been developed to identify metabolic constraints in the optimal solutions. The constraints of stoichiometry, energy, and redox have been determined in the conversions of glucose, glycerol, and acetate substrates into the biochemicals. Flux distributions corresponding to the maximal production of the biochemicals are presented. The goals of metabolic engineering are the optimal redirection of fluxes from generating biomass toward producing the desired biochemical. Optimal biomass generation is shown to decrease in a piecewise linear manner with increasing product formation. In some cases, synergy is observed between biochemical production and growth, leading to an increased overall carbon conversion. Balanced growth and product formation are important in a bioprocess, particularly for nonsecreted products.

Response of the central metabolism of Corynebacterium glutamicum to different flux burdens

To evaluate the importance of reactions within the central metabolism under different flux burdens the fluxes within the pentose phosphate pathway (PPP), as well as the other reactions of the central metabolism, were intensively analyzed and quantitated. For this purpose, Corynebacterium glutamicum was grown with [1-(13)C]glucose to metabolic and isotopic steady state and the fractional enrichments in precursor metabolites (e.g., pentose 5-phosphate) were quantified. Matrix calculus was used to express these data together with metabolite mass data. The detailed analysis of the dependence of (13)C enrichments on exchange fluxes enabled the transketolase-catalyzed exchange rate (2 pentose 5-phosphate <--> sedoheptulose 7-phosphate + glyceraldehyde 3-phosphate) to be quantified as 74.3% (molar metabolite flux) at a net flux of 10.3% and the exchange rate (pentose 5-phosphate + erythrose 4-phosphate <--> fructose 6-phosphate + glyceraldehyde 3-phosphate) to be quantified as 5.6% at a net flux of 8.1%. The flux entering the tricarboxylic acid cycle was 93.3%. The same comprehensive flux analysis as performed for the nonexcreting condition was done with the identical strain that had been forced to excrete L-glutamate. Because we had already quantified the fluxes for L-lysine excretion with an isogenic strain, three directly comparable flux situations are thus available. Consequently, this comparison permits a direct cause-and-effect relationship to be specified. In response to the different flux burdens of the cell, the PPP flux decreased from a maximum of 67% to 26%, with the glycolytic flux increasing accordingly. The carbon flux through isocitrate dehydrogenase increased from 20% to 36%. The bidirectional carbon flux between pyruvate and oxaloacetate decreased from 36% to 9%. Since the cause of the three different flux states was the allelic exchange in the final L-lysine assembling pathway or the glutamate export activity, respectively, the flexible response is the effect. This shows conclusively the enormous flexibility within the central metabolism of C. glutamicum to supply precursors upon their withdrawal for the synthesis of amino acids. (c) 1997 John Wiley & Sons, Inc. Biotechnol Bioeng 56: 168-180, 1997.

The biomass objective function

Flux balance analysis (FBA) is a mathematical approach for analyzing the flow of metabolites through a metabolic network. To computationally predict cell growth using FBA, one has to determine the biomass objective function that describes the rate at which all of the biomass precursors are made in the correct proportions. Here we review fundamental issues associated with its formulation and use to compute optimal growth states.

Stable isotope fingerprinting: a novel method for identifying plant, fungal, or bacterial origins of amino acids

Amino acids play an important role in ecology as essential nutrients for animals and as currencies in symbiotic associations. Here we present a new approach to tracing the origins of amino acids by identifying unique patterns of carbon isotope signatures generated by amino acid synthesis in plants, fungi, and bacteria ("13C fingerprints"). We measured amino acid delta 13C from 10 C3 plants, 13 fungi, and 10 bacteria collected and isolated from a boreal forest in interior Alaska, USA, using gas chromatography-combustion-isotope ratio mass spectrometry (GC-C-IRMS). Microorganisms were cultured under amino-acid-free conditions and identified based on DNA sequences. Bacteria, fungi, and plants generated consistent, unique 13C fingerprints based on the more complex amino acids (five or more biosynthetic steps) that are classified as essential for animals. Linear discriminant analysis classified all samples correctly with >99% certainty and correctly classified nearly all insect samples from a previous study by diet. Our results suggest that 13C fingerprints of amino acids could provide a powerful in situ assay of the biosynthetic sources of amino acids and a potential new tool for understanding nutritional linkages in food webs.

Genome-wide gene regulation of biosynthesis and energy generation by a novel transcriptional repressor in Geobacter species

Geobacter species play important roles in bioremediation of contaminated environments and in electricity production from waste organic matter in microbial fuel cells. To better understand physiology of Geobacter species, expression and function of citrate synthase, a key enzyme in the TCA cycle that is important for organic acid oxidation in Geobacter species, was investigated. Geobacter sulfurreducens did not require citrate synthase for growth with hydrogen as the electron donor and fumarate as the electron acceptor. Expression of the citrate synthase gene, gltA, was repressed by a transcription factor under this growth condition. Functional and comparative genomics approaches, coupled with genetic and biochemical assays, identified a novel transcription factor termed HgtR that acts as a repressor for gltA. Further analysis revealed that HgtR is a global regulator for genes involved in biosynthesis and energy generation in Geobacter species. The hgtR gene was essential for growth with hydrogen, during which hgtR expression was induced. These findings provide important new insights into the mechanisms by which Geobacter species regulate their central metabolism under different environmental conditions.

Disulfide bond formation in yeast NAD+-specific isocitrate dehydrogenase

The tricarboxylic acid cycle NAD+-specific isocitrate dehydrogenase (IDH) of Saccharomyces cerevisiae is an octameric enzyme composed of four heterodimers of regulatory IDH1 and catalytic IDH2 subunits. Recent structural analyses revealed the close proximity of Cys-150 residues from IDH2 in adjacent heterodimers, and features of the structure for the ligand-free enzyme suggested that formation of a disulfide bond between these residues might stabilize an inactive form of the enzyme. We constructed two mutant forms of IDH, one containing a C150S substitution in IDH2 and the other containing C56S/C242S substitutions in IDH2 leaving Cys-150 as the sole cysteine residue. Treatment of the affinity-purified enzymes with diamide resulted in the formation of disulfide bonds and in decreased activities for the wild-type and C56S/C242S enzymes. Both effects were reversible by the addition of dithiothreitol. Diamide had no effect on the C150S mutant enzyme, suggesting that Cys-150 is essential for the formation of a disulfide bond that inhibits IDH activity. Diamide-induced formation of the Cys-150 disulfide bond was also observed in vivo for yeast transformants expressing the wild-type or C56S/C242S enzymes but not for a transformant expressing the C150S enzyme. Finally, natural formation of the Cys-150 disulfide bond with a concomitant decrease in cellular IDH activity was observed during the stationary phase for the parental strain and for transformants expressing wild-type or C56S/C242S enzymes but not for a transformant expressing the C150S enzyme. A reduction in viability for the latter strain suggests that a decrease in IDH activity is important for metabolic changes in stationary phase cells.

In silico genome-scale reconstruction and validation of the Corynebacterium glutamicum metabolic network

A genome-scale metabolic model of the Gram-positive bacteria Corynebacterium glutamicum ATCC 13032 was constructed comprising 446 reactions and 411 metabolites, based on the annotated genome and available biochemical information. The network was analyzed using constraint based methods. The model was extensively validated against published flux data, and flux distribution values were found to correlate well between simulations and experiments. The split pathway of the lysine synthesis pathway of C. glutamicum was investigated, and it was found that the direct dehydrogenase variant gave a higher lysine yield than the alternative succinyl pathway at high lysine production rates. The NADPH demand of the network was not found to be critical for lysine production until lysine yields exceeded 55% (mmol lysine (mmol glucose)(-1)). The model was validated during growth on the organic acids acetate and lactate. Comparable flux values between in silico model and experimental values were seen, although some differences in the phenotypic behavior between the model and the experimental data were observed.

Allosteric motions in structures of yeast NAD+-specific isocitrate dehydrogenase

Mitochondrial NAD(+)-specific isocitrate dehydrogenases (IDHs) are key regulators of flux through biosynthetic and oxidative pathways in response to cellular energy levels. Here we present the first structures of a eukaryotic member of this enzyme family, the allosteric, hetero-octameric, NAD(+)-specific IDH from yeast in three forms: 1) without ligands, 2) with bound analog citrate, and 3) with bound citrate + AMP. The structures reveal the molecular basis for ligand binding to homologous but distinct regulatory and catalytic sites positioned at the interfaces between IDH1 and IDH2 subunits and define pathways of communication between heterodimers and heterotetramers in the hetero-octamer. Disulfide bonds observed at the heterotetrameric interfaces in the unliganded IDH hetero-octamer are reduced in the ligand-bound forms, suggesting a redox regulatory mechanism that may be analogous to the "on-off" regulation of non-allosteric bacterial IDHs via phosphorylation. The results strongly suggest that eukaryotic IDH enzymes are exquisitely tuned to ensure that allosteric activation occurs only when concentrations of isocitrate are elevated.

Glyoxylate and Pyruvate Are Antagonistic Effectors of the Escherichia coli IclR Transcriptional Regulator

The Escherichia coli isocitrate lyase regulator (IclR) regulates the expression of the glyoxylate bypass operon (aceBAK). Founding member of a large family of common fold transcriptional regulators, IclR comprises a DNA binding domain that interacts with the operator sequence and a C-terminal domain (C-IclR) that binds a hitherto unknown small molecule. We screened a chemical library of more than 150 metabolic scaffolds using a high-throughput protein stability assay to identify molecules that bind IclR and then tested the active compounds in in vitro assays of operator binding. Glyoxylate and pyruvate, identified by this method, bound the C-IclR domain with KD values of 0.9+/-0.2 and 156.2+/-7.9 microM, as defined by isothermal titration calorimetry. Both compounds altered IclR interactions with operator DNA in electrophoretic mobility shift assays but showed an antagonistic effect. Glyoxylate disrupted the formation of the IclR/operator complex in vitro by favoring the inactive dimeric state of the protein, whereas pyruvate increased the binding of IclR to the aceBAK promoter by stabilizing the active tetrameric form of the protein. Structures of the C-IclR domain alone and in complex with each effector were determined at 2.3 A, confirming the binding of both molecules in the effector recognition site previously characterized for the other representative of the family, the E. coli AllR regulator. Site-directed mutagenesis demonstrated the importance of hydrophobic patch formed by Met-146, Leu-154, Leu-220, and Leu-143 in interactions with effector molecules. In general, our strategy of combining chemical screens with functional assays and structural studies has uncovered two small molecules with antagonistic effects that regulate the IclR-dependent transcription of the aceBAK operon.

Role of phosphoenolpyruvate in the NADP-isocitrate dehydrogenase and isocitrate lyase reaction in Escherichia coli

Phosphoenolpyruvate inhibited Escherichia coli NADP-isocitrate dehydrogenase allosterically (Ki of 0.31 mM) and isocitrate lyase uncompetitively (Ki' of 0.893 mM). Phosphoenolpyruvate enhances the uncompetitive inhibition of isocitrate lyase by increasing isocitrate, which protects isocitrate dehydrogenase from the inhibition, and contributes to the control through the tricarboxylic acid cycle and glyoxylate shunt.

Effect of the global redox sensing/regulation networks on Escherichia coli and metabolic flux distribution based on C-13 labeling experiments

Escherichia coli has several elaborate sensing mechanisms for response to the availability of oxygen and the presence of other electron acceptors. Among them, the one component Fnr protein and the two-component Arc system coordinate the adaptive responses to oxygen availability. To systematically investigate the contribution of Arc- and Fnr-dependent regulation in catabolism, glucose-limited chemostat cultures were conducted on wild-type E. coli, an arcA mutant, an fnr mutant, and an arcAfnr double mutant strains under a well-defined semi-aerobic condition. The metabolic flux distributions of the cultures of these strains were estimated based on C-13 labeling experiments. It was shown that the oxidative pentose phosphate (PP) pathway was functioning at low level under semi-aerobic condition. The fluxes through pyruvate dehydrogenase (PDH) and tricarboxylic acid (TCA) cycle were found to be lower in the arcA mutant and the arcAfnr double mutant strains than that in the wild-type strain, although the expression of the genes involved in these pathways have been proved to be derepressed in the mutant strains ([Shalel-Levanon, S., San, K.Y., Bennett, G.N., 2005a. Effect of ArcA and FNR on the expression of genes related to the oxygen regulation and the glycolysis pathway in Escherichia coli under microaerobic growth conditions. Biotechnol. Bioeng. 92, 147-159; Shalel-Levanon, S., San, K.Y., Bennett, G.N., 2005c. Effect of oxygen, and ArcA and FNR regulators on the expression of genes related to the electron transfer chain and the TCA cycle in Escherichia coli. Metab. Eng. 7, 364-374]). The significantly higher lactate production in the arcAfnr double mutant strain was shown to be an indirect effect caused by the reduced pyruvate formate-lyase (PFL) and PDH fluxes as well as the intracellular redox state.

Kinetic modeling of tricarboxylic acid cycle and glyoxylate bypass in Mycobacterium tuberculosis, and its application to assessment of drug targets

Background

Targeting persistent tubercule bacilli has become an important challenge in the development of anti-tuberculous drugs. As the glyoxylate bypass is essential for persistent bacilli, interference with it holds the potential for designing new antibacterial drugs. We have developed kinetic models of the tricarboxylic acid cycle and glyoxylate bypass in Escherichia coli and Mycobacterium tuberculosis, and studied the effects of inhibition of various enzymes in the M. tuberculosis model.

Results

We used E. coli to validate the pathway-modeling protocol and showed that changes in metabolic flux can be estimated from gene expression data. The M. tuberculosis model reproduced the observation that deletion of one of the two isocitrate lyase genes has little effect on bacterial growth in macrophages, but deletion of both genes leads to the elimination of the bacilli from the lungs. It also substantiated the inhibition of isocitrate lyases by 3-nitropropionate. On the basis of our simulation studies, we propose that: (i) fractional inactivation of both isocitrate dehydrogenase 1 and isocitrate dehydrogenase 2 is required for a flux through the glyoxylate bypass in persistent mycobacteria; and (ii) increasing the amount of active isocitrate dehydrogenases can stop the flux through the glyoxylate bypass, so the kinase that inactivates isocitrate dehydrogenase 1 and/or the proposed inactivator of isocitrate dehydrogenase 2 is a potential target for drugs against persistent mycobacteria. In addition, competitive inhibition of isocitrate lyases along with a reduction in the inactivation of isocitrate dehydrogenases appears to be a feasible strategy for targeting persistent mycobacteria.

Conclusion

We used kinetic modeling of biochemical pathways to assess various potential anti-tuberculous drug targets that interfere with the glyoxylate bypass flux, and indicated the type of inhibition needed to eliminate the pathogen. The advantage of such an approach to the assessment of drug targets is that it facilitates the study of systemic effect(s) of the modulation of the target enzyme(s) in the cellular environment.

rRNA and poly-beta-hydroxybutyrate dynamics in bioreactors subjected to feast and famine cycles

Feast and famine cycles are common in activated sludge wastewater treatment systems, and they select for bacteria that accumulate storage compounds, such as poly-beta-hydroxybutyrate (PHB). Previous studies have shown that variations in influent substrate concentrations force bacteria to accumulate high levels of rRNA compared to the levels in bacteria grown in chemostats. Therefore, it can be hypothesized that bacteria accumulate more rRNA when they are subjected to feast and famine cycles. However, PHB-accumulating bacteria can form biomass (grow) throughout a feast and famine cycle and thus have a lower peak biomass formation rate during the cycle. Consequently, PHB-accumulating bacteria may accumulate less rRNA when they are subjected to feast and famine cycles than bacteria that are not capable of PHB accumulation. These hypotheses were tested with Wautersia eutropha H16 (wild type) and W. eutropha PHB-4 (a mutant not capable of accumulating PHB) grown in chemostat and semibatch reactors. For both strains, the cellular RNA level was higher when the organism was grown in semibatch reactors than when it was grown in chemostats, and the specific biomass formation rates during the feast phase were linearly related to the cellular RNA levels for cultures. Although the two strains exhibited maximum uptake rates when they were grown in semibatch reactors, the wild-type strain responded much more rapidly to the addition of fresh medium than the mutant responded. Furthermore, the chemostat-grown mutant culture was unable to exhibit maximum substrate uptake rates when it was subjected to pulse-wise addition of fresh medium. These data show that the ability to accumulate PHB does not prevent bacteria from accumulating high levels of rRNA when they are subjected to feast and famine cycles. Our results also demonstrate that the ability to accumulate PHB makes the bacteria more responsive to sudden increases in substrate concentrations, which explains their ecological advantage.

The selective cause of an ancient adaptation

Phylogenetic analysis reveals that the use of nicotinamide adenine dinucleotide phosphate (NADP) by prokaryotic isocitrate dehydrogenase (IDH) arose around the time eukaryotic mitochondria first appeared, about 3.5 billion years ago. We replaced the wild-type gene that encodes the NADP-dependent IDH of Escherichia coli with an engineered gene that possesses the ancestral NAD-dependent phenotype. The engineered enzyme is disfavored during competition for acetate. The selection intensifies in genetic backgrounds where other sources of reduced NADP have been removed. A survey of sequenced prokaryotic genomes reveals that those genomes that encode isocitrate lyase, which is essential for growth on acetate, always have an NADP-dependent IDH. Those with only an NAD-dependent IDH never have isocitrate lyase. Hence, the NADP dependence of prokaryotic IDH is an ancient adaptation to anabolic demand for reduced NADP during growth on acetate.

Effect of zwf gene knockout on the metabolism of Escherichia coli grown on glucose or acetate

The mutant deficient in glucose-6-phosphate dehydrogenase (G6PDH) was constructed by disrupting zwf gene by one-step inactivation protocol using polymerase chain reaction primers. The knockout of zwf gene was shown to have different influence on the metabolism of Escherichia coli grown on glucose or acetate. The decreased rates of substrate uptake and CO(2) production were found for the mutant grown on acetate, whereas these two rates were increased during the growth on glucose. The metabolic flux analysis based on (13)C-labeling experiments indicates that the metabolism of the mutant grown on glucose is related to the higher flux via tricorboxylic acid (TCA) cycle to generate anabolic reducing equivalents normally provided by the oxidative pentose phosphate pathway. However, the metabolism of the mutant grown on acetate shows a lower flux towards the TCA cycle as compared with the parent strain. The decreased flux through TCA cycle is associated with an increased flux via the glyoxylate shunt, by which the carbon source can bypass the two decarboxylative steps of TCA cycle in which CO(2) is released, thus conserving more carbon for biosynthesis in response to the decreased uptake rate of the carbon source.

A kinetic model of the branch-point between the methionine and threonine biosynthesis pathways in Arabidopsis thaliana

This work proposes a model of the metabolic branch-point between the methionine and threonine biosynthesis pathways in Arabidopsis thaliana which involves kinetic competition for phosphohomoserine between the allosteric enzyme threonine synthase and the two-substrate enzyme cystathionine gamma-synthase. Threonine synthase is activated by S-adenosylmethionine and inhibited by AMP. Cystathionine gamma-synthase condenses phosphohomoserine to cysteine via a ping-pong mechanism. Reactions are irreversible and inhibited by inorganic phosphate. The modelling procedure included an examination of the kinetic links, the determination of the operating conditions in chloroplasts and the establishment of a computer model using the enzyme rate equations. To test the model, the branch-point was reconstituted with purified enzymes. The computer model showed a partial agreement with the in vitro results. The model was subsequently improved and was then found consistent with flux partition in vitro and in vivo. Under near physiological conditions, S-adenosylmethionine, but not AMP, modulates the partition of a steady-state flux of phosphohomoserine. The computer model indicates a high sensitivity of cystathionine flux to enzyme and S-adenosylmethionine concentrations. Cystathionine flux is sensitive to modulation of threonine flux whereas the reverse is not true. The cystathionine gamma-synthase kinetic mechanism favours a low sensitivity of the fluxes to cysteine. Though sensitivity to inorganic phosphate is low, its concentration conditions the dynamics of the system. Threonine synthase and cystathionine gamma-synthase display similar kinetic efficiencies in the metabolic context considered and are first-order for the phosphohomoserine substrate. Under these conditions outflows are coordinated.

A novel metabolic cycle catalyzes glucose oxidation and anaplerosis in hungry Escherichia coli

Complete oxidation of carbohydrates to CO2 is considered to be the exclusive property of the ubiquitous tricarboxylic acid cycle, the central process in cellular energy metabolism of aerobic organisms. Based on metabolism-wide in vivo quantification of intracellular carbon fluxes, we describe here complete oxidation of carbohydrates via the novel P-enolpyruvate (PEP)-glyoxylate cycle, in which two PEP molecules are oxidized by means of acetyl coenzyme A, citrate, glyoxylate, and oxaloacetate to CO2, and one PEP is regenerated. Key reactions are the constituents of the glyoxylate shunt and PEP carboxykinase, whose conjoint operation in this bi-functional catabolic and anabolic cycle is in sharp contrast to their generally recognized functions in anaplerosis and gluconeogenesis, respectively. Parallel operation of the PEP-glyoxylate cycle and the tricarboxylic acid cycle was identified in the bacterium Escherichia coli under conditions of glucose hunger in a slow-growing continuous culture. Because the PEP-glyoxylate cycle was also active in glucose excess batch cultures of an NADPH-overproducing phosphoglucose isomerase mutant, one function of this new central pathway may be the decoupling of catabolism from NADPH formation that would otherwise occur in the tricarboxylic acid cycle.

Analysis of Escherichia coli anaplerotic metabolism and its regulation mechanisms from the metabolic responses to altered dilution rates and phosphoenolpyruvate carboxykinase knockout

The gluconeogenic phosphoenolpyruvate (PEP) carboxykinase is active in Escherichia coli during its growth on glucose. The present study investigated the influence of growth rates and PEP carboxykinase knockout on the anaplerotic fluxes in E. coli. The intracellular fluxes were determined using the complementary methods of flux ratio analysis and metabolic flux analysis based on [U-(13)C(6)]glucose labeling experiments and 2D nuclear magnetic resonance (NMR) spectroscopy of cellular amino acids and glycerol. Significant activity of PEP carboxykinase was identified in wild-type E. coli, and the ATP dissipation for the futile cycling via this reaction accounted for up to 8.2% of the total energy flux. Flux analysis of pck deletion mutant revealed that abolishment of PEP carboxykinase activity resulted in a remarkably reduced flux through the anaplerotic PEP carboxylase and the activation of the glyoxylate shunt, with 23% of isocitrate found being channeled in the glyoxylate shunt. The changes in intracellular metabolite concentrations and specific enzyme activities associated with different growth rates and pck deletion, were also determined. Combining the measurement data of in vivo fluxes, metabolite concentrations and enzyme activities, the in vivo regulations of PEP carboxykinase flux, PEP carboxylation, and glyoxylate shunt in E. coli are discussed.

Two phenotypically compensating isocitrate dehydrogenases inRalstonia eutropha

The tricarboxylic acid (TCA) cycle enzyme isocitrate dehydrogenase (IDH) and the glyoxylate bypass enzyme isocitrate lyase are involved in catabolism of isocitrate and play a key role in controlling the metabolic flux between the TCA cycle and the glyoxylate shunt. Two IDH genes icd1 and icd2 of Ralstonia eutropha HF39, encoding isocitrate dehydrogenase 1 (IDH1) and isocitrate dehydrogenase 2 (IDH2), were identified and characterized. Icd1 was functionally expressed in Escherichia coli, whereas icd2 was expressed in E. coli but no activity was obtained. Interposon-mutants of icd1 (HF39Deltaicd1) and icd2 (HF39Deltaicd2) of R. eutropha HF39 were constructed and their phenotypes were investigated. HF39Deltaicd1 retained 43% of IDH activity, which was not induced by acetate, and HF39Deltaicd2 expressed 74% of acetate-induced IDH activity. Both HF39Deltaicd1and HF39Deltaicd2 kept the same growth rate on acetate as the wild-type. These data suggested that IDH1 is induced by acetate. The interposon-mutants HF39Deltaicd1 and HF39Deltaicd2 accumulated the same amount of poly(3-hydroxybutyric acid) as the wild-type.