Elucidating the inhibiting mode of AHPBA derivatives against HIV-1 protease and building predictive 3D-QSAR models

The Lamarckian genetic algorithm of AutoDock 3.0 has been used to dock 27 3(S)-amino-2(S)-hydroxyl-4-phenylbutanoic acids (AHPBAs) into the active site of HIV-1 protease (HIVPR). The binding mode was demonstrated in the aspects of the inhibitor's conformation, subsite interaction, and hydrogen bonding. The data of geometrical parameters (tau(1), tau(2), and tau(3) listed in Table 2) and root mean square deviation values as compared with the known inhibitor, kni272,(28) show that both kinds of inhibitors interact with HIVPR in a very similar way. The r(2) value of 0.860 indicates that the calculated binding free energies correlate well with the inhibitory activities. The structural and energetic differences in inhibitory potencies of AHPBAs were reasonably explored. Using the binding conformations of AHPBAs, consistent and highly predictive 3D-QSAR models were developed by performing CoMFA, CoMSIA, and HQSAR analyses. The reasonable r(corss)(2) values were 0.613, 0.530, and 0.717 for CoMFA, CoMSIA,and HQSAR models, respectively. The predictive ability of these models was validated by kni272 and a set of nine compounds that were not included in the training set. Mapping these models back to the topology of the active site of HIVPR leads to a better understanding of vital AHPBA-HIVPR interactions. Structural-based investigations and the final 3D-QSAR results provide clear guidelines and accurate activity predictions for novel HIVPR inhibitors.