Correlating Protein Hot Spot Surface Analysis Using ProBiS with Simulated Free Energies of Protein-Protein Interfacial Residues

A protocol was developed for the computational determination of the contribution of interfacial amino acid residues to the free energy of protein-protein binding. Thermodynamic integration, based on molecular dynamics simulation in CHARMM, was used to determine the free energy associated with single point mutations to glycine in a protein-protein interface. The hot spot amino acids found in this way were then correlated to structural similarity scores detected by the ProBiS algorithm for local structural alignment. We find that amino acids with high structural similarity scores contribute on average -3.19 kcal/mol to the free energy of protein-protein binding and are thus correlated with hot spot residues, while residues with low similarity scores contribute on average only -0.43 kcal/mol. This suggests that the local structural alignment method provides a good approximation of the contribution of a residue to the free energy of binding and is particularly useful for detection of hot spots in proteins with known structures but undetermined protein-protein complexes.

[1]  Stefan Bruckner,et al.  Efficiency of alchemical free energy simulations. I. A practical comparison of the exponential formula, thermodynamic integration, and Bennett's acceptance ratio method , 2011, J. Comput. Chem..

[2]  Urban Bren,et al.  DNA duplex stability: the role of preorganized electrostatics. , 2010, The journal of physical chemistry. B.

[3]  Dusanka Janezic,et al.  Harmonic analysis of large systems. III. Comparison with molecular dynamics , 1995, J. Comput. Chem..

[4]  Peter Gruber,et al.  Modulators of protein-protein interactions: novel approaches in targeting protein kinases and other pharmaceutically relevant biomolecules. , 2011, Current topics in medicinal chemistry.

[5]  Ozlem Keskin,et al.  Identification of computational hot spots in protein interfaces: combining solvent accessibility and inter-residue potentials improves the accuracy , 2009, Bioinform..

[6]  B. Rost,et al.  Predicted protein–protein interaction sites from local sequence information , 2003, FEBS letters.

[7]  Stefan Bruckner,et al.  Unorthodox uses of Bennett's acceptance ratio method , 2009, J. Comput. Chem..

[8]  P. Kollman,et al.  Computational Alanine Scanning To Probe Protein−Protein Interactions: A Novel Approach To Evaluate Binding Free Energies , 1999 .

[9]  Rune Matthiesen,et al.  Methods, algorithms and tools in computational proteomics: A practical point of view , 2007, Proteomics.

[10]  Stefano Alcaro,et al.  The polymorphisms of DNA G-quadruplex investigated by docking experiments with telomestatin enantiomers. , 2012, Current pharmaceutical design.

[11]  Zsuzsanna Dosztányi,et al.  ANCHOR: web server for predicting protein binding regions in disordered proteins , 2009, Bioinform..

[12]  D. Cussac,et al.  Molecular and cellular analysis of Grb2 SH3 domain mutants: interaction with Sos and dynamin. , 1999, Journal of molecular biology.

[13]  Urban Bren,et al.  Free energy simulations of uncatalyzed DNA replication fidelity: structure and stability of T.G and dTTP.G terminal DNA mismatches flanked by a single dangling nucleotide. , 2006, The journal of physical chemistry. B.

[14]  Urban Bren,et al.  Do all pieces make a whole? Thiele cumulants and the free energy decomposition , 2007 .

[15]  Gerhard König,et al.  Non‐Boltzmann sampling and Bennett's acceptance ratio method: How to profit from bending the rules , 2011, J. Comput. Chem..

[16]  A. Warshel,et al.  Free energy of charges in solvated proteins: microscopic calculations using a reversible charging process. , 1986, Biochemistry.

[17]  R. Nussinov,et al.  Hot regions in protein--protein interactions: the organization and contribution of structurally conserved hot spot residues. , 2005, Journal of molecular biology.

[18]  B. Brooks,et al.  Self-guided Langevin dynamics simulation method , 2003 .

