Atomistic molecular dynamics simulations represent an established technique for investigation of RNA structural dynamics. Despite continuous development, contemporary RNA simulations still suffer from suboptimal accuracy of empirical potentials (force fields, ffs) and sampling limitations. Development of efficient enhanced sampling techniques is important for two reasons. First, they allow us to overcome the sampling limitations, and second, they can be used to quantify ff imbalances provided they reach a sufficient convergence. Here, we study two RNA tetraloops (TLs), namely the GAGA and UUCG motifs. We perform extensive folding simulations and calculate folding free energies (ΔGfold°) with the aim to compare different enhanced sampling techniques and to test several modifications of the nonbonded terms extending the AMBER OL3 RNA ff. We demonstrate that replica-exchange solute tempering (REST2) simulations with 12-16 replicas do not show any sign of convergence even when extended to a timescale of 120 μs per replica. However, the combination of REST2 with well-tempered metadynamics (ST-MetaD) achieves good convergence on a timescale of 5-10 μs per replica, improving the sampling efficiency by at least 2 orders of magnitude. Effects of ff modifications on ΔGfold° energies were initially explored by the reweighting approach and then validated by new simulations. We tested several manually prepared variants of the gHBfix potential which improve stability of the native state of both TLs by ∼2 kcal/mol. This is sufficient to conveniently stabilize the folded GAGA TL while the UUCG TL still remains under-stabilized. Appropriate adjustment of van der Waals parameters for C-H···O5' base-phosphate interaction may further stabilize the native states of both TLs by ∼0.6 kcal/mol.

Department of Physical Chemistry Faculty of Science Palacký University tř 17 listopadu 12 771 46 Olomouc Czech Republic

Institute of Biophysics of the Czech Academy of Sciences Královopolská 135 612 65 Brno Czech Republic

IT4Innovations VSB─Technical University of Ostrava 17 listopadu 2172 15 708 00 Ostrava Poruba Czech Republic

Regional Centre of Advanced Technologies and Materials Czech Advanced Technology and Research Institute Palacký University Olomouc Šlechtitelů 27 779 00 Olomouc Czech Republic

Scuola Internazionale Superiore di Studi Avanzati SISSA via Bonomea 265 34136 Trieste Italy

Citace poskytuje Crossref.org

Biophysics in the Czech Republic: 70th anniversary of the institute of biophysics

Sequence-Dependent Shape and Stiffness of DNA and RNA Double Helices: Hexanucleotide Scale and Beyond

RNA G-quadruplexes emerge from a compacted coil-like ensemble via multiple pathways

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

Computational Methods for Modeling Lipid-Mediated Active Pharmaceutical Ingredient Delivery

Computer Folding of Parallel DNA G-Quadruplex: Hitchhiker's Guide to the Conformational Space

Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations

Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations

Automatic Learning of Hydrogen-Bond Fixes in the AMBER RNA Force Field