Abstract
The Long Interspersed Element-1 (L1) retrotransposon is an ancient genetic parasite that comprises a significant part of the human genome. ORF2p is a multifunctional enzyme with endonuclease (EN) and reverse transcriptase (RT) activities that mediate target-primed reverse transcription of RNA into DNA. Structural studies of LINE-1 ORF2p consistently show a single Mg(2+) cation in the reverse transcriptase active site, conflicting with the common DNA polymerase mechanism which involves two divalent cations. We explored a reaction pathway of the DNA elongation based on the recent high-resolution ternary complex structure of the ORF2p. The combined quantum and molecular mechanics approach at the QM (PBE0-D3/6-31G**)/MM (CHARMM) level is employed for biased umbrella sampling molecular dynamics simulations followed by umbrella integration utilized to obtain the free energy profile. The nucleotidyl transfer reaction proceeds in a single step with a free energy barrier of 15.1 ± 0.8 kcal/mol, and 7.8 ± 1.2 kcal/mol product stabilization relative to reagents. Concerted nucleophilic attack by DNA O3' and proton transfer to Asp703 occur without a second catalytic metal ion. Estimated rate constant ∼60 s(-1) aligns with RT kinetics, while analysis of the Laplacian of the electron density along the cleaving P-O bond identifies a dissociative mechanism.