[19]  Thomas Simonson,et al.  Free Energy Simulations Come of Age: Protein—Ligand Recognition , 2002 .

[20]  Dusanka Janezic,et al.  ProBiS: a web server for detection of structurally similar protein binding sites , 2010, Nucleic Acids Res..

[21]  Irina S. Moreira,et al.  Protein–protein recognition: a computational mutagenesis study of the MDM2–P53 complex , 2008 .

[22]  Ozlem Keskin,et al.  A survey of available tools and web servers for analysis of protein-protein interactions and interfaces , 2008, Briefings Bioinform..

[23]  A. Bogan,et al.  Anatomy of hot spots in protein interfaces. , 1998, Journal of molecular biology.

[24]  Desmond J. Higham,et al.  Geometric De-noising of Protein-Protein Interaction Networks , 2009, PLoS Comput. Biol..

[25]  Bernard R. Brooks,et al.  CHARMMing: A New, Flexible Web Portal for CHARMM , 2008, J. Chem. Inf. Model..

[26]  Bert L de Groot,et al.  Protein thermostability calculations using alchemical free energy simulations. , 2010, Biophysical journal.

[27]  Xue-wen Chen,et al.  Sequence-based prediction of protein interaction sites with an integrative method , 2009, Bioinform..

[28]  Charles H. Bennett,et al.  Efficient estimation of free energy differences from Monte Carlo data , 1976 .

[29]  Eduardo Garcia Urdiales,et al.  Accurate Prediction of Peptide Binding Sites on Protein Surfaces , 2009, PLoS Comput. Biol..

[30]  Dusanka Janezic,et al.  Protein-Protein Binding-Sites Prediction by Protein Surface Structure Conservation , 2007, J. Chem. Inf. Model..

[31]  Wilfred F van Gunsteren,et al.  Enveloping distribution sampling: a method to calculate free energy differences from a single simulation. , 2007, The Journal of chemical physics.

[32]  Mihaly Mezei,et al.  The PT1-Ca2+ Gla domain binds to a membrane through two dipalmitoylphosphatidylserines. A computational study. , 2008, Biochemistry.

[33]  Dusanka Janezic,et al.  ProBiS algorithm for detection of structurally similar protein binding sites by local structural alignment , 2010, Bioinform..

[34]  Richard M. Jackson,et al.  Predicting protein interaction sites: binding hot-spots in protein-protein and protein-ligand interfaces , 2006, Bioinform..

[35]  Kenneth M Merz,et al.  Limits of Free Energy Computation for Protein-Ligand Interactions. , 2010, Journal of chemical theory and computation.

[36]  K. Dill,et al.  Predicting absolute ligand binding free energies to a simple model site. , 2007, Journal of molecular biology.

[37]  Woody Sherman,et al.  Computational approaches for fragment-based and de novo design. , 2010, Current topics in medicinal chemistry.

[38]  Alexandre M J J Bonvin,et al.  How proteins get in touch: interface prediction in the study of biomolecular complexes. , 2008, Current protein & peptide science.

[39]  Hoover,et al.  Canonical dynamics: Equilibrium phase-space distributions. , 1985, Physical review. A, General physics.

[40]  Dušanka Janežič,et al.  Binding-sites Prediction Assisting Protein-protein Docking. , 2011, Acta chimica Slovenica.

[41]  Urban Bren,et al.  Decomposition of the solvation free energies of deoxyribonucleoside triphosphates using the free energy perturbation method. , 2006, The journal of physical chemistry. B.

[42]  Dusanka Janezic,et al.  Protein Surface Conservation in Binding Sites , 2008, J. Chem. Inf. Model..

[43]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[44]  Song Liu,et al.  Protein binding site prediction using an empirical scoring function , 2006, Nucleic acids research.

[45]  H. Wolfson,et al.  Spatial chemical conservation of hot spot interactions in protein-protein complexes , 2007, BMC Biology.

[46]  T. Pawson,et al.  Assembly of Cell Regulatory Systems Through Protein Interaction Domains , 2003, Science.

