A protocol for generating a high-quality genome-scale metabolic reconstruction

From Physics to Pharmacology?

Over the last fifty years there has been an explosion of biological data, leading to the realization that to fully explain biological mechanisms it is necessary to interpret them as complex dynamical systems. The first stage of this interpretation is to determine which components (proteins, genes or metabolites) of the system interact. This is usually represented by a graph, or network. The behavior of this network can then be investigated using mathematical modeling. In vivo these biological networks show several remarkable (and seemingly paradoxical) properties including robustness, plasticity and sensitivity. Erroneous behavior of these networks is often associated with disease. Hence understanding the system-level properties can have important implications for the treatment of disease. Systems biology is an organized approach to quantitatively describe and elucidate the behavior of these complex networks. This review focuses on the progress and future challenges of a systems approach to biology.

Whole-genome metabolic network reconstruction and constraint-based modeling

With the advent of modern high-throughput genomics, there is a significant need for genome-scale analysis techniques that can assist in complex systems analysis. Metabolic genome-scale network reconstructions (GENREs) paired with constraint-based modeling are an efficient method to integrate genomics, transcriptomics, and proteomics to conduct organism-specific analysis. This text explains key steps in the GENRE construction process and several methods of constraint-based modeling that can help elucidate basic life processes and development of disease treatment, bioenergy solutions, and industrial bioproduction applications.

Improving the iMM904 S. cerevisiae metabolic model using essentiality and synthetic lethality data

BACKGROUND

Saccharomyces cerevisiae is the first eukaryotic organism for which a multi-compartment genome-scale metabolic model was constructed. Since then a sequence of improved metabolic reconstructions for yeast has been introduced. These metabolic models have been extensively used to elucidate the organizational principles of yeast metabolism and drive yeast strain engineering strategies for targeted overproductions. They have also served as a starting point and a benchmark for the reconstruction of genome-scale metabolic models for other eukaryotic organisms. In spite of the successive improvements in the details of the described metabolic processes, even the recent yeast model (i.e., iMM904) remains significantly less predictive than the latest E. coli model (i.e., iAF1260). This is manifested by its significantly lower specificity in predicting the outcome of grow/no grow experiments in comparison to the E. coli model.

RESULTS

In this paper we make use of the automated GrowMatch procedure for restoring consistency with single gene deletion experiments in yeast and extend the procedure to make use of synthetic lethality data using the genome-scale model iMM904 as a basis. We identified and vetted using literature sources 120 distinct model modifications including various regulatory constraints for minimal and YP media. The incorporation of the suggested modifications led to a substantial increase in the fraction of correctly predicted lethal knockouts (i.e., specificity) from 38.84% (87 out of 224) to 53.57% (120 out of 224) for the minimal medium and from 24.73% (45 out of 182) to 40.11% (73 out of 182) for the YP medium. Synthetic lethality predictions improved from 12.03% (16 out of 133) to 23.31% (31 out of 133) for the minimal medium and from 6.96% (8 out of 115) to 13.04% (15 out of 115) for the YP medium.

CONCLUSIONS

Overall, this study provides a roadmap for the computationally driven correction of multi-compartment genome-scale metabolic models and demonstrates the value of synthetic lethals as curation agents.

Reconstruction and flux-balance analysis of the Plasmodium falciparum metabolic network

Genome-scale metabolic reconstructions can serve as important tools for hypothesis generation and high-throughput data integration. Here, we present a metabolic network reconstruction and flux-balance analysis (FBA) of Plasmodium falciparum, the primary agent of malaria. The compartmentalized metabolic network accounts for 1001 reactions and 616 metabolites. Enzyme-gene associations were established for 366 genes and 75% of all enzymatic reactions. Compared with other microbes, the P. falciparum metabolic network contains a relatively high number of essential genes, suggesting little redundancy of the parasite metabolism. The model was able to reproduce phenotypes of experimental gene knockout and drug inhibition assays with up to 90% accuracy. Moreover, using constraints based on gene-expression data, the model was able to predict the direction of concentration changes for external metabolites with 70% accuracy. Using FBA of the reconstructed network, we identified 40 enzymatic drug targets (i.e. in silico essential genes), with no or very low sequence identity to human proteins. To demonstrate that the model can be used to make clinically relevant predictions, we experimentally tested one of the identified drug targets, nicotinate mononucleotide adenylyltransferase, using a recently discovered small-molecule inhibitor.

