Glucose catabolism of Escherichia coli strains with increased activity and altered regulation of key glycolytic enzymes

The Protein Interactome of Glycolysis in Escherichia coli

Glycolysis is regulated by numerous mechanisms including allosteric regulation, post-translational modification or protein-protein interactions (PPI). While glycolytic enzymes have been found to interact with hundreds of proteins, the impact of only some of these PPIs on glycolysis is well understood. Here we investigate which of these interactions may affect glycolysis in E. coli and possibly across numerous other bacteria, based on the stoichiometry of interacting protein pairs (from proteomic studies) and their conservation across bacteria. We present a list of 339 protein-protein interactions involving glycolytic enzymes but predict that ~70% of glycolytic interactors are not present in adequate amounts to have a significant impact on glycolysis. Finally, we identify a conserved but uncharacterized subset of interactions that are likely to affect glycolysis and deserve further study.

Identification of a rate-limiting step in a metabolic pathway using the kinetic model and in vitro experiment

Identification of the rate-limiting step in a metabolic pathway is an important challenge in metabolic engineering for enhancing pathway flow. Although specific enzyme activities (V_max) provide valuable clues for the identification, it is time-consuming and difficult to measure multiple enzymes in the pathway because different assay protocols are required for each enzyme. In the present study, we propose a method to simultaneously determine the V_max values of multiple enzymes using a kinetic model with a time course of the intermediate concentrations through an in vitro experiment. To demonstrate this method, nine glycolysis reactions for converting glucose-6-phosphate (G6P) to pyruvate in Escherichia coli were considered. In a reaction mixture containing G6P and cofactors, glycolysis was initiated by adding a crude cell extract obtained from stationary phase cells. The V_max values were optimized to minimize the difference between the measured and simulated time-courses using a kinetic model. Metabolic control analysis using the kinetic model with the estimated V_max values revealed that fructose bisphosphate aldolase (FBA) was the rate-limiting step in the upper part of glycolysis. The addition of FBA in the reaction mixture successfully increased the glycolytic flux in vitro. Furthermore, in vivo, the specific glucose consumption rate of an FBA overexpression strain was 1.4 times higher than that of the control strain during the stationary phase. These results confirmed that FBA was the rate-limiting step in glycolysis under the stationary phase. This approach provides V_max values of multiple enzymes in a pathway for metabolic control analysis with a kinetic model.

Synthesis, Molecular Docking Studies and Antibacterial Activities of Novel Monocationic Indole-benzimidazole Derivatives

BACKGROUND

Finding efficient therapy against hospital-acquired MRSA infections has become rather important in the last decade. To this end, inhibition of the enzyme pyruvate kinase (PK) is being investigated for antibacterial activity, since this enzyme controls energy generation and metabolic flux distribution. Our main scaffold consists of benzimidazole and indole rings fused together. Both rings are famous for antibacterial properties and promising anti-MRSA compounds including indole ring.

METHODS

Several 1-substituted-2-(1H-indol-3-yl)-N-substituted-1H-benzimidazole-5-carboxamidine analogues were developed, synthesized and their antibacterial activities were evaluated against Staphylococcus aureus (ATCC 25923), Methicillin resistant Staphylococcus aureus (MRSA) (ATCC 43300), and Staphylococcus epidermidis (ATCC 12228) by using tube dilution method. Molecular docking analysis with a characteristic protein called MRSA- Pyruvate Kinase has been conducted for the assessment of the activities of our compounds against Methicillinresistant S. aureus (MRSA).

RESULTS

Among all the tested compounds, the most potent compound 36 had MIC values as 3.12, 3.12 and 6.25 μg/mL against S. aureus, Methicillin-resistant S. aureus (MRSA), and S. epidermidis, respectively. This compound had much better docking energy value than standard ampicillin and also created the link between two residues in different monomers of PK.

DISCUSSION

This approach of using indol-amidine conjugate systems as anti-MRSA agents may include MRSA-PK as potential target. To further increase the affinity, some other H-bonding parts may be added. By doing so, another bridge with Ile361 residues on both sides can be created. Our compounds tend to violate log P limit of Lipinski, therefore some optimizations with formulation can be made.

CONCLUSION

This study mainly includes the design, synthesis and optimization of indolebenzimidazole- amidine derivatives. Docking studies confirmed our results, since our most potent hit compound 36 created the necessary interactions between two chains of MRSA-PK. Further optimization can be considered to increase drug ability.

Coupling xylitol dehydrogenase with NADH oxidase improves l-xylulose production in Escherichia coli culture

Escherichia coli expressing NAD-dependent xylitol-4-dehydrogenase (XDH) from Pantoea ananatis and growing on glucose or glycerol converts xylitol to the rare sugar l-xylulose. Although blocking potential l-xylulose consumption (l-xylulosekinase, lyxK) or co-expression of the glycerol facilitator (glpF) did not significantly affect l-xylulose formation, co-expressing XDH with water-forming NADH oxidase (NOX) from Streptococcus pneumoniae increased l-xylulose formation in shake flasks when glycerol was the carbon source. Controlled batch processes at the 1L scale demonstrated that the final equilibrium l-xylulose/xylitol ratio was correlated to the intracellular NAD⁺/NADH ratio, with 69% conversion of xylitol to l-xylulose and a yield of 0.88g l-xylulose/g xylitol consumed attained for MG1655/pZE12-xdh/pCS27-nox growing on glycerol. NADH oxidase was less effective at improving l-xylulose formation in the bioreactor than in shake flasks, likely as a result of an intrinsic maximum NAD⁺/NADH and l-xylulose/xylitol equilibrium ratio being attained. Intermittently feeding carbon source was ineffective at increasing the final l-xylulose concentration because introduction of carbon source was accompanied by a reduction in NAD⁺/NADH ratio. A batch process using 12g/L glycerol and 22g/L xylitol generated over 14g/L l-xylulose after 80h, corresponding to 65% conversion and a yield of 0.89g l-xylulose/g xylitol consumed.

