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.     

Network Reconstruction Using Nonparametric Additive ODE Models

Network representations of biological systems are widespread and reconstructing unknown networks from data is a focal problem for computational biologists. For example, the series of biochemical reactions in a metabolic pathway can be represented as a network, with nodes corresponding to metabolites and edges linking reactants to products. In a different context, regulatory relationships among genes are commonly represented as directed networks with edges pointing from influential genes to their targets. Reconstructing such networks from data is a challenging problem receiving much attention in the literature. There is a particular need for approaches tailored to time-series data and not reliant on direct intervention experiments, as the former are often more readily available. In this paper, we introduce an approach to reconstructing directed networks based on dynamic systems models. Our approach generalizes commonly used ODE models based on linear or nonlinear dynamics by extending the functional class for the functions involved from parametric to nonparametric models. Concomitantly we limit the complexity by imposing an additive structure on the estimated slope functions. Thus the submodel associated with each node is a sum of univariate functions. These univariate component functions form the basis for a novel coupling metric that we define in order to quantify the strength of proposed relationships and hence rank potential edges. We show the utility of the method by reconstructing networks using simulated data from computational models for the glycolytic pathway of Lactocaccus Lactis and a gene network regulating the pluripotency of mouse embryonic stem cells. For purposes of comparison, we also assess reconstruction performance using gene networks from the DREAM challenges. We compare our method to those that similarly rely on dynamic systems models and use the results to attempt to disentangle the distinct roles of linearity, sparsity, and derivative estimation.     

A Systems Approach to Predict Oncometabolites via Context-Specific Genome-Scale Metabolic Networks

Altered metabolism in cancer cells has been viewed as a passive response required for a malignant transformation. However, this view has changed through the recently described metabolic oncogenic factors: mutated isocitrate dehydrogenases (IDH), succinate dehydrogenase (SDH), and fumarate hydratase (FH) that produce oncometabolites that competitively inhibit epigenetic regulation. In this study, we demonstrate in silico predictions of oncometabolites that have the potential to dysregulate epigenetic controls in nine types of cancer by incorporating massive scale genetic mutation information (collected from more than 1,700 cancer genomes), expression profiling data, and deploying Recon 2 to reconstruct context-specific genome-scale metabolic models. Our analysis predicted 15 compounds and 24 substructures of potential oncometabolites that could result from the loss-of-function and gain-of-function mutations of metabolic enzymes, respectively. These results suggest a substantial potential for discovering unidentified oncometabolites in various forms of cancers.     

Network Thermodynamic Curation of Human and Yeast Genome-Scale Metabolic Models

Genome-scale models are used for an ever-widening range of applications. Although there has been much focus on specifying the stoichiometric matrix, the predictive power of genome-scale models equally depends on reaction directions. Two-thirds of reactions in the two eukaryotic reconstructions Homo sapiens Recon 1 and Yeast 5 are specified as irreversible. However, these specifications are mainly based on biochemical textbooks or on their similarity to other organisms and are rarely underpinned by detailed thermodynamic analysis. In this study, a to our knowledge new workflow combining network-embedded thermodynamic and flux variability analysis was used to evaluate existing irreversibility constraints in Recon 1 and Yeast 5 and to identify new ones. A total of 27 and 16 new irreversible reactions were identified in Recon 1 and Yeast 5, respectively, whereas only four reactions were found with directions incorrectly specified against thermodynamics (three in Yeast 5 and one in Recon 1). The workflow further identified for both models several isolated internal loops that require further curation. The framework also highlighted the need for substrate channeling (in human) and ATP hydrolysis (in yeast) for the essential reaction catalyzed by phosphoribosylaminoimidazole carboxylase in purine metabolism. Finally, the framework highlighted differences in proline metabolism between yeast (cytosolic anabolism and mitochondrial catabolism) and humans (exclusively mitochondrial metabolism). We conclude that network-embedded thermodynamics facilitates the specification and validation of irreversibility constraints in compartmentalized metabolic models, at the same time providing further insight into network properties.     

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.     

Context-Specific Mechanisms of Cell Polarity Regulation

Cell polarity is an essential process shared by almost all animal tissues. Moreover, cell polarity enables cells to sense and respond to the cues provided by the neighboring cells and the surrounding microenvironment. These responses play a critical role in regulating key physiological processes, including cell migration, proliferation, differentiation, vesicle trafficking and immune responses. The polarity protein complexes regulating these interactions are highly evolutionarily conserved between vertebrates and invertebrates. Interestingly, these polarity complexes interact with each other and key signaling pathways in a cell-polarity context-dependent manner. However, the exact mechanisms by which these interactions take place are poorly understood. In this review, we will focus on the roles of the key polarity complexes SCRIB, PAR and Crumbs in regulating different forms of cell polarity, including epithelial cell polarity, cell migration, asymmetric cell division and the T-cell immunological synapse assembly and signaling.     