Flux-based analysis of sulfur metabolism in desulfurizing strains of Rhodococcus erythropolis

Rhodococcus erythropolis has been studied widely for potential applications in biodesulfurization. Previous works have been largely experimental with an emphasis on the characterization and genetic engineering of desulfurizing strains for improved biocatalysis. A systems modeling approach that can complement these experimental efforts by providing useful insights into the complex interactions of desulfurization reactions with various other metabolic activities is absent in the literature. In this work, we report the first attempt at reconstructing a flux-based model to analyze sulfur utilization by R. erythropolis. The model includes the 4S pathway for dibenzothiophene (DBT) desulfurization. It predicts closely the growth rates reported by two independent experimental studies, and gives a clear and comprehensive picture of the pathways that assimilate the sulfur from DBT into biomass. In addition, it successfully elucidates that sulfate promotes higher cell growth than DBT and its presence in the medium reduces DBT desulfurization rates. A study using eight carbon sources suggests that ethanol and lactate yield higher cell growth and desulfurization rates than citrate, fructose, glucose, gluconate, glutamate, and glycerol.

In situ to in silico and back: elucidating the physiology and ecology of Geobacter spp. using genome-scale modelling

There is a wide diversity of unexplored metabolism encoded in the genomes of microorganisms that have an important environmental role. Genome-scale metabolic modelling enables the individual reactions that are encoded in annotated genomes to be organized into a coherent whole, which can then be used to predict metabolic fluxes that will optimize cell function under a range of conditions. In this Review, we summarize a series of studies in which genome-scale metabolic modelling of Geobacter spp. has resulted in an in-depth understanding of their central metabolism and ecology. A similar iterative modelling and experimental approach could accelerate elucidation of the physiology and ecology of other microorganisms inhabiting a diversity of environments, and could guide optimization of the practical applications of these species.

C4GEM, a genome-scale metabolic model to study C4 plant metabolism

Leaves of C(4) grasses (such as maize [Zea mays], sugarcane [Saccharum officinarum], and sorghum [Sorghum bicolor]) form a classical Kranz leaf anatomy. Unlike C(3) plants, where photosynthetic CO(2) fixation proceeds in the mesophyll (M), the fixation process in C(4) plants is distributed between two cell types, the M cell and the bundle sheath (BS) cell. Here, we develop a C(4) genome-scale model (C4GEM) for the investigation of flux distribution in M and BS cells during C(4) photosynthesis. C4GEM, to our knowledge, is the first large-scale metabolic model that encapsulates metabolic interactions between two different cell types. C4GEM is based on the Arabidopsis (Arabidopsis thaliana) model (AraGEM) but has been extended by adding reactions and transporters responsible to represent three different C(4) subtypes (NADP-ME [for malic enzyme], NAD-ME, and phosphoenolpyruvate carboxykinase). C4GEM has been validated for its ability to synthesize 47 biomass components and consists of 1,588 unique reactions, 1,755 metabolites, 83 interorganelle transporters, and 29 external transporters (including transport through plasmodesmata). Reactions in the common C(4) model have been associated with well-annotated C(4) species (NADP-ME subtypes): 3,557 genes in sorghum, 11,623 genes in maize, and 3,881 genes in sugarcane. The number of essential reactions not assigned to genes is 131, 135, and 156 in sorghum, maize, and sugarcane, respectively. Flux balance analysis was used to assess the metabolic activity in M and BS cells during C(4) photosynthesis. Our simulations were consistent with chloroplast proteomic studies, and C4GEM predicted the classical C(4) photosynthesis pathway and its major effect in organelle function in M and BS. The model also highlights differences in metabolic activities around photosystem I and photosystem II for three different C(4) subtypes. Effects of CO(2) leakage were also explored. C4GEM is a viable framework for in silico analysis of cell cooperation between M and BS cells during photosynthesis and can be used to explore C(4) plant metabolism.

