In silico predictions of Escherichia coli metabolic capabilities are consistent with experimental data

A significant goal in the post-genome era is to relate the annotated genome sequence to the physiological functions of a cell. Working from the annotated genome sequence, as well as biochemical and physiological information, it is possible to reconstruct complete metabolic networks. Furthermore, computational methods have been developed to interpret and predict the optimal performance of a metabolic network under a range of growth conditions. We have tested the hypothesis that Escherichia coli uses its metabolism to grow at a maximal rate using the E. coli MG1655 metabolic reconstruction. Based on this hypothesis, we formulated experiments that describe the quantitative relationship between a primary carbon source (acetate or succinate) uptake rate, oxygen uptake rate, and maximal cellular growth rate. We found that the experimental data were consistent with the stated hypothesis, namely that the E. coli metabolic network is optimized to maximize growth under the experimental conditions considered. This study thus demonstrates how the combination of in silico and experimental biology can be used to obtain a quantitative genotype-phenotype relationship for metabolism in bacterial cells.     

Development of a software tool for in silico simulation of Escherichia coli using a visual programming environment

This study describes the development of a software tool, EcoSim, to assist users in implementing quantitative in silico simulation easily. It consists of four parts: extracellular environment and constraints setting mode, table for optimal metabolic flux distribution and chart for changes of substrate concentration, dynamic flux distribution viewer and dynamic hierarchical regulatory network viewer. Representation of a hierarchical regulatory network was constructed with defined modeling symbols and weight in the central Escherichia coli metabolism. All programming procedures for EcoSim were accomplished in a visual programming environment (LabVIEW). To illustrate quantitative in silico simulation with EcoSim, this program was performed on E. coli using glucose and acetate as carbon sources. The simulation results were in agreement with the experimental data obtained from the literature. EcoSim can be used to assist biologists and engineers in predicting and interpreting dynamic behaviors of E. coli under a variety of environmental conditions.     

EcoCyc: an encyclopedia of Escherichia coli genes and metabolism.

The encyclopedia of Escherichia coli genes and metabolism (EcoCyc) is a database that combines information about the genome and the intermediary metabolism of E.coli. It describes 2034 genes, 306 enzymes encoded by these genes, 580 metabolic reactions that occur in E.coli and the organization of these reactions into 100 metabolic pathways. The EcoCyc graphical user interface allows query and exploration of 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 investigation of 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 silico model of E.coli that can be probed and analyzed through computational means.     

Description and interpretation of adaptive evolution of Escherichia coli K-12 MG1655 by using a genome-scale in silico metabolic model

Genome-scale in silico metabolic networks of Escherichia coli have been reconstructed. By using a constraint-based in silico model of a reconstructed network, the range of phenotypes exhibited by E. coli under different growth conditions can be computed, and optimal growth phenotypes can be predicted. We hypothesized that the end point of adaptive evolution of E. coli could be accurately described a priori by our in silico model since adaptive evolution should lead to an optimal phenotype. Adaptive evolution of E. coli during prolonged exponential growth was performed with M9 minimal medium supplemented with 2 g of alpha-ketoglutarate per liter, 2 g of lactate per liter, or 2 g of pyruvate per liter at both 30 and 37 degrees C, which produced seven distinct strains. The growth rates, substrate uptake rates, oxygen uptake rates, by-product secretion patterns, and growth rates on alternative substrates were measured for each strain as a function of evolutionary time. Three major conclusions were drawn from the experimental results. First, adaptive evolution leads to a phenotype characterized by maximized growth rates that may not correspond to the highest biomass yield. Second, metabolic phenotypes resulting from adaptive evolution can be described and predicted computationally. Third, adaptive evolution on a single substrate leads to changes in growth characteristics on other substrates that could signify parallel or opposing growth objectives. Together, the results show that genome-scale in silico metabolic models can describe the end point of adaptive evolution a priori and can be used to gain insight into the adaptive evolutionary process for E. coli.     

Characterizing Escherichia coli DH5alpha growth and metabolism in a complex medium using genome-scale flux analysis

