Abstract
Pyrophosphate ion (PPi) release after nucleotide incorporation is a necessary step for RNA polymerase II (pol II) to enter the next nucleotide addition cycle during transcription elongation. However, the role of pol II residues in PPi release and the mechanistic relationship between PPi release and the conformational change of the trigger loop remain unclear. In this study, we constructed a Markov state model (MSM) from extensive all-atom molecular dynamics (MD) simulations in the explicit solvent to simulate the PPi release process along the pol II secondary channel. Our results show that the trigger loop has significantly larger intrinsic motion after catalysis and formation of PPi, which in turn aids PPi release mainly through the hydrogen bonding between the trigger loop residue H1085 and the (Mg–PPi)2– group. Once PPi leaves the active site, it adopts a hopping model through several highly conserved positively charged residues such as K752 and K619 to release from the pol II pore region of the secondary channel. These positive hopping sites form favorable interactions with PPi and generate four kinetically metastable states as identified by our MSM. Furthermore, our single-mutant simulations suggest that H1085 and K752 aid PPi exit from the active site after catalysis, whereas K619 facilitates its passage through the secondary channel. Finally, we suggest that PPi release could help the opening motion of the trigger loop, even though PPi release precedes full opening of the trigger loop due to faster PPi dynamics. Our simulations provide predictions to guide future experimental tests.