Genome-scale reconstruction and analysis of the metabolic network in the hyperthermophilic archaeon sulfolobus solfataricus

We describe the reconstruction of a genome-scale metabolic model of the crenarchaeon Sulfolobus solfataricus, a hyperthermoacidophilic microorganism. It grows in terrestrial volcanic hot springs with growth occurring at pH 2-4 (optimum 3.5) and a temperature of 75-80°C (optimum 80°C). The genome of Sulfolobus solfataricus P2 contains 2,992,245 bp on a single circular chromosome and encodes 2,977 proteins and a number of RNAs. The network comprises 718 metabolic and 58 transport/exchange reactions and 705 unique metabolites, based on the annotated genome and available biochemical data. Using the model in conjunction with constraint-based methods, we simulated the metabolic fluxes induced by different environmental and genetic conditions. The predictions were compared to experimental measurements and phenotypes of S. solfataricus. Furthermore, the performance of the network for 35 different carbon sources known for S. solfataricus from the literature was simulated. Comparing the growth on different carbon sources revealed that glycerol is the carbon source with the highest biomass flux per imported carbon atom (75% higher than glucose). Experimental data was also used to fit the model to phenotypic observations. In addition to the commonly known heterotrophic growth of S. solfataricus, the crenarchaeon is also able to grow autotrophically using the hydroxypropionate-hydroxybutyrate cycle for bicarbonate fixation. We integrated this pathway into our model and compared bicarbonate fixation with growth on glucose as sole carbon source. Finally, we tested the robustness of the metabolism with respect to gene deletions using the method of Minimization of Metabolic Adjustment (MOMA), which predicted that 18% of all possible single gene deletions would be lethal for the organism.     

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.     

Dissecting the energy metabolism in Mycoplasma pneumoniae through genome‐scale metabolic modeling

Mycoplasma pneumoniae, a threatening pathogen with a minimal genome, is a model organism for bacterial systems biology for which substantial experimental information is available. With the goal of understanding the complex interactions underlying its metabolism, we analyzed and characterized the metabolic network of M. pneumoniae in great detail, integrating data from different omics analyses under a range of conditions into a constraint‐based model backbone. Iterating model predictions, hypothesis generation, experimental testing, and model refinement, we accurately curated the network and quantitatively explored the energy metabolism. In contrast to other bacteria, M. pneumoniae uses most of its energy for maintenance tasks instead of growth. We show that in highly linear networks the prediction of flux distributions for different growth times allows analysis of time‐dependent changes, albeit using a static model. By performing an in silico knock‐out study as well as analyzing flux distributions in single and double mutant phenotypes, we demonstrated that the model accurately represents the metabolism of M. pneumoniae. The experimentally validated model provides a solid basis for understanding its metabolic regulatory mechanisms.     

Predicting Growth Conditions from Internal Metabolic Fluxes in an In-Silico Model of E. coli

A widely studied problem in systems biology is to predict bacterial phenotype from growth conditions, using mechanistic models such as flux balance analysis (FBA). However, the inverse prediction of growth conditions from phenotype is rarely considered. Here we develop a computational framework to carry out this inverse prediction on a computational model of bacterial metabolism. We use FBA to calculate bacterial phenotypes from growth conditions in E. coli, and then we assess how accurately we can predict the original growth conditions from the phenotypes. Prediction is carried out via regularized multinomial regression. Our analysis provides several important physiological and statistical insights. First, we show that by analyzing metabolic end products we can consistently predict growth conditions. Second, prediction is reliable even in the presence of small amounts of impurities. Third, flux through a relatively small number of reactions per growth source (∼10) is sufficient for accurate prediction. Fourth, combining the predictions from two separate models, one trained only on carbon sources and one only on nitrogen sources, performs better than models trained to perform joint prediction. Finally, that separate predictions perform better than a more sophisticated joint prediction scheme suggests that carbon and nitrogen utilization pathways, despite jointly affecting cellular growth, may be fairly decoupled in terms of their dependence on specific assortments of molecular precursors.     

Microbial Community Metabolic Modeling: A Community Data‐Driven Network Reconstruction