von Bertalanffy 1.0: a COBRA toolbox extension to thermodynamically constrain metabolic models

MOTIVATION

In flux balance analysis of genome scale stoichiometric models of metabolism, the principal constraints are uptake or secretion rates, the steady state mass conservation assumption and reaction directionality. Here, we introduce an algorithmic pipeline for quantitative assignment of reaction directionality in multi-compartmental genome scale models based on an application of the second law of thermodynamics to each reaction. Given experimental or computationally estimated standard metabolite species Gibbs energy and metabolite concentrations, the algorithms bounds reaction Gibbs energy, which is transformed to in vivo pH, temperature, ionic strength and electrical potential.

RESULTS

This cross-platform MATLAB extension to the COnstraint-Based Reconstruction and Analysis (COBRA) toolbox is computationally efficient, extensively documented and open source.

AVAILABILITY

http://opencobra.sourceforge.net.

Large-scale in silico modeling of metabolic interactions between cell types in the human brain

Metabolic interactions between multiple cell types are difficult to model using existing approaches. Here we present a workflow that integrates gene expression data, proteomics data and literature-based manual curation to model human metabolism within and between different types of cells. Transport reactions are used to account for the transfer of metabolites between models of different cell types via the interstitial fluid. We apply the method to create models of brain energy metabolism that recapitulate metabolic interactions between astrocytes and various neuron types relevant to Alzheimer's disease. Analysis of the models identifies genes and pathways that may explain observed experimental phenomena, including the differential effects of the disease on cell types and regions of the brain. Constraint-based modeling can thus contribute to the study and analysis of multicellular metabolic processes in the human tissue microenvironment and provide detailed mechanistic insight into high-throughput data analysis.

Environments that induce synthetic microbial ecosystems

Interactions between microbial species are sometimes mediated by the exchange of small molecules, secreted by one species and metabolized by another. Both one-way (commensal) and two-way (mutualistic) interactions may contribute to complex networks of interdependencies. Understanding these interactions constitutes an open challenge in microbial ecology, with applications ranging from the human microbiome to environmental sustainability. In parallel to natural communities, it is possible to explore interactions in artificial microbial ecosystems, e.g. pairs of genetically engineered mutualistic strains. Here we computationally generate artificial microbial ecosystems without re-engineering the microbes themselves, but rather by predicting their growth on appropriately designed media. We use genome-scale stoichiometric models of metabolism to identify media that can sustain growth for a pair of species, but fail to do so for one or both individual species, thereby inducing putative symbiotic interactions. We first tested our approach on two previously studied mutualistic pairs, and on a pair of highly curated model organisms, showing that our algorithms successfully recapitulate known interactions, robustly predict new ones, and provide novel insight on exchanged molecules. We then applied our method to all possible pairs of seven microbial species, and found that it is always possible to identify putative media that induce commensalism or mutualism. Our analysis also suggests that symbiotic interactions may arise more readily through environmental fluctuations than genetic modifications. We envision that our approach will help generate microbe-microbe interaction maps useful for understanding microbial consortia dynamics and evolution, and for exploring the full potential of natural metabolic pathways for metabolic engineering applications.

Metabolic systems analysis to advance algal biotechnology

Algal fuel sources promise unsurpassed yields in a carbon neutral manner that minimizes resource competition between agriculture and fuel crops. Many challenges must be addressed before algal biofuels can be accepted as a component of the fossil fuel replacement strategy. One significant challenge is that the cost of algal fuel production must become competitive with existing fuel alternatives. Algal biofuel production presents the opportunity to fine-tune microbial metabolic machinery for an optimal blend of biomass constituents and desired fuel molecules. Genome-scale model-driven algal metabolic design promises to facilitate both goals by directing the utilization of metabolites in the complex, interconnected metabolic networks to optimize production of the compounds of interest. Network analysis can direct microbial development efforts towards successful strategies and enable quantitative fine-tuning of the network for optimal product yields while maintaining the robustness of the production microbe. Metabolic modeling yields insights into microbial function, guides experiments by generating testable hypotheses, and enables the refinement of knowledge on the specific organism. While the application of such analytical approaches to algal systems is limited to date, metabolic network analysis can improve understanding of algal metabolic systems and play an important role in expediting the adoption of new biofuel technologies.

