BLEEP—potential of mean force describing protein–ligand interactions: II. Calculation of binding energies and comparison with experimental data

We have developed BLEEP (biomolecular ligand energy evaluation protocol), an atomic level potential of mean force (PMF) describing protein–ligand interactions. Here, we present four tests designed to assess different attributes of BLEEP. Calculating the energy of a small hydrogen‐bonded complex allows us to compare BLEEP's description of this system with a quantum‐chemical description. The results suggest that BLEEP gives an adequate description of hydrogen bonding. A study of the relative energies of various heparin binding geometries for human basic fibroblast growth factor (bFGF) demonstrates that BLEEP performs excellently in identifying low‐energy binding modes from decoy conformations for a given protein–ligand complex. We also calculate binding energies for a set of 90 protein–ligand complexes, obtaining a correlation coefficient of 0.74 when compared with experiment. This shows that BLEEP can perform well in the difficult area of ranking the interaction energies of diverse complexes. We also study a set of nine serine proteinase–inhibitor complexes; BLEEP's good performance here illustrates its ability to determine the relative energies of a series of similar complexes. We find that a protocol for incorporating solvation does not improve correlation with experiment. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 1177–1185, 1999

[1]  J. Thornton,et al.  Satisfying hydrogen bonding potential in proteins. , 1994, Journal of molecular biology.

[2]  G M Verkhivker,et al.  A mean field model of ligand-protein interactions: implications for the structural assessment of human immunodeficiency virus type 1 protease complexes and receptor-specific binding. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[3]  W. Pitt,et al.  New methods for the analysis of the protein-solvent interface , 1995 .

[4]  G. Waksman,et al.  FGF binding and FGF receptor activation by synthetic heparan-derived di- and trisaccharides. , 1995, Science.

[5]  Anthony J. Stone,et al.  Distributed multipole analysis, or how to describe a molecular charge distribution , 1981 .

[6]  E. Katchalski‐Katzir,et al.  Molecular surface recognition: determination of geometric fit between proteins and their ligands by correlation techniques. , 1992, Proceedings of the National Academy of Sciences of the United States of America.

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

[8]  Hans-Joachim Böhm,et al.  The development of a simple empirical scoring function to estimate the binding constant for a protein-ligand complex of known three-dimensional structure , 1994, J. Comput. Aided Mol. Des..

[9]  K. Ng,et al.  A reliable semi-empirical approach for evaluating the isotropic intermolecular forces between closed-shell systems. , 1979 .

[10]  M. Sippl,et al.  Helmholtz free energies of atom pair interactions in proteins. , 1996, Folding & design.

[11]  Anthony J. Stone,et al.  An intermolecular perturbation theory for the region of moderate overlap , 1984 .

[12]  William R. Pitt,et al.  AQUARIUS2: Knowledge‐based modeling of solvent sites around proteins , 1993, J. Comput. Chem..

[13]  P. Kraulis A program to produce both detailed and schematic plots of protein structures , 1991 .

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

[15]  M. Sternberg,et al.  Modelling protein docking using shape complementarity, electrostatics and biochemical information. , 1997, Journal of molecular biology.

[16]  J. Thornton,et al.  Protein recognition of adenylate: an example of a fuzzy recognition template. , 1996, Journal of molecular biology.

[17]  D C Rees,et al.  Heparin Structure and Interactions with Basic Fibroblast Growth Factor , 1996, Science.

[18]  A. Itai,et al.  A new method for predicting binding free energy between receptor and ligand , 1998, Proteins.

[19]  John B. O. Mitchell,et al.  The nature of the N  H…︁OC hydrogen bond: An intermolecular perturbation theory study of the formamide/formaldehyde complex , 1990 .

[20]  Gennady M Verkhivker,et al.  Mean field analysis of FKBP12 complexes with FK506 and rapamycin: Implications for a role of crystallographic water molecules in molecular recognition and specificity , 1997, Proteins.

[21]  G J Williams,et al.  The Protein Data Bank: a computer-based archival file for macromolecular structures. , 1977, Journal of molecular biology.

[22]  John B. O. Mitchell,et al.  On the relative strengths of amide…amide and amide…water hydrogen bonds , 1991 .

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

[24]  C. DeLisi,et al.  Determination of atomic desolvation energies from the structures of crystallized proteins. , 1997, Journal of molecular biology.