Molecular Dynamics Simulations of Viral RNA Polymerases Link Conserved and Correlated Motions of Functional Elements to Fidelity

Abstract The viral RNA-dependent RNA polymerase (RdRp) is essential for multiplication of all RNA viruses. The sequence diversity of an RNA virus population contributes to its ability to infect the host. This diversity emanates from errors made by the RdRp during RNA synthesis. The physical basis for RdRp fidelity is unclear but is linked to conformational changes occurring during the nucleotide-addition cycle. To understand RdRp dynamics that might influence RdRp function, we have analyzed all-atom molecular dynamics simulations on the nanosecond timescale of four RdRps from the picornavirus family that exhibit 30–74% sequence identity. Principal component analysis showed that the major motions observed during the simulations derived from conserved structural motifs and regions of known function. The dynamics of residues participating in the same biochemical property, for example, RNA binding, nucleotide binding or catalysis, were correlated even when spatially distant on the RdRp structure. The conserved and correlated dynamics of functional structural elements suggest coevolution of dynamics with structure and function of the RdRp. Crystal structures of all picornavirus RdRps exhibit a template–nascent RNA duplex channel too small to fully accommodate duplex RNA. Simulations revealed opening and closing motions of the RNA and nucleoside triphosphate channels, which might be relevant to nucleoside triphosphate entry, inorganic pyrophosphate exit and translocation. A role for nanosecond timescale dynamics in RdRp fidelity is supported by the altered dynamics of the high-fidelity G64S derivative of PV RdRp relative to wild-type enzyme.