Alberto Noronha

Generation of genome-scale metabolic reconstructions for 773 members of the human gut microbiota

Genome-scale metabolic models derived from human gut metagenomic data can be used as a framework to elucidate how microbial communities modulate human metabolism and health. We present AGORA (assembly of gut organisms through reconstruction and analysis), a resource of genome-scale metabolic reconstructions semi-automatically generated for 773 human gut bacteria. Using this resource, we identified a defined growth medium for Bacteroides caccae ATCC 34185. We also showed that interactions among modeled species depend on both the metabolic potential of each species and the nutrients available. AGORA reconstructions can integrate either metagenomic or 16S rRNA sequencing data sets to infer the metabolic diversity of microbial communities. AGORA reconstructions could provide a starting point for the generation of high-quality, manually curated metabolic reconstructions. AGORA is fully compatible with Recon 2, a comprehensive metabolic reconstruction of human metabolism, which will facilitate studies of host–microbiome interactions.

Do genome-scale models need exact solvers or clearer standards?

Constraint‐based analysis of genome‐scale models (GEMs) arose shortly after the first genome sequences became available. As numerous reviews of the field show, this approach and methodology has proven to be successful in studying a wide range of biological phenomena (McCloskey et al, 2013; Bordbar et al, 2014). However, efforts to expand the user base are impeded by hurdles in correctly formulating these problems to obtain numerical solutions. In particular, in a study entitled “An exact arithmetic toolbox for a consistent and reproducible structural analysis of metabolic network models” (Chindelevitch et al, 2014), the authors apply an exact solver to 88 genome‐scale constraint‐based models of metabolism. The authors claim that COBRA calculations (Orth et al, 2010) are inconsistent with their results and that many published and actively used (Lee et al, 2007; McCloskey et al, 2013) genome‐scale models do support cellular growth in existing studies only because of numerical errors. They base these broad claims on two observations: (i) three reconstructions (iAF1260, iIT341, and iNJ661) compute feasibly in COBRA, but are infeasible when exact numerical algorithms are used by their software (entitled MONGOOSE); (ii) linear programs generated by MONGOOSE for iIT341 were submitted to the NEOS Server (a Web site that runs linear programs through various solvers) and gave inconsistent results. They further claim that a large percentage of these COBRA models are actually unable to produce biomass flux. Here, we demonstrate that the claims made by Chindelevitch et al (2014) stem from an incorrect parsing of models from files rather than actual problems with numerical error or COBRA computations.

Recon3D enables a three-dimensional view of gene variation in human metabolism

Genome-scale network reconstructions have helped uncover the molecular basis of metabolism. Here we present Recon3D, a computational resource that includes three-dimensional (3D) metabolite and protein structure data and enables integrated analyses of metabolic functions in humans. We use Recon3D to functionally characterize mutations associated with disease, and identify metabolic response signatures that are caused by exposure to certain drugs. Recon3D represents the most comprehensive human metabolic network model to date, accounting for 3,288 open reading frames (representing 17% of functionally annotated human genes), 13,543 metabolic reactions involving 4,140 unique metabolites, and 12,890 protein structures. These data provide a unique resource for investigating molecular mechanisms of human metabolism. Recon3D is available at http://vmh.life.

ReconMap: an interactive visualization of human metabolism.

Motivation: A genome‐scale reconstruction of human metabolism, Recon 2, is available but no interface exists to interactively visualize its content integrated with omics data and simulation results. Results: We manually drew a comprehensive map, ReconMap 2.0, that is consistent with the content of Recon 2. We present it within a web interface that allows content query, visualization of custom datasets and submission of feedback to manual curators. Availability and Implementation: ReconMap can be accessed via http://vmh.uni.lu, with network export in a Systems Biology Graphical Notation compliant format released under a Creative Commons Attribution‐NonCommercial‐NoDerivatives 4.0 International License. A Constraint‐Based Reconstruction and Analysis (COBRA) Toolbox extension to interact with ReconMap is available via https://github.com/opencobra/cobratoolbox. Contact: ronan.mt.fleming@gmail.com

Leigh map: A novel computational diagnostic resource for mitochondrial disease

Mitochondrial disorders are amongst the most severe metabolic disorders and are beset by genetic, biochemical, and clinical heterogeneity. Variation between individuals and poor understanding of disease pathophysiology pose significant diagnostic challenges. We present a novel interactive computational network, the Leigh Map, cataloguing >1700 gene-to-phenotype interactions in Leigh syndrome, the most common and genetically heterogeneous mitochondrial disorder. Blinded validation of the Leigh Map yielded an 80% success rate in correct identification of causative genes. We conclude that the Leigh Map is an efficacious resource that, in combination with whole-exome sequencing, can be utilized as a novel diagnostic resource for mitochondrial disease. This article is protected by copyright. All rights reserved.

An integrated network visualization framework towards metabolic engineering applications

