Improving the iMM904 S. cerevisiae metabolic model using essentiality and synthetic lethality data

BackgroundSaccharomyces cerevisiae is the first eukaryotic organism for which a multi-compartment genome-scale metabolic model was constructed. Since then a sequence of improved metabolic reconstructions for yeast has been introduced. These metabolic models have been extensively used to elucidate the organizational principles of yeast metabolism and drive yeast strain engineering strategies for targeted overproductions. They have also served as a starting point and a benchmark for the reconstruction of genome-scale metabolic models for other eukaryotic organisms. In spite of the successive improvements in the details of the described metabolic processes, even the recent yeast model (i.e., i MM904) remains significantly less predictive than the latest E. coli model (i.e., i AF1260). This is manifested by its significantly lower specificity in predicting the outcome of grow/no grow experiments in comparison to the E. coli model.ResultsIn this paper we make use of the automated GrowMatch procedure for restoring consistency with single gene deletion experiments in yeast and extend the procedure to make use of synthetic lethality data using the genome-scale model i MM904 as a basis. We identified and vetted using literature sources 120 distinct model modifications including various regulatory constraints for minimal and YP media. The incorporation of the suggested modifications led to a substantial increase in the fraction of correctly predicted lethal knockouts (i.e., specificity) from 38.84% (87 out of 224) to 53.57% (120 out of 224) for the minimal medium and from 24.73% (45 out of 182) to 40.11% (73 out of 182) for the YP medium. Synthetic lethality predictions improved from 12.03% (16 out of 133) to 23.31% (31 out of 133) for the minimal medium and from 6.96% (8 out of 115) to 13.04% (15 out of 115) for the YP medium.ConclusionsOverall, this study provides a roadmap for the computationally driven correction of multi-compartment genome-scale metabolic models and demonstrates the value of synthetic lethals as curation agents.

[1]  S C HARTMAN,et al.  The biosynthesis of the purines. , 1951, Ergebnisse der Physiologie, biologischen Chemie und experimentellen Pharmakologie.

[2]  M. McCammon,et al.  Genetic and biochemical interactions involving tricarboxylic acid cycle (TCA) function using a collection of mutants defective in all TCA cycle genes. , 1999, Genetics.

[3]  Jeffrey D. Orth,et al.  Systematizing the generation of missing metabolic knowledge , 2010, Biotechnology and bioengineering.

[4]  J. Bader,et al.  A robust toolkit for functional profiling of the yeast genome. , 2004, Molecular cell.

[5]  Adam M. Feist,et al.  A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information , 2007, Molecular systems biology.

[6]  B. Palsson,et al.  A protocol for generating a high-quality genome-scale metabolic reconstruction , 2010 .

[7]  André Goffeau,et al.  The yeast genome directory. , 1997, Nature.

[8]  M. Basson,et al.  Identifying mutations in duplicated functions in Saccharomyces cerevisiae: recessive mutations in HMG-CoA reductase genes. , 1987, Genetics.

[9]  Gary D Bader,et al.  Systematic Genetic Analysis with Ordered Arrays of Yeast Deletion Mutants , 2001, Science.

[10]  S. Henry,et al.  Saccharomyces cerevisiae cho2 mutants are deficient in phospholipid methylation and cross-pathway regulation of inositol synthesis. , 1988, Genetics.

[11]  B. Palsson,et al.  Regulation of gene expression in flux balance models of metabolism. , 2001, Journal of theoretical biology.

[12]  Ronald W. Davis,et al.  Functional profiling of the Saccharomyces cerevisiae genome , 2002, Nature.

[13]  T. Kodaki,et al.  Yeast phosphatidylethanolamine methylation pathway. Cloning and characterization of two distinct methyltransferase genes. , 1987, The Journal of biological chemistry.

[14]  Intawat Nookaew,et al.  The genome-scale metabolic model iIN800 of Saccharomyces cerevisiae and its validation: a scaffold to query lipid metabolism , 2008, BMC Syst. Biol..

[15]  B. Kaufman,et al.  Aconitase Couples Metabolic Regulation to Mitochondrial DNA Maintenance , 2005, Science.

[16]  Adilson E. Motter,et al.  Spontaneous Reaction Silencing in Metabolic Optimization , 2008, PLoS Comput. Biol..

