We have carried out an extended reference set of explicit solvent molecular dynamics simulations (63 simulations with 8.4 μs of simulation data) of canonical A-RNA duplexes. Most of the simulations were done using the latest variant of the Cornell et al. AMBER RNA force field bsc0χ(OL3), while several other RNA force fields have been tested. The calculations show that the A-RNA helix compactness, described mainly by geometrical parameters inclination, base pair roll, and helical rise, is sequence-dependent. In the calculated set of structures, the inclination varies from 10° to 24°. On the basis of simulations with modified bases (inosine and 2,6-diaminopurine), we suggest that the sequence-dependence of purely canonical A-RNA double helix is caused by the steric shape of the base pairs, i.e., the van der Waals interactions. The electrostatic part of stacking does not appear to affect the A-RNA shape. Especially visible is the role of the minor groove amino group of purines. This resembles the so-called Dickerson-Calladine mechanical rules suggested three decades ago for the DNA double helices. We did not identify any long-living backbone substate in A-RNA double helices that would resemble, for example, the B-DNA BI/BII dynamics. The variability of the A-RNA compactness is due to mutual movements of the consecutive base pairs coupled with modest change of the glycosidic χ torsion. The simulations further show that the A-RNA compactness is modestly affected by the water model used, while the effect of ionic conditions, investigated in the range from net-neutral condition to ~0.8 M monovalent ion excess salt, is smaller.

Institute of Biophysics Academy of Sciences of the Czech Republic Kralovopolska 135 612 65 Brno Czech Republic

Citace poskytuje Crossref.org

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

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes

Improving the Performance of the Amber RNA Force Field by Tuning the Hydrogen-Bonding Interactions

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

Structural study of the Fox-1 RRM protein hydration reveals a role for key water molecules in RRM-RNA recognition

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

The role of an active site Mg(2+) in HDV ribozyme self-cleavage: insights from QM/MM calculations

Isosteric and nonisosteric base pairs in RNA motifs: molecular dynamics and bioinformatics study of the sarcin-ricin internal loop

Toward Improved Description of DNA Backbone: Revisiting Epsilon and Zeta Torsion Force Field Parameters

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study