Effect of Escherichia coli biomass composition on central metabolic fluxes predicted by a stoichiometric model

Reconstruction and metabolic profiling of the genome-scale metabolic network model of Pseudomonas stutzeri A1501

Pseudomonas stutzeri A1501 is a non-fluorescent denitrifying bacteria that belongs to the gram-negative bacterial group. As a prominent strain in the fields of agriculture and bioengineering, there is still a lack of comprehensive understanding regarding its metabolic capabilities, specifically in terms of central metabolism and substrate utilization. Therefore, further exploration and extensive studies are required to gain a detailed insight into these aspects. This study reconstructed a genome-scale metabolic network model for P. stutzeri A1501 and conducted extensive curations, including correcting energy generation cycles, respiratory chains, and biomass composition. The final model, iQY1018, was successfully developed, covering more genes and reactions and having higher prediction accuracy compared with the previously published model iPB890. The substrate utilization ability of 71 carbon sources was investigated by BIOLOG experiment and was utilized to validate the model quality. The model prediction accuracy of substrate utilization for P. stutzeri A1501 reached 90 %. The model analysis revealed its new ability in central metabolism and predicted that the strain is a suitable chassis for the production of Acetyl CoA-derived products. This work provides an updated, high-quality model of P. stutzeri A1501for further research and will further enhance our understanding of the metabolic capabilities.

Mitigating biomass composition uncertainties in flux balance analysis using ensemble representations

The biomass equation is a critical component in genome-scale metabolic models (GEMs): it is used as the de facto objective function in flux balance analysis (FBA). This equation accounts for the quantities of all known biomass precursors that are required for cell growth based on the macromolecular and monomer compositions measured at certain conditions. However, it is often reported that the macromolecular composition of cells could change across different environmental conditions and thus the use of the same single biomass equation in FBA, under multiple conditions, is questionable. Herein, we first investigated the qualitative and quantitative variations of macromolecular compositions of three representative host organisms, Escherichia coli, Saccharomyces cerevisiae and Cricetulus griseus, across different environmental/genetic variations. While macromolecular building blocks such as RNA, protein, and lipid composition vary notably, changes in fundamental biomass monomer units such as nucleotides and amino acids are not appreciable. We also observed that flux predictions through FBA is quite sensitive to macromolecular compositions but not the monomer compositions. Based on these observations, we propose ensemble representations of biomass equation in FBA to account for the natural variation of cellular constituents. Such ensemble representations of biomass better predicted the flux through anabolic reactions as it allows for the flexibility in the biosynthetic demands of the cells. The current study clearly highlights that certain component of the biomass equation indeed vary across different conditions, and the ensemble representation of biomass equation in FBA by accounting for such natural variations could avoid inaccuracies that may arise from in silico simulations.

H+-Translocating Membrane-Bound Pyrophosphatase from Rhodospirillum rubrum Fuels Escherichia coli Cells via an Alternative Pathway for Energy Generation

Inorganic pyrophosphatases (PPases) catalyze an essential reaction, namely, the hydrolysis of PP_i, which is formed in large quantities as a side product of numerous cellular reactions. In the majority of living species, PP_i hydrolysis is carried out by soluble cytoplasmic PPase (S-PPases) with the released energy dissipated in the form of heat. In Rhodospirillum rubrum, part of this energy can be conserved by proton-pumping pyrophosphatase (H⁺-PPase^Rru) in the form of a proton electrochemical gradient for further ATP synthesis. Here, the codon-harmonized gene hppa^Rru encoding H⁺-PPase^Rru was expressed in the Escherichia coli chromosome. We demonstrate, for the first time, that H⁺-PPase^Rru complements the essential native S-PPase in E. coli cells. ¹³C-MFA confirmed that replacing native PPase to H⁺-PPase^Rru leads to the re-distribution of carbon fluxes; a statistically significant 36% decrease in tricarboxylic acid (TCA) cycle fluxes was found compared with wild-type E. coli MG1655. Such a flux re-distribution can indicate the presence of an additional method for energy generation (e.g., ATP), which can be useful for the microbiological production of a number of compounds, the biosynthesis of which requires the consumption of ATP.

Genome-Scale Reconstruction of Microbial Dynamic Phenotype: Successes and Challenges

This review is a part of the SI 'Genome-Scale Modeling of Microorganisms in the Real World'. The goal of GEM is the accurate prediction of the phenotype from its respective genotype under specified environmental conditions. This review focuses on the dynamic phenotype; prediction of the real-life behaviors of microorganisms, such as cell proliferation, dormancy, and mortality; balanced and unbalanced growth; steady-state and transient processes; primary and secondary metabolism; stress responses; etc. Constraint-based metabolic reconstructions were successfully started two decades ago as FBA, followed by more advanced models, but this review starts from the earlier nongenomic predecessors to show that some GEMs inherited the outdated biokinetic frameworks compromising their performances. The most essential deficiencies are: (i) an inadequate account of environmental conditions, such as various degrees of nutrients limitation and other factors shaping phenotypes; (ii) a failure to simulate the adaptive changes of MMCC (MacroMolecular Cell Composition) in response to the fluctuating environment; (iii) the misinterpretation of the SGR (Specific Growth Rate) as either a fixed constant parameter of the model or independent factor affecting the conditional expression of macromolecules; (iv) neglecting stress resistance as an important objective function; and (v) inefficient experimental verification of GEM against simple growth (constant MMCC and SGR) data. Finally, we propose several ways to improve GEMs, such as replacing the outdated Monod equation with the SCM (Synthetic Chemostat Model) that establishes the quantitative relationships between primary and secondary metabolism, growth rate and stress resistance, process kinetics, and cell composition.

Understanding and mathematical modelling of cellular resource allocation in microorganisms: a comparative synthesis

BACKGROUND

The rising consensus that the cell can dynamically allocate its resources provides an interesting angle for discovering the governing principles of cell growth and metabolism. Extensive efforts have been made in the past decade to elucidate the relationship between resource allocation and phenotypic patterns of microorganisms. Despite these exciting developments, there is still a lack of explicit comparison between potentially competing propositions and a lack of synthesis of inter-related proposals and findings.

RESULTS

In this work, we have reviewed resource allocation-derived principles, hypotheses and mathematical models to recapitulate important achievements in this area. In particular, the emergence of resource allocation phenomena is deciphered by the putative tug of war between the cellular objectives, demands and the supply capability. Competing hypotheses for explaining the most-studied phenomenon arising from resource allocation, i.e. the overflow metabolism, have been re-examined towards uncovering the potential physiological root cause. The possible link between proteome fractions and the partition of the ribosomal machinery has been analysed through mathematical derivations. Finally, open questions are highlighted and an outlook on the practical applications is provided. It is the authors' intention that this review contributes to a clearer understanding of the role of resource allocation in resolving bacterial growth strategies, one of the central questions in microbiology.

CONCLUSIONS

We have shown the importance of resource allocation in understanding various aspects of cellular systems. Several important questions such as the physiological root cause of overflow metabolism and the correct interpretation of 'protein costs' are shown to remain open. As the understanding of the mechanisms and utility of resource application in cellular systems further develops, we anticipate that mathematical modelling tools incorporating resource allocation will facilitate the circuit-host design in synthetic biology.

Computation of condition-dependent proteome allocation reveals variability in the macro and micro nutrient requirements for growth