Refactoring the Embden-Meyerhof-Parnas Pathway as a Whole of Portable GlucoBricks for Implantation of Glycolytic Modules in Gram-Negative Bacteria

The Embden–Meyerhof–Parnas (EMP) pathway is generally considered to be the biochemical standard for glucose catabolism. Alas, its native genomic organization and the control of gene expression in Escherichia coli are both very intricate, which limits the portability of the EMP pathway to other biotechnologically important bacterial hosts that lack the route. In this work, the genes encoding all the enzymes of the linear EMP route have been individually recruited from the genome of E. coli K-12, edited in silico to remove their endogenous regulatory signals, and synthesized de novo following a standard (GlucoBrick) that enables their grouping in the form of functional modules at the user’s will. After verifying their activity in several glycolytic mutants of E. coli, the versatility of these GlucoBricks was demonstrated in quantitative physiology tests and biochemical assays carried out in Pseudomonas putida KT2440 and P. aeruginosa PAO1 as the heterologous hosts. Specific configurations of GlucoBricks were also adopted to streamline the downward circulation of carbon from hexoses to pyruvate in E. coli recombinants, thereby resulting in a 3-fold increase of poly(3-hydroxybutyrate) synthesis from glucose. Refactoring whole metabolic blocks in the fashion described in this work thus eases the engineering of biochemical processes where the optimization of carbon traffic is facilitated by the operation of the EMP pathway—which yields more ATP than other glycolytic routes such as the Entner–Doudoroff pathway.

Functional role of pyruvate kinase from Lactobacillus bulgaricus in acid tolerance and identification of its transcription factor by bacterial one-hybrid

Lactobacillus delbrueckii subsp. bulgaricus develops acid tolerance response when subjected to acid stress conditions, such as the induction of enzymes associated with carbohydrate metabolism. In this study, pyk gene encoding pyruvate kinase was over-expressed in heterologous host Lactococcus lactis NZ9000, and SDS-PAGE analysis revealed the successful expression of this gene in NZ9000. The survival rate of Pyk-overproducing strain was 45-fold higher than the control under acid stress condition (pH 4.0). In order to determine the transcription factor (TF) which regulates the expression of pyk by bacterial one-hybrid, we constructed a TF library including 65 TFs of L. bulgaricus. Western blotting indicated that TFs in this library could be successfully expressed in host strains. Subsequently, the promoter of pfk-pyk operon in L. bulgaricus was identified by 5′-RACE PCR. The bait plasmid pH3U3-p01 carrying the deletion fragment of pfk-pyk promoter captured catabolite control protein A (CcpA) which could regulate the expression of pyk by binding to a putative catabolite-responsive element (5′-TGTAAGCCCTAACA-3′) upstream the -35 region. Real-time qPCR analysis revealed the transcription of pyk was positively regulated by CcpA. This is the first report about identifying the TF of pyk in L. bulgaricus, which will provide new insight into the regulatory network.

Dysregulation of glucose transport, glycolysis, TCA cycle and glutaminolysis by oncogenes and tumor suppressors in cancer cells