Metabolic network modeling of microbial communities provides an in‐depth understanding of community‐wide metabolic and regulatory processes. Compared to single organism analyses, community metabolic network modeling is more complex because it needs to account for interspecies interactions. To date, most approaches focus on reconstruction of high‐quality individual networks so that, when combined, they can predict community behaviors as a result of interspecies interactions. However, this conventional method becomes ineffective for communities whose members are not well characterized and cannot be experimentally interrogated in isolation. Here, we tested a new approach that uses community‐level data as a critical input for the network reconstruction process. This method focuses on directly predicting interspecies metabolic interactions in a community, when axenic information is insufficient. We validated our method through the case study of a bacterial photoautotroph–heterotroph consortium that was used to provide data needed for a community‐level metabolic network reconstruction. Resulting simulations provided experimentally validated predictions of how a photoautotrophic cyanobacterium supports the growth of an obligate heterotrophic species by providing organic carbon and nitrogen sources. J. Cell. Physiol. 231: 2339–2345, 2016. © 2016 The Authors. Journal of Cellular Physiology Published by Wiley Periodicals, Inc.     

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.     

Flux balance analysis of the ammonia-oxidizing bacterium Nitrosomonas europaea ATCC19718 unravels specific metabolic activities while degrading toxic compounds

The ammonia-oxidizing bacterium Nitrosomonas europaea has been widely recognized as an important player in the nitrogen cycle as well as one of the most abundant members in microbial communities for the treatment of industrial or sewage wastewater. Its natural metabolic versatility and extraordinary ability to degrade environmental pollutants (e.g., aromatic hydrocarbons such as benzene and toluene) enable it to thrive under various harsh environmental conditions. Constraint-based metabolic models constructed from genome sequences enable quantitative insight into the central and specialized metabolism within a target organism. These genome-scale models have been utilized to understand, optimize, and design new strategies for improved bioprocesses. Reduced modeling approaches have been used to elucidate Nitrosomonas europaea metabolism at a pathway level. However, genome-scale knowledge about the simultaneous oxidation of ammonia and pollutant metabolism of N. europaea remains limited. Here, we describe the reconstruction, manual curation, and validation of the genome-scale metabolic model for N. europaea, iGC535. This reconstruction is the most accurate metabolic model for a nitrifying organism to date, reaching an average prediction accuracy of over 90% under several growth conditions. The manually curated model can predict phenotypes under chemolithotrophic and chemolithoorganotrophic conditions while oxidating methane and wastewater pollutants. Calculated flux distributions under different trophic conditions show that several key pathways are affected by the type of carbon source available, including central carbon metabolism and energy production.     

High-quality genome-scale metabolic modelling of Pseudomonas putida highlights its broad metabolic capabilities

Genome-scale reconstructions of metabolism are computational species-specific knowledge bases able to compute systemic metabolic properties. We present a comprehensive and validated reconstruction of the biotechnologically relevant bacterium Pseudomonas putida KT2440 that greatly expands computable predictions of its metabolic states. The reconstruction represents a significant reactome expansion over available reconstructed bacterial metabolic networks. Specifically, iJN1462 (i) incorporates several hundred additional genes and associated reactions resulting in new predictive capabilities, including new nutrients supporting growth; (ii) was validated by in vivo growth screens that included previously untested carbon (48) and nitrogen (41) sources; (iii) yielded gene essentiality predictions showing large accuracy when compared with a knock-out library and Bar-seq data; and (iv) allowed mapping of its network to 82 P. putida sequenced strains revealing functional core that reflect the large metabolic versatility of this species, including aromatic compounds derived from lignin. Thus, this study provides a thoroughly updated metabolic reconstruction and new computable phenotypes for P. putida, which can be leveraged as a first step toward understanding the pan metabolic capabilities of Pseudomonas.     

Metabolic flux analysis of Cyanothece sp. ATCC 51142 under mixotrophic conditions

Cyanobacteria are a group of photosynthetic prokaryotes capable of utilizing solar energy to fix atmospheric carbon dioxide to biomass. Despite several “proof of principle” studies, low product yield is an impediment in commercialization of cyanobacteria-derived biofuels. Estimation of intracellular reaction rates by ¹³C metabolic flux analysis (¹³C-MFA) would be a step toward enhancing biofuel yield via metabolic engineering. We report ¹³C-MFA for Cyanothece sp. ATCC 51142, a unicellular nitrogen-fixing cyanobacterium, known for enhanced hydrogen yield under mixotrophic conditions. Rates of reactions in the central carbon metabolism under nitrogen-fixing and -non-fixing conditions were estimated by monitoring the competitive incorporation of ¹²C and ¹³C from unlabeled CO₂ and uniformly labeled glycerol, respectively, into terminal metabolites such as amino acids. The observed labeling patterns suggest mixotrophic growth under both the conditions, with a larger fraction of unlabeled carbon in nitrate-sufficient cultures asserting a greater contribution of carbon fixation by photosynthesis and an anaplerotic pathway. Indeed, flux analysis complements the higher growth observed under nitrate-sufficient conditions. On the other hand, the flux through the oxidative pentose phosphate pathway and tricarboxylic acid cycle was greater in nitrate-deficient conditions, possibly to supply the precursors and reducing equivalents needed for nitrogen fixation. In addition, an enhanced flux through fructose-6-phosphate phosphoketolase possibly suggests the organism’s preferred mode under nitrogen-fixing conditions. The ¹³C-MFA results complement the reported predictions by flux balance analysis and provide quantitative insight into the organism’s distinct metabolic features under nitrogen-fixing and -non-fixing conditions.     

