Abstract:The RNA hairpin loops represent important RNA topologies with indispensable biological functions in RNA folding and tertiary interactions. 5'-UNCG-3' and 5'-GNRA-3' RNA tetraloops are the most important classes of RNA hairpin loops. Both tetraloops are highly structured with characteristic signature three-dimensional features and are recurrently seen in functional RNAs and ribonucleoprotein particles. Explicit solvent molecular dynamics (MD) simulation is a computational technique which can efficiently complement the experimental data and provide unique structural dynamics information on the atomic scale. Nevertheless, the outcome of simulations is often compromised by imperfections in the parametrization of simplified pairwise additive empirical potentials referred to also as force fields. We have pointed out in several recent studies that a force field description of single-stranded hairpin segments of nucleic acids may be particularly challenging for the force fields. In this paper, we report a critical assessment of a broad set of MD simulations of UUCG, GAGA, and GAAA tetraloops using various force fields. First, we utilized the three widely used variants of Cornell et al. (AMBER) force fields known as ff94, ff99, and ff99bsc0. Some simulations were also carried out with CHARMM27. The simulations reveal several problems which show that these force fields are not able to retain all characteristic structural features (structural signature) of the studied tetraloops. Then we tested four recent reparameterizations of glycosidic torsion of the Cornell et al. force field (two of them being currently parametrized in our laboratories). We show that at least some of the new versions show an improved description of the tetraloops, mainly in the syn glycosidic torsion region of the UNCG tetraloop. The best performance is achieved in combination with the bsc0 parametrization of the α/γ angles. Another critically important region to properly describe RNA molecules is the anti/high-anti region of the glycosidic torsion, where there are significant differences among the tested force fields. The tetraloop simulations are complemented by simulations of short A-RNA stems, which are especially sensitive to an appropriate description of the anti/high-anti region. While excessive accessibility of the high-anti region converts the A-RNA into a senseless "ladder-like" geometry, excessive penalization of the high-anti region shifts the simulated structures away from typical A-RNA geometry to structures with a visibly underestimated inclination of base pairs with respect to the helical axis.

Severe axial bending of RNA induced by the U1A binding element present in the 3′ untranslated region of the U1A mRNA

Induced fit or conformational selection for RNA/U1A folding.

A Structural Determinant Required for Rna Editing

Myotonic dystrophy type 1 RNA crystal structures reveal heterogeneous 1 × 1 nucleotide UU internal loop conformations.

106 Severe Helical Distortion of the U1A Mrna 3′ Untranslated Region Induced by the U1A Protein Binding Site

Retained Stability of the RNA Structure in DNA Packaging Motor with a Single Mg2+ Ion Bound at the Double Mg-Clamp Structure.

NMR analysis of a loop-bulge structure of UUCGA pentaloop

Dynamic stacking of an expected branch point adenosine in duplexes containing pseudouridine-modified or unmodified U2 snRNA sites.

Structure and Mechanical Stabilities of the Three-Way Junction Motifs in Prohead RNA

Alternative Mode of E-Site tRNA Binding in the Presence of a Downstream mRNA Stem Loop at the Entrance Channel

Crystal Structure of Human UP1, the Domain of Hnrnp A1 That Contains Two RNA-recognition Motifs.

Urea destabilizes RNA by forming stacking interactions and multiple hydrogen bonds with nucleic acid bases

NMR structures of small molecules bound to a model of an RNA CUG repeat expansion

RNA adapts its flexibility to efficiently fold and resist unfolding

Visualizing RNA Conformational and Architectural Heterogeneity in Solution

Correlated Alternative Side Chain Conformations in the RNA-recognition Motif of Heterogeneous Nuclear Ribonucleoprotein A1

Shortening the HIV-1 TAR RNA Bulge by a Single Nucleotide Preserves Motional Modes over a Broad Range of Time Scales.

The Mechanism of UP1 Binding and Unfolding of Human Telomeric DNA G-quadruplex.

The origin of different bending stiffness between double-stranded RNA and DNA revealed by magnetic tweezers and simulations

A-minor tertiary interactions in RNA kink-turns. Molecular dynamics and quantum chemical analysis

Performance of Molecular Mechanics Force Fields for RNA Simulations: Stability of UUCG and GNRA Hairpins