Sustaining a robust metabolic network requires a balanced and fully functioning proteome. In addition to amino acids, many enzymes require cofactors (coenzymes and engrafted prosthetic groups) to function properly. Extensively validated resource allocation models, such as genome-scale models of metabolism and gene expression (ME-models), have the ability to compute an optimal proteome composition underlying a metabolic phenotype, including the provision of all required cofactors. Here we apply the ME-model for Escherichia coli K-12 MG1655 to computationally examine how environmental conditions change the proteome and its accompanying cofactor usage. We found that: (1) The cofactor requirements computed by the ME-model mostly agree with the standard biomass objective function used in models of metabolism alone (M-models); (2) ME-model computations reveal non-intuitive variability in cofactor use under different growth conditions; (3) An analysis of ME-model predicted protein use in aerobic and anaerobic conditions suggests an enrichment in the use of peroxyl scavenging acids in the proteins used to sustain aerobic growth; (4) The ME-model could describe how limitation in key protein components affect the metabolic state of E. coli. Genome-scale models have thus reached a level of sophistication where they reveal intricate properties of functional proteomes and how they support different E. coli lifestyles.

Escherichia coli is capable of growing in many environments, each of which requires a different collection of enzymes to metabolize the nutrients within that environment. Each individual enzyme requires its own set of amino acids and oftentimes cofactors, which are accessory molecules essential for the enzyme to function. Thus, the composition of the micronutrients (amino acids, cofactors, etc.) within a cell will differ depending on its metabolic needs. The presented work is the first effort to employ metabolic models to probe the connection between E. coli’s diverse growth environments and its biomass composition. We first show how differences in model-predicted enzyme use for aerobic or anaerobic growth results in distinct amino acid and cofactor usage. Alternatively, we show that the metabolic models can predict how modifying the cell’s biomass composition will affect growth. For example, by modeling the exposure of E. coli to trimethoprim or sulfamethoxazole—two antibiotics that target folate (vitamin B9) synthesis—we predicted how E. coli could adapt to grow under folate-limited conditions. This work demonstrates how models can be used to study antibiotic resistance of drugs that target amino acid or cofactor synthesis.

Comparative 13 C-metabolic flux analysis indicates elevation of ATP regeneration, carbon dioxide, and heat production in industrial Saccharomyces cerevisiae strains

BACKGROUND

Various industrial Saccharomyces cerevisiae strains are used for specific processes, such as sake, wine brewing and bread making. Understanding mechanisms underlying the fermentation performance of these strains would be useful for further engineering of the S. cerevisiae metabolism. However, the relationship between the fermentation performance, intra-cellular metabolic states, and other phenotypic characteristics of industrial yeasts is still unclear. In this study, ¹³ C-metabolic flux analysis of four diploid yeast strains-laboratory, sake, bread, and wine yeasts-was conducted.

RESULTS

While the Crabtree effect was observed for all strains, the metabolic flux level of glycolysis was elevated in bread and sake yeast. Furthermore, increased flux levels of the TCA cycle were commonly observed in the three industrial strains. The specific rates of CO₂ production, net ATP regeneration, and metabolic heat generation estimated from the metabolic flux distribution were two to three times greater than those of the laboratory strain. The elevation in metabolic heat generation was correlated with the tolerance to low-temperature stress.

CONCLUSION

These results indicate that the metabolic flux distribution of sake and bread yeast strains contributes to faster production of ethanol and CO₂ . It is also suggested that the generation of metabolic heat is preferable under the actual industrial fermentation conditions.

Proteome Regulation Patterns Determine Escherichia coli Wild-Type and Mutant Phenotypes

It is generally recognized that proteins constitute the key cellular component in shaping microbial phenotypes. Due to limited cellular resources and space, optimal allocation of proteins is crucial for microbes to facilitate maximum proliferation rates while allowing a flexible response to environmental changes. To account for the growth condition-dependent proteome in the constraint-based metabolic modeling of Escherichia coli, we consolidated a coarse-grained protein allocation approach with the explicit consideration of enzymatic constraints on reaction fluxes. Besides representing physiologically relevant wild-type phenotypes and flux distributions, the resulting protein allocation model (PAM) advances the predictability of the metabolic responses to genetic perturbations. A main driver of mutant phenotypes was ascribed to inherited regulation patterns in protein distribution among metabolic enzymes. Moreover, the PAM correctly reflected metabolic responses to an augmented protein burden imposed by the heterologous expression of green fluorescent protein. In summary, we were able to model the effects of important and frequently applied metabolic engineering approaches on microbial metabolism. Therefore, we want to promote the integration of protein allocation constraints into classical constraint-based models to foster their predictive capabilities and application for strain analysis and engineering purposes.

IMPORTANCE Predictive metabolic models are important, e.g., for generating biological knowledge and designing microbes with superior performance for target compound production. Yet today’s whole-cell models either show insufficient predictive capabilities or are computationally too expensive to be applied to metabolic engineering purposes. By linking the inherent genotype-phenotype relationship to a complete representation of the proteome, the PAM advances the accuracy of simulated phenotypes and intracellular flux distributions of E. coli. Being equally computationally lightweight as classical stoichiometric models and allowing for the application of established in silico tools, the PAM and related simulation approaches will foster the use of a model-driven metabolic research. Applications range from the investigation of mechanisms of microbial evolution to the determination of optimal strain design strategies in metabolic engineering, thus supporting basic scientists and engineers alike.

Addressing uncertainty in genome-scale metabolic model reconstruction and analysis

The reconstruction and analysis of genome-scale metabolic models constitutes a powerful systems biology approach, with applications ranging from basic understanding of genotype-phenotype mapping to solving biomedical and environmental problems. However, the biological insight obtained from these models is limited by multiple heterogeneous sources of uncertainty, which are often difficult to quantify. Here we review the major sources of uncertainty and survey existing approaches developed for representing and addressing them. A unified formal characterization of these uncertainties through probabilistic approaches and ensemble modeling will facilitate convergence towards consistent reconstruction pipelines, improved data integration algorithms, and more accurate assessment of predictive capacity.

Reconstructing organisms in silico: genome-scale models and their emerging applications

Escherichia coli is considered to be the best-known microorganism given the large number of published studies detailing its genes, its genome and the biochemical functions of its molecular components. This vast literature has been systematically assembled into a reconstruction of the biochemical reaction networks that underlie E. coli's functions, a process which is now being applied to an increasing number of microorganisms. Genome-scale reconstructed networks are organized and systematized knowledge bases that have multiple uses, including conversion into computational models that interpret and predict phenotypic states and the consequences of environmental and genetic perturbations. These genome-scale models (GEMs) now enable us to develop pan-genome analyses that provide mechanistic insights, detail the selection pressures on proteome allocation and address stress phenotypes. In this Review, we first discuss the overall development of GEMs and their applications. Next, we review the evolution of the most complete GEM that has been developed to date: the E. coli GEM. Finally, we explore three emerging areas in genome-scale modelling of microbial phenotypes: collections of strain-specific models, metabolic and macromolecular expression models, and simulation of stress responses.

Model metabolic strategy for heterotrophic bacteria in the cold ocean based on Colwellia psychrerythraea 34H

Colwellia psychrerythraea 34H is a model psychrophilic bacterium found in the cold ocean-polar sediments, sea ice, and the deep sea. Although the genomes of such psychrophiles have been sequenced, their metabolic strategies at low temperature have not been quantified. We measured the metabolic fluxes and gene expression of 34H at 4 °C (the mean global-ocean temperature and a normal-growth temperature for 34H), making comparative analyses at room temperature (above its upper-growth temperature of 18 °C) and with mesophilic Escherichia coli When grown at 4 °C, 34H utilized multiple carbon substrates without catabolite repression or overflow byproducts; its anaplerotic pathways increased flux network flexibility and enabled CO₂ fixation. In glucose-only medium, the Entner-Doudoroff (ED) pathway was the primary glycolytic route; in lactate-only medium, gluconeogenesis and the glyoxylate shunt became active. In comparison, E. coli, cold stressed at 4 °C, had rapid glycolytic fluxes but no biomass synthesis. At their respective normal-growth temperatures, intracellular concentrations of TCA cycle metabolites (α-ketoglutarate, succinate, malate) were 4-17 times higher in 34H than in E. coli, while levels of energy molecules (ATP, NADH, NADPH) were 10- to 100-fold lower. Experiments with E. coli mutants supported the thermodynamic advantage of the ED pathway at cold temperature. Heat-stressed 34H at room temperature (2 hours) revealed significant down-regulation of genes associated with glycolytic enzymes and flagella, while 24 hours at room temperature caused irreversible cellular damage. We suggest that marine heterotrophic bacteria in general may rely upon simplified metabolic strategies to overcome thermodynamic constraints and thrive in the cold ocean.