Glycolysis/gluconeogenesis specialization in microbes is driven by biochemical constraints of flux sensing

Central carbon metabolism is highly conserved across microbial species, but can catalyze very different pathways depending on the organism and their ecological niche. Here, we study the dynamic reorganization of central metabolism after switches between the two major opposing pathway configurations of central carbon metabolism, glycolysis, and gluconeogenesis in Escherichia coli, Pseudomonas aeruginosa, and Pseudomonas putida. We combined growth dynamics and dynamic changes in intracellular metabolite levels with a coarse‐grained model that integrates fluxes, regulation, protein synthesis, and growth and uncovered fundamental limitations of the regulatory network: After nutrient shifts, metabolite concentrations collapse to their equilibrium, rendering the cell unable to sense which direction the flux is supposed to flow through the metabolic network. The cell can partially alleviate this by picking a preferred direction of regulation at the expense of increasing lag times in the opposite direction. Moreover, decreasing both lag times simultaneously comes at the cost of reduced growth rate or higher futile cycling between metabolic enzymes. These three trade‐offs can explain why microorganisms specialize for either glycolytic or gluconeogenic substrates and can help elucidate the complex growth patterns exhibited by different microbial species.     

A Strategy to Calculate the Patterns of Nutrient Consumption by Microorganisms Applying a Two-Level Optimisation Principle to Reconstructed Metabolic Networks

Bacterial responses to environmental changes rely on a complex network of biochemical reactions. The properties of the metabolic network determining these responses can be divided into two groups: the stoichiometric properties, given by the stoichiometry matrix, and the kinetic/thermodynamic properties, given by the rate equations of the reaction steps. The stoichiometry matrix represents the maximal metabolic capabilities of the organism, and the regulatory mechanisms based on the rate laws could be considered as being responsible for the administration of these capabilities. Post-genomic reconstruction of metabolic networks provides us with the stoichiometry matrix of particular strains of microorganisms, but the kinetic aspects of in vivo rate laws are still largely unknown. Therefore, the validity of predictions of cellular responses requiring detailed knowledge of the rate equations is difficult to assert. In this paper, we show that by applying optimisation criteria to the core stoichiometric network of the metabolism of Escherichia coli, and including information about reversibility/irreversibility only of the reaction steps, it is possible to calculate bacterial responses to growth media with different amounts of glucose and galactose. The target was the minimisation of the number of active reactions (subject to attaining a growth rate higher than a lower limit) and subsequent maximisation of the growth rate (subject to the number of active reactions being equal to the minimum previously calculated). Using this two-level target, we were able to obtain by calculation four fundamental behaviours found experimentally: inhibition of respiration at high glucose concentrations in aerobic conditions, turning on of respiration when glucose decreases, induction of galactose utilisation when the system is depleted of glucose and simultaneous use of glucose and galactose as carbon sources when both sugars are present in low concentrations. Preliminary results of the coarse pattern of sugar utilisation were also obtained with a genome-scale E. coli reconstructed network, yielding similar qualitative results.     

Systems analysis of metabolism in the pathogenic trypanosomatid Leishmania major

Systems analyses have facilitated the characterization of metabolic networks of several organisms. We have reconstructed the metabolic network of Leishmania major, a poorly characterized organism that causes cutaneous leishmaniasis in mammalian hosts. This network reconstruction accounts for 560 genes, 1112 reactions, 1101 metabolites and 8 unique subcellular localizations. Using a systems‐based approach, we hypothesized a comprehensive set of lethal single and double gene deletions, some of which were validated using published data with approximately 70% accuracy. Additionally, we generated hypothetical annotations to dozens of previously uncharacterized genes in the L. major genome and proposed a minimal medium for growth. We further demonstrated the utility of a network reconstruction with two proof‐of‐concept examples that yielded insight into robustness of the network in the presence of enzymatic inhibitors and delineation of promastigote/amastigote stage‐specific metabolism. This reconstruction and the associated network analyses of L. major is the first of its kind for a protozoan. It can serve as a tool for clarifying discrepancies between data sources, generating hypotheses that can be experimentally validated and identifying ideal therapeutic targets.     

