Different approaches toward an automatic structural alignment of drug molecules: applications to sterol mimics, thrombin and thermolysin inhibitors.

A relative comparison of the binding properties of different drug molecules requires their mutual superposition with respect to various alignment criteria. In order to validate the results of different alignment methods, the crystallographically observed binding geometries of ligands in the pocket of a common protein receptor have been used. The alignment function in the program SEAL that calculates the mutual superposition of molecules has been optimized with respect to these references. Across the reference data set, alignments could be produced that show mean rms deviations of approximately 1 A compared to the experimental situation. For structures with obvious skeletal similarities a multiple-flexible fit, linking common pharmacophoric groups by virtual springs, has been incorporated into the molecular mechanics program MOMO. In order to combine conformational searching with comparative alignments, the optimized SEAL approach has been applied to sets of conformers generated by MIMUMBA, a program for conformational analysis. Multiple-flexible fits have been calculated for inhibitors of ergosterol biosynthesis. Sets of different thrombin and thermolysin inhibitors have been conformationally analyzed and subsequently aligned by a combined MIMUMBA/SEAL approach. Since for these examples crystallographic data on their mutual alignment are available, an objective assessment of the computed results could be performed. Among the generated conformers, one geometry could be selected for the thrombin and thermolysin inhibitors that approached reasonably well the experimentally observed alignment.