Population FBA predicts metabolic phenotypes in yeast

Using protein counts sampled from single cell proteomics distributions to constrain fluxes through a genome-scale model of metabolism, Population flux balance analysis (Population FBA) successfully described metabolic heterogeneity in a population of independent Escherichia coli cells growing in a defined medium. We extend the methodology to account for correlations in protein expression arising from the co-regulation of genes and apply it to study the growth of independent Saccharomyces cerevisiae cells in two different growth media. We find the partitioning of flux between fermentation and respiration predicted by our model agrees with recent ¹³C fluxomics experiments, and that our model largely recovers the Crabtree effect (the experimentally known bias among certain yeast species toward fermentation with the production of ethanol even in the presence of oxygen), while FBA without proteomics constraints predicts respirative metabolism almost exclusively. The comparisons to the ¹³C study showed improvement upon inclusion of the correlations and motivated a technique to systematically identify inconsistent kinetic parameters in the literature. The minor secretion fluxes for glycerol and acetate are underestimated by our method, which indicate a need for further refinements to the metabolic model. For yeast cells grown in synthetic defined (SD) medium, the calculated broad distribution of growth rates matches experimental observations from single cell studies, and we characterize several metabolic phenotypes within our modeled populations that make use of diverse pathways. Fast growing yeast cells are predicted to perform significant amount of respiration, use serine-glycine cycle and produce ethanol in mitochondria as opposed to slow growing cells. We use a genetic algorithm to determine the proteomics constraints necessary to reproduce the growth rate distributions seen experimentally. We find that a core set of 51 constraints are essential but that additional constraints are still necessary to recover the observed growth rate distribution in SD medium.

No two living cells are exactly the same. Even cells from a clonal population with identical genomes living in the same environment will express proteins in different numbers simply due to the random nature of the chemistry involved in gene expression. The consequences of this stochastic gene expression are complex and not well understood, especially at the level of large reaction networks like metabolism. Here we investigate how variability in the copy numbers of metabolic enzymes affects how individual cells extract nourishment from their environment and grow. We model 100,000 independent yeast cells, each with their own set of enzyme copy numbers sampled from experimental distributions, and use flux balance analysis (FBA) to compute the optimal way that each cell can use its metabolic pathways—an approach we dubbed Population FBA. We find that enzyme variability gives rise to a wide distribution of growth rates, and several metabolic phenotypes—subpopulations relying on diverse metabolic pathways. Most importantly, we compare the predicted fluxes through the different pathways to experimental values; we find that Population FBA is able to correctly predict Crabtree effect, while traditional FBA, which lacks the proteomics constraints our method imposes, differs both qualitatively and quantitatively from experiment.

lumpGEM: Systematic generation of subnetworks and elementally balanced lumped reactions for the biosynthesis of target metabolites

In the post-genomic era, Genome-scale metabolic networks (GEMs) have emerged as invaluable tools to understand metabolic capabilities of organisms. Different parts of these metabolic networks are defined as subsystems/pathways, which are sets of functional roles to implement a specific biological process or structural complex, such as glycolysis and TCA cycle. Subsystem/pathway definition is also employed to delineate the biosynthetic routes that produce biomass building blocks. In databases, such as MetaCyc and SEED, these representations are composed of linear routes from precursors to target biomass building blocks. However, this approach cannot capture the nested, complex nature of GEMs. Here we implemented an algorithm, lumpGEM, which generates biosynthetic subnetworks composed of reactions that can synthesize a target metabolite from a set of defined core precursor metabolites. lumpGEM captures balanced subnetworks, which account for the fate of all metabolites along the synthesis routes, thus encapsulating reactions from various subsystems/pathways to balance these metabolites in the metabolic network. Moreover, lumpGEM collapses these subnetworks into elementally balanced lumped reactions that specify the cost of all precursor metabolites and cofactors. It also generates alternative subnetworks and lumped reactions for the same metabolite, accounting for the flexibility of organisms. lumpGEM is applicable to any GEM and any target metabolite defined in the network. Lumped reactions generated by lumpGEM can be also used to generate properly balanced reduced core metabolic models.

Stoichiometric models have been used in the area of metabolic engineering and systems biology for many decades. The early examples of these models include simplified ad hoc built metabolic pathways, and biomass compositions. The development of genome scale models (GEMs) brought a standard to metabolic network modeling. However, the vast amount of detailed biochemistry in GEMs makes it necessary to develop methods to manage the complexity in them. In this study, we developed lumpGEM, a tool that can systematically identify subnetworks from metabolic networks that can perform certain tasks, such as biosynthesis of a biomass building block and any other target metabolite. By generating alternative subnetworks, lumpGEM also accounts for the redundancy in metabolic networks. We applied lumpGEM on latest E. coli GEM iJO1366 and identified subnetworks/lumped reactions for every biomass building block defined in its biomass formulation. We also compared the results from lumpGEM with existing knowledge in the literature. The lumped reactions generated by lumpGEM can be used to generate consistently reduced metabolic network models.

Improving Escherichia coli membrane integrity and fatty acid production by expression tuning of FadL and OmpF

Background

Construction of microbial biocatalysts for the production of biorenewables at economically viable yields and titers is frequently hampered by product toxicity. Membrane damage is often deemed as the principal mechanism of this toxicity, particularly in regards to decreased membrane integrity. Previous studies have attempted to engineer the membrane with the goal of increasing membrane integrity. However, most of these works focused on engineering of phospholipids and efforts to identify membrane proteins that can be targeted to improve fatty acid production have been unsuccessful.

Results

Here we show that deletion of outer membrane protein ompF significantly increased membrane integrity, fatty acid tolerance and fatty acid production, possibly due to prevention of re-entry of short chain fatty acids. In contrast, deletion of fadL resulted in significantly decreased membrane integrity and fatty acid production. Consistently, increased expression of fadL remarkably increased membrane integrity and fatty acid tolerance while also increasing the final fatty acid titer. This 34% increase in the final fatty acid titer was possibly due to increased membrane lipid biosynthesis. Tuning of fadL expression showed that there is a positive relationship between fadL abundance and fatty acid production. Combinatorial deletion of ompF and increased expression of fadL were found to have an additive role in increasing membrane integrity, and was associated with a 53% increase the fatty acid titer, to 2.3 g/L.

Conclusions

These results emphasize the importance of membrane proteins for maintaining membrane integrity and production of biorenewables, such as fatty acids, which expands the targets for membrane engineering.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-017-0650-8) contains supplementary material, which is available to authorized users.

Integration of Biomass Formulations of Genome-Scale Metabolic Models with Experimental Data Reveals Universally Essential Cofactors in Prokaryotes

The composition of a cell in terms of macromolecular building blocks and other organic molecules underlies the metabolic needs and capabilities of a species. Although some core biomass components such as nucleic acids and proteins are evident for most species, the essentiality of the pool of other organic molecules, especially cofactors and prosthetic groups, is yet unclear. Here we integrate biomass compositions from 71 manually curated genome-scale models, 33 large-scale gene essentiality datasets, enzyme-cofactor association data and a vast array of publications, revealing universally essential cofactors for prokaryotic metabolism and also others that are specific for phylogenetic branches or metabolic modes. Our results revise predictions of essential genes in Klebsiella pneumoniae and identify missing biosynthetic pathways in models of Mycobacterium tuberculosis. This work provides fundamental insights into the essentiality of organic cofactors and has implications for minimal cell studies as well as for modeling genotype-phenotype relations in prokaryotic metabolic networks.

