Stability of Metabolic Correlations under Changing Environmental Conditions in Escherichia coli – A Systems Approach

Background

Biological systems adapt to changing environments by reorganizing their cellular and physiological program with metabolites representing one important response level. Different stresses lead to both conserved and specific responses on the metabolite level which should be reflected in the underlying metabolic network.

Methodology/Principal Findings

Starting from experimental data obtained by a GC-MS based high-throughput metabolic profiling technology we here develop an approach that: (1) extracts network representations from metabolic condition-dependent data by using pairwise correlations, (2) determines the sets of stable and condition-dependent correlations based on a combination of statistical significance and homogeneity tests, and (3) can identify metabolites related to the stress response, which goes beyond simple observations about the changes of metabolic concentrations. The approach was tested with Escherichia coli as a model organism observed under four different environmental stress conditions (cold stress, heat stress, oxidative stress, lactose diauxie) and control unperturbed conditions. By constructing the stable network component, which displays a scale free topology and small-world characteristics, we demonstrated that: (1) metabolite hubs in this reconstructed correlation networks are significantly enriched for those contained in biochemical networks such as EcoCyc, (2) particular components of the stable network are enriched for functionally related biochemical pathways, and (3) independently of the response scale, based on their importance in the reorganization of the correlation network a set of metabolites can be identified which represent hypothetical candidates for adjusting to a stress-specific response.

Conclusions/Significance

Network-based tools allowed the identification of stress-dependent and general metabolic correlation networks. This correlation-network-based approach does not rely on major changes in concentration to identify metabolites important for stress adaptation, but rather on the changes in network properties with respect to metabolites. This should represent a useful complementary technique in addition to more classical approaches.     

The Escherichia coli Proteome: Past, Present, and Future Prospects

Proteomics has emerged as an indispensable methodology for large-scale protein analysis in functional genomics. The Escherichia coli proteome has been extensively studied and is well defined in terms of biochemical, biological, and biotechnological data. Even before the entire E. coli proteome was fully elucidated, the largest available data set had been integrated to decipher regulatory circuits and metabolic pathways, providing valuable insights into global cellular physiology and the development of metabolic and cellular engineering strategies. With the recent advent of advanced proteomic technologies, the E. coli proteome has been used for the validation of new technologies and methodologies such as sample prefractionation, protein enrichment, two-dimensional gel electrophoresis, protein detection, mass spectrometry (MS), combinatorial assays with n-dimensional chromatographies and MS, and image analysis software. These important technologies will not only provide a great amount of additional information on the E. coli proteome but also synergistically contribute to other proteomic studies. Here, we review the past development and current status of E. coli proteome research in terms of its biological, biotechnological, and methodological significance and suggest future prospects.     

Characterization of Unexpected Growth of Escherichia coli O157:H7 by Modeling

Modeling of batch kinetics in minimal synthetic medium was used to characterize Escherichia coli O157:H7 growth, which appeared to be different from the exponential growth expected in minimal synthetic medium and observed for E. coli K-12. The turbidimetric kinetics of 14 of the 15 O157:H7 strains tested (93%) were nonexponential, whereas 25 of the 36 other E. coli strains tested (70%) exhibited exponential kinetics. Moreover, the anomaly was almost corrected when the minimal medium was supplemented with methionine. These observations were confirmed with two reference strains by using plate count monitoring. In mixed cultures, E. coli K-12 had a positive effect on E. coli O157:H7 and corrected its growth anomaly. This demonstrated that commensalism occurred, as the growth curve for E. coli K-12 was not affected. The interaction could be explained by an exchange of methionine, as the effect of E. coli K-12 on E. coli O157:H7 appeared to be similar to the effect of methionine.     

A novel platform for automated high-throughput fluxome profiling of metabolic variants

Advances in metabolic engineering are enabling the creation of a large number of cell factories. However, high-throughput platforms do not yet exist for rapidly analyzing the metabolic network of the engineered cells. To fill the gap, we developed an integrated solution for fluxome profiling of large sets of biological systems and conditions. This platform combines a robotic system for ¹³C-labelling experiments and sampling of labelled material with NMR-based isotopic fingerprinting and automated data interpretation. As a proof-of-concept, this workflow was applied to discriminate between Escherichia coli mutants with gradual expression of the glucose-6-phosphate dehydrogenase. Metabolic variants were clearly discriminated while pathways that support metabolic flexibility towards modulation of a single enzyme were elucidating. By directly connecting the data flow between cell cultivation and flux quantification, considerable advances in throughput, robustness, release of resources and screening capacity were achieved. This will undoubtedly facilitate the development of efficient cell factories.     

