Predicting folding free energy changes upon single point mutations

MOTIVATION The folding free energy is an important characteristic of proteins stability and is directly related to protein's wild-type function. The changes of protein's stability due to naturally occurring mutations, missense mutations, are typically causing diseases. Single point mutations made in vitro are frequently used to assess the contribution of given amino acid to the stability of the protein. In both cases, it is desirable to predict the change of the folding free energy upon single point mutations in order to either provide insights of the molecular mechanism of the change or to design new experimental studies. RESULTS We report an approach that predicts the free energy change upon single point mutation by utilizing the 3D structure of the wild-type protein. It is based on variation of the molecular mechanics Generalized Born (MMGB) method, scaled with optimized parameters (sMMGB) and utilizing specific model of unfolded state. The corresponding mutations are built in silico and the predictions are tested against large dataset of 1109 mutations with experimentally measured changes of the folding free energy. Benchmarking resulted in root mean square deviation = 1.78 kcal/mol and slope of the linear regression fit between the experimental data and the calculations was 1.04. The sMMGB is compared with other leading methods of predicting folding free energy changes upon single mutations and results discussed with respect to various parameters. AVAILABILITY All the pdb files we used in this article can be downloaded from http://compbio.clemson.edu/downloadDir/mentaldisorders/sMMGB_pdb.rar. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

[1]  Molly Strom,et al.  Single Nucleotide Polymorphisms That Cause Structural Changes in the Cyclic AMP Receptor Protein Transcriptional Regulator of the Tuberculosis Vaccine Strain Mycobacterium bovis BCG Alter Global Gene Expression without Attenuating Growth , 2008, Infection and Immunity.

[2]  B. Matthews,et al.  Response of a protein structure to cavity-creating mutations and its relation to the hydrophobic effect. , 1992, Science.

[3]  E. Alexov,et al.  Approaches and resources for prediction of the effects of non-synonymous single nucleotide polymorphism on protein function and interactions. , 2008, Current pharmaceutical biotechnology.

[4]  Huan-Xiang Zhou,et al.  Direct test of the Gaussian-chain model for treating residual charge-charge interactions in the unfolded state of proteins. , 2003, Journal of the American Chemical Society.

[5]  T L Blundell,et al.  Prediction of the stability of protein mutants based on structural environment-dependent amino acid substitution and propensity tables. , 1997, Protein engineering.

[6]  Philip E. Bourne,et al.  The RCSB PDB information portal for structural genomics , 2005, Nucleic Acids Res..

[7]  Adam Godzik,et al.  Modeling and Analyzing Three-Dimensional Structures of Human Disease Proteins , 2005, Pacific Symposium on Biocomputing.

[8]  B. Shastry.  SNPs: impact on gene function and phenotype. , 2009, Methods in molecular biology.

[9]  Shuangye Yin,et al.  Eris: an automated estimator of protein stability , 2007, Nature Methods.

[10]  J. Moult,et al.  Identification and analysis of deleterious human SNPs. , 2006, Journal of molecular biology.

[11]  Pedro A Fernandes,et al.  Hot spots—A review of the protein–protein interface determinant amino‐acid residues , 2007, Proteins.

[12]  K Nishikawa,et al.  Knowledge-based potential defined for a rotamer library to design protein sequences. , 2001, Protein engineering.

[13]  C. Frenz,et al.  Neural network‐based prediction of mutation‐induced protein stability changes in Staphylococcal nuclease at 20 residue positions , 2005, Proteins.

[14]  Villegas,et al.  Stabilization of proteins by rational design of alpha-helix stability using helix/coil transition theory. , 1995, Folding & design.

[15]  Gennady Verkhivker,et al.  Computational modeling of structurally conserved cancer mutations in the RET and MET kinases: the impact on protein structure, dynamics, and stability. , 2009, Biophysical journal.

[16]  Z. Xiang,et al.  Extending the accuracy limits of prediction for side-chain conformations. , 2001, Journal of molecular biology.

