Conformational and Helicoidal Analysis of 30 PS of Molecular Dynamics on the d(CGCGAATTCGCG) Double Helix: “Curves”, Dials and Windows

Abstract

A new procedure for the analysis of the structure and molecular dynamics of duplex DNA is introduced, in which comprehensive visualization of results and pattern recognition is greatly facilitated. The method involves determining the values of the conformational and helicoidal parameters for each structure entering the analysis using the method “Curves” developed by Lavery and Sklenar, J. Biomol. Str. Dyn. 6, 63 (1988), followed by a novel computer graphic display of the results. The graphic display is organized systematically using conformation wheels, or “dials”, for each IUPAC torsional parameter and “windows” on the range of values assumed by the linear and angular helicoidal parameters, and is presented in a form isomorphous with the structure per se. The complete time evolution of the conformational and helicoidal parameters of a DNA double helix can then be depicted in a set of six composite figures. Dynamical aspects of helix bending are also subsumed in this analysis. The procedure is illustrated with an analysis of the structures of canonical A and B forms of DNA and the 300° K native dodecamer duplex d(CGCGAAITCGCG). The “dials and windows” are then used for a comprehensive analysis of 30 psec of molecular dynamics on the dodecamer in the vicinity of a canonical B-DNA energy minimum. This involves presentation of the time evolution of 206 conformational and 230 helicoidal parameters for the dodecamer. A number of interesting structural features can be recognized in the analysis, including crankshaft motions, B_I - B_II transitions, sugar repuckerings, and a description of spontaneous helix bending at what corresponds to the 1° and 2° “hinge points” indicated in the crystal structure. Our approach is expected to be directly useful for critical analysis of the effects of various assumptions about force field parameters, hydration and electrostatic effects and thus contribute to the development of reliable simulation protocols for nucleic acid systems. Extension of the method to present differential changes in conformational and helicoidal parameters is expected to be valuable for the analysis of structural and molecular dynamics studies of the reorganization and adaptation of DNA on complexation with various drags and regulatory proteins.