Sequence-dependent local variations of helical parameters, structure, and flexibility are crucial for molecular recognition processes involving B-DNA. A-tracts, i.e., stretches of several consecutive adenines in one strand that are in phase with the DNA helical repeat, mediate significant DNA bending. During the past few decades, there have been intense efforts to understand the sequence dependence of helical parameters in DNA. Molecular dynamics (MD) simulations can provide valuable insights into the molecular mechanism behind the relationship between sequence and structure. However, although recent improvements in empirical force fields have helped to capture many sequence-dependent B-DNA properties, several problems remain, such as underestimation of the helical twist and suspected underestimation of the propeller twist in A-tracts. Here, we employ reference quantum mechanical (QM) calculations, explicit solvent MD, and bioinformatics to analyze the underestimation of propeller twisting of A-tracts in simulations. Although we did not identify a straightforward explanation, we discovered two imbalances in the empirical force fields. The first was overestimation of stacking interactions accompanied by underestimation of base-pairing energy, which we attribute to anisotropic polarizabilities that are not reflected by the isotropic force fields. This may lead to overstacking with potentially important consequences for MD simulations of nucleic acids. The second observed imbalance was steric clash between A(N1) and T(N3) nitrogens of AT base pairs in force-field descriptions, resulting in overestimation of the AT pair stretch in MD simulations. We also substantially extend the available set of benchmark estimated CCSD(T)/CBS data for B-DNA base stacking and provide a code that allows the generation of diverse base-stacking geometries suitable for QM computations with predefined intra- and interbase pair parameters.

CEITEC Central European Institute of Technology Masaryk University Campus Bohunice 625 00 Brno Czech Republic

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

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

Laboratory of Informatics and Chemistry Faculty of Chemical Technology Institute of Chemical Technology Technicka 5 166 28 Prague Czech Republic

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

Citace poskytuje Crossref.org

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

MD simulations reveal the basis for dynamic assembly of Hfq-RNA complexes

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

Structural dynamics of propeller loop: towards folding of RNA G-quadruplex

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

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

Structure, Stiffness and Substates of the Dickerson-Drew Dodecamer

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

Molecular mechanism of preQ1 riboswitch action: a molecular dynamics study