Monte Carlo calculations on HIV-1 reverse transcriptase complexed with the non-nucleoside inhibitor 8-Cl TIBO: contribution of the L100I and Y181C variants to protein stability and biological activity

A computational model of the non-nucleoside inhibitor 8-Cl TIBO complexed with HIV-1 reverse transcriptase (RT) was constructed in order to determine the binding free energies. Using Monte Carlo simulations, both free energy perturbation and linear response calculations were carried out for the transformation of wild-type RT to two key mutants, Y181C and L100I. The newer linear response method estimates binding free energies based on changes in electrostatic and van der Waals energies and solvent-accessible surface areas. In addition, the change in stability of the protein between the folded and unfolded states was estimated for each of these mutations, which are known to emerge upon treatment with the inhibitor. Results from the calculations revealed that there is a large hydrophobic contribution to protein stability in the native, folded state. The calculated absolute free energies of binding from both the linear response, and also the more rigorous free energy perturbation method, gave excellent agreement with the experimental differences in activity. The success of the relatively rapid linear response method in predicting experimental activites holds promise for estimating the activity of the inhibitors not only against the wild-type RT, but also against key protein variants whose emergence undermines the efficacy of the drugs.