Comparative bending dynamics in DNA with and without regularly repeated adenine tracts

Normal-Mode Analysis of Circular DNA at the Base-Pair Level. 2. Large-Scale Configurational Transformation of a Naturally Curved Molecule

Fine structural and energetic details embedded in the DNA base sequence, such as intrinsic curvature, are important to the packaging and processing of the genetic material. Here we investigate the internal dynamics of a 200 bp closed circular molecule with natural curvature using a newly developed normal-mode treatment of DNA in terms of neighboring base-pair "step" parameters. The intrinsic curvature of the DNA is described by a 10 bp repeating pattern of bending distortions at successive base-pair steps. We vary the degree of intrinsic curvature and the superhelical stress on the molecule and consider the normal-mode fluctuations of both the circle and the stable figure-8 configuration under conditions where the energies of the two states are similar. To extract the properties due solely to curvature, we ignore other important features of the double helix, such as the extensibility of the chain, the anisotropy of local bending, and the coupling of step parameters. We compare the computed normal modes of the curved DNA model with the corresponding dynamical features of a covalently closed duplex of the same chain length constructed from naturally straight DNA and with the theoretically predicted dynamical properties of a naturally circular, inextensible elastic rod, i.e., an O-ring. The cyclic molecules with intrinsic curvature are found to be more deformable under superhelical stress than rings formed from naturally straight DNA. As superhelical stress is accumulated in the DNA, the frequency, i.e., energy, of the dominant bending mode decreases in value, and if the imposed stress is sufficiently large, a global configurational rearrangement of the circle to the figure-8 form takes place. We combine energy minimization with normal-mode calculations of the two states to decipher the configurational pathway between the two states. We also describe and make use of a general analytical treatment of the thermal fluctuations of an elastic rod to characterize the motions of the minicircle as a whole from knowledge of the full set of normal modes. The remarkable agreement between computed and theoretically predicted values of the average deviation and dispersion of the writhe of the circular configuration adds to the reliability in the computational approach. Application of the new formalism to the computed modes of the figure-8 provides insights into macromolecular motions which are beyond the scope of current theoretical treatments.

Physical rationale behind the nonlinear enthalpy-entropy compensation in DNA duplex stability

The physical-chemical sense of nonlinear entropy-enthalpy compensation based upon the standard thermodynamical parameters of high-temperature melting for doublet units in DNA duplexes has been considered. We are able to show that there are three, with no other constraints equally plausible, principal levels of DNA melting/hybridization description. First, DNA structure assembly/disassembly can be seen from the viewpoint of the conventional equilibrium thermodynamics without taking special care of the heat capacity DeltaC(p) value (by simply setting it equal to zero). Second, it is possible to assume that the DeltaC(p) is finite, but independent of temperature. At this approximation level the high-temperature DNA melting cannot be described, but only some special transition between metastable states of DNA duplexes in water solutions in the vicinity of ice melting point. Third, both the latter transition and the high-temperature DNA melting can be reproduced by one and the same approach, if the DeltaC(p) is assumed to be temperature dependent. These three approximation levels are equally justified from the nonlinear entropy-enthalpy compensation standpoint and by a generalized theory of temperature effects on themodynamical stability as is outlined here. Applicability of each of the approximation levels involved is discussed.

Wormlike chain theory and bending of short DNA

The probability distributions for bending angles in double helical DNA obtained in all-atom molecular dynamics simulations are compared with theoretical predictions. The computed distributions remarkably agree with the wormlike chain theory and qualitatively differ from predictions of the subelastic chain model. The computed data exhibit only small anomalies in the apparent flexibility of short DNA and cannot account for the recently reported AFM data. It is possible that the current atomistic DNA models miss some essential mechanisms of DNA bending on intermediate length scales. Analysis of bent DNA structures reveal, however, that the bending motion is structurally heterogeneous and directionally anisotropic on the length scales where the experimental anomalies were detected. These effects are essential for interpretation of the experimental data and they also can be responsible for the apparent discrepancy.

Evaluation of elastic properties of atomistic DNA models

A number of intriguing aspects in dynamics of double-helical DNA is related to the coupling between its macroscopic and microscopic states. A link between the elastic properties of long DNA chains and their atom-level dynamics can be established by comparing the worm-like chain model of polymer DNA with the conformational ensembles produced by molecular dynamics simulations. This problem is complicated by the complexity of the DNA structure, the small size of DNA fragments, and relatively short trajectory durations accessible in computer simulations of microscopic DNA dynamics. A careful study of all these aspects has been performed by using longer DNA fragments and increased durations of MD trajectories as compared to earlier such investigations. Special attention is paid to the necessary conditions and criteria of time convergence, and the possibility to increase the sampling by using constrained DNA models and simplified simulation conditions. It is found that dynamics of 25-mer duplexes with regular sequences agrees well with the worm-like chain theory and that accurate evaluation of DNA elastic parameters requires at least two turns of the double helix and approximately 20-ns duration of trajectories. Bond length and bond-angle constraints affect the estimates within numerical errors. In contrast, simplified treatment of solvation can strongly change the observed elastic parameters of DNA. The elastic parameters evaluated for AT- and GC-alternating duplexes reasonably agree with experimental data and suggest that, in different basepair sequences, the torsional and stretching elasticities vary stronger than the bending stiffness.

Structural and energetic origins of sequence-specific DNA bending: Monte Carlo simulations of papillomavirus E2-DNA binding sites