[17]  B. L. de Groot,et al.  Predicting free energy changes using structural ensembles. , 2009, Nature methods.

[18]  Nikolay V Dokholyan,et al.  FALS mutations in Cu, Zn superoxide dismutase destabilize the dimer and increase dimer dissociation propensity: A large-scale thermodynamic analysis , 2006, Amyloid : the international journal of experimental and clinical investigation : the official journal of the International Society of Amyloidosis.

[19]  D. Baker,et al.  Native protein sequences are close to optimal for their structures. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[20]  Chi-Chuan Tseng,et al.  Gut-enriched Krüppel-like Factor Represses Ornithine Decarboxylase Gene Expression and Functions as Checkpoint Regulator in Colonic Cancer Cells* , 2002, The Journal of Biological Chemistry.

[21]  M Karplus,et al.  Simulation analysis of the stability mutant R96H of T4 lysozyme. , 1991, Biochemistry.

[22]  A. Ben-Naim.  STATISTICAL POTENTIALS EXTRACTED FROM PROTEIN STRUCTURES : ARE THESE MEANINGFUL POTENTIALS? , 1997 .

[23]  A. Fersht,et al.  Crystal structural analysis of mutations in the hydrophobic cores of barnase. , 1994, Journal of molecular biology.

[24]  Piero Fariselli,et al.  I-Mutant2.0: predicting stability changes upon mutation from the protein sequence or structure , 2005, Nucleic Acids Res..

[25]  B. Matthews,et al.  The response of T4 lysozyme to large‐to‐small substitutions within the core and its relation to the hydrophobic effect , 1998, Protein science : a publication of the Protein Society.

[26]  Xiong Xing-Min.  Study of isochronal annealing behavior of neutron irradiated hydrogen FZ silicon by positron annihilation , 1986 .

[27]  Emil Alexov,et al.  Numerical calculations of the pH of maximal protein stability. The effect of the sequence composition and three-dimensional structure. , 2003, European journal of biochemistry.

[28]  Dietmar Schomburg,et al.  Prediction of protein thermostability with a direction‐ and distance‐dependent knowledge‐based potential , 2005, Protein science : a publication of the Protein Society.

[29]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[30]  Emil Alexov,et al.  Modeling effects of human single nucleotide polymorphisms on protein-protein interactions. , 2009, Biophysical journal.

[31]  D Gilis,et al.  PoPMuSiC, an algorithm for predicting protein mutant stability changes: application to prion proteins. , 2000, Protein engineering.

[32]  D Gilis,et al.  Predicting protein stability changes upon mutation using database-derived potentials: solvent accessibility determines the importance of local versus non-local interactions along the sequence. , 1997, Journal of molecular biology.

[33]  Feng Ding,et al.  Correction: Emergence of Protein Fold Families through Rational Design , 2006, PLoS Comput. Biol..

[34]  Tingjun Hou,et al.  Assessing the performance of the molecular mechanics/Poisson Boltzmann surface area and molecular mechanics/generalized Born surface area methods. II. The accuracy of ranking poses generated from docking , 2011, J. Comput. Chem..

[35]  R. Huber,et al.  Crystal structure analysis of oxidized Pseudomonas aeruginosa azurin at pH 5.5 and pH 9.0. A pH-induced conformational transition involves a peptide bond flip. , 1991, Journal of molecular biology.

[36]  M. Levitt,et al.  Accurate prediction of the stability and activity effects of site-directed mutagenesis on a protein core , 1991, Nature.

[37]  Petras J Kundrotas,et al.  Modeling of denatured state for calculation of the electrostatic contribution to protein stability , 2002, Protein science : a publication of the Protein Society.

[38]  Feng Ding,et al.  Modeling backbone flexibility improves protein stability estimation. , 2007, Structure.

[39]  P. Kollman,et al.  Pathways to a protein folding intermediate observed in a 1-microsecond simulation in aqueous solution. , 1998, Science.