Genome-scale flux analysis of Escherichia coli DH5alpha growth in a complex medium was performed to investigate the relationship between the uptake of various nutrients and their metabolic outcomes. During the exponential growth phase, we observed a sequential consumption order of serine, aspartate and glutamate in the complex medium as well as the complete consumption of key carbohydrate nutrients, glucose and trehalose. Based on the consumption and production rates of the measured metabolites, constraints-based flux analysis of a genome-scale E. coli model was then conducted to elucidate their utilization in the metabolism. The in silico analysis revealed that the cell exploited biosynthetic precursors taken up directly from the complex medium, through growth-related anabolic pathways. This suggests that the cell could be functioning in an energetically more efficient manner by reducing the energy needed to produce amino acids. The in silico simulation also allowed us to explain the observed rapid consumption of serine: excessively consumed external serine from the complex medium was mainly converted into pyruvate and glycine, which in turn, led to the acetate accumulation. The present work demonstrates the application of an in silico modeling approach to characterizing microbial metabolism under complex medium condition. This work further illustrates the use of in silico genome-scale analysis for developing better strategies related to improving microbial growth and enhancing the productivity of desirable metabolites.     

Elucidating and reprogramming Escherichia coli metabolisms for obligate anaerobic n-butanol and isobutanol production

Elementary mode (EM) analysis based on the constraint-based metabolic network modeling was applied to elucidate and compare complex fermentative metabolisms of Escherichia coli for obligate anaerobic production of n-butanol and isobutanol. The result shows that the n-butanol fermentative metabolism was NADH-deficient, while the isobutanol fermentative metabolism was NADH redundant. E. coli could grow and produce n-butanol anaerobically as the sole fermentative product but not achieve the maximum theoretical n-butanol yield. In contrast, for the isobutanol fermentative metabolism, E. coli was required to couple with either ethanol- or succinate-producing pathway to recycle NADH. To overcome these "defective" metabolisms, EM analysis was implemented to reprogram the native fermentative metabolism of E. coli for optimized anaerobic production of n-butanol and isobutanol through multiple gene deletion (~8-9 genes), addition (~6-7 genes), up- and downexpression (~6-7 genes), and cofactor engineering (e.g., NADH, NADPH). The designed strains were forced to couple both growth and anaerobic production of n-butanol and isobutanol, which is a useful characteristic to enhance biofuel production and tolerance through metabolic pathway evolution. Even though the n-butanol and isobutanol fermentative metabolisms were quite different, the designed strains could be engineered to have identical metabolic flux distribution in "core" metabolic pathways mainly supporting cell growth and maintenance. Finally, the model prediction in elucidating and reprogramming the native fermentative metabolism of E. coli for obligate anaerobic production of n-butanol and isobutanol was validated with published experimental data.     

Reconstruction of a genome-scale metabolic network of Rhodococcus erythropolis for desulfurization studies

The remarkable catabolic diversity of Rhodococcus erythropolis makes it an interesting organism for bioremediation and fuel desulfurization. However, a model that can describe and explain the combined influence of various intracellular metabolic activities on its desulfurizing capabilities is missing from the literature. Such a model can greatly aid the development of R. erythropolis as an effective desulfurizing biocatalyst. This work reports the reconstruction of the first genome-scale metabolic model for R. erythropolis using the available genomic, experimental, and biochemical information. We have validated our in silico model by successfully predicting cell growth results and explaining several experimental observations in the literature on biodesulfurization using dibenzothiophene. We report several in silico experiments and flux balance analyses to propose minimal media, determine gene and reaction essentiality, and compare effectiveness of carbon, nitrogen, and sulfur sources. We demonstrate the usefulness of our model by studying a few in silico mutants of R. erythropolis for improved biodesulfurization, and comparing the desulfurization abilities of R. erythropolis with an in silico mutant of E. coli.     

Genome-scale models of bacterial metabolism: reconstruction and applications