BackgroundOver the last years, several methods for the phenotype simulation of microorganisms, under specified genetic and environmental conditions have been proposed, in the context of Metabolic Engineering (ME). These methods provided insight on the functioning of microbial metabolism and played a key role in the design of genetic modifications that can lead to strains of industrial interest. On the other hand, in the context of Systems Biology research, biological network visualization has reinforced its role as a core tool in understanding biological processes. However, it has been scarcely used to foster ME related methods, in spite of the acknowledged potential.ResultsIn this work, an open-source software that aims to fill the gap between ME and metabolic network visualization is proposed, in the form of a plugin to the OptFlux ME platform. The framework is based on an abstract layer, where the network is represented as a bipartite graph containing minimal information about the underlying entities and their desired relative placement. The framework provides input/output support for networks specified in standard formats, such as XGMML, SBGN or SBML, providing a connection to genome-scale metabolic models. An user-interface makes it possible to edit, manipulate and query nodes in the network, providing tools to visualize diverse effects, including visual filters and aspect changing (e.g. colors, shapes and sizes). These tools are particularly interesting for ME, since they allow overlaying phenotype simulation results or elementary flux modes over the networks.ConclusionsThe framework and its source code are freely available, together with documentation and other resources, being illustrated with well documented case studies.

When metabolism meets physiology: Harvey and Harvetta

Precision medicine is an emerging paradigm that requires realistic, mechanistic models capturing the complexity of the human body. We present two comprehensive molecular to physiological-level, gender-specific whole-body metabolism (WBM) reconstructions, named Harvey, in recognition of William Harvey, and Harvetta. These validated, knowledge-based WBM reconstructions capture the metabolism of 20 organs, six sex organs, six blood cells, the gastrointestinal lumen, systemic blood circulation, and the blood-brain barrier. They represent 99% of the human body weight, when excluding the weight of the skeleton. Harvey and Harvetta can be parameterized based on physiological, dietary, and omics data. They correctly predict inter-organ metabolic cycles, basal metabolic rates, and energy use. We demonstrate the integration of microbiome data thereby allowing the assessment of individual-specific, organ-level modulation of host metabolism by the gut microbiota. The WBM reconstructions and the individual organ reconstructions are available under http://vmh.life. Harvey and Harvetta represent a pivotal step towards virtual physiological humans.

Comparative evaluation of atom mapping algorithms for balanced metabolic reactions: application to Recon 3D

The mechanism of each chemical reaction in a metabolic network can be represented as a set of atom mappings, each of which relates an atom in a substrate metabolite to an atom of the same element in a product metabolite. Genome-scale metabolic network reconstructions typically represent biochemistry at the level of reaction stoichiometry. However, a more detailed representation at the underlying level of atom mappings opens the possibility for a broader range of biological, biomedical and biotechnological applications than with stoichiometry alone. Complete manual acquisition of atom mapping data for a genome-scale metabolic network is a laborious process. However, many algorithms exist to predict atom mappings. How do their predictions compare to each other and to manually curated atom mappings? For more than four thousand metabolic reactions in the latest human metabolic reconstruction, Recon 3D, we compared the atom mappings predicted by six atom mapping algorithms. We also compared these predictions to those obtained by manual curation of atom mappings for over five hundred reactions distributed among all top level Enzyme Commission number classes. Five of the evaluated algorithms had similarly high prediction accuracy of over 91% when compared to manually curated atom mapped reactions. On average, the accuracy of the prediction was highest for reactions catalysed by oxidoreductases and lowest for reactions catalysed by ligases. In addition to prediction accuracy, the algorithms were evaluated on their accessibility, their advanced features, such as the ability to identify equivalent atoms, and their ability to map hydrogen atoms. In addition to prediction accuracy, we found that software accessibility and advanced features were fundamental to the selection of an atom mapping algorithm in practice.

Network Visualization Tools to Enhance Metabolic Engineering Platforms

In this work, we present a software platform for the visualization of metabolic models, which is implemented as a plug-in for the open-source metabolic engineering (ME) platform OptFlux. The tools provided by this plug-in allow the visualization of the models (or parts of the models) combined with the results from operations applied over these models, mainly regarding phenotype simulation, strain optimization and pathway analysis. The tool provides a generic input/ output framework that can import/ export layouts from different formats used by other tools, namely XGMML and SBML. Thus, this work provides a bridge between network visualization and ME.

The Virtual Metabolic Human database: integrating human and gut microbiome metabolism with nutrition and disease

A multitude of factors contribute to complex diseases and can be measured with “omics” methods. Databases facilitate data interpretation for underlying mechanisms. Here, we describe the Virtual Metabolic Human (VMH, http://vmh.life) database encapsulating current knowledge of human metabolism within five interlinked resources “Human metabolism”, “Gut microbiome”, “Disease”, “Nutrition”, and “ReconMaps”. The VMH captures 5,180 unique metabolites, 17,730 unique reactions, 3,288 human genes, 255 Mendelian diseases, 818 microbes, 632,685 microbial genes, and 8,790 food items. The VMH’s unique features are i) the hosting the metabolic reconstructions of human and gut microbes amenable for metabolic modeling; ii) seven human metabolic maps for data visualization; iii) a nutrition designer; iv) a user-friendly webpage and application-programming interface to access its content; and v) user feedback option for community engagement. We demonstrate with four examples the VMH’s utility. The VMH represents a novel, interdisciplinary database for data interpretation and hypothesis generation to the biomedical community.