Evaluation of a knowledge‐based potential of mean force for scoring docked protein–ligand complexes

The Biomolecular Ligand Energy Evaluation Protocol (BLEEP) is a knowledge‐based potential derived from high‐resolution X‐ray structures of protein–ligand complexes. The performance of this potential in ranking the hypothetical structures resulting from a docking study has been evaluated using fifteen protein–ligand complexes from the Protein Data Bank. In the majority of complexes BLEEP was successful in identifying the native (experimental) binding mode or an alternative of low rms deviation (from the native) as the lowest in energy. Overall BLEEP is slightly better than the DOCK energy function in discriminating native‐like modes. Even when alternative binding modes rank lower than the native structure, a reasonable energy is assigned to the latter. Breaking down the BLEEP scores into the atom–atom contributions reveals that this type of potential is grossly dominated by longer range interactions (>5 Å), which makes it relatively insensitive to small local variations in the binding site. However, despite this limitation, the lack, at present, of accurate protein–ligand potentials means that BLEEP is a promising approach to improve the filtering of structures resulting from docking programs. Moreover, BLEEP should improve with the continuously increasing number of complexes available in the PDB. © 2001 John Wiley & Sons, Inc. J Comput Chem 22: 673–688, 2001

[1]  Peter A. Kollman,et al.  FREE ENERGY CALCULATIONS : APPLICATIONS TO CHEMICAL AND BIOCHEMICAL PHENOMENA , 1993 .

[2]  M. Sippl.  Calculation of conformational ensembles from potentials of mean force. An approach to the knowledge-based prediction of local structures in globular proteins. , 1990, Journal of molecular biology.

[3]  M J Sippl,et al.  Helmholtz free energy of peptide hydrogen bonds in proteins. , 1996, Journal of molecular biology.

[4]  J. Tame,et al.  Crystallographic and calorimetric analysis of peptide binding to OppA protein. , 1999, Journal of molecular biology.

[5]  P. Koehl,et al.  Influence of protein structure databases on the predictive power of statistical pair potentials , 1998, Proteins.

[6]  Janet M. Thornton,et al.  BLEEP—potential of mean force describing protein–ligand interactions: II. Calculation of binding energies and comparison with experimental data , 1999 .

[7]  R. Cramer,et al.  Validation of the general purpose tripos 5.2 force field , 1989 .

[8]  John B. O. Mitchell,et al.  SATIS: Atom Typing from Chemical Connectivity , 1999, J. Chem. Inf. Comput. Sci..

[9]  Gennady M Verkhivker,et al.  Empirical free energy calculations of ligand-protein crystallographic complexes. I. Knowledge-based ligand-protein interaction potentials applied to the prediction of human immunodeficiency virus 1 protease binding affinity. , 1995, Protein engineering.

[10]  David S. Goodsell,et al.  Automated docking using a Lamarckian genetic algorithm and an empirical binding free energy function , 1998 .

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

[12]  R. Jernigan,et al.  Inter-residue potentials in globular proteins and the dominance of highly specific hydrophilic interactions at close separation. , 1997, Journal of molecular biology.

[13]  Paul A. Bartlett,et al.  Differential binding energy: a detailed evaluation of the influence of hydrogen-bonding and hydrophobic groups on the inhibition of thermolysin by phosphorus-containing inhibitors , 1991 .

[14]  David C. Jones,et al.  CATH--a hierarchic classification of protein domain structures. , 1997, Structure.

[15]  Ian A. Wilson,et al.  Molecular basis of crossreactivity and the limits of antibody–antigen complementarity , 1993, Nature.

[16]  T A Jones,et al.  Crystallographic refinement of human serum retinol binding protein at 2Å resolution , 1990, Proteins.

[17]  Todd J. A. Ewing,et al.  Critical evaluation of search algorithms for automated molecular docking and database screening , 1997 .

[18]  Y. Martin,et al.  A general and fast scoring function for protein-ligand interactions: a simplified potential approach. , 1999, Journal of medicinal chemistry.

[19]  Thomas Lengauer,et al.  A fast flexible docking method using an incremental construction algorithm. , 1996, Journal of molecular biology.

[20]  M. Noble,et al.  The adaptability of the active site of trypanosomal triosephosphate isomerase as observed in the crystal structures of three different complexes , 1991, Proteins.