Genome-scale metabolic models bridge the gap between genome-derived biochemical information and metabolic phenotypes in a principled manner, providing a solid interpretative framework for experimental data related to metabolic states, and enabling simple in silico experiments with whole-cell metabolism. Models have been reconstructed for almost 20 bacterial species, so far mainly through expert curation efforts integrating information from the literature with genome annotation. A wide variety of computational methods exploiting metabolic models have been developed and applied to bacteria, yielding valuable insights into bacterial metabolism and evolution, and providing a sound basis for computer-assisted design in metabolic engineering. Recent advances in computational systems biology and high-throughput experimental technologies pave the way for the systematic reconstruction of metabolic models from genomes of new species, and a corresponding expansion of the scope of their applications. In this review, we provide an introduction to the key ideas of metabolic modeling, survey the methods, and resources that enable model reconstruction and refinement, and chart applications to the investigation of global properties of metabolic systems, the interpretation of experimental results, and the re-engineering of their biochemical capabilities.     

Systematic analysis of conservation relations in Escherichia coli genome-scale metabolic network reveals novel growth media

A biochemical species is called producible in a constraints-based metabolic model if a feasible steady-state flux configuration exists that sustains its nonzero concentration during growth. Extreme semipositive conservation relations (ESCRs) are the simplest semipositive linear combinations of species concentrations that are invariant to all metabolic flux configurations. In this article, we outline a fundamental relationship between the ESCRs of a metabolic network and the producibility of a biochemical species under a nutrient media. We exploit this relationship in an algorithm that systematically enumerates all minimal nutrient sets that render an objective species weakly producible (i.e., producible in the absence of thermodynamic constraints) through a simple traversal of ESCRs. We apply our results to a recent genome scale model of Escherichia coli metabolism, in which we traverse the 51 anhydrous ESCRs of the metabolic network to determine all 928 minimal aqueous nutrient media that render biomass weakly producible. Applying irreversibility constraints, we find 287 of these 928 nutrient sets to be thermodynamically feasible. We also find that an additional 365 of these nutrient sets are thermodynamically feasible in the presence of oxygen. Since biomass producibility is commonly used as a surrogate for growth in genome scale metabolic models, our results represent testable hypotheses of alternate growth media derived from in silico analysis of the E. coli genome scale metabolic network.     

Integrating cybernetic modeling with pathway analysis provides a dynamic, systems-level description of metabolic control

Cybernetic modeling strives to uncover the inbuilt regulatory programs of biological systems and leverage them toward computational prediction of metabolic dynamics. Because of its focus on incorporating the global aims of metabolism, cybernetic modeling provides a systems-oriented approach for describing regulatory inputs and inferring the impact of regulation within biochemical networks. Combining cybernetic control laws with concepts from metabolic pathway analysis has culminated in a systematic strategy for constructing cybernetic models, which was previously lacking. The newly devised framework relies upon the simultaneous application of local controls that maximize the net flux through each elementary flux mode and global controls that modulate the activities of these modes to optimize the overall nutritional state of the cell. The modeling concepts are illustrated using a simple linear pathway and a larger network representing anaerobic E. coli central metabolism. The E. coli model successfully describes the metabolic shift that occurs upon deleting the pta-ackA operon that is responsible for fermentative acetate production. The model also furnishes predictions that are consistent with experimental results obtained from additional knockout strains as well as strains expressing heterologous genes. Because of the stabilizing influence of the included control variables, the resulting cybernetic models are more robust and reliable than their predecessors in simulating the network response to imposed genetic and environmental perturbations.     

Enzyme activity and fitness: Evolution in solution

Natural selection should be studied as an end in itself, and this requires rigorous experimental tests of theoretical models linking molecular phenotypes to differences in fitness. We describe the experimental verification of one such model and thereby demonstrate that the causal relations between genotype and fitness need not be as hopelessly complex as many have assumed. The model uses metabolic control theory to link enzyme activity to metabolic flux and then assumes that fitness is proportional to flux. The model was tested using the pathway for the uptake and metabolism of growth-rate-limiting concentrations of lactose in E. coli inhabiting chemostats. Many of the properties expected of natural selection are manifest in this system.     

Genome-scale analysis of Mannheimia succiniciproducens metabolism