[40]  P. A. Bash,et al.  Free energy calculations by computer simulation. , 1987, Science.

[41]  J A McCammon,et al.  Electrostatic contributions to the stability of halophilic proteins. , 1998, Journal of molecular biology.

[42]  Huan-Xiang Zhou,et al.  A Gaussian-chain model for treating residual charge–charge interactions in the unfolded state of proteins , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[43]  P. Bork,et al.  Human non-synonymous SNPs: server and survey. , 2002, Nucleic acids research.

[44]  James W Tunnell,et al.  Temperature-induced unfolding of epidermal growth factor (EGF): insight from molecular dynamics simulation. , 2010, Journal of molecular graphics & modelling.

[45]  Akinori Sarai,et al.  Thermodynamic database for proteins: features and applications. , 2010, Methods in molecular biology.

[46]  K. Dill,et al.  Statistical potentials extracted from protein structures: how accurate are they? , 1996, Journal of molecular biology.

[47]  L Serrano,et al.  Development of the multiple sequence approximation within the AGADIR model of alpha-helix formation: comparison with Zimm-Bragg and Lifson-Roig formalisms. , 1997, Biopolymers.

[48]  S J Wodak,et al.  Contribution of the hydrophobic effect to protein stability: analysis based on simulations of the Ile-96----Ala mutation in barnase. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[49]  Ruth Nussinov,et al.  Molecular dynamics simulations of the unfolding of beta(2)-microglobulin and its variants. , 2003, Protein engineering.

[50]  Jianpeng Ma,et al.  CHARMM: The biomolecular simulation program , 2009, J. Comput. Chem..

[51]  T. N. Bhat,et al.  The Protein Data Bank , 2000, Nucleic Acids Res..

[52]  Emil Alexov,et al.  A missense mutation in CLIC2 associated with intellectual disability is predicted by in silico modeling to affect protein stability and dynamics , 2011, Proteins.

[53]  F A Quiocho,et al.  Extensive features of tight oligosaccharide binding revealed in high-resolution structures of the maltodextrin transport/chemosensory receptor. , 1997, Structure.

[54]  Zhe Zhang,et al.  In Silico and In Vitro Investigations of the Mutability of Disease-Causing Missense Mutation Sites in Spermine Synthase , 2011, PloS one.

[55]  D Gilis,et al.  Stability changes upon mutation of solvent-accessible residues in proteins evaluated by database-derived potentials. , 1996, Journal of molecular biology.

[56]  P. Kollman,et al.  Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. , 2000, Accounts of chemical research.

[57]  A. Karshikoff,et al.  Model for calculation of electrostatic interactions in unfolded proteins. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[58]  D. Beveridge,et al.  Free energy via molecular simulation: applications to chemical and biomolecular systems. , 1989, Annual review of biophysics and biophysical chemistry.

[59]  J. Hermans,et al.  ES/IS: estimation of conformational free energy by combining dynamics simulations with explicit solvent with an implicit solvent continuum model. , 1999, Biophysical chemistry.

[60]  M. Permutt,et al.  Missense polymorphism in the human carboxypeptidase E gene alters enzymatic activity , 2001, Human mutation.

[61]  Iosif I. Vaisman,et al.  Accurate prediction of stability changes in protein mutants by combining machine learning with structure based computational mutagenesis , 2008, Bioinform..

[62]  Piero Fariselli,et al.  Predicting protein stability changes from sequences using support vector machines , 2005, ECCB/JBI.

[63]  Akinori Sarai,et al.  ProTherm, version 4.0: thermodynamic database for proteins and mutants , 2004, Nucleic Acids Res..

[64]  Emil Alexov,et al.  On the electrostatic component of protein-protein binding free energy , 2008, PMC biophysics.

[65]  Hongyi Zhou,et al.  Distance‐scaled, finite ideal‐gas reference state improves structure‐derived potentials of mean force for structure selection and stability prediction , 2002, Protein science : a publication of the Protein Society.