Solitary BioY Proteins Mediate Biotin Transport into Recombinant Escherichia coli

Energy-coupling factor (ECF) transporters form a large group of vitamin uptake systems in prokaryotes. They are composed of highly diverse, substrate-specific, transmembrane proteins (S units), a ubiquitous transmembrane protein (T unit), and homo- or hetero-oligomeric ABC ATPases. Biotin transporters represent a special case of ECF-type systems. The majority of the biotin-specific S units (BioY) is known or predicted to interact with T units and ABC ATPases. About one-third of BioY proteins, however, are encoded in organisms lacking any recognizable T unit. This finding raises the question of whether these BioYs function as transporters in a solitary state, a feature ascribed to certain BioYs in the past. To address this question in living cells, an Escherichia coli K-12 derivative deficient in biotin synthesis and devoid of its endogenous high-affinity biotin transporter was constructed as a reference strain. This organism is particularly suited for this purpose because components of ECF transporters do not naturally occur in E. coli K-12. The double mutant was viable in media containing either high levels of biotin or a precursor of the downstream biosynthetic path. Importantly, it was nonviable on trace levels of biotin. Eight solitary bioY genes of proteobacterial origin were individually expressed in the reference strain. Each of the BioYs conferred biotin uptake activity on the recombinants, which was inferred from uptake assays with [³H]biotin and growth of the cells on trace levels of biotin. The results underscore that solitary BioY transports biotin across the cytoplasmic membrane.     

Directed modification of Escherichia coli metabolism for the design of threonine-producing strains

The modification of Escherichia coli K-12 metabolism leading to threonine overproduction is the most studied system in synthetic biology that has been used to elaborate the majority of the currently known approaches to constructing microbial producers. They include optimization of biosynthesis through search for rate-limiting stages, modification of substrate and product transport, elimination of side metabolic pathways and degradation systems, reinforcement of the regeneration of coenzymes that are required for product biosynthesis, and exclusion of futile cycles and metabolic pathways with low energy efficiency. Extensive research in functional genomics made it possible to selectively remove the “unnecessary genes,” the functions of which are useless for producing a strain or adversely affect its properties. In total, using various approaches to designing threonine-producing strains, over 150 genome loci that affect more than 30% genes in E. coli were directly modified, thus providing interesting data for researchers in the field of microbial synthesis, as well as in related biological sciences. This review is dedicated to the assessment of genetic engineering modifications in E. coli metabolism (primarily, on the basis of modern patent literature) that ensure threonine overproduction.     

Physiological and Proteomic Analysis of Escherichia coli Iron-Limited Chemostat Growth

Iron bioavailability is a major limiter of bacterial growth in mammalian host tissue and thus represents an important area of study. Escherichia coli K-12 metabolism was studied at four levels of iron limitation in chemostats using physiological and proteomic analyses. The data documented an E. coli acclimation gradient where progressively more severe iron scarcity resulted in a larger percentage of substrate carbon being directed into an overflow metabolism accompanied by a decrease in biomass yield on glucose. Acetate was the primary secreted organic by-product for moderate levels of iron limitation, but as stress increased, the metabolism shifted to secrete primarily lactate (∼70% of catabolized glucose carbon). Proteomic analysis reinforced the physiological data and quantified relative increases in glycolysis enzyme abundance and decreases in tricarboxylic acid (TCA) cycle enzyme abundance with increasing iron limitation stress. The combined data indicated that E. coli responds to limiting iron by investing the scarce resource in essential enzymes, at the cost of catabolic efficiency (i.e., downregulating high-ATP-yielding pathways containing enzymes with large iron requirements, like the TCA cycle). Acclimation to iron-limited growth was contrasted experimentally with acclimation to glucose-limited growth to identify both general and nutrient-specific acclimation strategies. While the iron-limited cultures maximized biomass yields on iron and increased expression of iron acquisition strategies, the glucose-limited cultures maximized biomass yields on glucose and increased expression of carbon acquisition strategies. This study quantified ecologically competitive acclimations to nutrient limitations, yielding knowledge essential for understanding medically relevant bacterial responses to host and to developing intervention strategies.     

Genome-Wide Assessment of Outer Membrane Vesicle Production in Escherichia coli

