Flexible matching of test ligands to a 3D pharmacophore using a molecular superposition force field: comparison of predicted and experimental conformations of inhibitors of three enzymes.

A computer procedure TFIT, which uses a molecular superposition force field to flexibly match test compounds to a 3D pharmacophore, was evaluated to find out whether it could reliably predict the bioactive conformations of flexible ligands. The program superposition force field optimizes the overlap of those atoms of the test ligand and template that are of similar chemical type, by applying an attractive force between atoms of the test ligand and template which are close together and of similar type (hydrogen bonding, charge, hydrophobicity). A procedure involving Monte Carlo torsion perturbations, followed by torsional energy minimization, is used to find conformations of the test ligand which cominimize the internal energy of the ligand and the superposition energy of ligand and template. The procedure was tested by applying it to a series of flexible ligands for which the bioactive conformation was known experimentally. The 15 molecules tested were inhibitors of thermolysin, HIV-1 protease or endothiapepsin for which X-ray structures of the bioactive conformation were available. For each enzyme, one of the molecules served as a template and the others, after being conformationally randomized, were fitted. The fitted conformation was then compared to the known binding geometry. The matching procedure was successful in predicting the bioactive conformations of many of the structures tested. Significant deviation from experimental results was found only for parts of molecules where it was readily apparent that the template did not contain sufficient information to accurately determine the bioactive conformation.