HIV-1 reverse transcriptase variants: molecular modeling of Y181C, V106A, L100I, and K103N mutations with nonnucleoside inhibitors using Monte Carlo simulations in combination with a linear response method.

The energies and physical descriptors for the binding of 21 novel 1-(2,6-difluorobenzyl)-2-(2,6-difluorophenyl)-benzimidazole (BPBI) analogs to HIV-1 reverse transcriptase (RT) variants Y181C, L100I, V106A, and K103N have been determined using Monte Carlo (MC) simulations. The crystallographic structure of the lead compound, 4-methyl BPBI, was used as a starting point to model the inhibitors in both the mutant bound and the unbound states. The energy terms and physical descriptors obtained from the calculations were reasonably correlated with the respective experimental EC50 values for the inhibitors against the various mutant RTs. Using the linear response correlations from the calculations, 2 novel BPBI inhibitors have been designed and simulations have been carried out. The results show the computed deltaG(binding) values match the experimental data for the analogs. Given the ongoing problem with drug resistance, the ability to predict the activity of novel analogs against variants prior to synthesis is highly advantageous.