Harmonic elastic constants of 3-11 bp duplex DNA fragments were evaluated using four 5 ns unrestrained molecular dynamics simulation trajectories of 17 bp duplexes with explicit inclusion of solvent and counterions. All simulations were carried out with the Cornell et al. force-field and particle mesh Ewald method for long-range electrostatic interactions. The elastic constants including anisotropic bending and all coupling terms were derived by analyzing the correlations of fluctuations of structural properties along the trajectories. The following sequences have been considered: homopolymer d(ApA)(n) and d(GpG)(n), and alternating d(GPC)(n) and d(APT)(n). The calculated values of elastic constants are in very good overall agreement with experimental values for random sequences. The atomic-resolution molecular dynamics approach, however, reveals a pronounced sequence-dependence of the stretching and torsional rigidity of DNA, while sequence-dependence of the bending rigidity is smaller for the sequences considered. The earlier predicted twist-bend coupling emerged as the most important cross-term for fragments shorter than one helical turn. The calculated hydrodynamic relaxation times suggest that damping of bending motions may play a role in molecular dynamics simulations of long DNA fragments. A comparison of elasticity calculations using global and local helicoidal analyses is reported. The calculations reveal the importance of the fragment length definition. The present work shows that large-scale molecular dynamics simulations represent a unique source of data to study various aspects of DNA elasticity including its sequence-dependence.

J Heyrovský Institute of Physical Chemistry Dolejskova 3 Praha 8 18200 Czech Republic

References provided by Crossref.org

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

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

Long-timescale dynamics of the Drew-Dickerson dodecamer

Explaining the striking difference in twist-stretch coupling between DNA and RNA: A comparative molecular dynamics analysis

Mechanical properties of symmetric and asymmetric DNA A-tracts: implications for looping and nucleosome positioning

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

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

A measure of bending in nucleic acids structures applied to A-tract DNA

Critical effect of the N2 amino group on structure, dynamics, and elasticity of DNA polypurine tracts

Local conformational variations observed in B-DNA crystals do not improve base stacking: computational analysis of base stacking in a d(CATGGGCCCATG)(2) B<-->A intermediate crystal structure