DNA bending is an important structural feature for indirect readout in protein-DNA recognition. The binding of papillomavirus E2 transcription factors to their DNA binding sites is associated with DNA bending, providing an attractive model system to study the origins of sequence-specific DNA bending. The consensus E2 target is of the general form ACCGN(4)CGGT with a variable four base pair region. We applied a new all-atom Monte Carlo (MC) algorithm that combines effective sampling with fast conformational equilibration. The resulting MC ensembles resemble the corresponding high-resolution crystal structures very well. Distinct bending is observed for the E2-DNA binding site with a central AATT linker in contrast to an essentially straight DNA with a central ACGT linker. Contributions of specific base pair steps to the overall bending are shown in terms of local structural parameters. The analysis of conformational substates provides new insights into the energetic origins of intrinsic DNA bending.

GANN: genetic algorithm neural networks for the detection of conserved combinations of features in DNA

Background

The multitude of motif detection algorithms developed to date have largely focused on the detection of patterns in primary sequence. Since sequence-dependent DNA structure and flexibility may also play a role in protein-DNA interactions, the simultaneous exploration of sequence- and structure-based hypotheses about the composition of binding sites and the ordering of features in a regulatory region should be considered as well. The consideration of structural features requires the development of new detection tools that can deal with data types other than primary sequence.

Results

GANN (available at ) is a machine learning tool for the detection of conserved features in DNA. The software suite contains programs to extract different regions of genomic DNA from flat files and convert these sequences to indices that reflect sequence and structural composition or the presence of specific protein binding sites. The machine learning component allows the classification of different types of sequences based on subsamples of these indices, and can identify the best combinations of indices and machine learning architecture for sequence discrimination. Another key feature of GANN is the replicated splitting of data into training and test sets, and the implementation of negative controls. In validation experiments, GANN successfully merged important sequence and structural features to yield good predictive models for synthetic and real regulatory regions.

Conclusion

GANN is a flexible tool that can search through large sets of sequence and structural feature combinations to identify those that best characterize a set of sequences.

Compensating bends in a 16-base-pair DNA oligomer containing a T(3)A(3) segment: A NMR study of global DNA curvature

In-phase ligated DNA containing T(n)A(n) segments fail to exhibit the retarded polyacrylamide gel electrophoresis (PAGE) migration observed for in-phase ligated A(n)T(n) segments, a behavior thought to be correlated with macroscopic DNA curvature. The lack of macroscopic curvature in ligated T(n)A(n) segments is thought to be due to cancellation of bending in regions flanking the TpA steps. To address this issue, solution-state NMR, including residual dipolar coupling (RDC) restraints, was used to determine a high-resolution structure of [d(CGAGGTTTAAACCTCG)2], a DNA oligomer containing a T3A3 tract. The overall magnitude and direction of bending, including the regions flanking the central TpA step, was measured using a radius of curvature, Rc, analysis. The Rc for the overall molecule indicated a small magnitude of global bending (Rc = 138 +/- 23 nm) towards the major groove, whereas the Rc for the two halves (72 +/- 33 nm and 69 +/- 14 nm) indicated greater localized bending into the minor groove. The direction of bending in the regions flanking the TpA step is in partial opposition (109 degrees), contributing to cancellation of bending. The cancellation of bending did not correlate with a pattern of roll values at the TpA step, or at the 5' and 3' junctions, of the T3A3 segment, suggesting a simple junction/roll model is insufficient to predict cancellation of DNA bending in all T(n)A(n) junction sequence contexts. Importantly, Rc analysis of structures refined without RDC restraints lacked the precision and accuracy needed to reliably measure bending.

Single-stranded breaks relax intrinsic curvature in DNA?

The recent hypothesis of a compressed backbone state as the origin of the intrinsic curvature in DNA suggested that it could result from a geometric mismatch between the partial specific backbone length and optimal base stacking. It predicted that the long-known phenomenon of static curvature in A-tract repeats may be affected by single-stranded breaks (nicks) that should relax it in a position-dependent manner. To check the aforementioned prediction, a special series of nicked DNA fragments was prepared from two mother sequences, one including phased A-tract repeats and the other being random, and the curvature was probed experimentally by gel mobility assays. In agreement with earlier reports, single-stranded breaks produce virtually no effect upon the gel mobility of the random sequence DNA. In contrast, for nicked A-tract fragments, the curvature exhibits regular periodical behavior depending upon the position of the strand break with respect to the overall bend. The modulations are rather strong, with the maximal increase in gel mobility exceeding 30% of the initial difference with respect to the reference straight DNA. This effect has not been encountered before, and it is opposite the usual nonspecific retardation caused by single-stranded breaks. The amplitude of the observed modulation is increased for phosphorylated nicks and in the presence of Mg(2+) ions.

DNA dynamics in a water drop without counterions

Because of its polyionic character, the DNA double helix is stable and biologically active only in salty aqueous media where its charge is compensated by solvent counterions. Monovalent metal ions are ubiquitous in DNA environment, and they are usually considered as the possible driving force of sequence-dependent modulations of DNA structure that make it recognizable by proteins. In an effort to directly examine this hypothesis, MD simulations of DNA in a water drop surrounded by vacuum were carried out, which relieves the requirement of charge neutrality. Surprisingly, with zero concentration of counterions, a dodecamer DNA duplex appears metastable, and its structure remains similar to that observed in experiment, including the minor groove narrowing in the dodecamer d(CGCGAATTCGCG)(2) often considered as the most evident cation effect. It is suggested that the same computational approach will allow one to simulate dynamics of long DNA chains more efficiently than with periodical boundary conditions.