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.

Grand Canonical Monte Carlo and Deep Learning Assisted Enhanced Sampling to Characterize the Distribution of Mg2+ and Influence of the Drude Polarizable Force Field on the Stability of Folded States of the Twister Ribozyme

Learning Effective Molecular Models from Experimental Observables.

Efficient Chain Moves for Monte Carlo Simulations of A Wormlike Dna Model: Excluded Volume, Supercoils, Site Juxtapositions, Knots, and Comparisons with Random-Flight and Lattice Models

Unravelling Mg$^{2+}$-RNA binding with atomistic molecular dynamics

Watching ion-driven kinetics of ribozyme folding and misfolding caused by energetic and topological frustration one molecule at a time

Characterizing RNA oligomers using Stochastic Titration Constant-pH Metadynamics simulations

Magnesium ions mitigate metastable states in the regulatory landscape of mRNA elements

Can We Ever Develop an Ideal RNA Force Field? Lessons Learned from Simulations of UUCG RNA Tetraloop and Other Systems.

Polarizable AMOEBA Model for Simulating Mg2+·Protein·Nucleotide Complexes

Molecular Dynamics Simulations of RNA Motifs to Guide the Architectural Parameters and Design Principles of RNA Nanostructures

Molecular Dynamics Simulations of the Denaturation and Refolding of an RNA Tetraloop.

Competition among Li(+), Na(+), K(+), and Rb(+) monovalent ions for DNA in molecular dynamics simulations using the additive CHARMM36 and Drude polarizable force fields

Molecular simulations to investigate the impact of N6-methylation in RNA recognition: Improving accuracy and precision of binding free energy prediction

Optimized Magnesium Force Field Parameters for Biomolecular Simulations with Accurate Solvation, Ion-Binding, and Water-Exchange Properties

Grand-canonical molecular dynamics simulations powered by a hybrid 4D nonequilibrium MD/MC method: Implementation in LAMMPS and applications to electrolyte solutions

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

Non-specific binding of Na$^+$ and Mg$^{2+}$ to RNA determined by force spectroscopy methods

Computer Folding of RNA Tetraloops: Identification of Key Force Field Deficiencies

Adaptive Resolution Force Probe Simulations: Coarse Graining in the Ideal Gas Approximation

Atomistic Simulations Reveal Crucial Role of Metal Ions for Ligand Binding in Gd-I Riboswitch

Exploring Biomolecular Conformational Dynamics with Polarizable Force Field AMOEBA and Enhanced Sampling Method Milestoning