In silico genome-scale metabolic analysis of Pseudomonas putida KT2440 for polyhydroxyalkanoate synthesis, degradation of aromatics and anaerobic survival

Genome-scale metabolic models have been appearing with increasing frequency and have been employed in a wide range of biotechnological applications as well as in biological studies. With the metabolic model as a platform, engineering strategies have become more systematic and focused, unlike the random shotgun approach used in the past. Here we present the genome-scale metabolic model of the versatile Gram-negative bacterium Pseudomonas putida, which has gained widespread interest for various biotechnological applications. With the construction of the genome-scale metabolic model of P. putida KT2440, PpuMBEL1071, we investigated various characteristics of P. putida, such as its capacity for synthesizing polyhydroxyalkanoates (PHA) and degrading aromatics. Although P. putida has been characterized as a strict aerobic bacterium, the physiological characteristics required to achieve anaerobic survival were investigated. Through analysis of PpuMBEL1071, extended survival of P. putida under anaerobic stress was achieved by introducing the ackA gene from Pseudomonas aeruginosa and Escherichia coli.

The virus as metabolic engineer

Recent genome-wide screens of host genetic requirements for viral infection have reemphasized the critical role of host metabolism in enabling the production of viral particles. In this review, we highlight the metabolic aspects of viral infection found in these studies, and focus on the opportunities these requirements present for metabolic engineers. In particular, the objectives and approaches that metabolic engineers use are readily comparable to the behaviors exhibited by viruses during infection. As a result, metabolic engineers have a unique perspective that could lead to novel and effective methods to combat viral infection.

Towards a unifying, systems biology understanding of large-scale cellular death and destruction caused by poorly liganded iron: Parkinson's, Huntington's, Alzheimer's, prions, bactericides, chemical toxicology and others as examples

Exposure to a variety of toxins and/or infectious agents leads to disease, degeneration and death, often characterised by circumstances in which cells or tissues do not merely die and cease to function but may be more or less entirely obliterated. It is then legitimate to ask the question as to whether, despite the many kinds of agent involved, there may be at least some unifying mechanisms of such cell death and destruction. I summarise the evidence that in a great many cases, one underlying mechanism, providing major stresses of this type, entails continuing and autocatalytic production (based on positive feedback mechanisms) of hydroxyl radicals via Fenton chemistry involving poorly liganded iron, leading to cell death via apoptosis (probably including via pathways induced by changes in the NF-κB system). While every pathway is in some sense connected to every other one, I highlight the literature evidence suggesting that the degenerative effects of many diseases and toxicological insults converge on iron dysregulation. This highlights specifically the role of iron metabolism, and the detailed speciation of iron, in chemical and other toxicology, and has significant implications for the use of iron chelating substances (probably in partnership with appropriate anti-oxidants) as nutritional or therapeutic agents in inhibiting both the progression of these mainly degenerative diseases and the sequelae of both chronic and acute toxin exposure. The complexity of biochemical networks, especially those involving autocatalytic behaviour and positive feedbacks, means that multiple interventions (e.g. of iron chelators plus antioxidants) are likely to prove most effective. A variety of systems biology approaches, that I summarise, can predict both the mechanisms involved in these cell death pathways and the optimal sites of action for nutritional or pharmacological interventions.

A detailed genome-wide reconstruction of mouse metabolism based on human Recon 1

BACKGROUND

Well-curated and validated network reconstructions are extremely valuable tools in systems biology. Detailed metabolic reconstructions of mammals have recently emerged, including human reconstructions. They raise the question if the various successful applications of microbial reconstructions can be replicated in complex organisms.

RESULTS

