Abstract:The thermodynamic stability of double-stranded (ds)DNA depends on its sequence. It is influenced by the base pairing and stacking with neighboring bases along DNA molecules. Semiempirical schemes are available that allow us to predict the thermodynamic stability of DNA sequences based on empirically derived nearest-neighbor contributions of base pairs formed in the context of all possible nearest-neighbor base pairs. Current molecular dynamics (MD) simulations allow one to simulate the dynamics of DNA molecules in good agreement with experimentally obtained structures and available data on conformational flexibility. However, the suitability of current force field methods to reproduce dsDNA stability and its sequence dependence has been much less well tested. We have employed alchemical free-energy simulations of whole base pair transversions in dsDNA and in unbound single-stranded partner molecules. Such transversions change the sequence context but not the nucleotide content or base pairing in dsDNA and allow a direct comparison with the empirical nearest-neighbor dsDNA stability model. For the alchemical free-energy changes in the unbound single-stranded (ss)DNA partner molecules, we tested different setups assuming either complete unstacking or unrestrained simulations with partial stacking in the unbound ssDNA. The free-energy simulations predicted nearest-neighbor effects of similar magnitude, as observed experimentally but showed overall limited correlation with experimental data. An inaccurate description of stacking interactions and other possible reasons such as the neglect of electronic polarization effects are discussed. The results indicate the need to improve the realistic description of stacking interactions in current molecular mechanic force fields.

Toward Force Fields with Improved Base Stacking Descriptions

Tumuc1: A New Accurate DNA Force Field Consistent with High-Level Quantum Chemistry

The development of nucleic acids force fields: From an unchallenged past to a competitive future

Modeling DNA Structure, Elasticity and Deformations at the Base-pair Level

Tuning the stability of DNA tetrahedra with base-stacking interactions

Development of Force Field Parameters for the Simulation of Single- and Double-Stranded DNA Molecules and DNA–Protein Complexes

Simulations Meet Experiment to Reveal New Insights into DNA Intrinsic Mechanics

Structures and Dynamics of DNA Mini-Dumbbells Are Force Field Dependent

Single-molecule analysis of DNA base-stacking energetics using patterned DNA nanostructures

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

Z-DNA as a Touchstone for Additive Empirical Force Fields and a Refinement of the Alpha/Gamma DNA torsions for AMBER

Nearest-Neighbor dsDNA Stability Analysis Using Alchemical Free-Energy Simulations

Extended magnesium and calcium force field parameters for accurate ion-nucleic acid interactions in biomolecular simulations.

Extensive Evaluation of Force Fields for G-Quadruplexes

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

Bending and Base-Stacking Interactions in Double-Stranded Semiflexible Polymer

Simulations of the molecular dynamics of nucleic acids

Stacking correlation length in single-stranded DNA

Competition between Stacking and Divalent Cation-Mediated Electrostatic Interactions Determines the Conformations of Short DNA Sequences

Recent developments in empirical atomistic force fields for nucleic acids and applications to studies of folding and dynamics

Quantum Machine Learning Corrects Classical Force Fields: Stretching DNA Base Pairs in Explicit Solvent