A common set of functional characteristics of cancer cells is that cancer cells consume a large amount of glucose, maintain high rate of glycolysis and convert a majority of glucose into lactic acid even in the presence of oxygen compared to that of normal cells (Warburg's Effects). In addition, cancer cells exhibit substantial alterations in several energy metabolism pathways including glucose transport, tricarboxylic acid (TCA) cycle, glutaminolysis, mitochondrial respiratory chain oxidative phosphorylation and pentose phosphate pathway (PPP). In the present work, we focused on reviewing the current knowledge about the dysregulation of the proteins/enzymes involved in the key regulatory steps of glucose transport, glycolysis, TCA cycle and glutaminolysis by several oncogenes including c-Myc and hypoxia inducible factor-1 (HIF-1) and tumor suppressor, p53, in cancer cells. The dysregulation of glucose transport and energy metabolism pathways by oncogenes and lost functions of the tumor suppressors have been implicated as important biomarkers for cancer detection and as valuable targets for the development of new anticancer therapies.

Posttranslational modification of 6-phosphofructo-1-kinase as an important feature of cancer metabolism

Background

Human cancers consume larger amounts of glucose compared to normal tissues with most being converted and excreted as lactate despite abundant oxygen availability (Warburg effect). The underlying higher rate of glycolysis is therefore at the root of tumor formation and growth. Normal control of glycolytic allosteric enzymes appears impaired in tumors; however, the phenomenon has not been fully resolved.

Methodology/Principal Findings

In the present paper, we show evidence that the native 85-kDa 6-phosphofructo-1-kinase (PFK1), a key regulatory enzyme of glycolysis that is normally under the control of feedback inhibition, undergoes posttranslational modification. After proteolytic cleavage of the C-terminal portion of the enzyme, an active, shorter 47-kDa fragment was formed that was insensitive to citrate and ATP inhibition. In tumorigenic cell lines, only the short fragments but not the native 85-kDa PFK1 were detected by immunoblotting. Similar fragments were detected also in a tumor tissue that developed in mice after the subcutaneous infection with tumorigenic B16-F10 cells. Based on limited proteolytic digestion of the rabbit muscle PFK-M, an active citrate inhibition-resistant shorter form was obtained, indicating that a single posttranslational modification step was possible. The exact molecular masses of the active shorter PFK1 fragments were determined by inserting the truncated genes constructed from human muscle PFK1 cDNA into a pfk null E. coli strain. Two E. coli transformants encoding for the modified PFK1s of 45,551 Da and 47,835 Da grew in glucose medium. The insertion of modified truncated human pfkM genes also stimulated glucose consumption and lactate excretion in stable transfectants of non-tumorigenic human HEK cell, suggesting the important role of shorter PFK1 fragments in enhancing glycolytic flux.

Conclusions/Significance

Posttranslational modification of PFK1 enzyme might be the pivotal factor of deregulated glycolytic flux in tumors that in combination with altered signaling mechanisms essentially supports fast proliferation of cancer cells.

Switch between life history strategies due to changes in glycolytic enzyme gene dosage in Saccharomyces cerevisiae

Adaptation is the process whereby a population or species becomes better fitted to its habitat through modifications of various life history traits which can be positively or negatively correlated. The molecular factors underlying these covariations remain to be elucidated. Using Saccharomyces cerevisiae as a model system, we have investigated the effects on life history traits of varying the dosage of genes involved in the transformation of resources into energy. Changing gene dosage for each of three glycolytic enzyme genes (hexokinase 2, phosphoglucose isomerase, and fructose-1,6-bisphosphate aldolase) resulted in variation in enzyme activities, glucose consumption rate, and life history traits (growth rate, carrying capacity, and cell size). However, the range of effects depended on which enzyme was expressed differently. Most interestingly, these changes revealed a genetic trade-off between carrying capacity and cell size, supporting the discovery of two extreme life history strategies already described in yeast populations: the "ants," which have lower glycolytic gene dosage, take up glucose slowly, and have a small cell size but reach a high carrying capacity, and the "grasshoppers," which have higher glycolytic gene dosage, consume glucose more rapidly, and allocate it to a larger cell size but reach a lower carrying capacity. These results demonstrate antagonist pleiotropy for glycolytic genes and show that altered dosage of a single gene drives a switch between two life history strategies in yeast.

Redox biocatalysis and metabolism: molecular mechanisms and metabolic network analysis

Whole-cell biocatalysis utilizes native or recombinant enzymes produced by cellular metabolism to perform synthetically interesting reactions. Besides hydrolases, oxidoreductases represent the most applied enzyme class in industry. Oxidoreductases are attributed a high future potential, especially for applications in the chemical and pharmaceutical industries, as they enable highly interesting chemistry (e.g., the selective oxyfunctionalization of unactivated C-H bonds). Redox reactions are characterized by electron transfer steps that often depend on redox cofactors as additional substrates. Their regeneration typically is accomplished via the metabolism of whole-cell catalysts. Traditionally, studies towards productive redox biocatalysis focused on the biocatalytic enzyme, its activity, selectivity, and specificity, and several successful examples of such processes are running commercially. However, redox cofactor regeneration by host metabolism was hardly considered for the optimization of biocatalytic rate, yield, and/or titer. This article reviews molecular mechanisms of oxidoreductases with synthetic potential and the host redox metabolism that fuels biocatalytic reactions with redox equivalents. The tools discussed in this review for investigating redox metabolism provide the basis for studies aiming at a deeper understanding of the interplay between synthetically active enzymes and metabolic networks. The ultimate goal of rational whole-cell biocatalyst engineering and use for fine chemical production is discussed.

Enhanced citric acid biosynthesis in Pseudomonas fluorescens ATCC 13525 by overexpression of the Escherichia coli citrate synthase gene

Citric acid secretion by fluorescent pseudomonads has a distinct significance in microbial phosphate solubilization. The role of citrate synthase in citric acid biosynthesis and glucose catabolism in pseudomonads was investigated by overexpressing the Escherichia coli citrate synthase (gltA) gene in Pseudomonas fluorescens ATCC 13525. The resultant approximately 2-fold increase in citrate synthase activity in the gltA-overexpressing strain Pf(pAB7) enhanced the intracellular and extracellular citric acid yields during the stationary phase, by about 2- and 26-fold, respectively, as compared to the control, without affecting the growth rate, glucose depletion rate or biomass yield. Decreased glucose consumption was paralleled by increased gluconic acid production due to an increase in glucose dehydrogenase activity. While the extracellular acetic acid yield increased in Pf(pAB7), pyruvic acid secretion decreased, correlating with an increase in pyruvate carboxylase activity and suggesting an increased demand for the anabolic precursor oxaloacetate. Activities of two other key enzymes, glucose-6-phosphate dehydrogenase and isocitrate dehydrogenase, remained unaltered, and the contribution of phosphoenolpyruvate carboxylase and isocitrate lyase to glucose catabolism was negligible. Strain Pf(pAB7) demonstrated an enhanced phosphate-solubilizing ability compared to the control. Co-expression of the Synechococcus elongatus PCC 6301 phosphoenolpyruvate carboxylase and E. coli gltA genes in P. fluorescens ATCC 13525, so as to supplement oxaloacetate for citrate biosynthesis, neither significantly affected citrate biosynthesis nor caused any change in the other physiological and biochemical parameters measured, despite approximately 1.3- and 5-fold increases in citrate synthase and phosphoenolpyruvate carboxylase activities, respectively. Thus, our results demonstrate that citrate synthase is rate-limiting in enhancing citrate biosynthesis in P. fluorescens ATCC 13525. Significantly low extracellular citrate levels as compared to the intracellular levels in Pf(pAB7) suggested a probable limitation of efficient citrate transport.

Metabolic control analysis: a tool for designing strategies to manipulate metabolic pathways

The traditional experimental approaches used for changing the flux or the concentration of a particular metabolite of a metabolic pathway have been mostly based on the inhibition or over-expression of the presumed rate-limiting step. However, the attempts to manipulate a metabolic pathway by following such approach have proved to be unsuccessful. Metabolic Control Analysis (MCA) establishes how to determine, quantitatively, the degree of control that a given enzyme exerts on flux and on the concentration of metabolites, thus substituting the intuitive, qualitative concept of rate limiting step. Moreover, MCA helps to understand (i) the underlying mechanisms by which a given enzyme exerts high or low control and (ii) why the control of the pathway is shared by several pathway enzymes and transporters. By applying MCA it is possible to identify the steps that should be modified to achieve a successful alteration of flux or metabolite concentration in pathways of biotechnological (e.g., large scale metabolite production) or clinical relevance (e.g., drug therapy). The different MCA experimental approaches developed for the determination of the flux-control distribution in several pathways are described. Full understanding of the pathway properties when is working under a variety of conditions can help to attain a successful manipulation of flux and metabolite concentration.

Specific ethanol production rate in ethanologenic Escherichia coli strain KO11 Is limited by pyruvate decarboxylase

Modification of ethanol productivity and yield, using mineral medium supplemented with glucose or xylose as carbon sources, was studied in ethanologenic Escherichia coli KO11 by increasing the activity of five key carbon metabolism enzymes. KO11 efficiently converted glucose or xylose to ethanol with a yield close to 100% of the theoretical maximum when growing in rich medium. However, when KO11 ferments glucose or xylose in mineral medium, the ethanol yields decreased to only 70 and 60%, respectively. An increase in GALP(Ec) (permease of galactose-glucose-xylose) or PGK(Ec) (phosphoglycerate kinase) activities did not change xylose or glucose and ethanol flux. However, when PDC(Zm) (pyruvate decarboxylase from Zymomonas mobilis) activity was increased 7-fold, the yields of ethanol from glucose or xylose were increased to 85 and 75%, respectively, and organic acid formation rates were reduced. Furthermore, as a response to a reduction in acetate and ATP yield, and a limited PDC(Zm) activity, an increase in PFK(Ec) (phosphofructokinase) or PYK(Bs) (pyruvate kinase from Bacillus stearothermophilus) activity drastically reduced glucose or xylose consumption and ethanol formation flux. This experimental metabolic control analysis showed that ethanol flux in KO11 is negatively controlled by phosphofructokinase and pyruvate kinase, and positively influenced by the PDC(Zm) activity level.

Coutilization of glucose and glycerol enhances the production of aromatic compounds in an Escherichia coli strain lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system

Background

Escherichia coli strains lacking the phosphoenolpyruvate: carbohydrate phosphotransferase system (PTS) are capable of coutilizing glucose and other carbon sources due to the absence of catabolite repression by glucose. In these strains, the lack of this important regulatory and transport system allows the coexistence of glycolytic and gluconeogenic pathways. Strains lacking PTS have been constructed with the goal of canalizing part of the phosphoenolpyruvate (PEP) not consumed in glucose transport to the aromatic pathway. The deletion of the ptsHIcrr operon inactivates PTS causing poor growth on this sugar; nonetheless, fast growing mutants on glucose have been isolated (PB12 strain). However, there are no reported studies concerning the growth potential of a PTS^- strain in mixtures of different carbon sources to enhance the production of aromatics compounds.

Results

PB12 strain is capable of coutilizing mixtures of glucose-arabinose, glucose-gluconate and glucose-glycerol. This capacity increases its specific growth rate (μ) given that this strain metabolizes more moles of carbon source per unit time. The presence of plasmids pRW300aroG^fbr and pCLtktA reduces the μ of strain PB12 in all mixtures of carbon sources, but enhances the productivity and yield of aromatic compounds, especially in the glucose-glycerol mixture, as compared to glucose or glycerol cultures. No acetate was detected in the glycerol and the glucose-glycerol batch fermentations.

Conclusion

Due to the lack of catabolite repression, PB12 strain carrying multicopy plasmids containing tktA and aroG^fbr genes is capable of coutilizing glucose and other carbon sources; this capacity, reduces its μ but increases the production of aromatic compounds.

Minimizing acetate formation in E. coli fermentations

Escherichia coli remains the best-established production organism in industrial biotechnology. However, when aerobic fermentation runs at high growth rates, considerable amounts of acetate are accumulated as by-product. This by-product has negative effects on growth and protein production. Over the last 20 years, substantial research efforts have been expended on reducing acetate accumulation during aerobic growth of E. coli on glucose. From the onset it was clear that this quest would not be a simple or uncomplicated one. Simple deletion of the acetate pathway reduced the acetate accumulation, but other by-products were formed. This mini review gives a clear outline of these research efforts and their outcome, including bioprocess level approaches and genetic approaches. Recently, the latter seems to have some promising results.

Increased production of pyruvic acid by Escherichia coli RNase G mutants in combination with cra mutations

The Escherichia coli RNase G is known as an endoribonuclease responsible for the 5'-end maturation of 16S rRNA and degradation of several specific mRNAs such as adhE and eno mRNAs. In this study, we found that an RNase G mutant derived from the MC1061 strain did not grow on a glucose minimal medium. Genetic analysis revealed that simultaneous defects of cra and ilvIH, encoding a transcriptional regulator of glycolysis/gluconeogenesis and one of isozymes of acetohydroxy acid synthase, respectively, were required for this phenomenon to occur. The results of additional experiments presented here indicate that the RNase G mutation, in combination with cra mutation, caused the increased production of pyruvic acid from glucose, which was then preferentially converted to valine due to the ilvIH mutation, resulting in depletion of isoleucine. In fact, the rng cra double mutant produced increased amount of pyruvate in the medium. These results suggest that the RNase G mutation could be applied in the breeding of producer strains of pyruvate and its derivatives such as valine.

Metabolic profiling of Escherichia coli cultivations: evaluation of extraction and metabolite analysis procedures

The quantitative estimation of intracellular metabolite concentrations (metabolic profiling) is a prerequisite for a better understanding of biological processes and thus inevitable for the rational improvement of microbial production strains and process design. Since pool sizes of substrates regulate flux through different enzymes, the accurate determination of intracellular metabolite concentrations is necessary to understand in vivo reaction kinetics. Quantification of intracellular concentrations of glycolytic intermediates in Escherichia coli K12 was achieved by using a novel in situ rapid sampling and quenching procedure. A new extraction procedure using buffered hot water was established. By use of simultaneous multi-substrate feeding with various ratios of glucose, fructose and acetate during continuous cultivations several metabolic states were induced. Metabolic flux analysis and the newly developed metabolic profiling procedure were used to determine in vivo enzyme kinetics as exemplified for fructose 1,6-bisphosphate aldolase and citrate synthase.

Metabolic flux analysis at ultra short time scale: Isotopically non-stationary 13C labeling experiments

A novel approach to (13)C metabolic flux analysis (MFA) is presented using cytosolic metabolite pool sizes and their (13)C labeling data from an isotopically non-stationary (13)C labeling experiment (INST-CLE). The procedure is demonstrated with an E. coli wild type strain grown at fed batch conditions. The intra cellular labeling dynamics are excited by a sudden step increase of the (13)C portion in the substrate feed. Due to unchanged saturation of the substrate uptake system, the metabolic fluxes remain constant during the following sampling time period of only 16s, in which 20 samples are taken by an automated rapid sampling device immediately stopping metabolism by methanol quenching. Subsequent cell disruptive sample preparation and LC-MS/MS enabled simultaneous determination of pool sizes and mass isotopomers of intra cellular metabolites requiring detection limits in the nM range. Based on this data the new computational flux analysis tool 13CFLUX/INST is used to determine the intra cellular fluxes based on a complex carbon labeling network model. The measured data is in good agreement with the model predictions, thus proving the applicability of the new isotopically non-stationary (13)C metabolic flux analysis (INST-(13)C-MFA) concept. Moreover, it is shown that significant new information with respect to flux identifiability, non-measurable pool sizes, data consistency, or large storage pools can be taken from the novel kind of experimental data. This offers new insight into the biological operation of the metabolic network in vivo.

Putative regulatory sites unraveled by network-embedded thermodynamic analysis of metabolome data

As one of the most recent members of the omics family, large-scale quantitative metabolomics data are currently complementing our systems biology data pool and offer the chance to integrate the metabolite level into the functional analysis of cellular networks. Network-embedded thermodynamic analysis (NET analysis) is presented as a framework for mechanistic and model-based analysis of these data. By coupling the data to an operating metabolic network via the second law of thermodynamics and the metabolites' Gibbs energies of formation, NET analysis allows inferring functional principles from quantitative metabolite data; for example it identifies reactions that are subject to active allosteric or genetic regulation as exemplified with quantitative metabolite data from Escherichia coli and Saccharomyces cerevisiae. Moreover, the optimization framework of NET analysis was demonstrated to be a valuable tool to systematically investigate data sets for consistency, for the extension of sub-omic metabolome data sets and for resolving intracompartmental concentrations from cell-averaged metabolome data. Without requiring any kind of kinetic modeling, NET analysis represents a perfectly scalable and unbiased approach to uncover insights from quantitative metabolome data.

Glycolysis and Flux Control

Central metabolism of carbohydrates uses the Embden-Meyerhof-Parnas (EMP), pentose phosphate (PP), and Entner-Doudoroff (ED) pathways. This review reviews the biological roles of the enzymes and genes of these three pathways of E. coli. Glucose, pentoses, and gluconate are primarily discussed as the initial substrates of the three pathways, respectively. The genetic and allosteric regulatory mechanisms of glycolysis and the factors that affect metabolic flux through the pathways are considered here. Despite the fact that a lot of information on each of the reaction steps has been accumulated over the years for E. coli, surprisingly little quantitative information has been integrated to analyze glycolysis as a system. Therefore, the review presents a detailed description of each of the catalytic steps by a systemic approach. It considers both structural and kinetic aspects. Models that include kinetic information of the reaction steps will always contain the reaction stoichiometry and therefore follow the structural constraints, but in addition to these also kinetic rate laws must be fulfilled. The kinetic information obtained on isolated enzymes can be integrated using computer models to simulate behavior of the reaction network formed by these enzymes. Successful examples of such approaches are the modeling of glycolysis in S. cerevisiae, the parasite Trypanosoma brucei, and the red blood cell. With the rapid developments in the field of Systems Biology many new methods have been and will be developed, for experimental and theoretical approaches, and the authors expect that these will be applied to E. coli glycolysis in the near future.

A sensitive and robust method for quantification of intracellular short-chain coenzyme A esters

A procedure for the analysis of short-chain intracellular coenzyme A (CoA) esters and adenine nucleotide pools in microbial cells is described. The simultaneous isolation of bacterial cells from media, quenching of their metabolism, and extraction of metabolites was accomplished by centrifugation of cells through a layer of silicone oil into a denser solution of trichloroacetic acid. The acid was neutralized by extraction into Freon containing tri-n-octylamine to provide a salt-free solution of cell metabolites. After high-performance liquid chromatography separation, CoA, CoA esters, and adenine-containing nucleotides were derivatized by postcolumn reaction with bromoacetaldehyde to form the fluorescent 1,N6-ethenoadenine adducts which were analyzed by a fluorescence detector at picomolar levels.

Effect of pyruvate kinase overproduction on glucose metabolism of Lactococcus lactis

Lactococcus lactis strain NZ9000(pNZpyk), which overproduces pyruvate kinase (PK), was constructed. The pNZpyk plasmid carries the P(nisA)-pyk transcriptional fusion, and the overexpression of its pyk gene was accomplished by using the nisin-inducible expression system of the NZ9000 strain. In vivo (13)C- and (31)P-NMR spectroscopy was used to evaluate the effect of this modification on the metabolism of glucose in non-growing cells. A detailed description of the kinetics of glucose, end products, glycolytic intermediates, NAD(+) and NADH was obtained. A 15-fold increase in the level of PK did not increase the overall glycolytic flux, which, on the contrary, was slightly reduced. Significant differences were observed in (i) the level of 3-phosphoglycerate (3-PGA) and phosphoenolpyruvate (PEP), metabolites associated with starvation; (ii) the rate of fructose 1,6-bisphosphate (FBP) depletion upon glucose exhaustion; and (iii) the NAD(+)/NADH ratio during glucose catabolism. In the mutant, the rate of FBP consumption after glucose depletion was notably accelerated under anaerobic conditions, whereas 3-PGA and PEP decreased to undetectable levels. Furthermore, the level of NAD(+) decreased steadily during the utilization of glucose, probably due to the unanticipated reduction in the lactate dehydrogenase activity in comparison with the control strain, NZ9000(pNZ8020). The results show that PK is an important bottleneck to carbon flux only when glucose becomes limiting; in the overproducer this constriction was no longer present, as evidenced by the faster FBP consumption and lack of accumulation of 3-PGA and PEP in anaerobic as well as aerobic conditions. Despite these clear changes, the PK-overproducing strain showed typical homolactic metabolism under anaerobic conditions, as did the strain harbouring the vector plasmid without the pyk insert. However, under an oxygen atmosphere, there was increased channelling of carbon to the production of acetate and acetoin, to the detriment of lactate production.

The soluble and membrane-bound transhydrogenases UdhA and PntAB have divergent functions in NADPH metabolism of Escherichia coli

Pentose phosphate pathway and isocitrate dehydrogenase are generally considered to be the major sources of the anabolic reductant NADPH. As one of very few microbes, Escherichia coli contains two transhydrogenase isoforms with unknown physiological function that could potentially transfer electrons directly from NADH to NADP+ and vice versa. Using defined mutants and metabolic flux analysis, we identified the proton-translocating transhydrogenase PntAB as a major source of NADPH in E. coli. During standard aerobic batch growth on glucose, 35-45% of the NADPH that is required for biosynthesis was produced via PntAB, whereas pentose phosphate pathway and isocitrate dehydrogenase contributed 35-45% and 20-25%, respectively. The energy-independent transhydrogenase UdhA, in contrast, was essential for growth under metabolic conditions with excess NADPH formation, i.e. growth on acetate or in a phosphoglucose isomerase mutant that catabolized glucose through the pentose phosphate pathway. Thus, both isoforms have divergent physiological functions: energy-dependent reduction of NADP+ with NADH by PntAB and reoxidation of NADPH by UdhA. Expression appeared to be modulated by the redox state of cellular metabolism, because genetic and environmental manipulations that increased or decreased NADPH formation down-regulated pntA or udhA transcription, respectively. The two transhydrogenase isoforms provide E. coli primary metabolism with an extraordinary flexibility to cope with varying catabolic and anabolic demands, which raises two general questions: why do only a few bacteria contain both isoforms, and how do other organisms manage NADPH metabolism?

United States Department of Agriculture-Agricultural Research Service research programs in biological control of plant diseases

A number of USDA-ARS programs directed at overcoming impediments to the use of biocontrol agents on a commercial scale are described. These include improvements in screening techniques, taxonomic studies to identify beneficial strains more precisely, and studies on various aspects of the large-scale production of biocontrol agents. Another broad area of studies covers the ecological aspects of biocontrol agents-their interaction with the pathogen, with the plant and with other aspects of the environmental complex. Examples of these studies are given and their relevance to the further development and expansion of biocontrol agents is discussed.

Pyruvate kinase: current status of regulatory and functional properties

Pyruvate kinase (PK) is a key enzyme for the glycolytic pathway and carbon metabolism in general. On the basis of the relevance and enormous diverse properties of this enzyme, this paper describes the results of a current and extensive review that determines the sites of conservation and/or difference in PK sequences, and the differences in the functional and regulatory properties of the enzymes. An alignment and analysis of 50 PK sequences from different sources and a phylogenetic tree are presented. This analysis was performed with reference to crystallographically characterized PK principally from E. coli, cat and rabbit muscle. A number of attributes of the enzyme that make it of particular interest in biomedicine and industry are also discussed.

Mutation of rpiA in Enterobacter cloacae decreases seed and root colonization and biocontrol of damping-off caused by Pythium ultimum on cucumber

Strains of Enterobacter cloacae show promise as biocontrol agents for Pythium ultimum-induced damping-off on cucumber and other crops. E. cloacae A145 is a mini-Tn5 Km transposon mutant of strain 501R3 that was significantly reduced in suppression of damping-off on cucumber caused by P. ultimum. Strain A145 was deficient in colonization of cucumber, sunflower, and wheat seeds and significantly reduced in colonization of corn and cowpea seeds relative to strain 501R3. Populations of strain A145 were also significantly lower than those of strain 501R3 at all sampling times in cucumber, wheat, and sunflower rhizosphere. Populations of strain A145 were not detectable in any rhizosphere after 42 days, while populations of strain 501R3 remained at substantial levels throughout all experiments. Molecular characterization of strain A145 indicated mini-Tn5 Km was inserted in a region of the E. cloacae genome with a high degree of DNA and amino acid sequence similarity to rpiA, which encodes ribose-5-phosphate isomerase. In Escherichia coli, RpiA catalyzes the interconversion of ribose-5-phosphate and ribulose-5-phosphate and is a key enzyme in the pentose phosphate pathway. Ribose-5-phosphate isomerase activity in cell lysates from strain A145 was approximately 3.5% of that from strain 501R3. In addition, strain A145 was a ribose auxotroph, as expected for an rpiA mutant. Introduction of a 1.0-kb DNA fragment containing only the rpiA homologue into strain A145 restored ribose phosphate isomerase activity, prototrophy, seedling colonization, and disease suppression to levels similar to those associated with strain 501R3. Experiments reported here indicate a key role for rpiA and possibly the pentose phosphate pathway in suppression of damping-off and colonization of subterranean portions of plants by E. cloacae.

Evolutionary engineering of industrially important microbial phenotypes

The tremendous complexity of dynamic interactions in cellular systems often impedes practical applications of metabolic engineering that are largely based on available molecular or functional knowledge. In contrast, evolutionary engineering follows nature's 'engineering' principle by variation and selection. Thus, it is a complementary strategy that offers compelling scientific and applied advantages for strain development and process optimization, provided a desired phenotype is amenable to direct or indirect selection. In addition to simple empirical strain development by random mutation and direct selection on plates, evolutionary engineering also encompasses recombination and continuous evolution of large populations over many generations. Two distinct evolutionary engineering applications are likely to gain more relevance in the future: first, as an integral component in metabolic engineering of strains with improved phenotypes, and second, to elucidate the molecular basis of desired phenotypes for subsequent transfer to other hosts. The latter will profit from the broader availability of recently developed methodologies for global response analysis at the genetic and metabolic level. These methodologies facilitate identification of the molecular basis of evolved phenotypes. It is anticipated that, together with novel analytical techniques, bioinformatics, and computer modeling of cellular functions and activities, evolutionary engineering is likely to find its place in the metabolic engineer's toolbox for research and strain development. This review presents evolutionary engineering of whole cells as an emerging methodology that draws on the latest advances from a wide range of scientific and technical disciplines.

Metabolic flux responses to pyruvate kinase knockout in Escherichia coli

The intracellular carbon flux distribution in wild-type and pyruvate kinase-deficient Escherichia coli was estimated using biosynthetically directed fractional 13C labeling experiments with [U-13C6]glucose in glucose- or ammonia-limited chemostats, two-dimensional nuclear magnetic resonance (NMR) spectroscopy of cellular amino acids, and a comprehensive isotopomer model. The general response to disruption of both pyruvate kinase isoenzymes in E. coli was a local flux rerouting via the combined reactions of phosphoenolpyruvate (PEP) carboxylase and malic enzyme. Responses in the pentose phosphate pathway and the tricarboxylic acid cycle were strongly dependent on the environmental conditions. In addition, high futile cycling activity via the gluconeogenic PEP carboxykinase was identified at a low dilution rate in glucose-limited chemostat culture of pyruvate kinase-deficient E. coli, with a turnover that is comparable to the specific glucose uptake rate. Furthermore, flux analysis in mutant cultures indicates that glucose uptake in E. coli is not catalyzed exclusively by the phosphotransferase system in glucose-limited cultures at a low dilution rate. Reliability of the flux estimates thus obtained was verified by statistical error analysis and by comparison to intracellular carbon flux ratios that were independently calculated from the same NMR data by metabolic flux ratio analysis.

Dissection of central carbon metabolism of hemoglobin-expressing Escherichia coli by 13C nuclear magnetic resonance flux distribution analysis in microaerobic bioprocesses

Escherichia coli MG1655 cells expressing Vitreoscilla hemoglobin (VHb), Alcaligenes eutrophus flavohemoprotein (FHP), the N-terminal hemoglobin domain of FHP (FHPg), and a fusion protein which comprises VHb and the A. eutrophus C-terminal reductase domain (VHb-Red) were grown in a microaerobic bioreactor to study the effects of low oxygen concentrations on the central carbon metabolism, using fractional (13)C-labeling of the proteinogenic amino acids and two-dimensional [(13)C, (1)H]-correlation nuclear magnetic resonance (NMR) spectroscopy. The NMR data revealed differences in the intracellular carbon fluxes between E. coli cells expressing either VHb or VHb-Red and cells expressing A. eutrophus FHP or the truncated heme domain (FHPg). E. coli MG1655 cells expressing either VHb or VHb-Red were found to function with a branched tricarboxylic acid (TCA) cycle. Furthermore, cellular demands for ATP and reduction equivalents in VHb- and VHb-Red-expressing cells were met by an increased flux through glycolysis. In contrast, in E. coli cells expressing A. eutrophus hemeproteins, the TCA cycle is running cyclically, indicating a shift towards a more aerobic regulation. Consistently, E. coli cells displaying FHP and FHPg activity showed lower production of the typical anaerobic by-products formate, acetate, and D-lactate. The implications of these observations for biotechnological applications are discussed.

Metabolic networks: a signal-oriented approach to cellular models

Complete genomes, far advanced proteomes, and even 'metabolomes' are available for at least a few organisms, e.g., Escherichia coli. Systematic functional analyses of such complete data sets will produce a wealth of information and promise an understanding of the dynamics of complex biological networks and perhaps even of entire living organisms. Such complete and holistic descriptions of biological systems, however, will increasingly require a quantitative analysis and the help of mathematical models for simulating whole systems. In particular, new procedures are required that allow a meaningful reduction of the information derived from complex systems that will consequently be used in the modeling process. In this review the biological elements of such a modeling procedure will be described. In a first step, complex living systems must be structured into well-defined and clearly delimited functional units, the elements of which have a common physiological goal, belong to a single genetic unit, and respond to the signals of a signal transduction system that senses changes in physiological states of the organism. These functional units occur at each level of complexity and more complex units originate by grouping several lower level elements into a single, more complex unit. To each complexity level corresponds a global regulator that is epistatic over lower level regulators. After its structuring into modules (functional units), a biological system is converted in a second step into mathematical submodels that by progressive combination can also be assembled into more aggregated model structures. Such a simplification of a cell (an organism) reduces its complexity to a level amenable to present modeling capacities. The universal biochemistry, however, promises a set of rules valid for modeling biological systems, from unicellular microorganisms and cells, to multicellular organisms and to populations.

Altered regulation of pyruvate kinase or co-overexpression of phosphofructokinase increases glycolytic fluxes in resting Escherichia coli

Glycolytic fluxes in resting Escherichia coli were enhanced by overexpression of heterologous pyruvate kinases (Pyk) from Bacillus stearothermophilus and Zymomonas mobilis, but not homologous Pyk. Compared to the control, an increase of 10% in specific glucose consumption and of 15% in specific ethanol production rates was found in anaerobic resting cells, expressing the heterologous Pyks, that were harvested from exponentially growing aerobic cultures. A further increase in glycolytic flux was achieved by simultaneous overexpression of E. coli phosphofructokinase (Pfk) and Pyk with specific glucose consumption and ethanol production rates of 25% and 35% greater, respectively, than the control. Fluxes to lactate were not significantly affected, contrary to previous observations with resting cells harvested from anaerobically growing cultures. To correlate the physiology of resting cells with the physiology of cells prior to harvest, we determined the relevant growth parameters from aerobic and anaerobic precultures. We conclude that glycolytic fluxes in E. coli with submaximal (aerobic) metabolic activity can be increased by overexpression of pyruvate kinases which do not require allosteric activation or co-overexpression with Pfk. However, such improvements require more extensive engineering in E. coli with near maximal (anaerobic) metabolic activity.

Metabolic flux ratio analysis of genetic and environmental modulations of Escherichia coli central carbon metabolism

The response of Escherichia coli central carbon metabolism to genetic and environmental manipulation has been studied by use of a recently developed methodology for metabolic flux ratio (METAFoR) analysis; this methodology can also directly reveal active metabolic pathways. Generation of fluxome data arrays by use of the METAFoR approach is based on two-dimensional (13)C-(1)H correlation nuclear magnetic resonance spectroscopy with fractionally labeled biomass and, in contrast to metabolic flux analysis, does not require measurements of extracellular substrate and metabolite concentrations. METAFoR analyses of E. coli strains that moderately overexpress phosphofructokinase, pyruvate kinase, pyruvate decarboxylase, or alcohol dehydrogenase revealed that only a few flux ratios change in concert with the overexpression of these enzymes. Disruption of both pyruvate kinase isoenzymes resulted in altered flux ratios for reactions connecting the phosphoenolpyruvate (PEP) and pyruvate pools but did not significantly alter central metabolism. These data indicate remarkable robustness and rigidity in central carbon metabolism in the presence of genetic variation. More significant physiological changes and flux ratio differences were seen in response to altered environmental conditions. For example, in ammonia-limited chemostat cultures, compared to glucose-limited chemostat cultures, a reduced fraction of PEP molecules was derived through at least one transketolase reaction, and there was a higher relative contribution of anaplerotic PEP carboxylation than of the tricarboxylic acid (TCA) cycle for oxaloacetate synthesis. These two parameters also showed significant variation between aerobic and anaerobic batch cultures. Finally, two reactions catalyzed by PEP carboxykinase and malic enzyme were identified by METAFoR analysis; these had previously been considered absent in E. coli cells grown in glucose-containing media. Backward flux from the TCA cycle to glycolysis, as indicated by significant activity of PEP carboxykinase, was found only in glucose-limited chemostat culture, demonstrating that control of this futile cycle activity is relaxed under severe glucose limitation.