The production of outer membrane vesicles by Gram-negative bacteria has been well documented; however, the mechanism behind the biogenesis of these vesicles remains unclear. Here a high-throughput experimental method and systems-scale analysis was conducted to determine vesiculation values for the whole genome knockout library of Escherichia coli mutant strains (Keio collection). The resultant dataset quantitatively recapitulates previously observed phenotypes and implicates nearly 150 new genes in the process of vesiculation. Gene functional and biochemical pathway analyses suggest that mutations that truncate outer membrane structures such as lipopolysaccharide and enterobacterial common antigen lead to hypervesiculation, whereas mutants in oxidative stress response pathways result in lower levels. This study expands and refines the current knowledge regarding the cellular pathways required for outer membrane vesiculation in E. coli.     

mhpT Encodes an Active Transporter Involved in 3-(3-Hydroxyphenyl)Propionate Catabolism by Escherichia coli K-12

Escherichia coli K-12 utilizes 3-(3-hydroxyphenyl)propionate (3HPP) as a sole carbon and energy source. Among the genes in its catabolic cluster in the genome, mhpT was proposed to encode a hypothetical transporter. Since no transporter for 3HPP uptake has been identified, we investigated whether MhpT is responsible for 3HPP uptake. MhpT fused with green fluorescent protein was found to be located at the periphery of cells by confocal microscopy, consistent with localization to the cytoplasmic membrane. Gene knockout and complementation studies clearly indicated that mhpT is essential for 3HPP catabolism in E. coli K-12 W3110 at pH 8.2. Uptake assays with ¹⁴C-labeled substrates demonstrated that strain W3110 and strain W3110ΔmhpT containing recombinant MhpT specifically transported 3HPP but not benzoate, 3-hydroxybenzoate, or gentisate into cells. Energy dependence assays suggested that MhpT-mediated 3HPP transport was driven by the proton motive force. The change of Ala-272 of MhpT to a histidine, surprisingly, resulted in enhanced transport activity, and strain W3110ΔmhpT containing the MhpT A272H mutation had a slightly higher growth rate than the wild-type strain at pH 8.2. Hence, we demonstrated that MhpT is a specific 3HPP transporter and vital for E. coli K-12 W3110 growth on this substrate under basic conditions.     

Comparison of the function of L- and D-proline as compatible solute inEscherichia coli K-12 under high osmolarity

We compared the function of L- and D-proline as compatible solute inEscherichia coli K-12 cells under high osmolarity. Growth ofE. coli K-12 in a Davis minimal medium was inhibited at 0.5 M and 1 M NaCl, but it was recovered by the addition of L-proline. Glucose uptake was reduced with increase of external NaCl concentrations, but it was improved by the addition of L-proline. On the other hand, the addition of D-proline did not show the role of compatible solute although accumulated in cells. On the analysis ofE. coli proline transporter mutants, difference of the affinity of proline transporters for D-proline was observed at PutP and ProP. These results presumed that the functional disorder of D-proline as compatible solute was caused by its structural feature in cells.     

Eco Cyc: encyclopedia of Escherichia coli genes and metabolism.

The EcoCyc database describes the genome and gene products of Escherichia coli, its metabolic and signal-transduction pathways, and its tRNAs. The database describes 4391 genes of E.coli, 695 enzymes encoded by a subset of these genes, 904 metabolic reactions that occur in E.coli, and the organization of these reactions into 129 metabolic pathways. The EcoCyc graphical user interface allows scientists to query and explore the EcoCyc database using visualization tools such as genomic-map browsers and automatic layouts of metabolic pathways. EcoCyc has many references to the primary literature, and is a (qualitative) computational model of E. coli metabolism. EcoCyc is available at URL http://ecocyc. PangeaSystems.com/ecocyc/     

Complete Genome Sequence and Comparative Metabolic Profiling of the Prototypical Enteroaggregative Escherichia coli Strain 042

Background

Escherichia coli can experience a multifaceted life, in some cases acting as a commensal while in other cases causing intestinal and/or extraintestinal disease. Several studies suggest enteroaggregative E. coli are the predominant cause of E. coli-mediated diarrhea in the developed world and are second only to Campylobacter sp. as a cause of bacterial-mediated diarrhea. Furthermore, enteroaggregative E. coli are a predominant cause of persistent diarrhea in the developing world where infection has been associated with malnourishment and growth retardation.

Methods

In this study we determined the complete genomic sequence of E. coli 042, the prototypical member of the enteroaggregative E. coli, which has been shown to cause disease in volunteer studies. We performed genomic and phylogenetic comparisons with other E. coli strains revealing previously uncharacterised virulence factors including a variety of secreted proteins and a capsular polysaccharide biosynthetic locus. In addition, by using Biolog™ Phenotype Microarrays we have provided a full metabolic profiling of E. coli 042 and the non-pathogenic lab strain E. coli K-12. We have highlighted the genetic basis for many of the metabolic differences between E. coli 042 and E. coli K-12.

