The local minima method (LMM) of pharmacophore determination: a protocol for predicting the bioactive conformation of...

Software has been developed for potential energy surface analysis and the local minima method of pharmacophore determination. LMM is rigorous and systematic and employs multiple conformations which are the local minima from the potential energy surface of each compound in the data set. It produces a series of possible pharmacophores from a postulated set of pharmacophore elements. The best pharmacophore is then determined by performing a comparative molecular field analysis (CoMFA) on each one. The pharmacophore which produces the most self-consistent model is deemed the best. Local minima on the gas-phase potential energy surface are shown to be a reasonably close approximation to protein bound conformations, and these conformations can be found through systematic conformational searches followed by minimization of the local minima. LMM was used to develop a 3D-QSAR model for dopamine beta-hydroxylase (DBH) inhibitors which was highly predictive (predictive R2 = 0.71 and standard error of predictions = 0.41). The model predicted that the phenyl and thienyl series of inhibitors were acting as bioisosteres. Examination of compounds overlayed in the model indicated a possible hydrogen bond acceptor in the DBH active site. Three tyrosine residues previously labeled by mechanism based inhibitors may be acting as the acceptor and therefore represent excellent candidates for site-directed mutagenesis studies.