•

Seventy one biomass equations of manually curated genome-scale metabolic models are compared.

•

Eight classes of universally essential prokaryotic organic cofactors are proposed.

•

Conditionally essential organic cofactors are presented and discussed.

•

Gene essentiality predictions for Klebsiella pneumoniae are revised.

•

A missing essential pathway in models of Mycobacterium tuberculosis is predicted.

Genome-scale reconstruction of the metabolic network in Pseudomonas stutzeri A1501

Pseudomonas stutzeri A1501 is an endophytic bacterium capable of nitrogen fixation. This strain has been isolated from the rice rhizosphere and provides the plant with fixed nitrogen and phytohormones. These interesting features encouraged us to study the metabolism of this microorganism at the systems-level. In this work, we present the first genome-scale metabolic model (iPB890) for P. stutzeri, involving 890 genes, 1135 reactions, and 813 metabolites. A combination of automatic and manual approaches was used in the reconstruction process. Briefly, using the metabolic networks of Pseudomonas aeruginosa and Pseudomonas putida as templates, a draft metabolic network of P. stutzeri was reconstructed. Then, the draft network was driven through an iterative and curative process of gap filling. In the next step, the model was evaluated using different experimental data such as specific growth rate, Biolog substrate utilization data and other experimental observations. In most of the evaluation cases, the model was successful in correctly predicting the cellular phenotypes. Thus, we posit that the iPB890 model serves as a suitable platform to explore the metabolism of P. stutzeri.

Constrained Allocation Flux Balance Analysis

New experimental results on bacterial growth inspire a novel top-down approach to study cell metabolism, combining mass balance and proteomic constraints to extend and complement Flux Balance Analysis. We introduce here Constrained Allocation Flux Balance Analysis, CAFBA, in which the biosynthetic costs associated to growth are accounted for in an effective way through a single additional genome-wide constraint. Its roots lie in the experimentally observed pattern of proteome allocation for metabolic functions, allowing to bridge regulation and metabolism in a transparent way under the principle of growth-rate maximization. We provide a simple method to solve CAFBA efficiently and propose an “ensemble averaging” procedure to account for unknown protein costs. Applying this approach to modeling E. coli metabolism, we find that, as the growth rate increases, CAFBA solutions cross over from respiratory, growth-yield maximizing states (preferred at slow growth) to fermentative states with carbon overflow (preferred at fast growth). In addition, CAFBA allows for quantitatively accurate predictions on the rate of acetate excretion and growth yield based on only 3 parameters determined by empirical growth laws.

The intracellular protein levels of exponentially growing bacteria are known to vary strongly with growth conditions, as described by quantitative “growth laws”. This work introduces a computational genome-scale framework (Constrained Allocation Flux Balance Analysis, CAFBA) which incorporates growth laws into canonical Flux Balance Analysis. Upon introducing 3 parameters based on established growth laws for E. coli, CAFBA accurately reproduces empirical results on the growth-rate dependent rate of carbon overflow and growth yield, and generates testable predictions about cellular energetic strategies and protein expression levels. CAFBA therefore provides a simple, quantitative approach to balancing the trade-off between growth and its associated biosynthetic costs at genome-scale, without the burden of tuning many inaccessible parameters.

Mapping the landscape of metabolic goals of a cell

Genome-scale flux balance models of metabolism provide testable predictions of all metabolic rates in an organism, by assuming that the cell is optimizing a metabolic goal known as the objective function. We introduce an efficient inverse flux balance analysis (invFBA) approach, based on linear programming duality, to characterize the space of possible objective functions compatible with measured fluxes. After testing our algorithm on simulated E. coli data and time-dependent S. oneidensis fluxes inferred from gene expression, we apply our inverse approach to flux measurements in long-term evolved E. coli strains, revealing objective functions that provide insight into metabolic adaptation trajectories.

Flux Balance Analysis of Plant Metabolism: The Effect of Biomass Composition and Model Structure on Model Predictions

The biomass composition represented in constraint-based metabolic models is a key component for predicting cellular metabolism using flux balance analysis (FBA). Despite major advances in analytical technologies, it is often challenging to obtain a detailed composition of all major biomass components experimentally. Studies examining the influence of the biomass composition on the predictions of metabolic models have so far mostly been done on models of microorganisms. Little is known about the impact of varying biomass composition on flux prediction in FBA models of plants, whose metabolism is very versatile and complex because of the presence of multiple subcellular compartments. Also, the published metabolic models of plants differ in size and complexity. In this study, we examined the sensitivity of the predicted fluxes of plant metabolic models to biomass composition and model structure. These questions were addressed by evaluating the sensitivity of predictions of growth rates and central carbon metabolic fluxes to varying biomass compositions in three different genome-/large-scale metabolic models of Arabidopsis thaliana. Our results showed that fluxes through the central carbon metabolism were robust to changes in biomass composition. Nevertheless, comparisons between the predictions from three models using identical modeling constraints and objective function showed that model predictions were sensitive to the structure of the models, highlighting large discrepancies between the published models.

Proteome reallocation in Escherichia coli with increasing specific growth rate

Cells usually respond to changing growth conditions with a change in the specific growth rate (μ) and adjustment of their proteome to adapt and maintain metabolic efficiency. Description of the principles behind proteome resource allocation is important for understanding metabolic regulation in response to changing μ. Thus, we analysed the proteome resource allocation dynamics of Escherichia coli into different metabolic processes in response to changing μ. E. coli was grown on minimal and defined rich media in steady state continuous cultures at different μ and characterised combining two LC-MS/MS-based proteomics methods: stable isotope labelling by amino acids in cell culture (SILAC) and intensity based label-free absolute quantification. We detected slowly growing cells investing more proteome resources in energy generation and carbohydrate transport and metabolism whereas for achieving faster growth cells needed to devote most resources to translation and processes closely related to the protein synthesis pipeline. Furthermore, down-regulation of energy generation and carbohydrate metabolism proteins with faster growth displayed very similar expression dynamics with the global transcriptional regulator CRP (cyclic AMP receptor protein), pointing to a dominant protein resource allocating role of this protein. Our data also suggest that acetate overflow may be the result of global proteome resource optimisation as cells saved proteome resources by switching from fully respiratory to respiro-fermentative growth. The presented results give a quantitative overview of how E. coli adjusts its proteome to achieve faster growth and in future could contribute to the design of more efficient cell factories through proteome optimisation.

A Method to Constrain Genome-Scale Models with 13C Labeling Data

Current limitations in quantitatively predicting biological behavior hinder our efforts to engineer biological systems to produce biofuels and other desired chemicals. Here, we present a new method for calculating metabolic fluxes, key targets in metabolic engineering, that incorporates data from ¹³C labeling experiments and genome-scale models. The data from ¹³C labeling experiments provide strong flux constraints that eliminate the need to assume an evolutionary optimization principle such as the growth rate optimization assumption used in Flux Balance Analysis (FBA). This effective constraining is achieved by making the simple but biologically relevant assumption that flux flows from core to peripheral metabolism and does not flow back. The new method is significantly more robust than FBA with respect to errors in genome-scale model reconstruction. Furthermore, it can provide a comprehensive picture of metabolite balancing and predictions for unmeasured extracellular fluxes as constrained by ¹³C labeling data. A comparison shows that the results of this new method are similar to those found through ¹³C Metabolic Flux Analysis (¹³C MFA) for central carbon metabolism but, additionally, it provides flux estimates for peripheral metabolism. The extra validation gained by matching 48 relative labeling measurements is used to identify where and why several existing COnstraint Based Reconstruction and Analysis (COBRA) flux prediction algorithms fail. We demonstrate how to use this knowledge to refine these methods and improve their predictive capabilities. This method provides a reliable base upon which to improve the design of biological systems.