We mapped the published, detailed reconstruction of human metabolism (Recon 1) to other mammals. By searching for genes homologous to Recon 1 genes within mammalian genomes, we were able to create draft metabolic reconstructions of five mammals, including the mouse. Each draft reconstruction was created in compartmentalized and non-compartmentalized version via two different approaches. Using gap-filling algorithms, we were able to produce all cellular components with three out of four versions of the mouse metabolic reconstruction. We finalized a functional model by iterative testing until it passed a predefined set of 260 validation tests. The reconstruction is the largest, most comprehensive mouse reconstruction to-date, accounting for 1,415 genes coding for 2,212 gene-associated reactions and 1,514 non-gene-associated reactions.We tested the mouse model for phenotype prediction capabilities. The majority of predicted essential genes were also essential in vivo. However, our non-tissue specific model was unable to predict gene essentiality for many of the metabolic genes shown to be essential in vivo. Our knockout simulation of the lipoprotein lipase gene correlated well with experimental results, suggesting that softer phenotypes can also be simulated.

CONCLUSIONS

We have created a high-quality mouse genome-scale metabolic reconstruction, iMM1415 (Mus Musculus, 1415 genes). We demonstrate that the mouse model can be used to perform phenotype simulations, similar to models of microbe metabolism. Since the mouse is an important experimental organism, this model should become an essential tool for studying metabolic phenotypes in mice, including outcomes from drug screening.

Predicting novel pathways in genome-scale metabolic networks

Elementary-modes analysis has become a well-established theoretical tool in metabolic pathway analysis. It allows one to decompose complex metabolic networks into the smallest functional entities, which can be interpreted as biochemical pathways. This analysis has, in medium-size metabolic networks, led to the successful theoretical prediction of hitherto unknown pathways. For illustration, we discuss the example of the phosphoenolpyruvate-glyoxylate cycle in Escherichia coli. Elementary-modes analysis meets with the problem of combinatorial explosion in the number of pathways with increasing system size, which has hampered scaling it up to genome-wide models. We present a novel approach to overcoming this obstacle. That approach is based on elementary flux patterns, which are defined as sets of reactions representing the basic routes through a particular subsystem that are compatible with admissible fluxes in a (possibly) much larger metabolic network. The subsystem can be made up by reactions in which we are interested in, for example, reactions producing a certain metabolite. This allows one to predict novel metabolic pathways in genome-scale networks.

Functional characterization of alternate optimal solutions of Escherichia coli's transcriptional and translational machinery

The constraint-based reconstruction and analysis approach has recently been extended to describe Escherichia coli's transcriptional and translational machinery. Here, we introduce the concept of reaction coupling to represent the dependency between protein synthesis and utilization. These coupling constraints lead to a significant contraction of the feasible set of steady-state fluxes. The subset of alternate optimal solutions (AOS) consistent with maximal ribosome production was calculated. The majority of transcriptional and translational reactions were active for all of these AOS, showing that the network has a low degree of redundancy. Furthermore, all calculated AOS contained the qualitative expression of at least 92% of the known essential genes. Principal component analysis of AOS demonstrated that energy currencies (ATP, GTP, and phosphate) dominate the network's capability to produce ribosomes. Additionally, we identified regulatory control points of the network, which include the transcription reactions of sigma70 (RpoD) as well as that of a degradosome component (Rne) and of tRNA charging (ValS). These reactions contribute significant variance among AOS. These results show that constraint-based modeling can be applied to gain insight into the systemic properties of E. coli's transcriptional and translational machinery.

High-throughput generation, optimization and analysis of genome-scale metabolic models

Genome-scale metabolic models have proven to be valuable for predicting organism phenotypes from genotypes. Yet efforts to develop new models are failing to keep pace with genome sequencing. To address this problem, we introduce the Model SEED, a web-based resource for high-throughput generation, optimization and analysis of genome-scale metabolic models. The Model SEED integrates existing methods and introduces techniques to automate nearly every step of this process, taking approximately 48 h to reconstruct a metabolic model from an assembled genome sequence. We apply this resource to generate 130 genome-scale metabolic models representing a taxonomically diverse set of bacteria. Twenty-two of the models were validated against available gene essentiality and Biolog data, with the average model accuracy determined to be 66% before optimization and 87% after optimization.