[17]  Zhikang Yin,et al.  Multiple signalling pathways trigger the exquisite sensitivity of yeast gluconeogenic mRNAs to glucose , 1996, Molecular microbiology.

[18]  Yoshihiro Yamanishi,et al.  KEGG for linking genomes to life and the environment , 2007, Nucleic Acids Res..

[19]  B C Guha Investigations on vitamin B(2): The sources of vitamin B(2). II. The stability of vitamin B(2). III. The chemistry of vitamin B(2). , The Biochemical journal.

[20]  B. Palsson,et al.  Expanded Metabolic Reconstruction of Helicobacter pylori (iIT341 GSM/GPR): an In Silico Genome-Scale Characterization of Single- and Double-Deletion Mutants , 2005, Journal of bacteriology.

[21]  J. Boeke,et al.  DNA helicase gene interaction network defined using synthetic lethality analyzed by microarray , 2003, Nature Genetics.

[22]  B. Palsson,et al.  Genome-scale reconstruction of the Saccharomyces cerevisiae metabolic network. , 2003, Genome research.

[23]  A. T. Carter,et al.  Genetic Analysis and Enzyme Activity Suggest the Existence of More Than One Minimal Functional Unit Capable of Synthesizing Phosphoribosyl Pyrophosphate in Saccharomyces cerevisiae * , 1999, The Journal of Biological Chemistry.

[24]  Steven C Almo,et al.  Closed Site Complexes of Adenine Phosphoribosyltransferase fromGiardia lamblia Reveal a Mechanism of Ribosyl Migration* , 2002, The Journal of Biological Chemistry.

[25]  Bires Chandra Guha,et al.  Investigations on vitamin B2: The sources of vitamin B2. II. The stability of vitamin B2. III. The chemistry of vitamin B2 , 1931 .

[26]  Markus J. Herrgård,et al.  A consensus yeast metabolic network reconstruction obtained from a community approach to systems biology , 2008, Nature Biotechnology.

[27]  Rick L Stevens,et al.  iBsu1103: a new genome-scale metabolic model of Bacillus subtilis based on SEED annotations , 2009, Genome Biology.

[28]  S Haarasilta,et al.  On the activity and regulation of anaplerotic and gluconeogenetic enzymes during the growth process of baker's yeast. The biphasic growth. , 1975, European journal of biochemistry.

[29]  Markus J. Herrgård,et al.  Integrated analysis of regulatory and metabolic networks reveals novel regulatory mechanisms in Saccharomyces cerevisiae. , 2006, Genome research.

[30]  Bernard Henrissat,et al.  Localization of Synthesis of β1,6-Glucan inSaccharomyces cerevisiae , 1999, Journal of bacteriology.

[31]  Manal AbuOun,et al.  Genome Scale Reconstruction of a Salmonella Metabolic Model , 2009, The Journal of Biological Chemistry.

[32]  S. Schuster,et al.  Structural robustness of metabolic networks with respect to multiple knockouts. , 2008, Journal of theoretical biology.

[33]  B. Palsson,et al.  Escherichia coli K-12 undergoes adaptive evolution to achieve in silico predicted optimal growth , 2002, Nature.

[34]  Vinay Satish Kumar,et al.  Optimization based automated curation of metabolic reconstructions , 2007, BMC Bioinformatics.

[35]  S. Steavenson,et al.  The pentose phosphate pathway in the endoplasmic reticulum. , 1988, The Journal of biological chemistry.

[36]  Neil Swainston,et al.  Further developments towards a genome-scale metabolic model of yeast , 2010, BMC Systems Biology.

[37]  Isaac Meilijson,et al.  Can single knockouts accurately single out gene functions? , 2008, BMC Systems Biology.

[38]  E. O’Shea,et al.  Global analysis of protein localization in budding yeast , 2003, Nature.

[39]  F. Klis,et al.  The linkage of (1–3)‐β‐glucan to chitin during cell wall assembly in Saccharomyces cerevisiae , 1994, Yeast.

[40]  Matthew D. Jankowski,et al.  Genome-scale thermodynamic analysis of Escherichia coli metabolism. , 2006, Biophysical journal.

[41]  Ian G. Macreadie,et al.  Dihydropteroate Synthase Mutations in Pneumocystis jiroveci Can Affect Sulfamethoxazole Resistance in a Saccharomyces cerevisiae Model , 2004, Antimicrobial Agents and Chemotherapy.

