Z-DNA duplexes are a particularly complicated test case for current force fields. We performed a set of explicit solvent molecular dynamics (MD) simulations with various AMBER force field parametrizations including our recent refinements of the ε/ζ and glycosidic torsions. None of these force fields described the ZI/ZII and other backbone substates correctly, and all of them underpredicted the population of the important ZI substate. We show that this underprediction can be attributed to an inaccurate potential for the sugar-phosphate backbone torsion angle β. We suggest a refinement of this potential, β(OL1), which was derived using our recently introduced methodology that includes conformation-dependent solvation effects. The new potential significantly increases the stability of the dominant ZI backbone substate and improves the overall description of the Z-DNA backbone. It also has a positive (albeit small) impact on another important DNA form, the antiparallel guanine quadruplex (G-DNA), and improves the description of the canonical B-DNA backbone by increasing the population of BII backbone substates, providing a better agreement with experiment. We recommend using β(OL1) in combination with our previously introduced corrections, εζ(OL1) and χ(OL4), (the combination being named OL15) as a possible alternative to the current β torsion potential for more accurate modeling of nucleic acids.

Department of Medicinal Chemistry University of Utah 30 South 2000 East Skaggs 105 Salt Lake City Utah 84112 United States

Institute of Biophysics Academy of Sciences of the Czech Republic Královopolská 135 612 65 Brno 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

Atomistic Insights Into Interaction of Doxorubicin With DNA: From Duplex to Nucleosome

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

Complex Biophysical and Computational Analyses of G-Quadruplex Ligands: The Porphyrin Stacks Back

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

Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations

Temperature-dependent elasticity of DNA, RNA, and hybrid double helices

Are kuravirus capsid diameters quantized? The first all-atom genome tracing method for double-stranded DNA viruses

Temperature-Dependent Twist of Double-Stranded RNA Probed by Magnetic Tweezer Experiments and Molecular Dynamics Simulations

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

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations

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

The influence of Holliday junction sequence and dynamics on DNA crystal self-assembly

The beginning and the end: flanking nucleotides induce a parallel G-quadruplex topology

c-kit2 G-quadruplex stabilized via a covalent probe: exploring G-quartet asymmetry

2,6-diaminopurine promotes repair of DNA lesions under prebiotic conditions

Insight into formation propensity of pseudocircular DNA G-hairpins

Importance of base-pair opening for mismatch recognition

RuvC uses dynamic probing of the Holliday junction to achieve sequence specificity and efficient resolution

Parallel G-triplexes and G-hairpins as potential transitory ensembles in the folding of parallel-stranded DNA G-Quadruplexes