3D-quantitative structure-activity relationships of HEPT derivatives as HIV-1 reverse transcriptase inhibitors, based on ab initio calculations

Comparative molecular field analysis (CoMFA) has been applied to a large set of 1-[(2-hydroxyethoxy)-methyl]-6-(phenylthio)thymine (HEPT) analogues. The starting geometry of HEPT was obtained from crystallographic data of HEPT/HIV-1 reverse transcriptase (RT) complexes. The structures of 101 HEPT derivatives were considered and fully optimized by ab initio molecular orbital calculations at the HF/3-21G level. The best CoMFA model is satisfactory in both statistical significance and predictive ability. It shows excellent, high predictive ability as r(cv)(2) = 0.858. The derived model indicates the importance of steric contributions (64.4%) as well as electrostatic interactions for the HIV-1 RT inhibition. In addition, steric and electrostatic contour maps from this analysis agree well with the experimentally observed trend that there are steric interactions between the side chain of HEPT and an aromatic ring of Tyr181. It is concluded that a moderately sized group at C5 enhances contact with Tyr181 enough to push it into a position which renders the protein nonfunctional, but a smaller group has insufficient steric requirements to do this and a larger group renders the ligand too large for the cavity. The mutation-induced resistance of reverse transcriptase is explained by this analysis. The obtained results not only lead to a better understanding of structural requirements of this set of compounds for the inhibition but also enable the suggestions for new and more potent drugs.