A MOLECULAR MECHANICS GRID METHOD FOR EVALUATION OF LIGAND-RECEPTOR INTERACTIONS

We present a computational method for prediction of the conformation of a ligand when bound to a macromolecular receptor. The method is intended for use in systems in which the approximate location of the binding site is known and no large-scale rearrangements of the receptor are expected upon formation of the complex. The ligand is initially placed in the vicinity of the binding site and the atomic motions of the ligand and binding site are explicitly simulated, with solvent represented by an implicit solvation model and using a grid representation for the bulk of the receptor protein. These two approximations make the method computationally efficient and yet maintain accuracy close to that of an all-atom calculation. For the benzamidine/trypsin system, we ran 100 independent simulations, in many of which the ligand settled into the low-energy conformation observed in the crystal structure of the complex. The energy of these conformations was lower than and well-separated from that of others sampled. Extensions of this method are also discussed. (C) 1995 by John Wiley and Sons, Inc.