Docking flexible ligands to macromolecular receptors by molecular shape

We present a method to explore the interaction of flexible ligands with receptors of known geometry on the basis of molecular shape. This method is an extension of that described by Kuntz et al. (J. Mol. Biol. 1982, 161, 269). The shape of a binding site on a macromolecular receptor is represented as a set of overlapping spheres. Each ligand is divided into a small set of large rigid fragments that are docked separately into the binding site and then rejoined later in the calculation. The division of ligands into separate fragments allows a degree of flexibility at the position that joins them. The rejoined fragments are then energy minimized in the receptor site. We illustrate the method with two test cases: dihydrofolate reductase/methotrexate and prealbumin/thyroxine. For each test case, the method finds binding geometries for the ligand near that observed crystallographically as well as others that provide good steric fit with the receptor.