Abstract:A novel dynamical protocol for finding the low‐energy conformations of a protein‐ligand complex is described. The energy functions examined consist of an empirical force field with four different dielectric screening models; the generalized Born/surface area model also is examined. Application of the method to three complexes of known crystal structure provides insights into the energy functions used for selecting low‐energy docked conformations and into the structure of the binding‐energy surface. Evidence is presented that the local energy minima of a ligand in a binding site are arranged in a hierarchical fashion. This observation motivates the construction of a hierarchical docking algorithm that substantially enriches the population of ligand conformations close to the crystal conformation. The algorithm is also adapted to permit docking into a flexible binding site and preliminary tests of this method are presented. Proteins 33:475–495, 1998. © 1998 Wiley‐Liss, Inc.

A hierarchical method for generating low‐energy conformers of a protein‐ligand complex