Modeling core metabolism in cancer cells: surveying the topology underlying the Warburg effect

BACKGROUND

Alterations on glucose consumption and biosynthetic activity of amino acids, lipids and nucleotides are metabolic changes for sustaining cell proliferation in cancer cells. Irrevocable evidence of this fact is the Warburg effect which establishes that cancer cells prefers glycolysis over oxidative phosphorylation to generate ATP. Regulatory action over metabolic enzymes has opened a new window for designing more effective anti-cancer treatments. This enterprise is not trivial and the development of computational models that contribute to identifying potential enzymes for breaking the robustness of cancer cells is a priority.

METHODOLOGY/PRINCIPAL FINDINGS

This work presents a constraint-base modeling of the most experimentally studied metabolic pathways supporting cancer cells: glycolysis, TCA cycle, pentose phosphate, glutaminolysis and oxidative phosphorylation. To evaluate its predictive capacities, a growth kinetics study for Hela cell lines was accomplished and qualitatively compared with in silico predictions. Furthermore, based on pure computational criteria, we concluded that a set of enzymes (such as lactate dehydrogenase and pyruvate dehydrogenase) perform a pivotal role in cancer cell growth, findings supported by an experimental counterpart.

CONCLUSIONS/SIGNIFICANCE

Alterations on metabolic activity are crucial to initiate and sustain cancer phenotype. In this work, we analyzed the phenotype capacities emerged from a constructed metabolic network conformed by the most experimentally studied pathways sustaining cancer cell growth. Remarkably, in silico model was able to resemble the physiological conditions in cancer cells and successfully identified some enzymes currently studied by its therapeutic effect. Overall, we supplied evidence that constraint-based modeling constitutes a promising computational platform to: 1) integrate high throughput technology and establish a crosstalk between experimental validation and in silico prediction in cancer cell phenotype; 2) explore the fundamental metabolic mechanism that confers robustness in cancer; and 3) suggest new metabolic targets for anticancer treatments. All these issues being central to explore cancer cell metabolism from a systems biology perspective.

Integration of metabolic databases for the reconstruction of genome-scale metabolic networks

BACKGROUND

Genome-scale metabolic reconstructions have been recognised as a valuable tool for a variety of applications ranging from metabolic engineering to evolutionary studies. However, the reconstruction of such networks remains an arduous process requiring a high level of human intervention. This process is further complicated by occurrences of missing or conflicting information and the absence of common annotation standards between different data sources.

RESULTS

In this article, we report a semi-automated methodology aimed at streamlining the process of metabolic network reconstruction by enabling the integration of different genome-wide databases of metabolic reactions. We present results obtained by applying this methodology to the metabolic network of the plant Arabidopsis thaliana. A systematic comparison of compounds and reactions between two genome-wide databases allowed us to obtain a high-quality core consensus reconstruction, which was validated for stoichiometric consistency. A lower level of consensus led to a larger reconstruction, which has a lower quality standard but provides a baseline for further manual curation.

CONCLUSION

This semi-automated methodology may be applied to other organisms and help to streamline the process of genome-scale network reconstruction in order to accelerate the transfer of such models to applications.

BiGG: a Biochemical Genetic and Genomic knowledgebase of large scale metabolic reconstructions

BACKGROUND

Genome-scale metabolic reconstructions under the Constraint Based Reconstruction and Analysis (COBRA) framework are valuable tools for analyzing the metabolic capabilities of organisms and interpreting experimental data. As the number of such reconstructions and analysis methods increases, there is a greater need for data uniformity and ease of distribution and use.

DESCRIPTION

We describe BiGG, a knowledgebase of Biochemically, Genetically and Genomically structured genome-scale metabolic network reconstructions. BiGG integrates several published genome-scale metabolic networks into one resource with standard nomenclature which allows components to be compared across different organisms. BiGG can be used to browse model content, visualize metabolic pathway maps, and export SBML files of the models for further analysis by external software packages. Users may follow links from BiGG to several external databases to obtain additional information on genes, proteins, reactions, metabolites and citations of interest.

