Towards a MIP-based alignment and docking in computer-aided drug design.

Structural alignment of ligands in their biological conformation is a crucial step in the building of pharmacophoric models in structure-based drug design. In addition, docking algorithms are limited in some cases by the quality of the scoring functions and the limited flexibility of the environment that the different programs allow. On the other hand, GRID molecular interaction potentials (MIPs) have been used for a long time in 3D-QSAR studies. However, in most of these studies the alignment of the molecules is performed on the basis of geometrical or physico-chemical criteria that differ from the MIPs used in the partial least squares statistical analysis. We have previously developed a method to use the same scoring function for the molecular alignment and for 3D-QSAR studies. This methodology, based on the use of GRID potentials, consists in the weighted averaging of similarities of the relevant MIPs of the molecules to be aligned. Here we present a method to obtain the weights for the different GRID probes in the average based on the structural information on protein-ligand complexes for relevant systems. The method, implemented in MIPSIM, is shown to yield good accuracy in the prediction of the alignments for two systems: a set of three inhibitors of dihydrofolate reductase and a set of fifteen non-nucleoside HIV-1 reverse transcriptase inhibitors (NNRTIs). The smooth GRID potentials are shown to capture the flexible character of the active site, as opposed to traditional docking scoring energy functions.