Calculation and analysis of low frequency normal modes for DNA

Normal mode calculations for two alternating sequence dodecamers in A, B, and Z conformations have been performed in dihedral angle space extended to endocyclic valence angles to account for sugar ring flexibility. Normal modes are analyzed in terms of helicoidal and backbone parameter variations with special attention being paid to global deformations of the double helix such as bending, twisting, or stretching. Results show that the allomorphic form of DNA has the largest influence on the flexibility of the sugar-phosphate backbone. Amplitudes of these vibrations follow the order: B > Z > A. In contrast, the amplitudes of helicoidal parameter variations are much more dependent on the base sequence. Global deformations of the double helix occur with characteristic times in the range of 1 to 10 ps and can be of mixed character, the strongest bending mode being at the same the time strongest stretching mode.