Search

Abstract Conformational adaptation of non-nucleoside reverse transcriptase inhibitor (NNRTI) via torsional flexibility is found to be very significant for targeting human immunodeficiency virus type 1 reverse transcriptase (HIV-1 RT) mutants. Catechol diether derivative including flexible torsions is new potent NNRTI with picomolar activity. Moreover, this derivative also reveals the good solubility, low toxicity and potent inhibition for HIV-1 mutants. In this study, torsional flexibility of an undecorated catechol diether compound in the binding pocket of wild type and mutants (Y181C and K103N/Y181C) HIV-1 RT is investigated by using QM/MM calculations. From the results, the uracil ring is found to exhibit more flexibility in the NNIBP. On the contrary, potential energy surfaces show that high energy is encountered by changing of the corresponding torsion of the cyanovinyl aryl ring indicating the limitation for torsional flexibility. For pointing out the key interaction for the binding, the residual interaction energies are performed by means of QM calculations. Important attractive interactions through hydrogen bonds between the inhibitor and K102, K/N103, V106, and Y188 are observed. The catechol ring is proposed to be modified in order to strengthen interactions with surrounding amino acids. The results may help for the designing of new potent NNRTIs. Graphical abstract Image 1 Highlights • Conformational adaptation of NNRTI via torsional flexibility is very significant for targeting HIV-1 RT mutants. • Torsional flexibility of an undecorated catechol diether in the binding pocket is investigated by using QM/MM calculations. • The catechol ring is suggested to modify in order to gain more interactions with the surrounding amino acids. [ABSTRACT FROM AUTHOR]

The therapeutic effectiveness of the catechol diether analogs against both the wild-type and drug-resistant reverse transcriptase (RT) mutants of HIV strains are investigated by performing molecular docking and hybrid ONIOM calculations. The docking protocol has been used to predict the binding modes of the non-nucleoside inhibitors inside the active site cavity of the viral enzymes. For each enzyme-inhibitor adduct, the predicted docked poses are assessed by employing different scoring function based programs. However, the docking protocol fails to explain satisfactorily the antiviral activities of the drug molecules. Two-layered ONIOM calculations have been carried out to compute the relative binding affinities of the catechol diether derivatives to the binding pockets of RT variants. The binding efficacies of the inhibitors are significantly suppressed by the Y181C and K103N mutations, as revealed by the computed interaction energies at the ONIOM [B3LYP/6-31G(d,p):PM6] level of theory. Deformation energies for each bound ligand conformer are also estimated. The nature of interactions between the drug molecules and the active site residues are analyzed from the reduced density gradient (RDG) isosurfaces. The simulated ECD spectra support the conformational adaption upon inhibitor binding in the binding pockets of HIV strains. [ABSTRACT FROM AUTHOR]

Abstract: We have investigated the structure and dynamics of the HIV-1 reverse transcriptase (HIV-RT) active site, by modelling the active conformation of the HIV-1 RT/DNA/deoxythymidine triphosphate (dTTP) ternary complex. This has included molecular dynamics simulations with the CHARMM27 force field, and combined quantum mechanics/molecular mechanics (QM/MM) calculations. Three different ternary systems were studied to investigate the effects of different protonation states of the dTTP substrate (a deprotonated and two different mono-protonated triphosphate forms of dTTP at the active site), and the effects of different possible protonation state of potentially catalytic aspartate residues (Asp185 and Asp186) were tested. Several potentially important hydrogen-bonding interactions with amino acids and bound water molecules in the deoxyribonucleoside triphosphate (dNTP) binding pocket were examined. The model of the deprotonated form of the dTTP substrate with the two aspartates in their charged (basic) form seemed to be the most stable and its orientation was in good agreement with X-ray crystallographic structure. In addition, two different semiempirical (AM1/CHARMM and PM3/CHARMM) QM/MM methods were tested for the HIV-RT system, in structural optimizations. Both methods provided conformations of the triphosphate moiety in either fully deprotonated or mono-protonated forms, which agreed well with the experimental structure of dTTP. The only significant difference between the AM1/CHARMM and PM3/CHARMM minimized structures is that the PM3/CHARMM Pα–O3′ optimized distance (important for nucleotide addition) is longer by 0.66Å in the deprotonated system but shorter by 0.37Å in the mono-protonated triphosphate system as compared with those obtained from AM1/CHARMM minimized structure. The obtained results suggest that both of these QM/MM methods, and the stochastic boundary molecular dynamics approach applied in this work, can give reasonable results for modelling the catalytically active complex of this important enzyme. The results provide insight into the structure and interactions of the active site of this important enzyme, with implications for its mechanism, which may be useful in inhibitor design. [Copyright &y& Elsevier]