CONCLUSIONS

BiGG addresses a need in the systems biology community to have access to high quality curated metabolic models and reconstructions. It is freely available for academic use at http://bigg.ucsd.edu.

The biomass objective function

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

Use of genome-scale metabolic models for understanding microbial physiology

The exploitation of microorganisms in industrial, medical, food and environmental biotechnology requires a comprehensive understanding of their physiology. The availability of genome sequences and accumulation of high-throughput data allows gaining understanding of microbial physiology at the systems level, and genome-scale metabolic models represent a valuable framework for integrative analysis of metabolism of microorganisms. Genome-scale metabolic models are reconstructed based on a combination of genome sequence information and detailed biochemical information, and these reconstructed models can be used for analyzing and simulating the operation of metabolism in response to different stimuli. Here we discuss the requirement for having detailed physiological insight in order to exploit microorganisms for production of fuels, chemicals and pharmaceuticals. We further describe the reconstruction process of genome-scale metabolic models and different algorithms that can be used to apply these models to gain improved insight into microbial physiology.

Integrated stoichiometric, thermodynamic and kinetic modelling of steady state metabolism

The quantitative analysis of biochemical reactions and metabolites is at frontier of biological sciences. The recent availability of high-throughput technology data sets in biology has paved the way for new modelling approaches at various levels of complexity including the metabolome of a cell or an organism. Understanding the metabolism of a single cell and multi-cell organism will provide the knowledge for the rational design of growth conditions to produce commercially valuable reagents in biotechnology. Here, we demonstrate how equations representing steady state mass conservation, energy conservation, the second law of thermodynamics, and reversible enzyme kinetics can be formulated as a single system of linear equalities and inequalities, in addition to linear equalities on exponential variables. Even though the feasible set is non-convex, the reformulation is exact and amenable to large-scale numerical analysis, a prerequisite for computationally feasible genome scale modelling. Integrating flux, concentration and kinetic variables in a unified constraint-based formulation is aimed at increasing the quantitative predictive capacity of flux balance analysis. Incorporation of experimental and theoretical bounds on thermodynamic and kinetic variables ensures that the predicted steady state fluxes are both thermodynamically and biochemically feasible. The resulting in silico predictions are tested against fluxomic data for central metabolism in Escherichia coli and compare favourably with in silico prediction by flux balance analysis.

Metabolic systems biology

The first full genome sequences were established in the mid-1990s. Shortly thereafter, genome-scale metabolic network reconstructions appeared. Since that time, we have witnessed an exponential growth in their number and uses. Here I discuss, from a personal point of view, four topics: (1) the placement of metabolic systems biology in the context of broader scientific developments, (2) its foundational concepts, (3) some of its current uses, and (4) some of the expected future developments in the field.

Quantitative assignment of reaction directionality in constraint-based models of metabolism: application to Escherichia coli

Constraint-based modeling is an approach for quantitative prediction of net reaction flux in genome-scale biochemical networks. In vivo, the second law of thermodynamics requires that net macroscopic flux be forward, when the transformed reaction Gibbs energy is negative. We calculate the latter by using (i) group contribution estimates of metabolite species Gibbs energy, combined with (ii) experimentally measured equilibrium constants. In an application to a genome-scale stoichiometric model of Escherichia coli metabolism, iAF1260, we demonstrate that quantitative prediction of reaction directionality is increased in scope and accuracy by integration of both data sources, transformed appropriately to in vivo pH, temperature and ionic strength. Comparison of quantitative versus qualitative assignment of reaction directionality in iAF1260, assuming an accommodating reactant concentration range of 0.02-20mM, revealed that quantitative assignment leads to a low false positive, but high false negative, prediction of effectively irreversible reactions. The latter is partly due to the uncertainty associated with group contribution estimates. We also uncovered evidence that the high intracellular concentration of glutamate in E. coli may be essential to direct otherwise thermodynamically unfavorable essential reactions, such as the leucine transaminase reaction, in an anabolic direction.