ABSTRACT
We report here kinetic and thermodynamic studies on differential isomerization of intramolecular Pyr*Pur·Pyr triplexes in supercoiled plasmids. Two structural isomers of the triplex exist: one with the 3'-half of the Pyr strand as the third strand (H-y3 form) and the other with the 5'
-half as the third strand (H-y5 form). The relative populations of the two triplex isomers was determined using the chemical probe with diethyl pryrocarbonate as a function of incubation time. The results demonstrated that triplexes were formed rapidly after a pH change from pH 8.0 to 5.0 and that the initial population of the two isomers exponentially changed with incubation time to reach true thermodynamic equilibrium with a time constant of 0.6-10 h, depending on temperature and the presence of Mg2+. The results clearly demonstrated that interconversion occurs between the two isomers and that the presence of Mg2+ generally retarded the interconversion rates. Kinetic and thermodynamic analyses of the relative populations of the two isomers revealed that the apparent energy barrier for transition from duplex to the H-y3 form is higher than that to the H-y5 form, but H-y3 is more stable in enthalpy terms than H-y5. Therefore, H-y3 is kinetically inferior but thermodynamically favored at higher supercoil levels in plasmids. The presence of Mg2+ resulted in both a kinetic and a thermodynamic preference for H-y5 formation, relative to the H-y3 form.
A polypurine-polypyrimidine sequence with mirror repeat can adopt an intramolecular triplex structure (also called H-DNA) which consists of a triple-stranded stem in one half of the mirror sequence and a single purine strand (see 1 -3 for reviews). A major type of the triplex structure is composed of canonical base triplets via Hoogsteen hydrogen bonding, C*G·C and T*A·T, forming a Pyr*Pur·Pyr triplex (Pyr* strand as the third strand). Two structural isomers exist for an intramolecular Pyr*Pur·Pyr triplex (see 1 ,3 ,4 for reviews): one with the 3'-half of the Pyr strand as the third strand (H-y3 isomer) and the other with the 5'-half as the third strand (H-y5 isomer). Formation of H-y3 was greatly preferred over H-y5 in most of the sequences reported (5 -11 ).
Htun and Dahlberg (8 ) first demonstrated differential formation of the H-y3 and H-y5 isomers depending on the supercoil density, which was taken into account by an inequivalent topology of two isomers in the nucleation step of triplex formation. A model proposed by them could explain supercoil-dependent behavior of differentiation of triplex isomers in terms of the sense of pre-nucleation rotation at the center of Pur·Pyr tracts, but it was not sufficient to understand the other factors affecting triplex isomerization. Parniewski et al. (12 ) showed formation of a mixture of the two isomers for the sequence (AG)7ATCGATATATSTCG(AG)7. Kang et al. (13 ,14 ) demonstrated that the presence of divalent cations is one of the crucial factors in formation of H-y5 for the sequence (GAA)4TTCGC(GAA)4. We showed that a higher G+C content in the central region of Pur·Pyr tracts and divalent metal ions both caused a preference for the H-y5 isomer (15 ). All these results together led us to conclude that the size of the denaturation bubble or the melted region at the center of a Pur·Pyr tract is critical in determining differential preference for triplex isomers (15 ), i.e. larger opening of the base pairs is required to wrap the pyrimidine strand into the H-y3 form than into the H-y5 form. A recent kinetic study with oligonucleotide intramolecular triplexes by Roberts and Crothers (16 ) supported the above notions and, moreover, they explained the `isomer paradox of H-form DNA, i.e. a bias toward H-y3 formation' in terms of pre-nucleation geometry and accessibility of the major groove.
Htun and Dahlberg (8 ) suggested that the two isomers are not interconvertible but that their populations would be a kinetic consequence of the topological difference. Our previous study (15 ) also emphasized that the rate limiting step in intramolecular formation involves partial melting of the duplex near the center of Pur·Pyr tracts and that a higher G+C content at the center requires a higher supercoil density for triplex formation. On the other hand, because more negative supercoils are relaxed when H-y3 is formed (8 ,10 ,15 ), this isomer is thermodynamically more stable than H-y5 under negative supercoil stress. As pointed out by Glover et al. (10 ) and Shimizu et al. (15 ), the folded back loop in the pyrimidine strand of the H-y5 form is less accessible to solvent than that of H-y3, suggesting some differences in the base stacking or base pairing mode in the two isomers. Thus differential isomerization may be understood as a combination of thermodynamic stability and the kinetic terms in triplex formation.
We now ask whether the populations of the two isomers are determined kinetically or energetically and if interconversion could occur. To this end we determined the kinetic and thermodynamic parameters on relative preference and lifetime for the two isomers by chemical probing. We have shown that the two isomers are indeed interconvertible, with characteristic rates depending on temperature and concentration of Mg2+ cation, and that true thermal equilibrium can be attained within 1-20 h at room temperature or higher.
Construction of plasmid pMS215 was described previously (15 ). Plasmids were grown in Escherichia coli JM109 and prepared using Magic Max-prep (Promega).
Modification of plasmids by diethyl pyrocarbonate (DEPC; Sigma) was performed as described (9 ,17 ,18 ). Briefly, 3 µg plasmid in TE buffer (10 mM Tris-HCl, 0.1 mM EDTA, pH 8.0) were diluted into 100 µl 0.3 M Tris-acetate buffer, pH 5.0, in the presence or absence of 50 mM MgCl2 and incubated for the given periods of time. The reaction was initiated by addition of 10 µl DEPC and incubated with vortexing at the given temperature. After 30 min the reaction was stopped by precipitation of the DNA with ethanol. The DNA was digested with HindIII and SacI endonucleases and the purine strand of the HindIII-SacI fragments was radioactively labeled at the HindIII site by filling in the 5'-overhang with the Klenow fragment of DNA polymerase (New England Biolabs) in the presence of [[alpha]-32P]dATP (~3000 mCi/mmol; ICN), dGTP, dCTP and dTTP. The radioactively labeled inserts were isolated by gel electrophoresis and treated with piperidine at 90°C for 20 min. A sequencing gel was run and the DEPC modification sites were visualized by autoradiography.
We have previously shown that the sequence (GA)7TGGC(AG)7 in the supercoiled plasmid pMS215 used here predominantly forms the H-y3 isomer in the absence of Mg2+, while the H-y5 isomer is dominant in the presence of 50 mM Mg2+ (15 ). Therefore, this plasmid will provide a suitable example for both kinetic and thermodynamic studies of triplex isomerization. We used chemical probing with DEPC to quantify the relative populations of the two isomers in pMS215, since it preferentially reacts with purine residues (A > G) in the single-strand region of the 5'- and 3'-half of the purine strand in H-y3 and H-y5 respectively, but is much less reactive with the duplex form.
To elucidate dynamic isomerization of intramolecular triplex formation, fine mapping of chemical modification sites by DEPC in supercoiled plasmid pMS215 was performed as a function of incubation time under three different experimental procedures, as shown in Figure 1 . Supercoiled plasmid pMS215 dissolved in TE buffer (10 mM Tris-HCl, pH 8.0, 0.1 mM EDTA) was added to 0.3 M Tris-acetate buffer, pH 5.0, without Mg2+ at the given temperature and after appropriate time intervals chemical probing with DEPC was carried out for 30 min. This is denoted experiment A. In the second experiment, B, the supercoiled plasmid (pH 8.0) was added to 0.3 M Tris-acetate buffer, pH 5.0, in the presence of 50 mM Mg2+, followed by the same procedures as in experiment A. In the third experiment, C, the supercoiled plasmid was added to 0.3 M Tris-acetate buffer, pH 5.0, in the absence of Mg2+ and incubated for 30 min. After addition of MgCl2 solution to 50 mM, chemical probing was carried out at the given intervals of incubation time.
The mechanism of intramolecular triplex formation includes a denaturation bubble at the dyad axis of the Pur·Pyr tract and nucleation, followed by elongation of the triplex stem with strand separation of the duplex (8 ,10 ,15 ). In this work we have studied the kinetic and thermodynamic properties of isomerization in intramolecular triplex formation and demonstrated that interconversion between the two isomers occurs.
The time dependence of the relative populations of the two triplex isomers may arise from two causes. One of the isomers may be generated as a result of interconversion from the other or it could be newly generated from the duplex form remaining unchanged in the plasmid at the initial time. The latter case is unlikely for the following reasons. Transition probabilities for transition from the duplex to the two triplex isomers are considered to be constant, independent of time under a given experimental condition. Assuming that the isomers are generated from the duplex fraction and that no conversion between two isomers occurs, the relative populations of the two isomers should be time independent according to kinetic theory for a first order competitive reaction. However, this was not the case, as can be seen in Figure 3 a and b. Furthermore, as mentioned in Results, it seems clear that triplex formation in supercoiled plasmids was rapid, being complete within ~10 min (19 ,20 ). Hanvey et al. (20 ) demonstrated that triplexes were formed within a few minutes at 25°C after a pH jump from pH 8.0 to 5.0, as judged from S1 nuclease sensitivity of the Pur·Pyr insert in supercoiled plasmids. Thus all these results verify that the duplex is short lived and that interconversion indeed occurs between the two isomers.
Interconversion could reasonably be assumed to occur only through the duplex as a short-lived intermediate because the process will take the energetic minimum pathway. Unwrapping of the pyrimidine strand from the triplex stem initiated at the duplex-triplex junction at the 3'- or 5'-end of the pyrimidine strand would be followed by duplex formation with the complementary single-stranded purine loop in a concerted manner. Thus the free energy cost due to an increase in negative supercoiling in the unwrapping process will be compensated for by the stacking energy due to duplex formation with the single-stranded purine loop. The rate of conversion from H-y5 to H-y3 is faster than that of the reverse conversion (see the last two columns in Table 1 ), which may be accounted for by the asymmetric nature of the energy barriers for transition from duplex to two triplex isomers, as mentioned below.
The barrier for transition from duplex to the H-y5 form is significantly lower by 2.6 kcal/mol in the absence of Mg2+ and 3.7 kcal/mol in the presence of Mg2+ than the barrier for H-y3, consistent with our notion (15 ) that H-y3 formation requires a wider region of base pair opening than does H-y5. Such a disadvantage in H-y3 formation would be overcome by excess free energy due to the difference in relaxation of supercoils between H-y3 and H-y5 in the transition states. Therefore, H-y3 has a low transition probability from the duplex relative to H-y5 at lower supercoil densities, while it is predominant at high levels of supercoiling.
Three main factors affecting the relative populations of the two isomers can be considered at the equilibrium state (15 ): negative supercoil density, G+C content in the central region of the Pur·Pyr tract and divalent cations. Free energy change at equilibrium for conversion from H-y5 to H-y3 may contain the following supercoil-dependent and structural energy term
Here [Delta]Gsc° is the free energy change associated with supercoil density at equilibrium and [Delta]Gst° is a structural term that may be related to the hairpin loop in the central region of the Pur·Pyr tract of the triplex structure. The supercoil-dependent term is given as follows (21 -23 ),
Here, [Delta]L3r and [Delta]L5r are the amount of relaxation of supercoils due to H-y3 and H-y5 formation respectively, R is the gas constant, T is absolute temperature and n is the number of base pairs in the plasmid used (n = 2686 for pMS215).
Knowing that the amount of maximum relaxation of supercoil turns produced by triplex formation in pMS215 was [Delta]L3r = 3.9 ± 0.1 and [Delta]L5r = 3.0 ± 0.1 for the H-y3 and H-y5 isomers respectively (15 ), [Delta]Gsc° was calculated to be 4.74 kcal/mol from equation 5, using the fact that the supercoil density of the plasmid was [sigma] = -0.055 (corresponding to [Delta]Lk = -14). Substituting the observed value, -0.58 kcal/mol (Table 2 ), for free energy difference [Delta]G° in equation 4 gave an estimated free energy change, [Delta]Gst°, of ~4.16 kcal/mol at 35°C. This positive value of [Delta]Gst° indicates that the H-y5 form is thermodynamically more stable than H-y3 at low supercoil stress.
Thermodynamic data on oligonucleotide intramolecular triplexes are abundant. In contrast to the result obtained above, most of the data available indicate that the two triplex isomers have similar thermal stability (16 ,24 -26 ). There is some evidence, however, showing high stability of H-y5 relative to the H-y3 form of intramolecular triplex in supercoiled plasmids. First, the H-y5 form is dominant at lower supercoil densities (8 ,27 ) and it was formed even in the linearized plasmid at pH 4.0 (10 ). Second, as pointed out by Glover et al. (10 ) and Shimizu et al. (15 ), the hairpin loop region in the H-y5 form is much less accessible to solvent than that of H-y3, indicating a relatively ordered structure for the H-y5 loop. Such a relatively high reactivity of the H-y3 form can be seen in the autoradiograms shown in Figure 4 b: two guanines of the loop region are insensitive at t = 0, while they become sensitive at 12 h, when the H-y3 form is dominant. Less accessibility of the hairpin loop to solvent may arise from an extra Watson-Crick base pairing of the residue at the 3'-end of the folded back loop of the homopyrimidine strand: H-y5 would form such a Watson-Crick base pairing, while H-y3 would not (10 ). If this were the case, the residues near the 3'-end of the four-member loop in the center of the pyrimidine strand might play a key role in stabilizing the hairpin loop of the H-y5 form.
In conclusion, kinetic and thermodynamic studies revealed that a differential preference among intramolecular triplex isomers is initially determined by inequivalent kinetic barriers for the transition from the duplex to the triplex isomers and their relative populations change exponentially toward an increase in the H-y3 isomer to reach equilibrium. Such inequivalency of the kinetic barriers arises from differences in topology and size of a denaturation bubble in the transition states for triplex formation. The H-y3 form is thermodynamically less stable than H-y5 at lower supercoil level, which may arise from instability of the hairpin loop in H-y3 relative to that in H-y5. However, its disadvantage in triplex formation is overcome at higher negative supercoil level, because H-y3 formation results in relaxation of negative supercoils by one more turn than does H-y5. These will be general features for isomerization of intramolecular triplexes in supercoiled plasmids.
We thank Professor U.Matsumoto for his constant interest and support. This work was supported by a Grant-in Aid for Scientific Research from the Ministry of Education, Science, Sports and Culture of Japan to H.S. and M.S.
[Delta]G° = [Delta]Gsc° + [Delta]Gst°
4
[Delta]Gsc = 1100 RTn [([Delta]Lk + [Delta]L3r)2 - ([Delta]Lk + [Delta]L5r)2]
5
REFERENCES
![]()
CiteULike
Connotea
Del.icio.us What's this?
This Article ![]()
![]()
Abstract
![]()
Print PDF (114K)
![]()
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 (2)
![]()
Commercial Re-use Guidelines
for Open Access NAR Content
![]()
Google Scholar ![]()
![]()
Articles by Shindo, H.
![]()
Articles by Shimizu, M.
![]()
Search for Related Content
![]()
PubMed ![]()
![]()
PubMed Citation
![]()
Articles by Shindo, H.
![]()
Articles by Shimizu, M.
![]()
Social Bookmarking ![]()
![]()
What's this?