While metabolic fluxes constitute the most direct window into a cell’s metabolism, their accurate measurement is non trivial. The gold standard for flux measurement involves providing a labeled feed where some of the carbon atoms have been substituted by isotopes with higher atomic mass (¹³C instead of ¹²C). The ensuing labeling found in intracellular metabolites is then used to computationally infer the metabolic fluxes that produced the observed pattern. However, this procedure is typically performed with small metabolic models encompassing only central carbon metabolism. The genomic revolution has afforded us easily available genomes and, with them, comprehensive genome-scale models of cellular metabolism. It would be desirable to use the ¹³C labeling experimental data to constrain genome-scale models: these data constrain fluxes very effectively and provide in the labeling data fit an obvious proof that the underlying model correctly explains measured quantities. Here, we introduce a rigorous, self-consistent method that uses the full amount of information contained in ¹³C labeling data to constrain fluxes for a genome-scale model where underlying assumptions are explicitly stated.

Quantifying biomass composition by gas chromatography/mass spectrometry

We developed a set of methods for the quantification of four major components of microbial biomass using gas chromatography/mass spectrometry (GC/MS). Specifically, methods are described to quantify amino acids, RNA, fatty acids, and glycogen, which comprise an estimated 88% of the dry weight of Escherichia coli. Quantification is performed by isotope ratio analysis with fully (13)C-labeled biomass as internal standard, which is generated by growing E. coli on [U-(13)C]glucose. This convenient, reliable, and accurate single-platform (GC/MS) workflow for measuring biomass composition offers significant advantages over existing methods. We demonstrate the consistency, accuracy, precision, and utility of this procedure by applying it to three metabolically unique E. coli strains. The presented methods will have widespread applicability in systems microbiology and bioengineering.

Toward quantitative understanding on microbial community structure and functioning: a modeling-centered approach using degradation of marine oil spills as example

Molecular ecology approaches are rapidly advancing our insights into the microorganisms involved in the degradation of marine oil spills and their metabolic potentials. Yet, many questions remain open: how do oil-degrading microbial communities assemble in terms of functional diversity, species abundances and organization and what are the drivers? How do the functional properties of microorganisms scale to processes at the ecosystem level? How does mass flow among species, and which factors and species control and regulate fluxes, stability and other ecosystem functions? Can generic rules on oil-degradation be derived, and what drivers underlie these rules? How can we engineer oil-degrading microbial communities such that toxic polycyclic aromatic hydrocarbons are degraded faster? These types of questions apply to the field of microbial ecology in general. We outline how recent advances in single-species systems biology might be extended to help answer these questions. We argue that bottom-up mechanistic modeling allows deciphering the respective roles and interactions among microorganisms. In particular constraint-based, metagenome-derived community-scale flux balance analysis appears suited for this goal as it allows calculating degradation-related fluxes based on physiological constraints and growth strategies, without needing detailed kinetic information. We subsequently discuss what is required to make these approaches successful, and identify a need to better understand microbial physiology in order to advance microbial ecology. We advocate the development of databases containing microbial physiological data. Answering the posed questions is far from trivial. Oil-degrading communities are, however, an attractive setting to start testing systems biology-derived models and hypotheses as they are relatively simple in diversity and key activities, with several key players being isolated and a high availability of experimental data and approaches.

Constraint-based models predict metabolic and associated cellular functions

The prediction of cellular function from a genotype is a fundamental goal in biology. For metabolism, constraint-based modelling methods systematize biochemical, genetic and genomic knowledge into a mathematical framework that enables a mechanistic description of metabolic physiology. The use of constraint-based approaches has evolved over ~30 years, and an increasing number of studies have recently combined models with high-throughput data sets for prospective experimentation. These studies have led to validation of increasingly important and relevant biological predictions. As reviewed here, these recent successes have tangible implications in the fields of microbial evolution, interaction networks, genetic engineering and drug discovery.

Reconciling in vivo and in silico key biological parameters of Pseudomonas putida KT2440 during growth on glucose under carbon-limited condition

Background

Genome scale metabolic reconstructions are developed to efficiently engineer biocatalysts and bioprocesses based on a rational approach. However, in most reconstructions, due to the lack of appropriate measurements, experimentally determined growth parameters are simply taken from literature including other organisms, which reduces the usefulness and suitability of these models. Pseudomonas putida KT2440 is an outstanding biocatalyst given its versatile metabolism, its ability to generate sufficient energy and turnover of NADH and NAD. To apply this strain optimally in industrial production, a previously developed genome-scale metabolic model (iJP815) was experimentally assessed and streamlined to enable accurate predictions of the outcome of metabolic engineering approaches.

Results

To substantially improve the accuracy of the genome scale model (iJP815), continuous bioreactor cultures on a mineral medium with glucose as a sole carbon source were carried out at different dilution rates, which covered pulling analysis of the macromolecular composition of the biomass. Besides, the maximum biomass yield (on substrate) of 0.397 g_DCW · g_glc^-1, the maintenance coefficient of 0.037 g_glc · g_DCW^-1 · h^-1 and the maximum specific growth rate of 0.59 h^-1 were determined. Only the DNA fraction increased with the specific growth rate. This resulted in reliable estimation for the Growth-Associated Maintenance (GAM) of 85 mmol_ATP · g_DCW^-1 and the Non Growth-Associated Maintenance (NGAM) of 3.96 mmol_ATP · g_DCW^-1 · h^-1. Both values were found significantly different from previous assignment as a consequence of a lower yield and higher maintenance coefficient than originally assumed. Contrasting already published ¹³C flux measurements and the improved model allowed for constraining the solution space, by eliminating futile cycles. Furthermore, the model predictions were compared with transcriptomic data at overall good consistency, which helped to identify missing links.

Conclusions

By careful interpretation of growth stoichiometry and kinetics when grown in the presence of glucose, this work reports on an accurate genome scale metabolic model of Pseudomonas putida, providing a solid basis for its use in designing superior strains for biocatalysis. By consideration of substrate specific variation in stoichiometry and kinetics, it can be extended to other substrates and new mutants.

Relative amino acid composition signatures of organisms and environments

BACKGROUND

Identifying organism-environment interactions at the molecular level is crucial to understanding how organisms adapt to and change the chemical and molecular landscape of their habitats. In this work we investigated whether relative amino acid compositions could be used as a molecular signature of an environment and whether such a signature could also be observed at the level of the cellular amino acid composition of the microorganisms that inhabit that environment.

METHODOLOGIES/PRINCIPAL FINDINGS

To address these questions we collected and analyzed environmental amino acid determinations from the literature, and estimated from complete genomic sequences the global relative amino acid abundances of organisms that are cognate to the different types of environment. Environmental relative amino acid abundances clustered into broad groups (ocean waters, host-associated environments, grass land environments, sandy soils and sediments, and forest soils), indicating the presence of amino acid signatures specific for each environment. These signatures correlate to those found in organisms. Nevertheless, relative amino acid abundance of organisms was more influenced by GC content than habitat or phylogeny.

CONCLUSIONS

Our results suggest that relative amino acid composition can be used as a signature of an environment. In addition, we observed that the relative amino acid composition of organisms is not highly determined by environment, reinforcing previous studies that find GC content to be the major factor correlating to amino acid composition in living organisms.

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

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

Predictive modeling of biomass component tradeoffs in Brassica napus developing oilseeds based on in silico manipulation of storage metabolism

