A quantitative measure of DNA curvature enabling the comparison of predicted structures

Evaluation of DNA bending models in their capacity to predict electrophoretic migration anomalies of satellite DNA sequences

DNA containing a sequence that generates a local curvature exhibits a pronounced retardation in electrophoretic mobility. Various theoretical models have been proposed to explain relationship between DNA structural features and migration anomaly. Here, we studied the capacity of 15 static wedge-bending models to predict electrophoretic behavior of 69 satellite monomers derived from four divergent families. All monomers exhibited retarded mobility in PAGE corresponding to retardation factors ranging 1.02-1.54. The curvature varied both within and across the groups and correlated with the number, position, and lengths of A-tracts. Two dinucleotide models provided strong correlation between gel mobility and curvature prediction; two trinucleotide models were satisfactory while remaining dinucleotide models provided intermediate results with reliable prediction for subsets of sequences only. In some cases, similarly shaped molecules exhibited relatively large differences in mobility and vice versa. Generally less accurate predictions were obtained in groups containing less homogeneous sequences possessing distinct structural features. In conclusion, relatively universal theoretical models were identified suitable for the analysis of natural sequences known to harbor relatively moderate curvature. These models could be potentially applied to genome wide studies. However, in silico predictions should be viewed in context of experimental measurement of intrinsic DNA curvature.

Static and statistical bending of DNA evaluated by Monte Carlo simulations

To investigate the influence of thermal fluctuations on DNA curvature the Metropolis procedure at 300 K was applied to B-DNA decamers containing A5.T5 and A4.T4 blocks. Monte Carlo simulations have confirmed the DNA bending anisotropy: B-DNA bends most easily in a groove direction (roll). The A5.T5 block is more rigid than the other sequences; the pyrimidine-purine dimers are found to be the most flexible. For A5TCTCT, A5CTCTC, and A5GAGAG, the average bend angle per decamer is 20-25 degrees in a direction toward the minor groove in the center of the A5.T5 tract, which is consistent with both the "junction" and "wedge AA" models. However, in A5T5, A4T4CG, and T4A4GC, bending is directed into the grooves at the 5' and 3' ends of purine tracts. Thus, directionality of bending caused by An.Tn blocks strongly depends on their neighboring sequences. These calculations demonstrate that the sequence-dependent variation of the minor-groove width mimics the observed hydroxyl radical cleavage pattern. To estimate the effect of fluctuations on the overall shape of curved DNA fragments, longer pieces of DNA (up to 200 base pairs) were generated. For sequences with strong curvature (A5X5 and A4T4CG), the static model and Monte Carlo ensemble give similar results but, for moderately and slightly curved sequences (A5T5 or T4A4GC), the static model predicts a much smaller degree of bending than does the statistical representation. Considering fluctuations is important for quantitative interpretation of the gel electrophoresis measurements of DNA curvature, where both the static and statistical bends are operative.

Anomalous gel migration of DNA oligomers containing multiple conformational junctions

We have previously shown that a short 16 base pair DNA oligomer can accommodate a B-Z conformational junction [Sheardy, R. D., & Winkle, S. A. (1989) Biochemistry 28, 720-725]. Results from 1H NMR studies indicated that only three base pairs were involved in the junction and that one of these base pairs was highly distorted. Being interested in the nature of this distortion, we constructed DNA oligomers which have the potential to contain multiple B-Z junctions for polyacrylamide electrophoretic studies. We report that the mobilities displayed by these molecules through acrylamide gels in the absence and presence of cobalt suggest that these molecules run shorter than they actually are. This anomalous migration may be due to structural/dynamic properties of the DNA helix manifested by the periodic distortions of the potential B-Z junctions.

A comparison of six DNA bending models

The predictions of six DNA bending models were compared with experimental relative mobility data. The study showed that all the models are reasonably accurate in predicting bending in synthetic sequences and in a natural sequence. The least accurate of these models is the Calladine-Dickerson model. The most consistent model is the ApA Wedge, possibly because it distributes the bends into base-roll and base-tilt components.

Computer modelling of DNA structures involved in chromosome maintenance

Sequence-dependent DNA bending of synthetic and natural molecules was studied by computer analysis. Modelling of synthetic oligonucleotides and of 107 kb of natural sequences gave results which closely resembled published electrophoretic data, demonstrating the powerful predictive capacity of the procedure. The analysis was extended to the study of DNA structures involved in chromosome maintenance. Centromeric DNAs from yeast were found to have sequences in their functional elements which cause them to be unusually straight. Autonomous replicating sequences were found to have two structural domains, one consisting of unusually straight sequences surrounding the consensus and the other of bending elements in flanking DNA. In addition to a structural homology, centromeric and autonomous replicating sequences share common sequence elements. These observations show that computer modelling of natural sequences is a viable approach to the study of the biological implications of alternative DNA structures.