[42]  Jason A. Papin,et al.  Genome-Scale Reconstruction and Analysis of the Pseudomonas putida KT2440 Metabolic Network Facilitates Applications in Biotechnology , 2008, PLoS Comput. Biol..

[43]  G W Smithers,et al.  Hypoxanthine phosphoribosyltransferase: radiochemical assay procedures for the forward and reverse reactions. , 1985, Analytical biochemistry.

[44]  J M BUCHANAN,et al.  Biosynthesis of the purines. XVI. The synthesis of adenosine 5'-phosphate and 5-amino-4-imidazolecarboxamide ribotide by a nucleotide pyrophosphorylase. , 1957, The Journal of biological chemistry.

[45]  A. Barabasi,et al.  Predicting synthetic rescues in metabolic networks , 2008, Molecular systems biology.

[46]  Andrew R. Joyce,et al.  Experimental and Computational Assessment of Conditionally Essential Genes in Escherichia coli , 2006, Journal of bacteriology.

[47]  Sean R. Collins,et al.  A tool-kit for high-throughput, quantitative analyses of genetic interactions in E. coli , 2008, Nature Methods.

[48]  S. Henry,et al.  The phospholipid methyltransferases in yeast. , 1997, Biochimica et biophysica acta.

[49]  Kara Dolinski,et al.  Fungal BLAST and Model Organism BLASTP Best Hits: new comparison resources at the Saccharomyces Genome Database (SGD) , 2004, Nucleic Acids Res..

[50]  Markus J. Herrgård,et al.  Reconstruction and validation of Saccharomyces cerevisiae iND750, a fully compartmentalized genome-scale metabolic model. , 2004, Genome research.

[51]  Michael R. Green,et al.  Dissecting the Regulatory Circuitry of a Eukaryotic Genome , 1998, Cell.

[52]  Anurag Tiwari,et al.  Molecular characterization reveals that YMR278w encoded protein is environmental stress response homologue of Saccharomyces cerevisiae PGM2. , 2008, Biochemical and biophysical research communications.

[53]  Vinay Satish Kumar,et al.  A Genome-Scale Metabolic Reconstruction of Mycoplasma genitalium, iPS189 , 2009, PLoS Comput. Biol..

[54]  Albert Sickmann,et al.  The proteome of Saccharomyces cerevisiae mitochondria , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[55]  J. François,et al.  A new method for quantitative determination of polysaccharides in the yeast cell wall. Application to the cell wall defective mutants of Saccharomyces cerevisiae , 1998, Yeast.

[56]  Vinay Satish Kumar,et al.  GrowMatch: An Automated Method for Reconciling In Silico/In Vivo Growth Predictions , 2009, PLoS Comput. Biol..

[57]  C. Thomas Caskey,et al.  The HPRT locus , 1979, Cell.

[58]  B. Palsson,et al.  Systems approach to refining genome annotation , 2006, Proceedings of the National Academy of Sciences.

[59]  T. Kodaki,et al.  Characterization of the methyltransferases in the yeast phosphatidylethanolamine methylation pathway by selective gene disruption. , 1989, European journal of biochemistry.

[60]  S. Oliver,et al.  Plasticity of genetic interactions in metabolic networks of yeast , 2007, Proceedings of the National Academy of Sciences.

[61]  Ali R. Zomorrodi,et al.  Genome-scale gene/reaction essentiality and synthetic lethality analysis , 2009, Molecular systems biology.

[62]  Lee Irvin Smith,et al.  The chemistry of vitamin E. , 1940 .

[63]  Bernhard O. Palsson,et al.  Connecting Extracellular Metabolomic Measurements to Intracellular Flux States in Yeast , 2022 .

[64]  C. Hollenberg,et al.  Concurrent knock‐out of at least 20 transporter genes is required to block uptake of hexoses in Saccharomyces cerevisiae , 1999, FEBS letters.

[65]  Jennifer L. Reed,et al.  OptORF: Optimal metabolic and regulatory perturbations for metabolic engineering of microbial strains , 2010, BMC Systems Biology.

[66]  B. Palsson,et al.  Large-scale evaluation of in silico gene deletions in Saccharomyces cerevisiae. , 2003, Omics : a journal of integrative biology.

[67]  U. Sauer,et al.  Metabolic functions of duplicate genes in Saccharomyces cerevisiae. , 2005, Genome research.