Terminal base pairs of DNA and RNA molecules in solution are known to undergo frequent transient opening events (fraying). Accurate modeling of this process is important because of its involvement in nucleic acid end recognition and enzymatic catalysis. In this article, we describe fraying in molecular dynamics simulations with the ff99bsc0, ff99bsc0χOL3, and ff99bsc0χOL4 force fields, both for DNA and RNA molecules. Comparison with the experiment showed that while some features of fraying are consistent with the available data, others indicate potential problems with the force field description. In particular, multiple noncanonical structures are formed at the ends of the DNA and RNA duplexes. Among them are tWC/sugar edge pair, C-H edge/Watson-Crick pair, and stacked geometries, in which the terminal bases are stacked above each other. These structures usually appear within the first tens to hundreds of nanoseconds and substantially limit the usefulness of the remaining part of the simulation due to geometry distortions that are transferred to several neighboring base pairs ("end effects"). We show that stability of the noncanonical structures in ff99bsc0 may be partly linked to inaccurate glycosidic (χ) torsion potentials that overstabilize the syn region and allow for rapid anti to syn transitions. The RNA refined glycosidic torsion potential χOL3 provides an improved description and substantially more stable MD simulations of RNA molecules. In the case of DNA, the χOL4 correction gives only partial improvement. None of the tested force fields provide a satisfactory description of the terminal regions, indicating that further improvement is needed to achieve realistic modeling of fraying in DNA and RNA molecules.

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

Institute of Organic Chemistry and Biochemistry Academy of Sciences of the Czech Republic Flemingovo nam 2 166 10 Prague Czech Republic

Regional Centre of Advanced Technologies and Materials Department of Physical Chemistry Faculty of Science Palacky University 17 listopadu 12 77146 Olomouc Czech Republic

References provided by Crossref.org

Refinement of the Sugar Puckering Torsion Potential in the AMBER DNA Force Field

The interplay between peptides and RNA is critical for protoribosome compartmentalization and stability

Comprehensive Assessment of Force-Field Performance in Molecular Dynamics Simulations of DNA/RNA Hybrid Duplexes

Assessing the Current State of Amber Force Field Modifications for DNA─2023 Edition

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

RNA Structural Dynamics As Captured by Molecular Simulations: A Comprehensive Overview

Assessing the Current State of Amber Force Field Modifications for DNA

Long-timescale dynamics of the Drew-Dickerson dodecamer

Hairpins participating in folding of human telomeric sequence quadruplexes studied by standard and T-REMD simulations

Extended molecular dynamics of a c-kit promoter quadruplex