We present an extensive molecular dynamics study (0.6 micros in total) on three A-RNA duplexes. The dependence of the A-RNA geometry on force fields (Parm99 and Parmbsc0) and salt strength conditions (approximately 0.18 M net-neutralizing Na(+) and approximately 0.3 M KCl) was investigated. The Parmbsc0 force field makes the A-RNA duplex more compact in comparison to the Parm99 by preventing temporary alpha/gamma t/t flips common in Parm99 simulations. Nevertheless, since the alpha/gamma t/t sub-state occurs to certain extent in experimental A-RNA structures, we consider both force fields as viable. The stabilization of the A-RNA helices caused by the Parmbsc0 force field includes visible reduction of the major groove width, increase of the base pair roll, larger helical inclination and small increases of twist. Therefore, the Parmbsc0 shifts the simulated duplexes more deeply into the A-form. Further narrowing of the deep major groove is observed in excess salt simulations, again accompanied by larger roll, inclination and twist. The cumulative difference between Parm99/lower-salt and Parmbsc0/higher-salt simulations is approximately 4-8 A for the average PP distances, and -0.7 to -2.5 degrees, -2.0 to -5.4 degrees, -2.6 to -8.6 degrees and 1.7 to 7.0 degrees for the twist, roll, inclination and propeller, respectively. The effects of the force field and salt condition are sequence-dependent. Thus, the compactness of A-RNA is sensitive to the sequence and the salt strength which may, for example, modulate the end-to-end distance of the A-RNA helix. The simulations neatly reproduce the known base pair roll re-distribution in alternating purine-pyrimidine A-RNA helices.

Institute of Biophysics Academy of Sciences of the Czech Republic Brno Czech Republic

Citace poskytuje Crossref.org

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

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

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study

Reference simulations of noncanonical nucleic acids with different χ variants of the AMBER force field: quadruplex DNA, quadruplex RNA and Z-DNA

Structure and mechanical properties of the ribosomal L1 stalk three-way junction

Molecular dynamics simulations of G-DNA and perspectives on the simulation of nucleic acid structures

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

Quantum chemical studies of nucleic acids: can we construct a bridge to the RNA structural biology and bioinformatics communities?

Molecular dynamics simulations suggest that RNA three-way junctions can act as flexible RNA structural elements in the ribosome

Protonation states of the key active site residues and structural dynamics of the glmS riboswitch as revealed by molecular dynamics

Extensive molecular dynamics simulations showing that canonical G8 and protonated A38H+ forms are most consistent with crystal structures of hairpin ribozyme

An RNA molecular switch: Intrinsic flexibility of 23S rRNA Helices 40 and 68 5'-UAA/5'-GAN internal loops studied by molecular dynamics methods