[66]  K Nishikawa,et al.  Experimental verification of the 'stability profile of mutant protein' (SPMP) data using mutant human lysozymes. , 1999, Protein engineering.

[67]  John Moult,et al.  Three‐dimensional structural location and molecular functional effects of missense SNPs in the T cell receptor Vβ domain , 2003, Proteins.

[68]  P. Kollman,et al.  Exhaustive mutagenesis in silico: Multicoordinate free energy calculations on proteins and peptides , 2000, Proteins.

[69]  Hongmei Wang,et al.  Tumor Suppressor von Hippel-Lindau (VHL) Stabilization of Jade-1 Protein Occurs through Plant Homeodomains and Is VHL Mutation Dependent , 2004, Cancer Research.

[70]  L. Serrano,et al.  Predicting changes in the stability of proteins and protein complexes: a study of more than 1000 mutations. , 2002, Journal of molecular biology.

[71]  Ruth Nussinov,et al.  Molecular dynamics simulations of the unfolding of β2‐microglobulin and its variants , 2003 .

[72]  J. Moult,et al.  SNPs, protein structure, and disease , 2001, Human mutation.

[73]  François Stricher,et al.  The FoldX web server: an online force field , 2005, Nucleic Acids Res..

[74]  Peng Yue,et al.  SNPs3D: Candidate gene and SNP selection for association studies , 2006, BMC Bioinformatics.

[75]  R. Abagyan,et al.  Large‐scale prediction of protein geometry and stability changes for arbitrary single point mutations , 2004, Proteins.

[76]  M. Gromiha,et al.  Prediction of protein stability upon point mutations. , 2007, Biochemical Society transactions.

[77]  Lee.  Testing homology modeling on mutant proteins: predicting structural and thermodynamic effects in the Ala98-->Val mutants of T4 lysozyme. , 1995, Folding & design.

[78]  Piero Fariselli,et al.  Predicting Free Energy Contribution to the Conformational Stability of Folded Proteins From the Residue Sequence with Radial Basis Function Networks , 1995, ISMB.

[79]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[80]  R. Nussinov,et al.  Conformational ensembles, signal transduction and residue hot spots: application to drug discovery. , 2010, Current opinion in drug discovery & development.

[81]  Piero Fariselli,et al.  A neural-network-based method for predicting protein stability changes upon single point mutations , 2004, ISMB/ECCB.

[82]  R L Jernigan,et al.  Protein stability for single substitution mutants and the extent of local compactness in the denatured state. , 1994, Protein engineering.

[83]  W. C. Still,et al.  Semianalytical treatment of solvation for molecular mechanics and dynamics , 1990 .

[84]  David Haussler,et al.  LS-SNP: large-scale annotation of coding non-synonymous SNPs based on multiple information sources , 2005, Bioinform..

[85]  Thorsten Joachims,et al.  Learning to classify text using support vector machines - methods, theory and algorithms , 2002, The Kluwer international series in engineering and computer science.

[86]  Emil Alexov,et al.  Structural and functional consequences of single amino acid substitutions in the pyrimidine base binding pocket of Escherichia coli CMP kinase , 2007, The FEBS journal.

[87]  Mauno Vihinen,et al.  Performance of protein stability predictors , 2010, Human mutation.

[88]  H Oschkinat,et al.  Improving the refolding yield of interleukin-4 through the optimization of local interactions. , 2000, Journal of biotechnology.

[89]  T. Hynes,et al.  The crystal structure of staphylococcal nuclease refined at 1.7 Å resolution , 1991, Proteins.

[90]  Zhe Zhang,et al.  Computational analysis of missense mutations causing Snyder‐Robinson syndrome , 2010, Human mutation.

[91]  Marianne Rooman,et al.  Prediction of stability changes upon single-site mutations using database-derived potentials , 1999 .