From Enzyme Kinetics to Metabolic Network Modeling – Visualization Tool for Enhanced Kinetic Analysis of Biochemical Network Models

Model‐based analysis of enzyme kinetics allows the determination of optimal conditions for their use in biocatalysis. For biotransformations or fermentative approaches the modeling of metabolic pathways or complex metabolic networks is necessary to obtain model‐based predictions of steps which limit product formation within the network. To set up adequate kinetic models, relevant mechanistic information about enzyme properties is required and can be taken from in vitro studies with isolated enzymes or from in vivo investigations using stimulus‐response experiments which provide a lot of kinetic information about the metabolic network. But with increasing number of reaction steps and regulatory interdependencies in the network structure the amount of simulation data dramatically increases and the simulation results from the dynamic models become difficult to analyze and interpret. Demonstrated for an Escherichia coli model of the central carbon metabolism, methods for visualization and animation of simulation data were applied and extended to facilitate model analysis and biological interpretation. The dynamic metabolite pool and metabolic flux changes were visualized simultaneously by a software tool. In addition, a new quantification method for enzyme activation/inhibition was proposed, and this information was implemented in the metabolic visualization.     

Global Metabolic Reconstruction and Metabolic Gene Evolution in the Cattle Genome

The sequence of cattle genome provided a valuable opportunity to systematically link genetic and metabolic traits of cattle. The objectives of this study were 1) to reconstruct genome-scale cattle-specific metabolic pathways based on the most recent and updated cattle genome build and 2) to identify duplicated metabolic genes in the cattle genome for better understanding of metabolic adaptations in cattle. A bioinformatic pipeline of an organism for amalgamating genomic annotations from multiple sources was updated. Using this, an amalgamated cattle genome database based on UMD_3.1, was created. The amalgamated cattle genome database is composed of a total of 33,292 genes: 19,123 consensus genes between NCBI and Ensembl databases, 8,410 and 5,493 genes only found in NCBI or Ensembl, respectively, and 266 genes from NCBI scaffolds. A metabolic reconstruction of the cattle genome and cattle pathway genome database (PGDB) was also developed using Pathway Tools, followed by an intensive manual curation. The manual curation filled or revised 68 pathway holes, deleted 36 metabolic pathways, and added 23 metabolic pathways. Consequently, the curated cattle PGDB contains 304 metabolic pathways, 2,460 reactions including 2,371 enzymatic reactions, and 4,012 enzymes. Furthermore, this study identified eight duplicated genes in 12 metabolic pathways in the cattle genome compared to human and mouse. Some of these duplicated genes are related with specific hormone biosynthesis and detoxifications. The updated genome-scale metabolic reconstruction is a useful tool for understanding biology and metabolic characteristics in cattle. There has been significant improvements in the quality of cattle genome annotations and the MetaCyc database. The duplicated metabolic genes in the cattle genome compared to human and mouse implies evolutionary changes in the cattle genome and provides a useful information for further research on understanding metabolic adaptations of cattle.     

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.     

乙醇合成途径的计算机重构

目的:用计算机重构乙醇合成途径,为合成生物燃料乙醇提供理论依据。方法:利用KEGG反应、化合物数据提取反应等式,过滤掉42个通用代谢物参与的反应,然后利用剩下的反应构建反应矩阵;利用广度优先搜索算法在反应矩阵中搜索生成乙醇的代谢途径。结果:计算机重构了23 108条乙醇合成途径,以大肠杆菌作为产乙醇基因工程菌为例,通过限制改构菌整合的关键酶数目,分别得到了78条以酒精O-乙酰基转移酶为关键酶的乙醇合成通路和89条以丙酮酸脱羧酶和乙醇脱氢酶为关键酶的乙醇合成通路,并构建了相应的乙醇合成网络图,标注每个反应的酶及编码该酶的基因。结论:通过计算机方法重构了多种乙醇合成途径,可以为利用微生物工业化生产乙醇提供理论依据。     

代谢网络刚性与代谢工程

本文介绍了有关代谢网络,网络刚性和代谢工程的基本观点,特别是关于网络刚性和载流途径的工作假设,还讨论了代谢工程的应用与面临的难题。     

EGFR Signal-Network Reconstruction Demonstrates Metabolic Crosstalk in EMT

