merlin,an improved framework for the reconstruction of high-quality genome-scale metabolic models

Genome-scale metabolic models have been recognised as useful tools for better understanding living organisms’ metabolism. merlin (https://www.merlin-sysbio.org/) is an open-source and user-friendly resource that hastens the models’ reconstruction process, conjugating manual and automatic procedures, while leveraging the user''s expertise with a curation-oriented graphical interface. An updated and redesigned version of merlin is herein presented. Since 2015, several features have been implemented in merlin, along with deep changes in the software architecture, operational flow, and graphical interface. The current version (4.0) includes the implementation of novel algorithms and third-party tools for genome functional annotation, draft assembly, model refinement, and curation. Such updates increased the user base, resulting in multiple published works, including genome metabolic (re-)annotations and model reconstructions of multiple (lower and higher) eukaryotes and prokaryotes. merlin version 4.0 is the only tool able to perform template based and de novo draft reconstructions, while achieving competitive performance compared to state-of-the art tools both for well and less-studied organisms.     

Semi-automated Curation of Metabolic Models via Flux Balance Analysis: A Case Study with Mycoplasma gallisepticum

Primarily used for metabolic engineering and synthetic biology, genome-scale metabolic modeling shows tremendous potential as a tool for fundamental research and curation of metabolism. Through a novel integration of flux balance analysis and genetic algorithms, a strategy to curate metabolic networks and facilitate identification of metabolic pathways that may not be directly inferable solely from genome annotation was developed. Specifically, metabolites involved in unknown reactions can be determined, and potentially erroneous pathways can be identified. The procedure developed allows for new fundamental insight into metabolism, as well as acting as a semi-automated curation methodology for genome-scale metabolic modeling. To validate the methodology, a genome-scale metabolic model for the bacterium Mycoplasma gallisepticum was created. Several reactions not predicted by the genome annotation were postulated and validated via the literature. The model predicted an average growth rate of 0.358±0.12, closely matching the experimentally determined growth rate of M. gallisepticum of 0.244±0.03. This work presents a powerful algorithm for facilitating the identification and curation of previously known and new metabolic pathways, as well as presenting the first genome-scale reconstruction of M. gallisepticum.     

Genome-scale reconstruction of metabolic network in Bacillus subtilis based on high-throughput phenotyping and gene essentiality data

In this report, a genome-scale reconstruction of Bacillus subtilis metabolism and its iterative development based on the combination of genomic, biochemical, and physiological information and high-throughput phenotyping experiments is presented. The initial reconstruction was converted into an in silico model and expanded in a four-step iterative fashion. First, network gap analysis was used to identify 48 missing reactions that are needed for growth but were not found in the genome annotation. Second, the computed growth rates under aerobic conditions were compared with high-throughput phenotypic screen data, and the initial in silico model could predict the outcomes qualitatively in 140 of 271 cases considered. Detailed analysis of the incorrect predictions resulted in the addition of 75 reactions to the initial reconstruction, and 200 of 271 cases were correctly computed. Third, in silico computations of the growth phenotypes of knock-out strains were found to be consistent with experimental observations in 720 of 766 cases evaluated. Fourth, the integrated analysis of the large-scale substrate utilization and gene essentiality data with the genome-scale metabolic model revealed the requirement of 80 specific enzymes (transport, 53; intracellular reactions, 27) that were not in the genome annotation. Subsequent sequence analysis resulted in the identification of genes that could be putatively assigned to 13 intracellular enzymes. The final reconstruction accounted for 844 open reading frames and consisted of 1020 metabolic reactions and 988 metabolites. Hence, the in silico model can be used to obtain experimentally verifiable hypothesis on the metabolic functions of various genes.     

Reconstruction of Danio rerio Metabolic Model Accounting for Subcellular Compartmentalisation

Plant and microbial metabolic engineering is commonly used in the production of functional foods and quality trait improvement. Computational model-based approaches have been used in this important endeavour. However, to date, fish metabolic models have only been scarcely and partially developed, in marked contrast to their prominent success in metabolic engineering. In this study we present the reconstruction of fully compartmentalised models of the Danio rerio (zebrafish) on a global scale. This reconstruction involves extraction of known biochemical reactions in D. rerio for both primary and secondary metabolism and the implementation of methods for determining subcellular localisation and assignment of enzymes. The reconstructed model (ZebraGEM) is amenable for constraint-based modelling analysis, and accounts for 4,988 genes coding for 2,406 gene-associated reactions and only 418 non-gene-associated reactions. A set of computational validations (i.e., simulations of known metabolic functionalities and experimental data) strongly testifies to the predictive ability of the model. Overall, the reconstructed model is expected to lay down the foundations for computational-based rational design of fish metabolic engineering in aquaculture.     

Quantifying cumulative phenotypic and genomic evidence for procedural generation of metabolic network reconstructions

Genome-scale metabolic network reconstructions (GENREs) are valuable tools for understanding microbial metabolism. The process of automatically generating GENREs includes identifying metabolic reactions supported by sufficient genomic evidence to generate a draft metabolic network. The draft GENRE is then gapfilled with additional reactions in order to recapitulate specific growth phenotypes as indicated with associated experimental data. Previous methods have implemented absolute mapping thresholds for the reactions automatically included in draft GENREs; however, there is growing evidence that integrating annotation evidence in a continuous form can improve model accuracy. There is a need for flexibility in the structure of GENREs to better account for uncertainty in biological data, unknown regulatory mechanisms, and context-specificity associated with data inputs. To address this issue, we present a novel method that provides a framework for quantifying combined genomic, biochemical, and phenotypic evidence for each biochemical reaction during automated GENRE construction. Our method, Constraint-based Analysis Yielding reaction Usage across metabolic Networks (CANYUNs), generates accurate GENREs with a quantitative metric for the cumulative evidence for each reaction included in the network. The structuring of CANYUNs allows for the simultaneous integration of three data inputs while maintaining all supporting evidence for biochemical reactions that may be active in an organism. CANYUNs is designed to maximize the utility of experimental and annotation datasets and to ultimately assist in the curation of the reference datasets used for the automatic construction of metabolic networks. We validated CANYUNs by generating an E. coli K-12 model and compared it to the manually curated reconstruction iML1515. Finally, we demonstrated the use of CANYUNs to build a model by generating an E. coli Nissle CANYUNs model using novel phenotypic data that we collected. This method may address key challenges for the procedural construction of metabolic networks by leveraging uncertainty and redundancy in biological data.     

IdentiCS – Identification of coding sequence and <Emphasis Type="Italic">in silico</Emphasis> reconstruction of the metabolic network directly from unannotated low-coverage bacterial genome sequence

Background  

A necessary step for a genome level analysis of the cellular metabolism is the in silico reconstruction of the metabolic network from genome sequences. The available methods are mainly based on the annotation of genome sequences including two successive steps, the prediction of coding sequences (CDS) and their function assignment. The annotation process takes time. The available methods often encounter difficulties when dealing with unfinished error-containing genomic sequence.     

Evaluation of a Genome-Scale In Silico Metabolic Model for Geobacter metallireducens by Using Proteomic Data from a Field Biostimulation Experiment

Accurately predicting the interactions between microbial metabolism and the physical subsurface environment is necessary to enhance subsurface energy development, soil and groundwater cleanup, and carbon management. This study was an initial attempt to confirm the metabolic functional roles within an in silico model using environmental proteomic data collected during field experiments. Shotgun global proteomics data collected during a subsurface biostimulation experiment were used to validate a genome-scale metabolic model of Geobacter metallireducens—specifically, the ability of the metabolic model to predict metal reduction, biomass yield, and growth rate under dynamic field conditions. The constraint-based in silico model of G. metallireducens relates an annotated genome sequence to the physiological functions with 697 reactions controlled by 747 enzyme-coding genes. Proteomic analysis showed that 180 of the 637 G. metallireducens proteins detected during the 2008 experiment were associated with specific metabolic reactions in the in silico model. When the field-calibrated Fe(III) terminal electron acceptor process reaction in a reactive transport model for the field experiments was replaced with the genome-scale model, the model predicted that the largest metabolic fluxes through the in silico model reactions generally correspond to the highest abundances of proteins that catalyze those reactions. Central metabolism predicted by the model agrees well with protein abundance profiles inferred from proteomic analysis. Model discrepancies with the proteomic data, such as the relatively low abundances of proteins associated with amino acid transport and metabolism, revealed pathways or flux constraints in the in silico model that could be updated to more accurately predict metabolic processes that occur in the subsurface environment.     

Functional Metabolic Map of Faecalibacterium prausnitzii,a Beneficial Human Gut Microbe

The human gut microbiota plays a central role in human well-being and disease. In this study, we present an integrated, iterative approach of computational modeling, in vitro experiments, metabolomics, and genomic analysis to accelerate the identification of metabolic capabilities for poorly characterized (anaerobic) microorganisms. We demonstrate this approach for the beneficial human gut microbe Faecalibacterium prausnitzii strain A2-165. We generated an automated draft reconstruction, which we curated against the limited biochemical data. This reconstruction modeling was used to develop in silico and in vitro a chemically defined medium (CDM), which was validated experimentally. Subsequent metabolomic analysis of the spent medium for growth on CDM was performed. We refined our metabolic reconstruction according to in vitro observed metabolite consumption and secretion and propose improvements to the current genome annotation of F. prausnitzii A2-165. We then used the reconstruction to systematically characterize its metabolic properties. Novel carbon source utilization capabilities and inabilities were predicted based on metabolic modeling and validated experimentally. This study resulted in a functional metabolic map of F. prausnitzii, which is available for further applications. The presented workflow can be readily extended to other poorly characterized and uncharacterized organisms to yield novel biochemical insights about the target organism.     

Reconstruction of the central carbon metabolism of Aspergillus niger.

The topology of central carbon metabolism of Aspergillus niger was identified and the metabolic network reconstructed, by integrating genomic, biochemical and physiological information available for this microorganism and other related fungi. The reconstructed network may serve as a valuable database for annotation of genes identified in future genome sequencing projects on aspergilli. Based on the metabolic reconstruction, a stoichiometric model was set up that includes 284 metabolites and 335 reactions, of which 268 represent biochemical conversions and 67 represent transport processes between the different intracellular compartments and between the cell and the extracellular medium. The stoichiometry of the metabolic reactions was used in combination with biosynthetic requirements for growth and pseudo-steady state mass balances over intracellular metabolites for the quantification of metabolic fluxes using metabolite balancing. This framework was employed to perform an in silico characterisation of the phenotypic behaviour of A. niger grown on different carbon sources. The effects on growth of single reaction deletions were assessed and essential biochemical reactions were identified for different carbon sources. Furthermore, application of the stoichiometric model for assessing the metabolic capabilities of A. niger to produce metabolites was evaluated by using succinate production as a case study.     

Interaction between sulphur mobilisation proteins SufB and SufC: evidence for an iron-sulphur cluster biogenesis pathway in the apicoplast of Plasmodium falciparum

The plastid of Plasmodium falciparum, the apicoplast, performs metabolic functions essential to the parasite. Various reactions in the plastid require the assembly of [Fe-S] prosthetic groups on participating proteins as well as the reductant activity of ferredoxin that is converted from its apo-form by the assembly of [Fe-S] clusters inside the apicoplast. The [Fe-S] assembly pathway involving sulphur mobilising Suf proteins has been predicted to function in the apicoplast with one component (PfSufB) encoded by the plastid genome itself. We demonstrate the ATPase activity of recombinant P. falciparum nuclear-encoded SufC and its localisation in the apicoplast. Further, an internal region of apicoplast SufB was used to detect PfSufB-PfSufC interaction in vitro; co-elution of SufB from parasite lysate with recombinant PfSufC on an affinity column also indicated an interaction of the two proteins. As a departure from bacterial SufB and similar to reported plant plastid SufB, apicoplast SufB exhibited ATPase activity, suggesting the evolution of specialised functions in the plastid counterparts. Our results provide experimental evidence for an active Suf pathway in the Plasmodium apicoplast.     

Shared Metabolic Pathways in a Coevolved Insect-Bacterial Symbiosis

The symbiotic bacterium Buchnera aphidicola lacks key genes in the biosynthesis of five essential amino acids (EAAs), and yet its animal hosts (aphids) depend on the symbiosis for the synthesis of these EAAs (isoleucine, leucine, methionine, phenylalanine, and valine). We tested the hypothesis, derived from genome annotation, that the missing Buchnera reactions are mediated by host enzymes, with the exchange of metabolic intermediates between the partners. The specialized host cells bearing Buchnera were separated into a Buchnera fraction and a Buchnera-free host cell fraction (HF). Addition of HF to isolated Buchnera preparations significantly increased the production of leucine and phenylalanine, and recombinant enzymes mediating the final reactions in branched-chain amino acid and phenylalanine synthesis rescued the production of these EAAs by Buchnera preparations without HF. The likely precursors for the missing proximal reactions in isoleucine and methionine synthesis were identified, and they differed from predictions based on genome annotations: synthesis of 2-oxobutanoate, the aphid-derived precursor of isoleucine synthesis, was stimulated by homoserine and not threonine via threonine dehydratase, and production of the homocysteine precursor of methionine was driven by cystathionine, not cysteine, via reversal of the transsulfuration pathway. The evolution of shared metabolic pathways in this symbiosis can be attributed to host compensation for genomic deterioration in the symbiont, involving changes in host gene expression networks to recruit specific enzymes to the host cell.     

Genomically and biochemically accurate metabolic reconstruction of Methanosarcina barkeri Fusaro,iMG746

Methanosarcina barkeri is an Archaeon that produces methane anaerobically as the primary byproduct of its metabolism. M. barkeri can utilize several substrates for ATP and biomass production including methanol, acetate, methyl amines, and a combination of hydrogen and carbon dioxide. In 2006, a metabolic reconstruction of M. barkeri, iAF692, was generated based on a draft genome annotation. The iAF692 reconstruction enabled the first genome-Scale simulations for Archaea. Since the publication of the first metabolic reconstruction of M. barkeri, additional genomic, biochemical, and phenotypic data have clarified several metabolic pathways. We have used this newly available data to improve the M. barkeri metabolic reconstruction. Modeling simulations using the updated model, iMG746, have led to increased accuracy in predicting gene knockout phenotypes and simulations of batch growth behavior. We used the model to examine knockout lethality data and make predictions about metabolic regulation under different growth conditions. Thus, the updated metabolic reconstruction of M. barkeri metabolism is a useful tool for predicting cellular behavior, studying the methanogenic lifestyle, guiding experimental studies, and making predictions relevant to metabolic engineering applications.     

Reconstruction and In Silico Analysis of Metabolic Network for an Oleaginous Yeast,Yarrowia lipolytica

With the emergence of energy scarcity, the use of renewable energy sources such as biodiesel is becoming increasingly necessary. Recently, many researchers have focused their minds on Yarrowia lipolytica, a model oleaginous yeast, which can be employed to accumulate large amounts of lipids that could be further converted to biodiesel. In order to understand the metabolic characteristics of Y. lipolytica at a systems level and to examine the potential for enhanced lipid production, a genome-scale compartmentalized metabolic network was reconstructed based on a combination of genome annotation and the detailed biochemical knowledge from multiple databases such as KEGG, ENZYME and BIGG. The information about protein and reaction associations of all the organisms in KEGG and Expasy-ENZYME database was arranged into an EXCEL file that can then be regarded as a new useful database to generate other reconstructions. The generated model iYL619_PCP accounts for 619 genes, 843 metabolites and 1,142 reactions including 236 transport reactions, 125 exchange reactions and 13 spontaneous reactions. The in silico model successfully predicted the minimal media and the growing abilities on different substrates. With flux balance analysis, single gene knockouts were also simulated to predict the essential genes and partially essential genes. In addition, flux variability analysis was applied to design new mutant strains that will redirect fluxes through the network and may enhance the production of lipid. This genome-scale metabolic model of Y. lipolytica can facilitate system-level metabolic analysis as well as strain development for improving the production of biodiesels and other valuable products by Y. lipolytica and other closely related oleaginous yeasts.     

Comparative Analysis of Yeast Metabolic Network Models Highlights Progress,Opportunities for Metabolic Reconstruction

We have compared 12 genome-scale models of the Saccharomyces cerevisiae metabolic network published since 2003 to evaluate progress in reconstruction of the yeast metabolic network. We compared the genomic coverage, overlap of annotated metabolites, predictive ability for single gene essentiality with a selection of model parameters, and biomass production predictions in simulated nutrient-limited conditions. We have also compared pairwise gene knockout essentiality predictions for 10 of these models. We found that varying approaches to model scope and annotation reflected the involvement of multiple research groups in model development; that single-gene essentiality predictions were affected by simulated medium, objective function, and the reference list of essential genes; and that predictive ability for single-gene essentiality did not correlate well with predictive ability for our reference list of synthetic lethal gene interactions (R = 0.159). We conclude that the reconstruction of the yeast metabolic network is indeed gradually improving through the iterative process of model development, and there remains great opportunity for advancing our understanding of biology through continued efforts to reconstruct the full biochemical reaction network that constitutes yeast metabolism. Additionally, we suggest that there is opportunity for refining the process of deriving a metabolic model from a metabolic network reconstruction to facilitate mechanistic investigation and discovery. This comparative study lays the groundwork for developing improved tools and formalized methods to quantitatively assess metabolic network reconstructions independently of any particular model application, which will facilitate ongoing efforts to advance our understanding of the relationship between genotype and cellular phenotype.     

MitoCore: a curated constraint-based model for simulating human central metabolism

Background

The complexity of metabolic networks can make the origin and impact of changes in central metabolism occurring during diseases difficult to understand. Computer simulations can help unravel this complexity, and progress has advanced in genome-scale metabolic models. However, many models produce unrealistic results when challenged to simulate abnormal metabolism as they include incorrect specification and localisation of reactions and transport steps, incorrect reaction parameters, and confounding of prosthetic groups and free metabolites in reactions. Other common drawbacks are due to their scale, making them difficult to parameterise and simulation results hard to interpret. Therefore, it remains important to develop smaller, manually curated models.

Results

We present MitoCore, a manually curated constraint-based computer model of human metabolism that incorporates the complexity of central metabolism and simulates this metabolism successfully under normal and abnormal physiological conditions, including hypoxia and mitochondrial diseases. MitoCore describes 324 metabolic reactions, 83 transport steps between mitochondrion and cytosol, and 74 metabolite inputs and outputs through the plasma membrane, to produce a model of manageable scale for easy interpretation of results. Its key innovations include a more accurate partitioning of metabolism between cytosol and mitochondrial matrix; better modelling of connecting transport steps; differentiation of prosthetic groups and free co-factors in reactions; and a new representation of the respiratory chain and the proton motive force. MitoCore’s default parameters simulate normal cardiomyocyte metabolism, and to improve usability and allow comparison with other models and types of analysis, its reactions and metabolites have extensive annotation, and cross-reference identifiers from Virtual Metabolic Human database and KEGG. These innovations—including over 100 reactions absent or modified from Recon 2—are necessary to model central metabolism more accurately.

Conclusion

We anticipate MitoCore as a research tool for scientists, from experimentalists looking to interpret their data and test hypotheses, to experienced modellers predicting the consequences of disease or using computationally intensive methods that are infeasible with larger models, as well as a teaching tool for those new to modelling and needing a small, manageable model on which to learn and experiment.

     

Genome-wide metabolic re-annotation of Ashbya gossypii: new insights into its metabolism through a comparative analysis with Saccharomyces cerevisiae and Kluyveromyces lactis

Background

Ashbya gossypii is an industrially relevant microorganism traditionally used for riboflavin production. Despite the high gene homology and gene order conservation comparatively with Saccharomyces cerevisiae, it presents a lower level of genomic complexity. Its type of growth, placing it among filamentous fungi, questions how close it really is from the budding yeast, namely in terms of metabolism, therefore raising the need for an extensive and thorough study of its entire metabolism. This work reports the first manual enzymatic genome-wide re-annotation of A. gossypii as well as the first annotation of membrane transport proteins.

Results

After applying a developed enzymatic re-annotation pipeline, 847 genes were assigned with metabolic functions. Comparatively to KEGG’s annotation, these data corrected the function for 14% of the common genes and increased the information for 52 genes, either completing existing partial EC numbers or adding new ones. Furthermore, 22 unreported enzymatic functions were found, corresponding to a significant increase in the knowledge of the metabolism of this organism. The information retrieved from the metabolic re-annotation and transport annotation was used for a comprehensive analysis of A. gossypii’s metabolism in comparison to the one of S. cerevisiae (post-WGD – whole genome duplication) and Kluyveromyces lactis (pre-WGD), suggesting some relevant differences in several parts of their metabolism, with the majority being found for the metabolism of purines, pyrimidines, nitrogen and lipids. A considerable number of enzymes were found exclusively in A. gossypii comparatively with K. lactis (90) and S. cerevisiae (13). In a similar way, 176 and 123 enzymatic functions were absent on A. gossypii comparatively to K. lactis and S. cerevisiae, respectively, confirming some of the well-known phenotypes of this organism.

Conclusions

This high quality metabolic re-annotation, together with the first membrane transporters annotation and the metabolic comparative analysis, represents a new important tool for the study and better understanding of A. gossypii’s metabolism.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-810) contains supplementary material, which is available to authorized users.     

