Non-complementary DNA helical structure induced by positive torsional stress
Non-complementary DNA helical structure induced by positive torsional stressAlexander V. Vologodskii, Xiaoping Yang and Nadrian C. Seeman*
Department of Chemistry, New York University, New York, NY 10003, USA
Received October 27, 1997;Revised and Accepted January 22, 1998
ABSTRACT
We have induced a local conformational transition by positive torsional stress in small synthetic circular DNA molecules containing cruciforms with immobile or tetramobile branched junctions. The immobile species correspond to the extruded and intruded extrema of the tetramobile junction. Under normal conditions the sequences of all the branched species prevent them from being re-absorbed into the circle. We have induced positive stress by addition of ethidium to the circle, in a low ionic strength medium. Alterations in gel electrophoretic mobility under increasing concentrations of ethidium suggest that the cruciforms undergo a transition under torsional stress. The product of this transition contains mispaired nucleotides, but interwound backbones. By comparing the electrophoretic mobilities of circles containing these structures with that of a completely complementary circle of the same length, we conclude that the twist in the mispairing region is similiar to that of completely paired species.
Since the discovery of cruciform formation in negatively supercoiled DNA (1-3) a great deal of effort has been devoted to the study of local conformational transitions in DNA secondary structure induced by the action of negative supercoiling. These transitions, formation of open regions, cruciform structures and the Z and H forms of DNA (reviewed in 4-9), usually occur at specific nucleotide sequences. Despite great diversity, the non-canonical structures share one common feature: interwinding of the complementary DNA strands is much lower in these structures than in canonical B-form DNA. Consequently, the formation of altered structures reduces the free energy of negative supercoiling. Indeed, negative torsional stress promotes formation of any structure that is underwound relative to B-form DNA. For the same reason positive torsional stress would promote formation of structures that are overwound relative to B-DNA. Such overwound structures are not known, however, and their occurrence is unlikely. Thus it would seem reasonable to believe that no local transitions in DNA can occur under the action of positive torsional stress.
The last statement suggests that the minimum energy state of torsionally relaxed DNA has all of its nucleotides paired in the regular B structure. Although this is usually true, exceptions are possible. For example, under specific solution conditions DNA segments with d(CG)n sequences adopt the Z conformation without any torsional stress (see reviews in 10,11); the action of positive torsional stress is likely to cause these segments to adopt the B conformation. A more interesting example is DNA molecules that contain stable cruciform structures derived from immobile or partially mobile 4-arm DNA branched junctions (12); these structures cannot be converted into regular linear double helices because the base pairs on the `extruded' arms lack complementarity (13). In theory, however, the cruciform structures could be destroyed if the positive torsional stress were large enough; this is true because the opposite strands could form interwound conformations, thereby reducing the torsional stress (Fig. 1a). From the Nucleic Acid Database (14) we only know of DNA structures containing a small number of contiguous mismatched nucleotides and detailed reports of their energetics are few (e.g. 15). Hence, it was not clear that the transition could be observed in practice. Nevertheless, under the action of positive torsional stress we have detected the transition of stable cruciforms into an interwound conformation containing non-complementary nucleotide pairs. Our analysis of this interwound conformation suggests that its overall twist is similar to the twist of DNA with perfect complementarity.a
Figure 1. Cruciform molecules and their transitions. (a) The transition described here. The structure on the left is a DNA minicircle containing a cruciform with an immobile (or partially mobile) branched junction whose opposite hairpin arms are not complementary to each other. The application of sufficient positive superhelical torsion by means of ethidium intercalation appears capable of inducing a conformational transition in the hairpin regions. The resulting structure appears to be an interwound conformation lacking in Watson-Crick complementarity, as indicated by the jagged segments. (b) The cruciform-containing molecular species used in this work. Three molecules containing a total of 298 nt in each strand are shown. At the lower left is the extruded immobile cruciform and at the lower right is the intruded immobile cruciform. Both of these molecules contain branch points fixed by sequence asymmetry. The differences are indicated by the 4 nucleotide pairs flanking the branch point. In the extruded immobile cruciform there are 262 nucleotide pairs in the circle to the base of the branch pont and 16 immobile pairs in the hairpins, which are capped by a dT4 loop in each strand. In the intruded immobile cruciform there are 270 nucleotide pairs in the circle and 12 immobile pairs in the hairpins. In the tetramobile cruciform shown at the top of the drawing the 4 nucleotide pairs that differentiate the two immobile cruciforms have been symmetrized, so that the molecule can undergo four steps of branch migration, as illustrated (13).
All DNA molecules in this study were synthesized on an Applied Biosystems 380B automatic DNA synthesizer, removed from the support and deprotected using standard phosphoramidite procedures (16). DNA strands have been purified by denaturing gel electrophoresis. Each of the final circular strands was synthesized from three different synthetic fragments, one of which contained the cruciform in its entirety. The portion of each strand containing the cruciform has been made in three different versions, to generate extruded immobile, intruded immobile and tetramobile cruciforms, as illustrated in Figure 1b; each strand contains 298 nt. The sequences of these strands and protocols for their phosphorylation, annealing and ligation have been published previously (13). A fourth strand was prepared complementary to strand 1 of the junction region (13), permitting assembly of a completely paired double-helical circle. In addition to a nominally relaxed circle, negatively supercoiled topoisomers of the completely paired circle were prepared by closing it in the presence of varying amounts of ethidium. These molecules were purified on two-dimensional denaturing gels. Nicked versions of each of these molecules were prepared by leaving one of the constituent strands unphosphorylated prior to ligation; nicked species were purified by non-denaturing gel electrophoresis.
All gels contained a ratio of 19:1 acrylamide:bisacrylamide, 7.0 M urea and a buffer consisting of 89 mM boric acid, 89 mM Tris, pH 8.3 (at 25°C), 2 mM Na2EDTA and 30-40% formamide. The sample was run on a 4% polyacrylamide gel. Each lane was then excised and run on a 7% denaturing polyacrylamide gel in a second dimension. Prior to loading onto gels the samples were dissolved in a denaturing loading dye containing 90% formamide and 5 mM EDTA and heated at 90°C for 3 min. Gels were run at 60°C at 31 V/cm.
Gels contained 5% polyacrylamide (19:1 acrylamide:bisacrylamide) and a buffer consisting of 89 mM boric acid, 89 mM Tris, pH 8.3 (at 25°C), and 2 mM Na2EDTA. The sample buffer was the same, but contained 5% glycerol and 0.02% each of Bromophenol blue and Xylene cyanol FF tracking dyes. Gels were run on a Hoefer SE-200 gel electrophoresis unit at 20 V/cm for 3 h at 25°C, dried onto Whatman 3MM paper and exposed to X-ray film for up to 15 h; preparative gels were not dried.
Radioactively labeled target molecules were first purified as a mixture of topoisomers using a two-dimensional (4 then 7%) polyacrylamide (19:1 acrylamide:bisacrylamide) denaturing gel. Topoisomers were then resolved on a 5% non-denaturing polyacrylamide gel (19:1 acrylamide:bisacrylamide) containing varied concentrations of ethidium bromide. Ethidium bromide was first added to the TBE running buffer and then introduced electrophoretically into the gel. These gels were run at room temperature. The constituents were then visualized by autoradiography, as described above.
Three different circular DNA molecules containing cruciforms were used in this work. The first DNA molecule consisted of a closed double helix 262 bp in length and a cruciform with restricted mobility (Fig. 1b). This is a tetramobile junction, because 4 bp at the base of each hairpin of the molecule are homologous, so they can be transformed into 8 bp of the circle by branch migration. Two other cruciforms were constructed with immobile junctions, corresponding to the extrema of the excursions of the tetramobile junction (also shown in schematic form in Fig. 1b). We will refer to these species as tetramobile, extruded immobile and intruded immobile cruciforms. In addition, a circular DNA molecule with perfect strand complementarity and the same total length was prepared.
A set of closed circular DNA molecules with different linking number (Lk) values usually appeared after cyclization of linear molecules. In our case of very short DNA circles the set could contain only very few topoisomers (17,18). We detected only one topoisomer for the extruded immobile cruciform and two topoisomers for the tetramobile and intruded immobile cruciforms. Analysis of denaturing gel mobilities combined with theoretical treatment showed that the Lk values were 25 and 26 for the molecules with two topoisomers and 25 for the only topoisomer of the extruded immobile cruciform (13). The mobility of the topoisomers was nearly identical on non-denaturing polyacrylamide gels; Figure 2a shows that they moved as a single band, except for the intruded species, which does resolve. Interestingly, in the absence of ethidium the molecules containing cruciforms migrated more rapidly than a simple DNA circle of the same length also shown on the gel.
We have observed a conformational transition in synthetic DNA minicircles containing stable cruciform structures. This transition was caused by positive torsional stress induced by the binding of ethidium to the minicircles. This means that under torsional stress the minicircles adopted conformations with higher `equilibrium' twist than they had had initially and that this change occurred in the initially branched portions of the molecules. Indeed, we have shown that a perfectly paired DNA circle does not demonstrate this transition (lane D in Fig. 2). Denaturation of the hairpins and successive interwinding of their single-stranded components would certainly decrease the torsional stress in the molecules. As a result, the molecular writhe must also diminish and therefore the mobility of the molecules on the gels decreases. It is worth noting that the contribution of the hairpins to the linking number of the DNA strands is zero.
We have used different concentrations of EtBr to introduce a positive torsional stress into the minicircles. It is interesting to estimate the superhelix density needed to cause the transition. To do this we need to calculate the amount of EtBr bound to our DNA molecules. Binding of EtBr to supercoiled DNA is well described by the excluded site model, which accounts for the change in supercoiling free energy as a result of ligand binding (25). Statistical mechanical treatment of the model gives the following equation for the number of bound EtBr molecules per base pair, [nu]
1
where K is the equilibrium binding constant of EtBr with linear DNA, c is the concentration of free EtBr in solution, A is a constant specifying the free energy of supercoiling, which depends on the DNA length, and [sigma] is the DNA superhelix density (25). The parameter [alpha] specifies the fractional unwinding (relative to the average twist) of the double helix upon ligand binding
2
where [Phi] = 26° is the angle of unwinding of the double helix resulting from binding of an EtBr molecule and [gamma] = 10.5 is the DNA helical repeat (26,27). Equation 1 can be solved numerically for any particular value of c.
There are two parameters in equation 1, K and A, which are known only approximately in our case. The binding constant K is known for a solution containing sodium ions (see 28; for a review see 29)
3
In the TBE buffer that we used in this work the sodium ions are a small fraction of all cations, present only as the counterions of EDTA. If we take into account only 4 mM Na+ in the buffer, the value of K should be ~107/M. If we suggest that the Tris cations in the buffer system might behave similarly to sodium cations we find from equation 3 the lower limit of K, 105/M. This limit is too low, probably because Tris ions are not as effective in screening DNA charges as sodium ions (30), and a lower limit of 106/M for K seems a more reasonable estimate.
Although the value of A for this DNA size has been found to be equal to 40 (17,18), the presence of the cruciform is likely to reduce the rigidity of the minicircles and the value of A along with it. We found earlier that A equals 18 for the minicircles with stable cruciforms (13). However, our analysis required additional assumptions and 18 should be considered as a lower limit for the value of A.
Thus we have two sets of K and A, corresponding to lower and upper limits of [nu] for any particular concentration of EtBr. Using these sets and solving equation 1 numerically we have estimated the effective superhelix density ([sigma] + [alpha][nu]) of both topoisomers of the molecules with mobile and intruded immobile junctions in the range of the transition (Table 1). We conclude from these data that the transition corresponds to a positive superhelix density in the range 0.04-0.1.
Calculated effective superhelix density of the minicircles at different concentrations of EtBr
[EtBr] (µM)
Effective superhelix density ([sigma] + [alpha][nu])
Lk = 25
Lk = 26
Lower limit
Higher limit
Lower limit
Higher limit
0.0
-0.020
-0.020
0.190
0.190
0.1
0.012
0.061
0.030
0.076
0.5
0.030
0.096
0.043
0.109
1.5
0.043
0.121
0.055
0.133
5.0
0.058
0.148
0.068
0.160
The lower and higher limits of [sigma] + [alpha][nu] for two topoisomers of the molecules with mobile and intruded immobile junctions were obtained by numerical solution of equation 1. The lower limit values were obtained for K = 106/M and A = 40, the higher limit for K = 107/M and A = 20. Superhelix density values in bold in the table correspond to EtBr concentrations in the range of the transition.
The electrophoretic mobility data in Figure 3 have enabled us to directly estimate interwinding of the opposite DNA strands in the post-transition structures. We found that the mobilities of all molecules, including perfectly complementary circles, depend only on their linking numbers, not on their sequences. This finding implies that the twist in the region where the sequences differ is close to that of B-DNA, regardless of the extent of its complementarity. More data are required to form conclusions on the extent of regularity in the interwound structure formed by non-complementary base pairs under positive torsional stress.
This research has been supported by grants GM-54215 (A.V.V.) and GM-29554 (N.C.S.) from the National Institute of General Medical Sciences, N000149810093 from the Office of Naval Research (N.C.S.), NSF-CCR-97-25021 from the National Science Foundation (N.C.S.) and by a Margaret and Herman Sokol Fellowship (X.Y.). We are grateful to Dr J.B.Chaires for valuable discussions about this work.
1. Gellert,M., Mizuuchi,K., O'Dea,M.H., Ohmori,H. and Tomizawa,J. (1979) Cold Spring Harbor Symp. Quant. Biol., 43, 35-40.MEDLINE Abstract
2. Lilley,D.M.J. (1980) Proc. Natl. Acad. Sci. USA, 77, 6468-6472.
3. Panayotatos,N. and Wells,R.D. (1981) Nature, 289, 466-470.MEDLINE Abstract
4. Wang,J.C. (1986) In Semlyen,J.A. (ed.), Cyclic Polymers. Elsevier, Essex, UK, pp. 225-260.
5. Wells,R.D. (1988) J. Biol. Chem., 263, 1095-1098.MEDLINE Abstract
6. Frank-Kamenetskii,M.D. (1990) In Cozzarelli,N.R. and Wang,J.C. (eds), DNA Topology and its Biological Effects. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY, pp. 185-216.
7. Murchie,A.I.H. and Lilley,D.M.J. (1992) Methods Enzymol., 211, 158-180.