Seed oil content is a key agronomical trait, while the control of carbon allocation into different seed storage compounds is still poorly understood and hard to manipulate. Using bna572, a large-scale model of cellular metabolism in developing embryos of rapeseed (Brassica napus) oilseeds, we present an in silico approach for the analysis of carbon allocation into seed storage products. Optimal metabolic flux states were obtained by flux variability analysis based on minimization of the uptakes of substrates in the natural environment of the embryo. For a typical embryo biomass composition, flux sensitivities to changes in different storage components were derived. Upper and lower flux bounds of each reaction were categorized as oil or protein responsive. Among the most oil-responsive reactions were glycolytic reactions, while reactions related to mitochondrial ATP production were most protein responsive. To assess different biomass compositions, a tradeoff between the fractions of oil and protein was simulated. Based on flux-bound discontinuities and shadow prices along the tradeoff, three main metabolic phases with distinct pathway usage were identified. Transitions between the phases can be related to changing modes of the tricarboxylic acid cycle, reorganizing the usage of organic carbon and nitrogen sources for protein synthesis and acetyl-coenzyme A for cytosol-localized fatty acid elongation. The phase close to equal oil and protein fractions included an unexpected pathway bypassing α-ketoglutarate-oxidizing steps in the tricarboxylic acid cycle. The in vivo relevance of the findings is discussed based on literature on seed storage metabolism.

Extensive exometabolome analysis reveals extended overflow metabolism in various microorganisms

Overflow metabolism is well known for yeast, bacteria and mammalian cells. It typically occurs under glucose excess conditions and is characterized by excretions of by-products such as ethanol, acetate or lactate. This phenomenon, also denoted the short-term Crabtree effect, has been extensively studied over the past few decades, however, its basic regulatory mechanism and functional role in metabolism is still unknown. Here we present a comprehensive quantitative and time-dependent analysis of the exometabolome of Escherichia coli, Corynebacterium glutamicum, Bacillus licheniformis, and Saccharomyces cerevisiae during well-controlled bioreactor cultivations. Most surprisingly, in all cases a great diversity of central metabolic intermediates and amino acids is found in the culture medium with extracellular concentrations varying in the micromolar range. Different hypotheses for these observations are formulated and experimentally tested. As a result, the intermediates in the culture medium during batch growth must originate from passive or active transportation due to a new phenomenon termed “extended” overflow metabolism. Moreover, we provide broad evidence that this could be a common feature of all microorganism species when cultivated under conditions of carbon excess and non-inhibited carbon uptake. In turn, this finding has consequences for metabolite balancing and, particularly, for intracellular metabolite quantification and ¹³C-metabolic flux analysis.

Applications of system-level models of metabolism for analysis of bacterial physiology and identification of new drug targets

For nearly all of the 20th century, biologists gained considerable insights into the fundamental principles of cellular dynamics by examining select modules of biochemical processes. This form of analysis provides detailed information about the workings of the examined pathways. However, any attempt to alter the normal function of bacteria (perhaps for industrial or medicinal goals) requires a detailed global understanding of cellular mechanisms. The reductionist mode of analysis cannot provide the required information for developing the needed perspective on the complex interactions of biochemical pathways. Thankfully, the increasing availability of microbial genomic, transcriptomic, proteomic and other high-throughput data permits system-level analyses of microbiology. During the past two decades, systems biologists have developed constraint-based genome-scale models (GSM) of metabolism for a variety of pathogens. These models are important tools for assessing the metabolic capabilities of various genotypes. Simultaneously, new computational methods have been developed that use these network reconstructions to answer an array of important immunological questions. The objective of this article is to briefly review some of the uses of GSMs for studying bacterial metabolism under different conditions and to discuss how the calculated solutions can be used for rational design of drugs.

Metabolic modeling of Chlamydomonas reinhardtii: energy requirements for photoautotrophic growth and maintenance

In this study, a metabolic network describing the primary metabolism of Chlamydomonas reinhardtii was constructed. By performing chemostat experiments at different growth rates, energy parameters for maintenance and biomass formation were determined. The chemostats were run at low irradiances resulting in a high biomass yield on light of 1.25 g  mol(-1). The ATP requirement for biomass formation from biopolymers (K(x)) was determined to be 109 mmol g(-1) (18.9 mol mol(-1)) and the maintenance requirement (m(ATP)) was determined to be 2.85 mmol g(-1) h(-1). With these energy requirements included in the metabolic network, the network accurately describes the primary metabolism of C. reinhardtii and can be used for modeling of C. reinhardtii growth and metabolism. Simulations confirmed that cultivating microalgae at low growth rates is unfavorable because of the high maintenance requirements which result in low biomass yields. At high light supply rates, biomass yields will decrease due to light saturation effects. Thus, to optimize biomass yield on light energy in photobioreactors, an optimum between low and high light supply rates should be found. These simulations show that metabolic flux analysis can be used as a tool to gain insight into the metabolism of algae and ultimately can be used for the maximization of algal biomass and product yield. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1007/s10811-011-9674-3) contains supplementary material, which is available to authorized users.

Selective utilization of exogenous amino acids by Dehalococcoides ethenogenes strain 195 and its effects on growth and dechlorination activity

Bacteria of the genus Dehalococcoides are important members of bioremediation communities because of their ability to detoxify chloroethenes to the benign end product ethene. Genome-enabled studies conducted with Dehalococcoides ethenogenes 195 have revealed that two ATP-binding cassette (ABC)-type amino acid transporters are expressed during its exponential growth stages. In light of previous findings that Casamino Acids enhanced its dechlorination activity, we hypothesized that strain 195 is capable of importing amino acids from its environment to facilitate dechlorination and growth. To test this hypothesis, we applied isotopomer-based dilution analysis with (13)C-labeled acetate to differentiate the amino acids that were taken up by strain 195 from those synthesized de novo and to determine the physiological changes caused by the significantly incorporated amino acids. Our results showed that glutamate/glutamine and aspartate/asparagine were almost exclusively synthesized by strain 195, even when provided in excess in the medium. In contrast, phenylalanine, isoleucine, leucine, and methionine were identified as the four most highly incorporated amino acids, at levels >30% of respective proteinogenic amino acids. When either phenylalanine or all four highly incorporated amino acids were added to the defined mineral medium, the growth rates, dechlorination activities, and yields of strain 195 were enhanced to levels similar to those observed with supplementation with 20 amino acids. However, genes for the putative ABC-type amino acids transporters and phenylalanine biosynthesis exhibited insignificant regulation in response to the imported amino acids. This study also demonstrates that using isotopomer-based metabolite analysis can be an efficient strategy for optimizing nutritional conditions for slow-growing microorganisms.

Decrease of energy spilling in Escherichia coli continuous cultures with rising specific growth rate and carbon wasting

Background

Growth substrates, aerobic/anaerobic conditions, specific growth rate (μ) etc. strongly influence Escherichia coli cell physiology in terms of cell size, biomass composition, gene and protein expression. To understand the regulation behind these different phenotype properties, it is useful to know carbon flux patterns in the metabolic network which are generally calculated by metabolic flux analysis (MFA). However, rarely is biomass composition determined and carbon balance carefully measured in the same experiments which could possibly lead to distorted MFA results and questionable conclusions. Therefore, we carried out both detailed carbon balance and biomass composition analysis in the same experiments for more accurate quantitative analysis of metabolism and MFA.

Results

We applied advanced continuous cultivation methods (A-stat and D-stat) to continuously monitor E. coli K-12 MG1655 flux and energy metabolism dynamic responses to change of μ and glucose-acetate co-utilisation. Surprisingly, a 36% reduction of ATP spilling was detected with increasing μ and carbon wasting to non-CO₂ by-products under constant biomass yield. The apparent discrepancy between constant biomass yield and decline of ATP spilling could be explained by the rise of carbon wasting from 3 to 11% in the carbon balance which was revealed by the discovered novel excretion profile of E. coli pyrimidine pathway intermediates carbamoyl-phosphate, dihydroorotate and orotate. We found that carbon wasting patterns are dependent not only on μ, but also on glucose-acetate co-utilisation capability. Accumulation of these compounds was coupled to the two-phase acetate accumulation profile. Acetate overflow was observed in parallel with the reduction of TCA cycle and glycolysis fluxes, and induction of pentose phosphate pathway.

Conclusions

