Molecular dynamics of HIV-1 reverse transcriptase indicates increased flexibility upon DNA binding

HIV-1 reverse transcriptase (RT) is one of the main targets for drugs used in the treatment of AIDS, among them, the non-nucleoside RT inhibitors (NNRTIs). The flexibility of RT unliganded and complexed to double-stranded DNA (RT/dsDNA), in water, has been studied by means of molecular dynamics. The simulations show that RT flexibility depends on its ligation state. The RT/dsDNA trajectories show larger fluctuations in the atomic positions than uncomplexed RT, particularly at the tips of the p66 fingers and thumb subdomains. This increased flexibility is consistent with the ability of the p66 fingers of the RT/dsDNA complex to close down after the binding of a deoxynucleoside triphosphate (dNTP) molecule, as observed in the crystal structures of RT/dsDNA bound to dNTP. The two complexation states present different patterns of concerted motions, indicating that the bound dsDNA alters RT flexibility. The motions of amino acid residues that form the non-nucleoside RT inhibitor binding pocket upon complexation with a NNRTI are anticorrelated with the p66 fingers (in RT/dsDNA) and correlated to the RNase H subdomain (unliganded RT). These concerted motions indicate that binding of a NNRTI could alter the flexibility of the subdomains whose motions are correlated to those of the binding pocket.