Project description:Caenorhabditis elegans is a powerful model to study metabolism and how it relates to nutrition, gene expression, and life history traits. However, while numerous experimental techniques that enable perturbation of its diet and gene function are available, a high-quality metabolic network model has been lacking. Here, we reconstruct an initial version of the C. elegans metabolic network. This network model contains 1,273 genes, 623 enzymes, and 1,985 metabolic reactions and is referred to as iCEL1273. Using flux balance analysis, we show that iCEL1273 is capable of representing the conversion of bacterial biomass into C. elegans biomass during growth and enables the predictions of gene essentiality and other phenotypes. In addition, we demonstrate that gene expression data can be integrated with the model by comparing metabolic rewiring in dauer animals versus growing larvae. iCEL1273 is available at a dedicated website (wormflux.umassmed.edu) and will enable the unraveling of the mechanisms by which different macro- and micronutrients contribute to the animal's physiology.
Project description:Genome-scale metabolic models are widely used to enhance our understanding of metabolic features of organisms, host-pathogen interactions and to identify therapeutics for diseases. Here we present iTMU798, the genome-scale metabolic model of the mouse whipworm Trichuris muris. The model demonstrates the metabolic features of T. muris and allows the prediction of metabolic steps essential for its survival. Specifically, that Thioredoxin Reductase (TrxR) enzyme is essential, a prediction we validate in vitro with the drug auranofin. Furthermore, our observation that the T. muris genome lacks gsr-1 encoding Glutathione Reductase (GR) but has GR activity that can be inhibited by auranofin indicates a mechanism for the reduction of glutathione by the TrxR enzyme in T. muris. In addition, iTMU798 predicts seven essential amino acids that cannot be synthesised by T. muris, a prediction we validate for the amino acid tryptophan. Overall, iTMU798 is as a powerful tool to study not only the T. muris metabolism but also other Trichuris spp. in understanding host parasite interactions and the rationale design of new intervention strategies.
Project description:Genome scale metabolic models (GSMs) are a representation of the current knowledge on the metabolism of a given organism or superorganism. They group metabolites, genes, enzymes and reactions together to form a mathematical model and representation that can be used to analyze metabolic networks in silico or used for analysis of omics data. Beside correct mass and charge balance, correct structural annotation of metabolites represents an important factor for analysis of these metabolic networks. However, several metabolites in different GSMs have no or only partial structural information associated with them. Here, a new systematic nomenclature for acyl-based metabolites such as fatty acids, acyl-carnitines, acyl-coenzymes A or acyl-carrier proteins is presented. This nomenclature enables one to encode structural details in the metabolite identifiers and improves human readability of reactions. As proof of principle, it was applied to the fatty acid biosynthesis and degradation in the Caenorhabditis elegans consensus model WormJam.
Project description:BACKGROUND: Fermentation of xylose, the major component in hemicellulose, is essential for economic conversion of lignocellulosic biomass to fuels and chemicals. The yeast Scheffersomyces stipitis (formerly known as Pichia stipitis) has the highest known native capacity for xylose fermentation and possesses several genes for lignocellulose bioconversion in its genome. Understanding the metabolism of this yeast at a global scale, by reconstructing the genome scale metabolic model, is essential for manipulating its metabolic capabilities and for successful transfer of its capabilities to other industrial microbes. RESULTS: We present a genome-scale metabolic model for Scheffersomyces stipitis, a native xylose utilizing yeast. The model was reconstructed based on genome sequence annotation, detailed experimental investigation and known yeast physiology. Macromolecular composition of Scheffersomyces stipitis biomass was estimated experimentally and its ability to grow on different carbon, nitrogen, sulphur and phosphorus sources was determined by phenotype microarrays. The compartmentalized model, developed based on an iterative procedure, accounted for 814 genes, 1371 reactions, and 971 metabolites. In silico computed growth rates were compared with high-throughput phenotyping data and the model could predict the qualitative outcomes in 74% of substrates investigated. Model simulations were used to identify the biosynthetic requirements for anaerobic growth of Scheffersomyces stipitis on glucose and the results were validated with published literature. The bottlenecks in Scheffersomyces stipitis metabolic network for xylose uptake and nucleotide cofactor recycling were identified by in silico flux variability analysis. The scope of the model in enhancing the mechanistic understanding of microbial metabolism is demonstrated by identifying a mechanism for mitochondrial respiration and oxidative phosphorylation. CONCLUSION: The genome-scale metabolic model developed for Scheffersomyces stipitis successfully predicted substrate utilization and anaerobic growth requirements. Useful insights were drawn on xylose metabolism, cofactor recycling and mechanism of mitochondrial respiration from model simulations. These insights can be applied for efficient xylose utilization and cofactor recycling in other industrial microorganisms. The developed model forms a basis for rational analysis and design of Scheffersomyces stipitis metabolic network for the production of fuels and chemicals from lignocellulosic biomass.
Project description:BackgroundSpirulina (Arthrospira) platensis is a well-known filamentous cyanobacterium used in the production of many industrial products, including high value compounds, healthy food supplements, animal feeds, pharmaceuticals and cosmetics, for example. It has been increasingly studied around the world for scientific purposes, especially for its genome, biology, physiology, and also for the analysis of its small-scale metabolic network. However, the overall description of the metabolic and biotechnological capabilities of S. platensis requires the development of a whole cellular metabolism model. Recently, the S. platensis C1 (Arthrospira sp. PCC9438) genome sequence has become available, allowing systems-level studies of this commercial cyanobacterium.ResultsIn this work, we present the genome-scale metabolic network analysis of S. platensis C1, iAK692, its topological properties, and its metabolic capabilities and functions. The network was reconstructed from the S. platensis C1 annotated genomic sequence using Pathway Tools software to generate a preliminary network. Then, manual curation was performed based on a collective knowledge base and a combination of genomic, biochemical, and physiological information. The genome-scale metabolic model consists of 692 genes, 837 metabolites, and 875 reactions. We validated iAK692 by conducting fermentation experiments and simulating the model under autotrophic, heterotrophic, and mixotrophic growth conditions using COBRA toolbox. The model predictions under these growth conditions were consistent with the experimental results. The iAK692 model was further used to predict the unique active reactions and essential genes for each growth condition. Additionally, the metabolic states of iAK692 during autotrophic and mixotrophic growths were described by phenotypic phase plane (PhPP) analysis.ConclusionsThis study proposes the first genome-scale model of S. platensis C1, iAK692, which is a predictive metabolic platform for a global understanding of physiological behaviors and metabolic engineering. This platform could accelerate the integrative analysis of various "-omics" data, leading to strain improvement towards a diverse range of desired industrial products from Spirulina.
Project description:BackgroundTo date, several genome-scale network reconstructions have been used to describe the metabolism of the yeast Saccharomyces cerevisiae, each differing in scope and content. The recent community-driven reconstruction, while rigorously evidenced and well annotated, under-represented metabolite transport, lipid metabolism and other pathways, and was not amenable to constraint-based analyses because of lack of pathway connectivity.ResultsWe have expanded the yeast network reconstruction to incorporate many new reactions from the literature and represented these in a well-annotated and standards-compliant manner. The new reconstruction comprises 1102 unique metabolic reactions involving 924 unique metabolites--significantly larger in scope than any previous reconstruction. The representation of lipid metabolism in particular has improved, with 234 out of 268 enzymes linked to lipid metabolism now present in at least one reaction. Connectivity is emphatically improved, with more than 90% of metabolites now reachable from the growth medium constituents. The present updates allow constraint-based analyses to be performed; viability predictions of single knockouts are comparable to results from in vivo experiments and to those of previous reconstructions.ConclusionsWe report the development of the most complete reconstruction of yeast metabolism to date that is based upon reliable literature evidence and richly annotated according to MIRIAM standards. The reconstruction is available in the Systems Biology Markup Language (SBML) and via a publicly accessible database http://www.comp-sys-bio.org/yeastnet/.
Project description:The Gram-negative bacterium Bordetella pertussis is the causative agent of whooping cough, a serious respiratory infection causing hundreds of thousands of deaths annually worldwide. There are effective vaccines, but their production requires growing large quantities of B. pertussis. Unfortunately, B. pertussis has relatively slow growth in culture, with low biomass yields and variable growth characteristics. B. pertussis also requires a relatively expensive growth medium. We present a new, curated flux balance analysis-based model of B. pertussis metabolism. We enhance the model with an experimentally-determined biomass objective function, and we perform extensive manual curation. We test the model's predictions with a genome-wide screen for essential genes using a transposon-directed insertional sequencing (TraDIS) approach. We test its predictions of growth for different carbon sources in the medium. The model predicts essentiality with an accuracy of 83% and correctly predicts improvements in growth under increased glutamate:fumarate ratios. We provide the model in SBML format, along with gene essentiality predictions.
Project description:BackgroundMethylocella silvestris is a facultative aerobic methanotrophic bacterium which uses not only methane, but also other alkanes such as ethane and propane, as carbon and energy sources. Its high metabolic versatility, together with the availability of tools for its genetic engineering, make it a very promising platform for metabolic engineering and industrial biotechnology using natural gas as substrate.ResultsThe first Genome Scale Metabolic Model for M. silvestris is presented. The model has been used to predict the ability of M. silvestris to grow on 12 different substrates, the growth phenotype of two deletion mutants (ΔICL and ΔMS), and biomass yield on methane and ethanol. The model, together with phenotypic characterization of the deletion mutants, revealed that M. silvestris uses the glyoxylate shuttle for the assimilation of C1 and C2 substrates, which is unique in contrast to published reports of other methanotrophs. Two alternative pathways for propane metabolism have been identified and validated experimentally using enzyme activity tests and constructing a deletion mutant (Δ1641), which enabled the identification of acetol as one of the intermediates of propane assimilation via 2-propanol. The model was also used to integrate proteomic data and to identify key enzymes responsible for the adaptation of M. silvestris to different substrates.ConclusionsThe model has been used to elucidate key metabolic features of M. silvestris, such as its use of the glyoxylate shuttle for the assimilation of one and two carbon compounds and the existence of two parallel metabolic pathways for propane assimilation. This model, together with the fact that tools for its genetic engineering already exist, paves the way for the use of M. silvestris as a platform for metabolic engineering and industrial exploitation of methanotrophs.
Project description: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.
Project description:Since the first in silico generation of a genome-scale metabolic (GSM) model for Haemophilus influenzae in 1999, the GSM models have been reconstructed for various organisms including human and mouse. There are two important strategies for generating a GSM model: in the bottom-up approach, individual genomic and biochemical components are integrated to build a GSM model. Alternatively, the orthology-based strategy uses a previously reconstructed model of a reference organism to infer a GSM model of a target organism. Following the update and development of the metabolic network of reference organism, the model of the target organism can also be updated to eliminate defects. Here, we presented iMM1865 model as an orthology-based reconstruction of a GSM model for Mus musculus based on the last flux-consistent version of the human metabolic network, Recon3D. We proposed two versions of the new mouse model, iMM1865 and min-iMM1865, with the same number of gene-associated reactions but different subsets of non-gene-associated reactions. A third extended but flux-inconsistent model (iMM3254) was also created based on the extended version of Recon3D. Compared to the previously published mouse models, both versions of iMM1865 include more comprehensive annotations of metabolites and reactions with no dead-end metabolites and blocked reactions. We evaluated functionality of the models using 431 metabolic objective functions. iMM1865 and min-iMM1865 passed 93% and 87% of the tests, respectively, while iMM1415 and MMR (another available mouse GSM) passed 80% and 84% of the tests, respectively. Three versions of tissue-specific embryo heart models were also reconstructed from each of iMM1865 and min-iMM1865 using mCADRE algorithm with different thresholds on expression-based scores. The ability of corresponding GSM and embryo heart models to predict essential genes was assessed across experimentally derived lethal and viable gene sets. Our analysis revealed that tissue-specific models render much better predictions than GSM models.