[47]  Bernard R. Brooks,et al.  Predicting binding affinities of host-guest systems in the SAMPL3 blind challenge: the performance of relative free energy calculations , 2011, Journal of Computer-Aided Molecular Design.

[48]  Benoît Roux,et al.  Calculation of the standard binding free energy of sparsomycin to the ribosomal peptidyl‐transferase P‐site using molecular dynamics simulations with restraining potentials , 2010, Journal of molecular recognition : JMR.

[49]  Alexander D. MacKerell,et al.  All‐atom empirical force field for nucleic acids: II. Application to molecular dynamics simulations of DNA and RNA in solution , 2000 .

[50]  Stefan Bruckner,et al.  Avoiding the van der Waals endpoint problem using serial atomic insertion , 2011, J. Comput. Chem..

[51]  Thomas Simonson,et al.  Free energy simulations come of age: protein-ligand recognition. , 2002, Accounts of chemical research.

[52]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[53]  E. Lattman,et al.  High apparent dielectric constants in the interior of a protein reflect water penetration. , 2000, Biophysical journal.

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

[55]  M Karplus,et al.  Polar hydrogen positions in proteins: Empirical energy placement and neutron diffraction comparison , 1988, Proteins.

[56]  Benjamin A. Shoemaker,et al.  Inferred Biomolecular Interaction Server—a web server to analyze and predict protein interacting partners and binding sites , 2009, Nucleic Acids Res..

[57]  Stefan Bruckner,et al.  Efficiency of alchemical free energy simulations. II. Improvements for thermodynamic integration , 2011, J. Comput. Chem..

[58]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[59]  T. Straatsma,et al.  Separation‐shifted scaling, a new scaling method for Lennard‐Jones interactions in thermodynamic integration , 1994 .

[60]  W. L. Jorgensen,et al.  Comparison of simple potential functions for simulating liquid water , 1983 .

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

[62]  Dusanka Janezic,et al.  ProBiS-2012: web server and web services for detection of structurally similar binding sites in proteins , 2012, Nucleic Acids Res..

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

[64]  Jason A. Papin,et al.  Reconstruction of cellular signalling networks and analysis of their properties , 2005, Nature Reviews Molecular Cell Biology.

[65]  J. Aqvist,et al.  A new method for predicting binding affinity in computer-aided drug design. , 1994, Protein engineering.

[66]  Chris Oostenbrink,et al.  Cytochrome P450 3A4 Inhibition by Ketoconazole: Tackling the Problem of Ligand Cooperativity Using Molecular Dynamics Simulations and Free-Energy Calculations , 2012, J. Chem. Inf. Model..

[67]  Wilfred F van Gunsteren,et al.  Free energies of ligand binding for structurally diverse compounds. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[68]  Sandor Vajda,et al.  Identification of hot spots within druggable binding regions by computational solvent mapping of proteins. , 2007, Journal of medicinal chemistry.

[69]  Ozlem Keskin,et al.  HotSprint: database of computational hot spots in protein interfaces , 2007, Nucleic Acids Res..

[70]  Walter Thiel,et al.  QM/MM Free-Energy Perturbation Compared to Thermodynamic Integration and Umbrella Sampling:  Application to an Enzymatic Reaction. , 2006, Journal of chemical theory and computation.

[71]  Rachelle J Bienstock,et al.  Computational drug design targeting protein-protein interactions. , 2012, Current pharmaceutical design.

[72]  Dusanka Janezic,et al.  Protein-Protein Binding Site Prediction by Local Structural Alignment , 2010, J. Chem. Inf. Model..

[73]  A. Warshel,et al.  Calculations of antibody-antigen interactions: microscopic and semi-microscopic evaluation of the free energies of binding of phosphorylcholine analogs to McPC603. , 1992, Protein engineering.

[74]  Jennifer L. Knight,et al.  k-Dynamics Free Energy Simulation Methods , 2009 .