Conclusion

This study provides a genetic context for the vast amount of experimental and epidemiological data published thus far and provides a template for future diagnostic and intervention strategies.     

Improving cellular malonyl-CoA level in Escherichia coli via metabolic engineering

Escherichia coli only maintains a small amount of cellular malonyl-CoA, impeding its utility for overproducing natural products such as polyketides and flavonoids. Here, we report the use of various metabolic engineering strategies to redirect the carbon flux inside E. coli to pathways responsible for the generation of malonyl-CoA. Overexpression of acetyl-CoA carboxylase (Acc) resulted in 3-fold increase in cellular malonyl-CoA concentration. More importantly, overexpression of Acc showed a synergistic effect with increased acetyl-CoA availability, which was achieved by deletion of competing pathways leading to the byproducts acetate and ethanol as well as overexpression of an acetate assimilation enzyme. These engineering efforts led to the creation of an E. coli strain with 15-fold elevated cellular malonyl-CoA level. To demonstrate its utility, this engineered E. coli strain was used to produce an important polyketide, phloroglucinol, and showed near 4-fold higher titer compared with wild-type E. coli, despite the toxicity of phloroglucinol to cell growth. This engineered E. coli strain with elevated cellular malonyl-CoA level should be highly useful for improved production of important natural products where the cellular malonyl-CoA level is rate-limiting.     

EcoCyc: Enyclopedia of Escherichia coli Genes and Metabolism.

The Encyclopedia of Genes and Metabolism (EcoCyc) is a database that combines information about the genome and the intermediary metabolism of Escherichia coli. It describes 2970 genes of E.coli, 547 enzymes encoded by these genes, 702 metabolic reactions that occur in E.coli and the organization of these reactions into 107 metabolic pathways. The EcoCyc graphical user interface allows scientists to query and explore the EcoCyc database using visualization tools such as genomic-map browsers and automatic layouts of metabolic pathways. EcoCyc spans the space from sequence to function to allow scientists to investigate an unusually broad range of questions. EcoCyc can be thought of as both an electronic review article because of its copious references to the primary literature, and as an in silicio model of E.coli metabolism that can be probed and analyzed through computational means.     

Proteome response of Escherichia coli fed-batch culture to temperature downshift

During fed-batch cultivation of Escherichia coli K-12, the proteomic response to a temperature downshift from 37 to 20°C was quantitatively monitored and analyzed by using two-dimensional electrophoresis. When the temperature of exponentially growing E. coli K-12 culture was downshifted to 20°C, the synthesis level of 57 intracellular proteins showed significant changes for a prolonged period of time, compared to the fed-batch culture controlled at 37°C. Thus, these proteins are regarded as important stress proteins responsive to cold shock, which were analyzed by using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry and identified using the E. coli SWISS-2DPAGE database. Most of the identified proteins were shown to be involved in energy metabolism, several cellular molecule biosynthetic pathways and catabolism, cell processes, flagellar biosynthesis and motility, and protein translation and folding. The systematic approach to the monitoring of proteomic responses and the detailed analysis results reported in this article would be useful in understanding the metabolic adaptation to lowered culture temperature and designing efficient fermentation strategies for the production of recombinant proteins and metabolites using E. coli strains.     

Reconstruction of the Saccharopolyspora erythraea genome-scale model and its use for enhancing erythromycin production

Genome-scale metabolic reconstructions are routinely used for the analysis and design of metabolic engineering strategies for production of primary metabolites. The use of such reconstructions for metabolic engineering of antibiotic production is not common due to the lack of simple design algorithms in the absence of a cellular growth objective function. Here, we present the metabolic network reconstruction for the erythromycin producer Saccharopolyspora erythraea NRRL23338. The model was manually curated for primary and secondary metabolism pathways and consists of 1,482 reactions (2,075 genes) and 1,646 metabolites. As part of the model validation, we explored the potential benefits of supplying amino acids and identified five amino acids “compatible” with erythromycin production, whereby if glucose is supplemented with this amino acid on a carbon mole basis, the in silico model predicts that high erythromycin yield is possible without lowering biomass yield. Increased erythromycin titre was confirmed for four of the five amino acids, namely valine, isoleucine, threonine and proline. In bioreactor experiments, supplementation with 2.5?% carbon mole of valine increased the growth rate by 20?% and simultaneously the erythromycin yield on biomass by 50?%. The model presented here can be used as a framework for the future integration of high-throughput biological data sets in S. erythraea and ultimately to realise strain designs capable of increasing erythromycin production closer to the theoretical yield.