HIV-1 RT catalyses the reverse transcription of viral genetic material (RNA) into double-stranded DNA, and is an important target of antiviral therapy in the treatment of AIDS. Better understanding of the structure, mechanism and functional role of residues involved in the resistance of HIV-1 RT against nucleoside-analog drugs may assist in the development of improved inhibitors, and also in understanding the effects of genetic variation on RT specificity and activity. In this study, firstly, molecular dynamics simulations (with CHARMM27) have been used to investigate binding interactions at the active site and the conformational behavior of the enzyme, then, mechanisms of deprotonation and DNA polymerization reactions have been modelled by the QM/MM method. A combined quantum mechanical and molecular mechanical (QM/MM) method (AM1/CHARMM) has been used to study the triphosphate substrate and the active site of HIV-1 reverse transcriptase complex structure, a virally-encoded enzyme. Free energy profiles for the reaction are also calculated. The obtained results provide important insight into the mechanistic activity of HIV-1 RT.