[21]  U Heinemann,et al.  RNase T1 mutant Glu46Gln binds the inhibitors 2'GMP and 2'AMP at the 3' subsite. , 1992, Journal of molecular biology.

[22]  P. Hajduk,et al.  Evaluation of PMF scoring in docking weak ligands to the FK506 binding protein. , 1999, Journal of medicinal chemistry.

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

[24]  M J Sippl,et al.  Knowledge-based potentials for proteins. , 1995, Current opinion in structural biology.

[25]  P. Mark Rodger,et al.  Effect of available volumes on radial distribution functions , 1998 .

[26]  E. Shakhnovich,et al.  Analysis of knowledge‐based protein‐ligand potentials using a self‐consistent method , 2008, Protein science : a publication of the Protein Society.

[27]  G. Schulz,et al.  Structure of the complex between adenylate kinase from Escherichia coli and the inhibitor Ap5A refined at 1.9 A resolution. A model for a catalytic transition state. , 1992, Journal of molecular biology.

[28]  H M Holden,et al.  Structures of two thermolysin-inhibitor complexes that differ by a single hydrogen bond. , 1987, Science.

[29]  Janet M. Thornton,et al.  BLEEP—potential of mean force describing protein–ligand interactions: I. Generating potential , 1999 .

[30]  E. F. Schumacher,et al.  A Guide for the Perplexed , 1977 .

[31]  R J Fletterick,et al.  Relocating a negative charge in the binding pocket of trypsin. , 1994, Journal of molecular biology.

[32]  J C Sacchettini,et al.  Escherichia coli-derived rat intestinal fatty acid binding protein with bound myristate at 1.5 A resolution and I-FABPArg106-->Gln with bound oleate at 1.74 A resolution. , 1994, The Journal of biological chemistry.

[33]  A. Finkelstein,et al.  Why do protein architectures have boltzmann‐like statistics? , 1995, Proteins.

[34]  H. Scheraga,et al.  Medium- and long-range interaction parameters between amino acids for predicting three-dimensional structures of proteins. , 1976, Macromolecules.

[35]  J. Abrahams,et al.  The anticoagulant activation of antithrombin by heparin. , 1997, Proceedings of the National Academy of Sciences of the United States of America.

[36]  A Rojnuckarin,et al.  Knowledge‐based interaction potentials for proteins , 1999, Proteins.

[37]  William L. Jorgensen,et al.  OPLS potential functions for nucleotide bases. Relative association constants of hydrogen-bonded base pairs in chloroform , 1991 .

[38]  P. Kollman,et al.  An all atom force field for simulations of proteins and nucleic acids , 1986, Journal of computational chemistry.

[39]  P Willett,et al.  Development and validation of a genetic algorithm for flexible docking. , 1997, Journal of molecular biology.

[40]  Chao Zhang,et al.  Extracting contact energies from protein structures: A study using a simplified model , 1998, Proteins.

[41]  G. Casari,et al.  Identification of native protein folds amongst a large number of incorrect models. The calculation of low energy conformations from potentials of mean force. , 1990, Journal of molecular biology.

[42]  A M Hassell,et al.  Hydroxyethylene isostere inhibitors of human immunodeficiency virus-1 protease: structure-activity analysis using enzyme kinetics, X-ray crystallography, and infected T-cell assays. , 1992, Biochemistry.

[43]  J L Sussman,et al.  Refined crystal structure of dogfish M4 apo-lactate dehydrogenase. , 1989, Journal of molecular biology.

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

[45]  G. Klebe,et al.  Knowledge-based scoring function to predict protein-ligand interactions. , 2000, Journal of molecular biology.

[46]  F. Melo,et al.  Novel knowledge-based mean force potential at atomic level. , 1997, Journal of molecular biology.

[47]  S E Ealick,et al.  Application of crystallographic and modeling methods in the design of purine nucleoside phosphorylase inhibitors. , 1993, Proceedings of the National Academy of Sciences of the United States of America.

[48]  J. Dunitz,et al.  Can statistical analysis of structural parameters from different crystal environments lead to quantitative energy relationships , 1988 .

[49]  I A Vakser.  Long-distance potentials: an approach to the multiple-minima problem in ligand-receptor interaction. , 1996, Protein engineering.