A scalable algorithm to explore the Gibbs energy landscape of genome-scale metabolic networks

The integration of various types of genomic data into predictive models of biological networks is one of the main challenges currently faced by computational biology. Constraint-based models in particular play a key role in the attempt to obtain a quantitative understanding of cellular metabolism at genome scale. In essence, their goal is to frame the metabolic capabilities of an organism based on minimal assumptions that describe the steady states of the underlying reaction network via suitable stoichiometric constraints, specifically mass balance and energy balance (i.e. thermodynamic feasibility). The implementation of these requirements to generate viable configurations of reaction fluxes and/or to test given flux profiles for thermodynamic feasibility can however prove to be computationally intensive. We propose here a fast and scalable stoichiometry-based method to explore the Gibbs energy landscape of a biochemical network at steady state. The method is applied to the problem of reconstructing the Gibbs energy landscape underlying metabolic activity in the human red blood cell, and to that of identifying and removing thermodynamically infeasible reaction cycles in the Escherichia coli metabolic network (iAF1260). In the former case, we produce consistent predictions for chemical potentials (or log-concentrations) of intracellular metabolites; in the latter, we identify a restricted set of loops (23 in total) in the periplasmic and cytoplasmic core as the origin of thermodynamic infeasibility in a large sample (10(6)) of flux configurations generated randomly and compatibly with the prior information available on reaction reversibility.     

Bioenergetic strategy of microalgae for the biodegradation of phenolic compounds: exogenously supplied energy and carbon sources adjust the level of biodegradation

The biodegradation of phenolic compounds by microalgae seems to be not a simple feature of a particular organism, but mostly a bioenergetic process depending on the growth conditions, especially on the exogenously supplied energy (carbon and light) sources. By using chlorophyll fluorescence induction measurements to estimate the molecular structure and function of the photosynthetic apparatus and therefore the tolerance/sensitivity of microalgae incubated with phenols, it can be assumed that, at least in low concentrations, phenol have no toxic effects on the cultures and can be used as alternative carbon source in them. Halophenols (chlorophenols, bromophenols and iodophenols) are quite toxic for the microalgal cultures. In halophenols the first step of the biodegradation is the split of the halogen substituent (dehalogenation). This is strongly determined by the bond dissociation energy of the corresponding substituent and therefore the energetic requirement for the biodegradation of halophenols increases following the sequence: iodophenol相似文献   

Carbon and Energy Sources for the Nitrifying Autotroph Nitrobacter.

Delwiche, C. C. (University of California, Davis), and M. S. Finstein. Carbon and energy sources for the nitrifying autotroph Nitrobacter. J. Bacteriol. 90:102-107. 1965.-The effect of various organic compounds on the growth and metabolism of the obligatively autotrophic nitrifying organism Nitrobacter was studied. A slight stimulation of both nitrification and growth was obtainable with a number of organic amendments, including yeast extract, Vitamin Free Casamino Acids, and some amino acids. Depending upon culture conditions, a strong stimulation of growth was obtained with acetate as an amendment to the culture solution. Several compounds, including valine, hydroxyproline, and threonine, were inhibitory at concentrations of 10(-3)m. The incorporation of carbon from isotopically labeled organic compounds was demonstrated. Acetate and glycine were particularly strong contributors to cell carbons. These could not substitute for carbon dioxide as a sole carbon source for growth, however, nor could any other of the carbon compounds that were tried.     

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.     

Optimization of 13C isotopic tracers for metabolic flux analysis in mammalian cells

Mammalian cells consume and metabolize various substrates from their surroundings for energy generation and biomass synthesis. Glucose and glutamine, in particular, are the primary carbon sources for proliferating cancer cells. While this combination of substrates generates static labeling patterns for use in (13)C metabolic flux analysis (MFA), the inability of single tracers to effectively label all pathways poses an obstacle for comprehensive flux determination within a given experiment. To address this issue we applied a genetic algorithm to optimize mixtures of (13)C-labeled glucose and glutamine for use in MFA. We identified tracer combinations that minimized confidence intervals in an experimentally determined flux network describing central carbon metabolism in tumor cells. Additional simulations were used to determine the robustness of the [1,2-(13)C(2)]glucose/[U-(13)C(5)]glutamine tracer combination with respect to perturbations in the network. Finally, we experimentally validated the improved performance of this tracer set relative to glucose tracers alone in a cancer cell line. This versatile method allows researchers to determine the optimal tracer combination to use for a specific metabolic network, and our findings applied to cancer cells significantly enhance the ability of MFA experiments to precisely quantify fluxes in higher organisms.