Mannheimia succiniciproducens MBEL55E isolated from bovine rumen is a capnophilic gram-negative bacterium that efficiently produces succinic acid, an industrially important four carbon dicarboxylic acid. In order to design a metabolically engineered strain which is capable of producing succinic acid with high yield and productivity, it is essential to optimize the whole metabolism at the systems level. Consequently, in silico modeling and simulation of the genome-scale metabolic network was employed for genome-scale analysis and efficient design of metabolic engineering experiments. The genome-scale metabolic network of M. succiniciproducens consisting of 686 reactions and 519 metabolites was constructed based on reannotation and validation experiments. With the reconstructed model, the network structure and key metabolic characteristics allowing highly efficient production of succinic acid were deciphered; these include strong PEP carboxylation, branched TCA cycle, relative weak pyruvate formation, the lack of glyoxylate shunt, and non-PTS for glucose uptake. Constraints-based flux analyses were then carried out under various environmental and genetic conditions to validate the genome-scale metabolic model and to decipher the altered metabolic characteristics. Predictions based on constraints-based flux analysis were mostly in excellent agreement with the experimental data. In silico knockout studies allowed prediction of new metabolic engineering strategies for the enhanced production of succinic acid. This genome-scale in silico model can serve as a platform for the systematic prediction of physiological responses of M. succiniciproducens to various environmental and genetic perturbations and consequently for designing rational strategies for strain improvement.     

Genome-scale metabolic model of Helicobacter pylori 26695

A genome-scale metabolic model of Helicobacter pylori 26695 was constructed from genome sequence annotation, biochemical, and physiological data. This represents an in silico model largely derived from genomic information for an organism for which there is substantially less biochemical information available relative to previously modeled organisms such as Escherichia coli. The reconstructed metabolic network contains 388 enzymatic and transport reactions and accounts for 291 open reading frames. Within the paradigm of constraint-based modeling, extreme-pathway analysis and flux balance analysis were used to explore the metabolic capabilities of the in silico model. General network properties were analyzed and compared to similar results previously generated for Haemophilus influenzae. A minimal medium required by the model to generate required biomass constituents was calculated, indicating the requirement of eight amino acids, six of which correspond to essential human amino acids. In addition a list of potential substrates capable of fulfilling the bulk carbon requirements of H. pylori were identified. A deletion study was performed wherein reactions and associated genes in central metabolism were deleted and their effects were simulated under a variety of substrate availability conditions, yielding a number of reactions that are deemed essential. Deletion results were compared to recently published in vitro essentiality determinations for 17 genes. The in silico model accurately predicted 10 of 17 deletion cases, with partial support for additional cases. Collectively, the results presented herein suggest an effective strategy of combining in silico modeling with experimental technologies to enhance biological discovery for less characterized organisms and their genomes.     

Systems biotechnology for strain improvement

Various high-throughput experimental techniques are routinely used for generating large amounts of omics data. In parallel, in silico modelling and simulation approaches are being developed for quantitatively analyzing cellular metabolism at the systems level. Thus informative high-throughput analysis and predictive computational modelling or simulation can be combined to generate new knowledge through iterative modification of an in silico model and experimental design. On the basis of such global cellular information we can design cells that have improved metabolic properties for industrial applications. This article highlights the recent developments in these systems approaches, which we call systems biotechnology, and discusses future prospects.     

Quantitative analysis of metabolic fluxes in Escherichia coli, using two-dimensional NMR spectroscopy and complete isotopomer models.

Complete isotopomer models that simulate distribution of label in 13C tracer experiments are applied to the quantification of metabolic fluxes in the primary carbon metabolism of E. coli under aerobic and anaerobic conditions. The concept of isotopomer mapping matrices (IMMs) is used to simplify the formulation of isotopomer mass balances by expressing all isotopomer mass balances of a metabolite pool in a single matrix equation. A numerically stable method to calculate the steady-state isotopomer distribution in metabolic networks in introduced. Net values of intracellular fluxes and the degree of reversibility of enzymatic steps are estimated by minimization of the deviations between experimental and simulated measurements. The metabolic model applied includes the Embden-Meyerhof-Parnas and the pentose phosphate pathway, the tricarboxylic acid cycle, anaplerotic reaction sequences and pathways involved in amino acid synthesis. The study clearly demonstrates the value of complete isotopomer models for maximizing the information obtainable from 13C tracer experiments. The approach applied here offers a completely general and comprehensive analysis of carbon tracer experiments where any set of experimental data on the labeling state and extracellular fluxes can be used for the quantification of metabolic fluxes in complex metabolic networks.