It can be concluded that acetate metabolism is one of the major regulating factors of central carbon metabolism. More importantly, our model calculations with actual biomass composition and detailed carbon balance analysis in steady state conditions with -omics data comparison demonstrate the importance of a comprehensive systems biology approach for more advanced understanding of metabolism and carbon re-routing mechanisms potentially leading to more successful metabolic engineering.

Expanding a dynamic flux balance model of yeast fermentation to genome-scale

Background

Yeast is considered to be a workhorse of the biotechnology industry for the production of many value-added chemicals, alcoholic beverages and biofuels. Optimization of the fermentation is a challenging task that greatly benefits from dynamic models able to accurately describe and predict the fermentation profile and resulting products under different genetic and environmental conditions. In this article, we developed and validated a genome-scale dynamic flux balance model, using experimentally determined kinetic constraints.

Results

Appropriate equations for maintenance, biomass composition, anaerobic metabolism and nutrient uptake are key to improve model performance, especially for predicting glycerol and ethanol synthesis. Prediction profiles of synthesis and consumption of the main metabolites involved in alcoholic fermentation closely agreed with experimental data obtained from numerous lab and industrial fermentations under different environmental conditions. Finally, fermentation simulations of genetically engineered yeasts closely reproduced previously reported experimental results regarding final concentrations of the main fermentation products such as ethanol and glycerol.

Conclusion

A useful tool to describe, understand and predict metabolite production in batch yeast cultures was developed. The resulting model, if used wisely, could help to search for new metabolic engineering strategies to manage ethanol content in batch fermentations.

Effect of iclR and arcA knockouts on biomass formation and metabolic fluxes in Escherichia coli K12 and its implications on understanding the metabolism of Escherichia coli BL21 (DE3)

Background

Gene expression is regulated through a complex interplay of different transcription factors (TFs) which can enhance or inhibit gene transcription. ArcA is a global regulator that regulates genes involved in different metabolic pathways, while IclR as a local regulator, controls the transcription of the glyoxylate pathway genes of the aceBAK operon. This study investigates the physiological and metabolic consequences of arcA and iclR deletions on E. coli K12 MG1655 under glucose abundant and limiting conditions and compares the results with the metabolic characteristics of E. coli BL21 (DE3).

Results

The deletion of arcA and iclR results in an increase in the biomass yield both under glucose abundant and limiting conditions, approaching the maximum theoretical yield of 0.65 c-mole/c-mole glucose under glucose abundant conditions. This can be explained by the lower flux through several CO₂ producing pathways in the E. coli K12 ΔarcAΔiclR double knockout strain. Due to iclR gene deletion, the glyoxylate pathway is activated resulting in a redirection of 30% of the isocitrate molecules directly to succinate and malate without CO₂ production. Furthermore, a higher flux at the entrance of the TCA was noticed due to arcA gene deletion, resulting in a reduced production of acetate and less carbon loss. Under glucose limiting conditions the flux through the glyoxylate pathway is further increased in the ΔiclR knockout strain, but this effect was not observed in the double knockout strain. Also a striking correlation between the glyoxylate flux data and the isocitrate lyase activity was observed for almost all strains and under both growth conditions, illustrating the transcriptional control of this pathway. Finally, similar central metabolic fluxes were observed in E. coli K12 ΔarcA ΔiclR compared to the industrially relevant E. coli BL21 (DE3), especially with respect to the pentose pathway, the glyoxylate pathway, and the TCA fluxes. In addition, a comparison of the genome sequences of the two strains showed that BL21 possesses two mutations in the promoter region of iclR and rare codons are present in arcA implying a lower tRNA acceptance. Both phenomena presumably result in a reduced ArcA and IclR synthesis in BL21, which contributes to the similar physiology as observed in E. coli K12 ΔarcAΔiclR.

Conclusions

The deletion of arcA results in a decrease of repression on transcription of TCA cycle genes under glucose abundant conditions, without significantly affecting the glyoxylate pathway activity. IclR clearly represses transcription of glyoxylate pathway genes under glucose abundance, a condition in which Crp activation is absent. Under glucose limitation, Crp is responsible for the high glyoxylate flux, but IclR still represses transcription. Finally, in E. coli BL21 (DE3), ArcA and IclR are poorly expressed, explaining the similar fluxes observed compared to the ΔarcAΔiclR strain.

Flux balance analysis of ammonia assimilation network in E. coli predicts preferred regulation point

Nitrogen assimilation is a critical biological process for the synthesis of biomolecules in Escherichia coli. The central ammonium assimilation network in E. coli converts carbon skeleton α-ketoglutarate and ammonium into glutamate and glutamine, which further serve as nitrogen donors for nitrogen metabolism in the cell. This reaction network involves three enzymes: glutamate dehydrogenase (GDH), glutamine synthetase (GS) and glutamate synthase (GOGAT). In minimal media, E. coli tries to maintain an optimal growth rate by regulating the activity of the enzymes to match the availability of the external ammonia. The molecular mechanism and the strategy of the regulation in this network have been the research topics for many investigators. In this paper, we develop a flux balance model for the nitrogen metabolism, taking into account of the cellular composition and biosynthetic requirements for nitrogen. The model agrees well with known experimental results. Specifically, it reproduces all the (15)N isotope labeling experiments in the wild type and the two mutant (ΔGDH and ΔGOGAT) strains of E. coli. Furthermore, the predicted catalytic activities of GDH, GS and GOGAT in different ammonium concentrations and growth rates for the wild type, ΔGDH and ΔGOGAT strains agree well with the enzyme concentrations obtained from western blots. Based on this flux balance model, we show that GS is the preferred regulation point among the three enzymes in the nitrogen assimilation network. Our analysis reveals the pattern of regulation in this central and highly regulated network, thus providing insights into the regulation strategy adopted by the bacteria. Our model and methods may also be useful in future investigations in this and other networks.

Genome-scale metabolic networks

During the last decade, models have been developed to characterize cellular metabolism at the level of an entire metabolic network. The main concept that underlies whole-network metabolic modeling is the identification and mathematical definition of constraints. Here, we review large-scale metabolic network modeling, in particular, stoichiometric- and constraint-based approaches. Although many such models have been reconstructed, few networks have been extensively validated and tested experimentally, and we focus on these. We describe how metabolic networks can be represented using stoichiometric matrices and well-defined constraints on metabolic fluxes. We then discuss relatively successful approaches, including flux balance analysis (FBA), pathway analysis, and common extensions or modifications to these approaches. Finally, we describe techniques for integrating these approaches with models of other biological processes.

Synergy between (13)C-metabolic flux analysis and flux balance analysis for understanding metabolic adaptation to anaerobiosis in E. coli

Genome-based Flux Balance Analysis (FBA) and steady-state isotopic-labeling-based Metabolic Flux Analysis (MFA) are complimentary approaches to predicting and measuring the operation and regulation of metabolic networks. Here, genome-derived models of Escherichia coli (E. coli) metabolism were used for FBA and ¹³C-MFA analyses of aerobic and anaerobic growths of wild-type E. coli (K-12 MG1655) cells. Validated MFA flux maps reveal that the fraction of maintenance ATP consumption in total ATP production is about 14% higher under anaerobic (51.1%) than aerobic conditions (37.2%). FBA revealed that an increased ATP utilization is consumed by ATP synthase to secrete protons from fermentation. The TCA cycle is shown to be incomplete in aerobically growing cells and submaximal growth is due to limited oxidative phosphorylation. An FBA was successful in predicting product secretion rates in aerobic culture if both glucose and oxygen uptake measurement were constrained, but the most-frequently predicted values of internal fluxes yielded from sampling the feasible space differ substantially from MFA-derived fluxes.

Genome-derived minimal metabolic models for Escherichia coli MG1655 with estimated in vivo respiratory ATP stoichiometry

