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.     

Genome‐scale metabolic reconstruction and constraint‐based modelling of the Antarctic bacterium Pseudoalteromonas haloplanktis TAC125

The Antarctic strain Pseudoalteromonas haloplanktis TAC125 is one of the model organisms of cold‐adapted bacteria and is currently exploited as a new alternative expression host for numerous biotechnological applications. Here, we investigated several metabolic features of this strain through in silico modelling and functional integration of –omics data. A genome‐scale metabolic model of P. haloplanktis TAC125 was reconstructed, encompassing information on 721 genes, 1133 metabolites and 1322 reactions. The predictive potential of this model was validated against a set of experimentally determined growth rates and a large dataset of growth phenotypic data. Furthermore, evidence synthesis from proteomics, phenomics, physiology and metabolic modelling data revealed possible drawbacks of cold‐dependent changes in gene expression on the overall metabolic network of P. haloplanktis TAC125. These included, for example, variations in its central metabolism, amino acid degradation and fatty acid biosynthesis. The genome‐scale metabolic model described here is the first one reconstructed so far for an Antarctic microbial strain. It allowed a system‐level investigation of variations in cellular metabolic fluxes following a temperature downshift. It represents a valuable platform for further investigations on P. haloplanktis TAC125 cellular functional states and for the design of more focused strategies for its possible biotechnological exploitation.     

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.     

The genome‐scale metabolic network of Ectocarpus siliculosus (EctoGEM): a resource to study brown algal physiology and beyond