Epithelial to mesenchymal transition (EMT) is an important event during development and cancer metastasis. There is limited understanding of the metabolic alterations that give rise to and take place during EMT. Dysregulation of signalling pathways that impact metabolism, including epidermal growth factor receptor (EGFR), are however a hallmark of EMT and metastasis. In this study, we report the investigation into EGFR signalling and metabolic crosstalk of EMT through constraint-based modelling and analysis of the breast epithelial EMT cell model D492 and its mesenchymal counterpart D492M. We built an EGFR signalling network for EMT based on stoichiometric coefficients and constrained the network with gene expression data to build epithelial (EGFR_E) and mesenchymal (EGFR_M) networks. Metabolic alterations arising from differential expression of EGFR genes was derived from a literature review of AKT regulated metabolic genes. Signaling flux differences between EGFR_E and EGFR_M models subsequently allowed metabolism in D492 and D492M cells to be assessed. Higher flux within AKT pathway in the D492 cells compared to D492M suggested higher glycolytic activity in D492 that we confirmed experimentally through measurements of glucose uptake and lactate secretion rates. The signaling genes from the AKT, RAS/MAPK and CaM pathways were predicted to revert D492M to D492 phenotype. Follow-up analysis of EGFR signaling metabolic crosstalk in three additional breast epithelial cell lines highlighted variability in in vitro cell models of EMT. This study shows that the metabolic phenotype may be predicted by in silico analyses of gene expression data of EGFR signaling genes, but this phenomenon is cell-specific and does not follow a simple trend.     

Models of Quantitative Variation of Flux in Metabolic Pathways

As a model of variation in a quantitative character, enzyme activity variation segregating in a population is assumed to affect the flux in simple metabolic pathways. The genetic variation of flux is partitioned into additive and nonadditive components. An interaction component of flux variance is present because the effect of an allelic substitution is modified by other substitutions which change the concentrations of shared metabolites. In a haploid population, the the proportion of interaction variance is a function of the gene frequencies at the loci contributing to the flux variation, enzyme activities of mutant and wild type at variable loci and activities at nonvariable loci. The proportion of interaction variance is inversely related to the ratio of mutant to wild-type activities at the loci controlling the enzyme activities. The interaction component as a function of gene frequencies is at a maximum with high mutant allele frequencies. In contrast, the dominance component which would apply to a diploid population is maximal as a proportion of the total when mutant alleles are at low frequencies. Unless there are many loci with large differences in activity between the alleles, the interaction component is a small proportion of the total variance. Data on enzyme activity variation from natural and artificial populations suggest that such variation generates little nonadditive variance despite the highly interactive nature of the underlying biochemical system.     

Semantics of Multimodal Network Models

A multimodal network (MMN) is a novel graph-theoretic formalism designed to capture the structure of biological networks and to represent relationships derived from multiple biological databases. MMNs generalize the standard notions of graphs and hypergraphs, which are the bases of current diagrammatic representations of biological phenomena, and incorporate the concept of mode. Each vertex of an MMN is a biological entity, a biot, while each modal hyperedge is a typed relationship, where the type is given by the mode of the hyperedge. The semantics of each modal hyperedge e is given through denotational semantics, where a valuation function f_{e} defines the relationship among the values of the vertices incident on e. The meaning of an MMN is denoted in terms of the semantics of a hyperedge sequence. A companion paper defines MMNs and concentrates on the structural aspects of MMNs. This paper develops MMN denotational semantics when used as a representation of the semantics of biological networks and discusses applications of MMNs in managing complex biological data.     

Systematic Construction of Kinetic Models from Genome-Scale Metabolic Networks

The quantitative effects of environmental and genetic perturbations on metabolism can be studied in silico using kinetic models. We present a strategy for large-scale model construction based on a logical layering of data such as reaction fluxes, metabolite concentrations, and kinetic constants. The resulting models contain realistic standard rate laws and plausible parameters, adhere to the laws of thermodynamics, and reproduce a predefined steady state. These features have not been simultaneously achieved by previous workflows. We demonstrate the advantages and limitations of the workflow by translating the yeast consensus metabolic network into a kinetic model. Despite crudely selected data, the model shows realistic control behaviour, a stable dynamic, and realistic response to perturbations in extracellular glucose concentrations. The paper concludes by outlining how new data can continuously be fed into the workflow and how iterative model building can assist in directing experiments.     

基于复杂网络理论的代谢网络结构研究进展

后基因组时代研究的一个重要内容就是在系统生物学的基础上对多种分子和基因相互作用网络进行分析,理解生物系统是如何从单个构造模块的基础上组织起来的。近几年,复杂网络理论得到了迅速发展,其理论方法在生物网络的研究中得到了广泛应用。一些学者运用该理论方法研究了大量有关代谢网络的结构组成以及网络中子集团的层次组成关系,并获得了一些有意义的结果。这些结果对生物功能的预测具有一定的指导作用。对近几年来有关这方面的进展作简要综述。