Systematic assignment of thermodynamic constraints in metabolic network models

Background  

The availability of genome sequences for many organisms enabled the reconstruction of several genome-scale metabolic network models. Currently, significant efforts are put into the automated reconstruction of such models. For this, several computational tools have been developed that particularly assist in identifying and compiling the organism-specific lists of metabolic reactions. In contrast, the last step of the model reconstruction process, which is the definition of the thermodynamic constraints in terms of reaction directionalities, still needs to be done manually. No computational method exists that allows for an automated and systematic assignment of reaction directions in genome-scale models.     

metaSHARK: software for automated metabolic network prediction from DNA sequence and its application to the genomes of Plasmodium falciparum and Eimeria tenella

The metabolic SearcH And Reconstruction Kit (metaSHARK) is a new fully automated software package for the detection of enzyme-encoding genes within unannotated genome data and their visualization in the context of the surrounding metabolic network. The gene detection package (SHARKhunt) runs on a Linux system and requires only a set of raw DNA sequences (genomic, expressed sequence tag and/or genome survey sequence) as input. Its output may be uploaded to our web-based visualization tool (SHARKview) for exploring and comparing data from different organisms. We first demonstrate the utility of the software by comparing its results for the raw Plasmodium falciparum genome with the manual annotations available at the PlasmoDB and PlasmoCyc websites. We then apply SHARKhunt to the unannotated genome sequences of the coccidian parasite Eimeria tenella and observe that, at an E-value cut-off of 10⁻²⁰, our software makes 142 additional assertions of enzymatic function compared with a recent annotation package working with translated open reading frame sequences. The ability of the software to cope with low levels of sequence coverage is investigated by analyzing assemblies of the E.tenella genome at estimated coverages from 0.5× to 7.5×. Lastly, as an example of how metaSHARK can be used to evaluate the genomic evidence for specific metabolic pathways, we present a study of coenzyme A biosynthesis in P.falciparum and E.tenella.