Metabolic network models describing growth of Escherichia coli on glucose, glycerol and acetate were derived from a genome scale model of E. coli. One of the uncertainties in the metabolic networks is the exact stoichiometry of energy generating and consuming processes. Accurate estimation of biomass and product yields requires correct information on the ATP stoichiometry. The unknown ATP stoichiometry parameters of the constructed E. coli network were estimated from experimental data of eight different aerobic chemostat experiments carried out with E. coli MG1655, grown at different dilution rates (0.025, 0.05, 0.1, and 0.3 h(-1)) and on different carbon substrates (glucose, glycerol, and acetate). Proper estimation of the ATP stoichiometry requires proper information on the biomass composition of the organism as well as accurate assessment of net conversion rates under well-defined conditions. For this purpose a growth rate dependent biomass composition was derived, based on measurements and literature data. After incorporation of the growth rate dependent biomass composition in a metabolic network model, an effective P/O ratio of 1.49 +/- 0.26 mol of ATP/mol of O, K(X) (growth dependent maintenance) of 0.46 +/- 0.27 mol of ATP/C-mol of biomass and m(ATP) (growth independent maintenance) of 0.075 +/- 0.015 mol of ATP/C-mol of biomass/h were estimated using a newly developed Comprehensive Data Reconciliation (CDR) method, assuming that the three energetic parameters were independent of the growth rate and the used substrate. The resulting metabolic network model only requires the specific rate of growth, micro, as an input in order to accurately predict all other fluxes and yields.

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.

UPLC/MS based method for quantitative determination of fatty acid composition in Gram-negative and Gram-positive bacteria

Quantitative fatty acid composition of microorganisms at various growth space points is required for understanding membrane associated processes of cells, but the majority of the relevant publications still restrict to the relative compositions. In the current study, a simple and reliable method for quantitative measurement of fatty acid content in bacterial biomass without prior derivatization using ultra performance liquid chromatography-electrospray ionization mass spectrometry was developed. The method was applied for investigating the influence of specific growth rate and pH on the fatty acid profiles of two biotechnologically important microorganisms - Gram-negative bacteria Escherichia coli and Gram-positive bacteria Lactococcus lactis grown in controlled physiological states. It was found that the membranes of slowly growing cells are more rigid and that the fatty acid fraction of the cells of L. lactis diminishes considerably with increasing growth rate.

Constraint-based model of Shewanella oneidensis MR-1 metabolism: a tool for data analysis and hypothesis generation

Shewanellae are gram-negative facultatively anaerobic metal-reducing bacteria commonly found in chemically (i.e., redox) stratified environments. Occupying such niches requires the ability to rapidly acclimate to changes in electron donor/acceptor type and availability; hence, the ability to compete and thrive in such environments must ultimately be reflected in the organization and utilization of electron transfer networks, as well as central and peripheral carbon metabolism. To understand how Shewanella oneidensis MR-1 utilizes its resources, the metabolic network was reconstructed. The resulting network consists of 774 reactions, 783 genes, and 634 unique metabolites and contains biosynthesis pathways for all cell constituents. Using constraint-based modeling, we investigated aerobic growth of S. oneidensis MR-1 on numerous carbon sources. To achieve this, we (i) used experimental data to formulate a biomass equation and estimate cellular ATP requirements, (ii) developed an approach to identify cycles (such as futile cycles and circulations), (iii) classified how reaction usage affects cellular growth, (iv) predicted cellular biomass yields on different carbon sources and compared model predictions to experimental measurements, and (v) used experimental results to refine metabolic fluxes for growth on lactate. The results revealed that aerobic lactate-grown cells of S. oneidensis MR-1 used less efficient enzymes to couple electron transport to proton motive force generation, and possibly operated at least one futile cycle involving malic enzymes. Several examples are provided whereby model predictions were validated by experimental data, in particular the role of serine hydroxymethyltransferase and glycine cleavage system in the metabolism of one-carbon units, and growth on different sources of carbon and energy. This work illustrates how integration of computational and experimental efforts facilitates the understanding of microbial metabolism at a systems level.

Flux balance analysis accounting for metabolite dilution

Flux balance analysis is a common method for predicting steady-state flux distributions within metabolic networks, accounting for the growth demand for the synthesis of a predefined set of essential biomass precursors. Ignoring the growth demand for the synthesis of intermediate metabolites required for balancing their dilution leads flux balance analysis to false predictions in some cases. Here, we present metabolite dilution flux balance analysis, which addresses this problem, resulting in improved metabolic phenotype predictions.

(13)C-based metabolic flux analysis

Stable isotope, and in particular (13)C-based flux analysis, is the exclusive approach to experimentally quantify the integrated responses of metabolic networks. Here we describe a protocol that is based on growing microbes on (13)C-labeled glucose and subsequent gas chromatography mass spectrometric detection of (13)C-patterns in protein-bound amino acids. Relying on publicly available software packages, we then describe two complementary mathematical approaches to estimate either local ratios of converging fluxes or absolute fluxes through different pathways. As amino acids in cell protein are abundant and stable, this protocol requires a minimum of equipment and analytical expertise. Most other flux methods are variants of the principles presented here. A true alternative is the analytically more demanding dynamic flux analysis that relies on (13)C-pattern in free intracellular metabolites. The presented protocols take 5-10 d, have been used extensively in the past decade and are exemplified here for the central metabolism of Escherichia coli.

Reconstruction of biochemical networks in microorganisms

Systems analysis of metabolic and growth functions in microbial organisms is rapidly developing and maturing. Such studies are enabled by reconstruction, at the genomic scale, of the biochemical reaction networks that underlie cellular processes. The network reconstruction process is organism specific and is based on an annotated genome sequence, high-throughput network-wide data sets and bibliomic data on the detailed properties of individual network components. Here we describe the process that is currently used to achieve comprehensive network reconstructions and discuss how these reconstructions are curated and validated. This review should aid the growing number of researchers who are carrying out reconstructions for particular target organisms.

Modeling Neisseria meningitidis B metabolism at different specific growth rates

Neisseria meningitidis is a human pathogen that can infect diverse sites within the human host. The major diseases caused by N. meningitidis are responsible for death and disability, especially in young infants. At the Netherlands Vaccine Institute (NVI) a vaccine against serogroup B organisms is currently being developed. This study describes the influence of the growth rate of N. meningitidis on its macro-molecular composition and its metabolic activity and was determined in chemostat cultures. In the applied range of growth rates, no significant changes in RNA content and protein content with growth rate were observed in N. meningitidis. The DNA content in N. meningitidis was somewhat higher at the highest applied growth rate. The phospholipid and lipopolysaccharide content in N. meningitidis changed with growth rate but no specific trends were observed. The cellular fatty acid composition and the amino acid composition did not change significantly with growth rate. Additionally, it was found that the PorA content in outer membrane vesicles was significantly lower at the highest growth rate. The metabolic fluxes at various growth rates were calculated using flux balance analysis. Errors in fluxes were calculated using Monte Carlo Simulation and the reliability of the calculated flux distribution could be indicated, which has not been reported for this type of analysis. The yield of biomass on substrate (Y(x/s)) and the maintenance coefficient (m(s)) were determined as 0.44 (+/-0.04) g g(-1) and 0.04 (+/-0.02) g g(-1) h(-1), respectively. The growth associated energy requirement (Y(x/ATP)) and the non-growth associated ATP requirement for maintenance (m(ATP)) were estimated as 0.13 (+/-0.04) mol mol(-1) and 0.43 (+/-0.14) mol mol(-1) h(-1), respectively. It was found that the split ratio between the Entner-Doudoroff and the pentose phosphate pathway, the sole glucose utilizing pathways in N. meningitidis, had a minor effect on ATP formation rate but a major effect on the fluxes going through for instance the citric-acid cycle. For this reason, we presented flux ranges for underdetermined parts of metabolic network rather than presenting single flux values, which is more commonly done in literature.