Brown algae (stramenopiles) are key players in intertidal ecosystems, and represent a source of biomass with several industrial applications. Ectocarpus siliculosus is a model to study the biology of these organisms. Its genome has been sequenced and a number of post‐genomic tools have been implemented. Based on this knowledge, we report the reconstruction and analysis of a genome‐scale metabolic network for E. siliculosus, EctoGEM ( http://ectogem.irisa.fr ). This atlas of metabolic pathways consists of 1866 reactions and 2020 metabolites, and its construction was performed by means of an integrative computational approach for identifying metabolic pathways, gap filling and manual refinement. The capability of the network to produce biomass was validated by flux balance analysis. EctoGEM enabled the reannotation of 56 genes within the E. siliculosus genome, and shed light on the evolution of metabolic processes. For example, E. siliculosus has the potential to produce phenylalanine and tyrosine from prephenate and arogenate, but does not possess a phenylalanine hydroxylase, as is found in other stramenopiles. It also possesses the complete eukaryote molybdenum co‐factor biosynthesis pathway, as well as a second molybdopterin synthase that was most likely acquired via horizontal gene transfer from cyanobacteria by a common ancestor of stramenopiles. EctoGEM represents an evolving community resource to gain deeper understanding of the biology of brown algae and the diversification of physiological processes. The integrative computational method applied for its reconstruction will be valuable to set up similar approaches for other organisms distant from biological benchmark models.     

A regulatory network modeled from wild-type gene expression data guides functional predictions in Caenorhabditis elegans development

Background

The iJO1366 reconstruction of the metabolic network of Escherichia coli is one of the most complete and accurate metabolic reconstructions available for any organism. Still, because our knowledge of even well-studied model organisms such as this one is incomplete, this network reconstruction contains gaps and possible errors. There are a total of 208 blocked metabolites in iJO1366, representing gaps in the network.

Results

A new model improvement workflow was developed to compare model based phenotypic predictions to experimental data to fill gaps and correct errors. A Keio Collection based dataset of E. coli gene essentiality was obtained from literature data and compared to model predictions. The SMILEY algorithm was then used to predict the most likely missing reactions in the reconstructed network, adding reactions from a KEGG based universal set of metabolic reactions. The feasibility of these putative reactions was determined by comparing updated versions of the model to the experimental dataset, and genes were predicted for the most feasible reactions.

Conclusions

Numerous improvements to the iJO1366 metabolic reconstruction were suggested by these analyses. Experiments were performed to verify several computational predictions, including a new mechanism for growth on myo-inositol. The other predictions made in this study should be experimentally verifiable by similar means. Validating all of the predictions made here represents a substantial but important undertaking.     

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.     

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.     

Fast Reconstruction of Compact Context-Specific Metabolic Network Models

Systemic approaches to the study of a biological cell or tissue rely increasingly on the use of context-specific metabolic network models. The reconstruction of such a model from high-throughput data can routinely involve large numbers of tests under different conditions and extensive parameter tuning, which calls for fast algorithms. We present fastcore, a generic algorithm for reconstructing context-specific metabolic network models from global genome-wide metabolic network models such as Recon X. fastcore takes as input a core set of reactions that are known to be active in the context of interest (e.g., cell or tissue), and it searches for a flux consistent subnetwork of the global network that contains all reactions from the core set and a minimal set of additional reactions. Our key observation is that a minimal consistent reconstruction can be defined via a set of sparse modes of the global network, and fastcore iteratively computes such a set via a series of linear programs. Experiments on liver data demonstrate speedups of several orders of magnitude, and significantly more compact reconstructions, over a rival method. Given its simplicity and its excellent performance, fastcore can form the backbone of many future metabolic network reconstruction algorithms.     

Metabolic networks are NP-hard to reconstruct

High-throughput data from various omics and sequencing techniques have rendered the automated metabolic network reconstruction a highly relevant problem. Our approach reflects the inherent probabilistic nature of the steps involved in metabolic network reconstruction. Here, the goal is to arrive at networks which combine probabilistic information with the possibility to obtain a small number of disconnected network constituents by reduction of a given preliminary probabilistic metabolic network. We define automated metabolic network reconstruction as an optimization problem on four-partite graph (nodes representing genes, enzymes, reactions, and metabolites) which integrates: (1) probabilistic information obtained from the existing process for metabolic reconstruction from a given genome, (2) connectedness of the raw metabolic network, and (3) clustering of components in the reconstructed metabolic network. The practical implications of our theoretical analysis refer to the quality of reconstructed metabolic networks and shed light on the problem of finding more efficient and effective methods for automated reconstruction. Our main contributions include: a completeness result for the defined problem, polynomial-time approximation algorithm, and an optimal polynomial-time algorithm for trees. Moreover, we exemplify our approach by the reconstruction of the sucrose biosynthesis pathway in Chlamydomonas reinhardtii.     

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.     

Genome-scale reconstruction and <Emphasis Type="Italic">in silico</Emphasis> analysis of the <Emphasis Type="Italic">Clostridium acetobutylicum</Emphasis> ATCC 824 metabolic network

To understand the metabolic characteristics of Clostridium acetobutylicum and to examine the potential for enhanced butanol production, we reconstructed the genome-scale metabolic network from its annotated genomic sequence and analyzed strategies to improve its butanol production. The generated reconstructed network consists of 502 reactions and 479 metabolites and was used as the basis for an in silico model that could compute metabolic and growth performance for comparison with fermentation data. The in silico model successfully predicted metabolic fluxes during the acidogenic phase using classical flux balance analysis. Nonlinear programming was used to predict metabolic fluxes during the solventogenic phase. In addition, essential genes were predicted via single gene deletion studies. This genome-scale in silico metabolic model of C. acetobutylicum should be useful for genome-wide metabolic analysis as well as strain development for improving production of biochemicals, including butanol. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. J. L. and H. Y. equally contributed to this work.     

Stoichiometric model and metabolic flux analysis for Leptospirillum ferrooxidans

A metabolic model for Leptospirillum ferrooxidans was developed based on the genomic information of an analogous iron oxidizing bacteria and on the pathways of ferrous iron oxidation, nitrogen and CO₂ assimilation based on experimental evidence for L. ferrooxidans found in the literature. From this metabolic reconstruction, a stoichiometric model was built, which includes 86 reactions describing the main catabolic and anabolic aspects of its metabolism. The model obtained has 2 degrees of freedom, so two external fluxes were estimated to achieve a determined and observable system. By using the external oxygen consumption rate and the generation flux biomass as input data, a metabolic flux map with a distribution of internal fluxes was obtained. The results obtained were verified with experimental data from the literature, achieving a very good prediction of the metabolic behavior of this bacterium at steady state. Biotechnol. Bioeng. 2010;107:696–706. © 2010 Wiley Periodicals, Inc.     

Metabolic network reconstruction and flux variability analysis of storage synthesis in developing oilseed rape (Brassica napus L.) embryos

Computational simulation of large‐scale biochemical networks can be used to analyze and predict the metabolic behavior of an organism, such as a developing seed. Based on the biochemical literature, pathways databases and decision rules defining reaction directionality we reconstructed bna572, a stoichiometric metabolic network model representing Brassica napus seed storage metabolism. In the highly compartmentalized network about 25% of the 572 reactions are transport reactions interconnecting nine subcellular compartments and the environment. According to known physiological capabilities of developing B. napus embryos, four nutritional conditions were defined to simulate heterotrophy or photoheterotrophy, each in combination with the availability of inorganic nitrogen (ammonia, nitrate) or amino acids as nitrogen sources. Based on mathematical linear optimization the optimal solution space was comprehensively explored by flux variability analysis, thereby identifying for each reaction the range of flux values allowable under optimality. The range and variability of flux values was then categorized into flux variability types. Across the four nutritional conditions, approximately 13% of the reactions have variable flux values and 10–11% are substitutable (can be inactive), both indicating metabolic redundancy given, for example, by isoenzymes, subcellular compartmentalization or the presence of alternative pathways. About one‐third of the reactions are never used and are associated with pathways that are suboptimal for storage synthesis. Fifty‐seven reactions change flux variability type among the different nutritional conditions, indicating their function in metabolic adjustments. This predictive modeling framework allows analysis and quantitative exploration of storage metabolism of a developing B. napus oilseed.     

Integrative genome‐scale metabolic analysis of Vibrio vulnificus for drug targeting and discovery

Although the genomes of many microbial pathogens have been studied to help identify effective drug targets and novel drugs, such efforts have not yet reached full fruition. In this study, we report a systems biological approach that efficiently utilizes genomic information for drug targeting and discovery, and apply this approach to the opportunistic pathogen Vibrio vulnificus CMCP6. First, we partially re‐sequenced and fully re‐annotated the V. vulnificus CMCP6 genome, and accordingly reconstructed its genome‐scale metabolic network, VvuMBEL943. The validated network model was employed to systematically predict drug targets using the concept of metabolite essentiality, along with additional filtering criteria. Target genes encoding enzymes that interact with the five essential metabolites finally selected were experimentally validated. These five essential metabolites are critical to the survival of the cell, and hence were used to guide the cost‐effective selection of chemical analogs, which were then screened for antimicrobial activity in a whole‐cell assay. This approach is expected to help fill the existing gap between genomics and drug discovery.     

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.     

冷冻电镜单颗粒三维重构密度掩模的自动生成方法研究

冷冻电镜单颗粒三维重构技术是用来解析生物大分子三维结构的常用方法.然而目前在单颗粒三维重构过程中,溶剂平滑操作还存在一定缺陷:没有一款主流的单颗粒三维重构程序能够自动寻找掩模(mask)三维密度图,使得三维重构过程难免受到噪音统计学模型计算偏差的干扰.为解决这一问题,本研究借鉴X射线晶体学中解析优化相位所广泛采用的溶剂平滑方法,采用高斯滤波、坎尼边缘检测、最小误差阈值处理等方法处理重构所得三维密度图,优化溶剂平滑操作,发展在单颗粒三维重构过程中自动寻找mask三维密度图的方法.运用三维密度图傅里叶壳层相关系数(fourier shell correlation,FSC)曲线图、模拟颗粒数据重构角度误差散点图等指标评估此方法的效果.结果表明,自动寻找mask密度图的方法能够较好地找到涵盖分子结构信号区域的mask密度图,较为明显提高三维重构所得密度图分辨率.     

Dissecting Leishmania infantum Energy Metabolism - A Systems Perspective

Leishmania infantum, causative agent of visceral leishmaniasis in humans, illustrates a complex lifecycle pertaining to two extreme environments, namely, the gut of the sandfly vector and human macrophages. Leishmania is capable of dynamically adapting and tactically switching between these critically hostile situations. The possible metabolic routes ventured by the parasite to achieve this exceptional adaptation to its varying environments are still poorly understood. In this study, we present an extensively reconstructed energy metabolism network of Leishmania infantum as an attempt to identify certain strategic metabolic routes preferred by the parasite to optimize its survival in such dynamic environments. The reconstructed network consists of 142 genes encoding for enzymes performing 237 reactions distributed across five distinct model compartments. We annotated the subcellular locations of different enzymes and their reactions on the basis of strong literature evidence and sequence-based detection of cellular localization signal within a protein sequence. To explore the diverse features of parasite metabolism the metabolic network was implemented and analyzed as a constraint-based model. Using a systems-based approach, we also put forth an extensive set of lethal reaction knockouts; some of which were validated using published data on Leishmania species. Performing a robustness analysis, the model was rigorously validated and tested for the secretion of overflow metabolites specific to Leishmania under varying extracellular oxygen uptake rate. Further, the fate of important non-essential amino acids in L. infantum metabolism was investigated. Stage-specific scenarios of L. infantum energy metabolism were incorporated in the model and key metabolic differences were outlined. Analysis of the model revealed the essentiality of glucose uptake, succinate fermentation, glutamate biosynthesis and an active TCA cycle as driving forces for parasite energy metabolism and its optimal growth. Finally, through our in silico knockout analysis, we could identify possible therapeutic targets that provide experimentally testable hypotheses.