| Nucleic Acids Research | Pages |
Unconventional helical phasing of repetitive DNA motifs reveals their relative bending contributions
Introduction
Materials And Methods
Purification, labeling and ligation of oligonucleotides and markers
Two-dimensional gel electrophoresis
Computer modeling
Results And Discussion
Conclusions
Acknowledgements
References
Unconventional helical phasing of repetitive DNA motifs reveals their relative bending contributions
ABSTRACT
INTRODUCTION
DNA bending was initially noticed with repeated adenine tracts (A-tracts) because such molecules migrate more slowly on polyacrylamide gels than does DNA of the same length but composed of random sequence (1,2). This observation provided the experimental support for an earlier hypothesis stating that local deflections within DNA, when repeated in phase with the helical turn, can add up and lead to macroscopic DNA curvature (3). Further attempts to characterize DNA curvature involved gel mobility experiments on repetitive synthetic DNA sequences (4-6). The idea behind this approach was that local structural variations of DNA (bends), if separated by an integral number of helical turns, will retain the same global direction and add coherently into a defined macroscopic shape (curve). Conversely, local DNA bends with the same direction will cancel each other when separated by half a helical turn (6). Although anomalous gel mobility of curved DNA is not yet fully understood theoretically (7), this method has often provided unique information about the direction and magnitude of intrinsic and protein-induced DNA bending (2,8). It was also used in developing an important predictive model for DNA structure (9). However, a serious problem with gel mobility studies was that they were usually conducted in low ionic strength buffers, which are not necessarily representative of in vivo conditions. Another weakness in this approach has been its relative inability to differentiate flexible versus bent DNA sequences (deformable versus deformed) (10-12). In addition, it has been shown that physiological ionic conditions give rise to DNA curvature in a very different way from that deduced by gel experiments without divalent cations (11-13). Finally, the sequence context in which DNA sequences of interest are located is also a major determinant of DNA structural properties (12). To incorporate these phenomena fully, higher order predictive models are required.
One of the most serious problems in defining the origin of DNA curvature has not been the limitations of gel mobility studies but rather the fact that results from crystal and solution studies disagree on the question of A-tract-associated DNA bending (14-16). Crystallographic results to date uniformly show essentially straight A-tracts (summarized in 17), whereas earlier studies in solution strongly suggest that these sequence elements are bent (5,18,19). To reconcile this dichotomy, it has been proposed that dehydrating agents used to facilitate DNA crystallization might induce structural changes in DNA (20), either directly or through their effect on water activity. This hypothesis has been investigated with gel mobility experiments and shown to be feasible (21). It is also supported by many independent experiments using different dehydrating agents and their effects on DNA structure (22-28). Although the dilemma concerning the origin of DNA bending is not yet fully resolved (12,17,29), several possible explanations have been offered (30,31). Thus, evidence currently available strongly suggests that A-tracts may be polymorphic and their specific structural details may depend on their local hydration environment.
In the present work, we have developed a multiple phasing analysis using AAAAA(A), GGGCCC and GAGAG motifs, which have been previously characterized by us in various combinations (12,31). These sequence motifs were assembled differently into several oligomers of identical base composition (Fig.
MATERIALS AND METHODS
Purification, labeling and ligation of oligonucleotides and markers
Oligonucleotides were obtained from commercial sources. Their purification, labeling and ligation have been described elsewhere (11). A BamHI linker and a HpaII digest of pBR322 were purchased from New England Biolabs and used as molecular weight markers. The BamHI linker was kinased at its 5[prime]-ends with [[gamma]-32P]ATP. Single strands were then annealed and ligated at room temperature for 4 h. The HpaII digest of pBR322 was dephosphorylated using alkaline phosphatase and subsequently labeled with [[gamma]-32P]ATP. Gel mobility experiments, both with and without 10 mM MgCl2, were performed according to a published protocol (12).
Two-dimensional gel electrophoresis
Products of ligation were analyzed by two-dimensional (2D) gel electrophoresis using methods described elsewhere (12,31). Briefly, the first dimension was on a 5% polyacrylamide gel with 25% glycerol. The gel was run at 10 V/cm until xylene cyanol migrated 9 cm on a 14 cm long gel. The gel lane was excised and placed horizontally between the slab gel plates to be used for the second dimension separation. Chloroquine phosphate was added to a final concentration of 50 µg/ml to both the running buffer and the 8% polyacrylamide mixture used for the second dimension gel. Glycerol was added to a final concentration of 10%. After running for 12-18 h at 7 V/cm, X-ray film was placed over the wet gel at room temperature and kept overnight. Sizes of linear and circular DNA fragments were determined by comparison with the original ligation products on 8% denaturing polyacrylamide gels. BamHI linker ligation products were used as additional markers.
Computer modeling
Spatial trajectories in Figures
RESULTS AND DISCUSSION
Figure
Figure 1. The sequences were designed so that one of the three used sequence motifs [AAAAA(A), GAGAG and GGGCCC] is out of phase, while the other two are in phase. For example, in AA-out the distance between centers of AAAAA(A) motifs is ~1.5 helical turns, which means that their bending contributions cancel each other out. At the same time, the GAGAG and GGGCCC motifs are in phase because their centers are separated by one or two helical turns. Similar logic applies to the other two sequences. In GC-out, the GGGCCC motif is out of phase, while the GAGAG and AAAAA(A) motifs are separated by either one or two helical turns. In GA-out, the GAGAG motif is out of phase, while both the AAAAA(A) and GGGCCC motifs are in phase. All sequences have identical base composition. Complementary strands were designed to produce 3 bp 5[prime]-end overhangs after annealing. Numbers indicate the interval between centers of identical motifs in terms of helical turns. Gel mobility experiments under standard conditions (low ionic strength buffer, no divalent cations) were used initially to determine bending properties of the oligomers from Figure Figure 2. (A) Gel retardation is AA-out < GC-out < GA-out without Mg2+ in the gel and buffer. (B) The general trend of retardation values remains, but retardation increases in the presence of 10 mM Mg2+. As noted earlier (38), the addition of divalent cations usually increases DNA curvature. However, this increase in retardation anomaly was generally smaller in sequences with A-tracts and larger for GC-rich sequences when only one of these motifs was present in the sequence (11-13). We observed an increase in retardation when 10 mM Mg2+ was added to the gel (Fig. Figure 3. Gel mobility of ligated oligomers (in 10 mM MgCl2) against the ligated BamHI linker (M1) and HpaII digest of pBR322 (M2). Lanes T1, T2 and T3 correspond to multimers of AA-out, GC-out and GA-out, respectively. It was shown earlier that cyclization experiments sometimes produce very different results from gel mobility studies, attributed to the presence of Mg2+ in cyclization experiments and its absence in gel electrophoresis experiments (11,12). The results of cyclization and gel experiments without Mg2+ are usually more similar in a qualitative sense for DNA sequences containing A-tracts (11). The pattern observed here in cyclization experiments is consistent with that seen in gel retardation analysis. Sequence AA-out failed to produce circles after ligation (Fig. Figure 4. (A) 2D cyclization gel of AA-out. (B) 2D cyclization gel of GA-out. With regard to the orientation of the picture, the first dimension is from left to right and the second is from top to bottom. Linear fragments are in the lower diagonal, as they migrate faster. Circles, when present, are located above linear molecules. Numbers show how many original oligomers were ligated in each DNA fragment. Cyclization parallels gel mobility experiments, as sequences with low retardation anomalies, e.g. AA-out, do not produce small circles, which in this assay is a hallmark of bent DNA sequences. Sequence GA-out shows the highest gel mobility anomaly (Fig. 2) and forms the smallest circles. Figure 5. Modeling five repeats of sequences from Figure 1 using the dinucleotide models of: (A) Bolshoy et al. (9); (B) Gorin et al. (34); (C) Ulyanov and James (35). Sequences AA-out, GC-out and GA-out correspond to rows I, II and III, respectively. Plots are shown in the plane of largest curvature. Figure 6. Modeling five repeats of sequences from Figure 1 using the trinucleotide models of: (A) Goodsell and Dickerson (15) based on experimental results by Satchwell et al. (43) and (B) Brukner et al. (36). Sequences AA-out, GC-out and GA-out correspond to rows I, II and III, respectively. Plots are shown in the plane of largest curvature. Figures 5 and 6 were produced using DIAMOD (32). Comparative gel mobility experiments and cyclization assays register only global DNA structural features (15,16). For that reason, our results can be explained in either of two ways: (i) the GAGAG motif is not bent, while A-tracts and the GGGCCC motifs produce bends of similar magnitude but of opposite direction; (ii) A-tracts are not bent, while GAGAG and GGGCCC motifs both bend in the same direction, with the magnitude of the latter being roughly twice that of the former (11,12). The experiments presented here cannot clearly distinguish between these possibilities. However, we can determine with certainty the relative bending contributions of sequence motifs used in this study. We will make the well-justified assumption, based on both theory and experiments, that only roll contributes significantly to bending (14,39,40). In this case, the bending can be described as an absolute difference between the roll contributions in the sequence motifs used here (16,41). For example, the GAGAG motifs cancel in GA-out and, for that reason, bending in this sequence can be described by means of a difference in the roll angles only between the AAAAA(A) and GGGCCC motifs. Since GA-out shows the largest bending, it can be inferred that roll angles within the AAAAA(A) and GGGCCC motifs have the greatest net difference (31,42). This is consistent with the notion that the roll angles of AAA and GGC are on completely opposite sides of the bending scale (15,36,43). Applying the same argument to the other two sequences, we observe that the roll-dominated bending in the GAGAG motif is between AAAAA(A) and GGGCCC. However, it must be closer to the latter because the net difference between bends in the GGGCCC and GAGAG motifs in AA-out is smaller than the difference between the AAAAA(A) and GAGAG motifs in GC-out. This result is slightly different from our earlier results showing that bending in the GAGAG motif is central between the AAAAA(A) and GGGCCC motifs (12). One explanation for this discrepancy comes from the fact that stable DNA bending seems to require phased A-tracts (11), regardless of whether A-tracts are themselves bent or straight. Since AA-out has A-tracts out of phase, it is feasible that bending of this sequence might be underestimated. It is also possible that different sequence contexts of the GAGAG motif in this analysis, compared with that in the previous study (12), contribute partially to the observed variance. As remarked above, different experimental approaches (e.g. gel mobility versus cyclization assays) and different conditions (e.g. low ionic strength buffers versus buffers with divalent cations) often produce variable results even on the same DNA sequence (11,31). In that regard, the results are not always easy to bring under a common denominator and analyze by objective theoretical methods. The sequences used in the present study overcome that problem as they show consistent results by several techniques. In addition, the monotonic increase in both gel retardation and cyclization from AA-out through GC-out to GA-out offers a good way of objectively testing the predictive models for DNA bending. To carry out this test, we generated pentamer oligos for each sequence and determined their bending properties. This size was chosen because pentamers of GA-out already produce circles (Fig. Several statistical parameters are presented in Table 1. The ratio of end-to-end distance (EED) to the contour length of a DNA molecule is one indicator of bending. For straight DNA this ratio equals one and it approaches zero as bending increases (bending brings the ends of the molecule closer together). The ratio of largest to smallest principle moments of inertia (RMI) is roughly comparable with the length/diameter ratio of the molecule. This value is large for an extended molecule, as its length is always bigger than any other dimension, but it converges to one as the molecule becomes more tightly packed by bending. Finally, the radius of gyration (RG) indicates the compactness of the molecule. When comparing DNA molecules of identical sizes, it will be smaller for curved and larger for extended molecules. For the purposes of our evaluation, a DNA molecule is straight when EED is close to 1 and RMI assumes a large double digit value (e.g. the Brukner model prediction for AA-out in Table 1). Conversely, DNA is bent when EED is close to 0 and RMI approaches 1 (e.g. the Bolshoy model prediction for GC-out in Table 1). Almost all models are successful in predicting that AA-out is essentially straight (EED close to 1 and large RMI), which is reflected in Figures Table 1. The statistics for dinucleotide predictive models by Bolshoy et al. (9), De Santis et al. (33), Gorin et al. (34) and Ulyanov and James (35) and the models by Goodsell and Dickerson (15) and Brukner et al. (36) based on trinucleotides Having said that the trinucleotide models exhibit better predictive power, we should point out that they cannot be used to distinguish between models for DNA curvature. Although the trinucleotide models reported to date all have roll values scaled from 0-10° (15,36), these models cannot support straight A-tracts and bent general sequence DNA per se, since the choice of scaling was completely arbitrary. For example, the experimental data used for these models (43,45) could have been scaled from -5° to 5° (the absolute difference between the smallest and largest roll values would also be 10°), in which case they would support bent A-tracts while still producing identical global DNA trajectories (M.Dlakic and R.E.Harrington, in preparation).
CONCLUSIONS
We have reported a novel, multiple phasing analysis to facilitate interpretation of DNA curvature studies. By combining three sequence motifs at a time and placing them into different phasing arrangements, we can obtain new information about their relative bending contributions. Extending these studies to a larger number of sequences is straightforward and could be used to construct a scale of relative bending propensities, similar to those already described (15,36).
The design of our sequences and the nature of their bending have provided a way of testing several predictive models of DNA structure. Confirming earlier results (15,44), we show that trinucleotide models are more successful in explaining our data than dinucleotide models. A potential problem with trinucleotide models is that they utilize sequence-independent values of twist, which is not true from a structural point of view. However, this incorrect assumption seems to be offset by three correct postulates that are embedded in trinucleotide models: (i) roll contributes mostly to DNA bending; (ii) tilt is insignificant compared with roll; (iii) sequence-dependent differential bendability, as opposed to fixed angular parameters used in dinucleotide models, can be used to capture structural features of DNA (all these points are reviewed in 46).
ACKNOWLEDGEMENTS
We thank Professor Steve Harvey for many constructive discussions. This work was supported by grants CA70274 and GM53517 from the NIH and from USDA Hatch Project NEV032D (R.E.H.) and by a supplementary grant from the Open Society Institute (M.D.).
REFERENCES
This article has been cited by other articles:
This page is run by Oxford University Press, Great Clarendon Street, Oxford OX2 6DP, as part of the OUP Journals
Comments and feedback: www-admin{at}oup.co.uk
Last modification: 29 Aug 1998
Copyright©Oxford University Press, 1998.
![]()
CiteULike
Connotea
Del.icio.us What's this?
![]()
![]()

![]()
![]()
![]()
A. Kanhere and M. Bansal
An assessment of three dinucleotide parameters to predict DNA curvature by quantitative comparison with experimental data
Nucleic Acids Res.,
May 15, 2003;
31(10):
2647 - 2658.
[Abstract]
[Full Text]
[PDF]
![]()
This Article ![]()
![]()
Abstract
![]()
Print PDF (120K)
![]()
Alert me when this article is cited
![]()
Alert me if a correction is posted
![]()
Services ![]()
![]()
Email this article to a friend
![]()
Similar articles in this journal
![]()
Similar articles in ISI Web of Science
![]()
Similar articles in PubMed
![]()
Alert me to new issues of the journal
![]()
Add to My Personal Archive
![]()
Download to citation manager
![]()
Search for citing articles in:
ISI Web of Science (18)
![]()
Request Permissions ![]()
Commercial Re-use Guidelines
for Open Access NAR Content
![]()
Google Scholar ![]()
![]()
Articles by Dlakic, M.
![]()
Articles by Harrington, R. E.
![]()
Search for Related Content
![]()
PubMed ![]()
![]()
PubMed Citation
![]()
Articles by Dlakic, M.
![]()
Articles by Harrington, R. E.
![]()
Social Bookmarking ![]()
![]()
What's this?