Studies of DNA dumbbells. III. Theoretical analysis of optical melting curves of dumbbells with a 16 base-pair duplex stem and Tn end loops (n = 2, 3, 4, 6, 8, 10, 14).

Optical melting curves of seven DNA dumbbells with the 16 base-pair duplex sequence 5'G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C3' linked on both ends by Tn (n = 2, 3, 4, 6, 8, 10, and 14) loops measured in 30, 70, and 120 mM Na+ are analyzed in terms of the numerically exact statistical thermodynamic model of DNA melting. The construction and characterization of these molecules were described in the previous paper (Amaratunga et al., 1992). As was recently reported for hairpins (T. M. Paner, M. Amaratunga, M. J. Doktycz, and A. S. Benight, 1990, Biopolymers, Vol. 29, pp. 1715-1734) theoretically calculated melting curves were fitted to experimental curves by simultaneously adjusting the parameters representing loop and circle formation to optimize the fits. The systematically determined empirical parameters provide evaluations of the free energies of hairpin loop formation delta Gloop (n) and single-strand circles delta Gcircle (N), as a function of end loop size, n = 2-14, and circle size, N = 32 + 2n. The dependence of these quantities on solvent ionic strength over the range from 30 to 120 mM Na+ was evaluated. An approximately analytical expression for the partition function Q(T) of the dumbbells was formulated that allowed a means for determining the transition enthalpy delta H degrees and entropy delta S degrees for every dumbbell, revealing the dependence of these quantities on loop size. In this multistate approach a manifold of partially melted intermediate microstates are considered and therefore no assumptions regarding the nature of the melting transitions (that they are two-state) are required. The transition thermodynamic parameters were also determined from a van't Hoff analysis of the melting curves. Comparisons between the results of the multistate analysis and the two-state van't Hoff analysis revealed significant differences for the dumbbells with larger end loops, indicating that the melting transitions of the larger looped dumbbells deviate considerably from two-state behavior. Results are then compared with published melting studies of much larger DNA dumbbells (D. B. Naritsin and Y. L. Lyubchenko, 1990, Journal of Biomolecular Structure and Dynamics, Vol. 8, pp. 1-13), of small hairpins (Paner et al., 1990; M. J. Doktycz, T. M. Paner, M. Amaratunga and A. S. Benight, 1990, Biopolymers, Vol. 30, pp. 829-845) and another dumbbell (A. S. Benight, J. M. Schurr, P. F. Flynn, B. R. Reid, and D. E. Wemmer, 1988) Journal of Molecular Biology, Vol. 200, pp. 377-399).(ABSTRACT TRUNCATED AT 400 WORDS)     

Analysis of melting transitions of the DNA hairpins formed from the oligomer sequences d[GGATAC(X)4GTATCC] (X = A, T, G, C)

Optical melting transitions of the short DNA hairpins formed from the self-complementary DNA oligomers d[GGATACX4GTATCC] where X = A, T, G, or C measured in 100 mM NaCl are presented. A significant dependence of the melting transitions on loop sequence is observed and transition temperatures, tm, of the hairpins vary from 58.3 degrees C for the T4 loop hairpin to 55.3 degrees C for the A4 loop. A nearest-neighbor sequence-dependent theoretical algorithm for calculating melting curves of DNA hairpins is presented and employed to analyze the experimental melting transitions. Experimental melting curves were fit by adjustment of a single theoretical parameter, Fend(n), the weighting function for a hairpin loop comprised of n single-strand bases. Empirically determined values of Fend(n) provide an evaluation of the free-energy of hairpin loop formation and stability. Effects of heterogeneous nearest-neighbor sequence interactions in the duplex stem on hairpin loop formation were investigated by evaluating Fend(n) in individual fitting procedures using two of the published sets of nearest-neighbor stacking interactions in DNA evaluated in 100 mM NaCl and given by Wartell and Benight, 1985. In all cases, evaluated values of Fend(n) were obtained that provided exact theoretical predictions of the experimental transitions. Results of the evaluations indicate: (1) Evaluated free-energies of hairpin loop formation are only slightly dependent on loop sequences examined. At the transition temperature, Tm, the free-energy of forming a loop of four bases is approximately equal for T4, G4, or C4 loops and varies from 3.9 to 4.8 kcal/mole depending on the set of nearest-neighbor interactions employed in the evaluations. This result suggests, in light of the observed differences in stability between the T4, G4, and C4 loop hairpins, that sequence-dependent interactions between base residues of the loop are most likely not the source of the enhanced stability of a T4 loop.(ABSTRACT TRUNCATED AT 400 WORDS)     

Studies of DNA dumbbells. VI. Analysis of optical melting curves of dumbbells with a sixteen-base pair duplex stem and end-loops of variable size and sequence

Optical melting curves of 22 DNA dumbbells with the 16-base pair duplex sequence 5′-G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C-3′ linked on both ends by single-strand loops of A_t or C_t sequences (˛ = 2, 3, 4, 6, 8, 10, 14), T_t sequences (˛ = 2, 3, 4, 6, 8, 10), and G_t sequences (t = 2, 4) were measured in phosphate buffered solvents containing 30, 70, and 120 mM Na⁺. For dumbbells with loops comprised of at least three nucleotides, stability is inversely proportional to end-loop size. Dumbbells with loops comprised of only two nucleotide bases generally have lower stabilities than dumbbells with three base nucleotide loops. Experimental melting curves were analyzed in terms of the numerically exact (multistate) statistical thermodynamic model of DNA dumbbell melting previously described (T. M. Paner, M. Amaratunga & A. S. Benight (1992), Biopolymers 32, 881). Theoretically calculated melting curves were fitted to experimental curves by simultaneously adjusting model parameters representing statistical weights of intramolecular hairpin loop and single-strand circle states. The systematically determined empirical parameters provided evaluations of the energetics of hairpin loop formation as a function of loop size, sequence, and salt environment. Values of the free energies of hairpin loop formation ΔG_loop(n > t) and single-strand circles, ΔG_cir(N) as a function of end-loop size, t = 2-14, circle size, N = 32 + 2_t, and loop sequence were obtained. These quantities were found to depend on end-loop size but not loop sequence. Their empirically determined values also varied with solvent ionic strength. Analytical expressions for the partition function Q(T) of the dumbbells were evaluated using the empirically evaluated best-fit loop parameters. From Q(T), the melting transition enthalpy ΔH, entropy ΔS, and free energy ΔG, were evaluated for the dumbbells as a function of end-loop size, sequence, and [Na⁺]. Since the multistate analysis is based on the numerically exact model, and considers a statistically significant number of theoretically possible partially melted states, it does not require prior assumptions regarding the nature of the melting transition, i.e., whether or not it occurs in a two-state manner. For comparison with the multistate analysis, thermodynamic transition parameters were also evaluated directly from experimental melting curves assuming a two-state transition and using the graphical van't Hoff analysis. Comparisons between results of the multistate and two-state analyses suggested dumbbells with loops comprised of six or fewer residues melted in a two-state manner, while the melting processes for dumbbells with larger end-loops deviate from two-state behavior.Dependence of thermodynamic parameters on[Na⁺] as a function of loop size suggests single-strand end-loops have different counterion binding properties than the melted circle. Results are compared with those obtained in an earlier study of dumbbells with the slightly different stem sequence 5'-G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C-3' linked on the ends by T loops (˛ = 2,3,4,6,8,10,14).© 1996 John Wiley &Sons, Inc.     

Thermal stability, structural features, and B-to-Z transition in DNA tetraloop hairpins as determined by optical spectroscopy in d(CG)(3)T(4)(CG)(3) and d(CG)(3)A(4)(CG)(3) oligodeoxynucleotides

NMR and CD data have previously shown the formation of the T(4) tetraloop hairpin in aqueous solutions, as well as the possibility of the B-to-Z transition in its stem in high salt concentration conditions. It has been shown that the stem B-to-Z transition in T(4) hairpins leads to S (south)- to N (north)-type conformational changes in the loop sugars, as well as anti to syn orientations in the loop bases. In this article, we have compared by means of UV absorption, CD, Raman, and Fourier transform infrared (FTIR), the thermodynamic and structural properties of the T(4) and A(4) tetraloop hairpins formed in 5'-d(CGCGCG-TTTT-CGCGCG)-3' and 5'-d(CGCGCG-AAAA-CGCGCG)-3', respectively. In presence of 5M NaClO(4), a complete B-to-Z transition of the stems is first proved by CD spectra. UV melting profiles are consistent with a higher thermal stability of the T(4) hairpin compared to the A(4) hairpin. Order-to-disorder transition of both hairpins has also been analyzed by means of Raman spectra recorded as a function of temperature. A clear Z-to-B transition of the stem has been confirmed in the T(4) hairpin, and not in the A(4) hairpin. With a right-handed stem, Raman and FTIR spectra have confirmed the C2'-endo/anti conformation for all the T(4) loop nucleosides. With a left-handed stem, a part of the T(4) loop sugars adopt a N-type (C3'-endo) conformation, and the C3'-endo/syn conformation seems to be the preferred one for the dA residues involved in the A(4) tetraloop.     

B to Z transitions of the short DNA hairpins formed from the oligomer sequences: d[(CG)3X4(CG)3] (X = A, T, G, C).

Circular Dichroism (CD) spectra were collected as a function of sodium perchlorate concentration [NaClO4] for the set of DNA hairpins formed from the oligomer sequences d[(CG)3X4(CG)3] where X = A, T, G or C. Over the range in salt concentration from 0 to 4.0 M NaClO4, the CD spectra invert in a manner characteristic of the B to Z transition. A factor analysis routine is described and employed to determine the least number of basis spectra required to fit the measured spectra of each hairpin over the entire salt range examined. In every case, linear combinations of only two sub-spectra fit the experimental spectra of the hairpins with greater than 98% accuracy, indicating the spectrally monitored structural transitions are two-state. From the relative weights of the individual sub-spectra, B-Z transition curves are constructed. The transitions are analyzed in terms of a simple two-state equilibrium model which yields an evaluation of the transition free-energy, delta GB-Z, as a function of NaClO4 concentration. At 1.0 M NaClO4 and 21 degrees C, delta GB-Z = 5.4, 4.9, 3.6 and 2.3 kcal/mole for the G4, T4, A4 and C4 loop hairpins, respectively.     

Extra stable 2',5'-linked RNA loops

We prepared hairpins that differ in the connectivity of phosphodiester linkages in the loop (RNA vs 2', 5'-RNA). We find that the stability of the extra stable RNA hairpin 5'-rGGAC(UUCG)GUCC-3' is the same as that observed for the hairpin containing a 2',5'RNA loop, i.e. 5'-rGGAC(UUCG)GUCC-3' (where UUCG = U2'p5'U2'p5' C2'p5'G2'p5'). Also significant is the finding that when the stem is duplex DNA, duplex 2',5'-RNA, or DNA:2',5'-RNA, hairpins with the UUCG loop are more stable than those with UUCG loop.     

Right- and left-handed (Z) helical conformations of the hairpin d(C-G)5T4(C-G)5 monomer and dimer

The partial self-complementary 24-mer oligodeoxynucleotide d(C-G)5T4(C-G)5 forms a hairpin which can be enzymatically dimerized to a dumbbell structure. The blunt-ended nature of the hairpin is indicated by its ability to inhibit the T4 DNA ligase catalyzed joining of phi X174 HaeIII fragments. The hairpin monomer and dimer (dumbbell) undergo a reversible B to Z transition as shown by ultraviolet, circular dichroism, and 31P NMR spectroscopy. The Z form of the hairpin monomer and dimer is supported by monovalent ions (Na+), divalent ions (Mg2+ but not Mn2+), and dehydrating (ethanol) conditions. The conformational transition of d(C-G)5T4(C-G)5 monomer requires higher ionic or dehydrating conditions than those necessary for the corresponding linear oligomer d(C-G)5. The contribution of the loop (-T4-) of the hairpin to the apparent free energy change for the B to Z conformational transition at the midpoint was calculated to be 3.8 kJ mol-1.     

Characterization of the structure and melting of DNAs containing backbone nicks and gaps

A DNA molecule containing a gap (a missing phosphate) has been examined and compared to two other molecules of the same sequence, one containing a nick (a phosphorylated gap) and the other a normal duplex containing no break in the backbone. A second gapped sequence was also compared to a normal duplex of the same sequence. The molecules containing nicks or gaps were generated as dumbbell molecules, short helices closed by a loop at each end. The dumbbells were formed by the association of two hairpins with self-complementary dangling 5'-ends. Nuclear magnetic resonance was used to monitor the melting transition and to probe structural differences between molecules. Under the conditions used here no change in stability was observed upon phosphorylation of the gap. Structural changes upon phosphorylation of a gap or closure of a nick were minimal and were localized to the region immediately around the gap or nick. Two transitions can be observed as a gapped or nicked molecule melts, although the resolution of the two transitions varies with the salt concentration. At moderate to high salt (greater than or equal to 30 mM) the molecule melts essentially all at once. At low salt the two transitions occur at temperatures that differ by as much as 15 degrees C. In addition, comparison with other NMR melting studies indicates that the duplex formed by the overlap of the dangling ends of the hairpins is stabilized relative to a free duplex of the same sequence, probably by stacking onto the hairpin stem.     

Stability of 3' double nucleotide overhangs that model the 3' ends of siRNA

Thermodynamic parameters are reported for duplex formation in 1 M NaCl for 16 RNA sequences, each containing a core tetramer duplex, GGCC, and a 3' overhang consisting of two bases. The results indicate additional double-helical stability is conferred by the double 3' terminal overhang relative to the single 3' terminal overhang. A nearest-neighbor analysis of the data indicates that the free energy contribution at 37 degrees C of the second base in the double 3' terminal overhang varies from 0 to 0.7 kcal/mol. The second base in the 3' double overhang can contribute nearly the same stability to a duplex as a base pair or a 3' dangling overhang. Stability contribution of a dangling base, two nucleotides removed from the 3' end of a duplex, is dependent upon both the identity of the base as well as that of the dangling base that it neighbors. A second dangling base only increases the stability of the duplex when it is neighboring a 3' purine dangling nucleotide. Furthermore, a second dangling pyrimidine provides a greater contribution to duplex stability than a purine. A nearest-neighbor model was developed to predict the influence of 3' double overhang on the stability of duplex formation. The model improves the prediction of free energy and melting temperature when tested against six sequences with different core duplexes.     

Stability of single-nucleotide bulge loops embedded in a GAAA RNA hairpin stem

Forty-six RNA hairpins containing combinations of 3' or 5' bulge loops and a 3' or 5' fluorescein label were optically melted in 1 M NaCl, and the thermodynamic parameters ΔH°, ΔS°, ΔG°(37), and T(M) for each hairpin were determined. The bulge loops were of the group I variety, in which the identity of the bulge is known, and the group II variety, in which the bulged nucleotide is identical to one of its nearest neighbors, leading to ambiguity as to the exact position of the bulge. The fluorescein label at either the 3' end or 5' end of the hairpin did not significantly influence the stability of the hairpin. As observed with bulge loops inserted into a duplex motif, the insertion of a bulge loop into the stem of a hairpin loop was destabilizing. The model developed to predict the influence of bulge loops on the stability of duplex formation was extended to predict the influence of bulge loops on hairpin stability. Specifically, the influence of the bulge is related to the stability of the hairpin stem distal from the hairpin loop.     

Studies of DNA dumbbells. IV. Preparation and melting of a DNA dumbbell with the 16 base-pair sequence 5′G-T-A-T-C-C-C-T-C-T-G-G-A-T-A-C3′ linked on the ends by dodecyl chains

The preparation and melting of a 16 base-pair duplex DNA linked on both ends by C¹²H²⁴ (dodecyl) chains is described. Absorbance vs temperature curves (optical melting curves) were measured for the dodecyl-linked molecule and the same duplex molecule linked on the ends instead by T₄ loops. Optical melting curves of both molecules were measured in 25, 55, and 85 mM Na⁺ and revealed, regardless of [Na ⁺], the duplex linked by dodecyl loops is more stable by at least 6°C than the same duplex linked by T₄ loops. Experimental curves in each salt environment were analyzed in terms of the two-state and multistate theoretical models. In the two-state, or van't Hoff analysis, the melting transition is assumed to occur in an all-or-none manner. Thus, the only possible states accessible to the molecule throughout the melting transition are the completely intact duplex and the completely melted duplex or minicircle. In the multistate analysis no assumptions regarding the melting transition are required and the statistical occurrence of every possible partially melted state of the duplex is explicitly considered. Results of the analysis revealed the melting transitions of both the dodecyl-linked molecule and the dumbbell with T₄ end loops are essentially two state in 25 and 55 mM Na⁺. In contrast, significant deviations from two-state behavior were observed in 85 m MNa⁺. From our previously published melting data of DNA dumbbells with T_n end loops where n = 2, 3, 4, 6, 8, 10, 14 [T. M. Paner, M. Amaratunga, and A. S. Benight, (1992) Biopolymers, Vol. 32, pp. 881–892] and the dumbbell with T₄ end loops of this study, a plot of d(T_m)/d ln [Na⁺] was constructed. Extrapolation of this data to n = 1 intersects with the value of d (T_m)/d ln [Na⁺] obtained for the alkyl-linked dumbbell, suggesting the salt-dependent stability of the alkyl-linked molecule behaves as though the duplex of this molecule were linked by end loops comprised of a single T residue. © 1993 John Wiley & Sons, Inc.     

Melting of Oligodeoxynucleotides with Various Structures

Abstract

Effects of DNA fragments end structures on their melting profiles were studied experimentally and theoretically. We examined melting of hairpins and dumbbells obtained from 62- bp-long linear DNA duplex which is a perfect palindromic sequence. To fit theoretical melting profile to experimental ones additional theoretical parameters were incorporated into the standard statistical mechanical helix-coil transition theory. From comparison theoretical and experimental melting profiles theoretical parameters connected with end- structure effects were evaluated. Analysis revealed the stabilization effect of the hairpin loops and helix ends with respect to DNA duplex melting. Both type of ends make melting these oligodeoxynucleotides more cooperative than predicted by the standard helix-coil transition theory. At low ionic strength ([Na⁺] < 0.04 M) this effect becomes so pronounced that melting of the DNA duplexes 30–40 bp-long conforms to the two state model.

From the analysis experimental data obtained for dumbbell structures loop-weighting factor for single-stranded loop consisting of 132 nucleotides was determined. This parameter decreases 10 times with the ionic strength decreasing by an order of magnitude from 0.2 to 0.02 M Na⁺.     

Melting studies of short DNA hairpins containing the universal base 5-nitroindole.

Effects of the universal base 5-nitroindole on the thermodynamic stability of DNA hairpins having a 6 bp stem and four base loops were investigated by optical absorbance and differential scanning calorimetry techniques. Melting studies were conducted in buffer containing 115 mM Na(+). Five different modified versions of DNA hairpins containing a 5-nitroindole base or bases substituted at different positions in the stem and loop regions were examined. Thermo-dynamic parameters of the melting transitions estimated from a two-state analysis of optical melting curves and measured directly by calorimetry revealed that the presence of 5-nitroindole bases in the duplex stem or loop regions of short DNA hairpins significantly affects both their enthalpic and entropic melting components in a compensating manner, while the transition free energy varies linearly with the transition temperature. The calorimetrically determined enthalpy and entropy values of the modified hairpins were considerably smaller (43-53%) than the two-state optical parameters, suggesting that solvent effects may be significant in the melting processes of these hairpins. Results of circular dichroism measurements also revealed slight differences between the modified hairpins and the control in both the duplex and melted states, suggesting subtle structural differences between the control and DNA hairpins containing a 5-nitroindole base or bases.     

Effect of base-pair sequence on the conformations and thermally induced transitions in oligodeoxyribonucleotides containing only AT base pairs

Tm curves, CD spectra, and kinetics results of the self-complementary DNA dodecamers d(A6T6), d(A3T3A3T3), d(A2T2A2T2A2T2), d(ATATATATATAT), and d(T6A6) demonstrate that the thermal transitions of these oligomers at low salt concentration involve a hairpin intermediate. At high salt concentrations (greater than 0.1 M Na+) only a duplex to denatured-strand transition appears to occur. The temperature and salt-concentration regions of the transitions are very sequence dependent. Alternating-type AT sequences have a lower duplex stability and a greater tendency to form hairpins than sequences containing more nonalternating AT base pairs. Of the two nonalternating sequences, d(T6A6) is significantly less stable than d(A6T6). Both oligomers have CD curves that are very similar to the unusual CD spectrum of poly(dA).poly(dT). The Raman spectra of these two oligomers are also quite similar, but at low temperature, small intensity differences in two backbone modes and three nucleoside vibrations are obtained. The hairpin to duplex transition for the AT dodecamers was examined by salt-jump kinetics measurements. The transition is faster than transitions for palindromic-sequence oligomers containing terminal GC base pairs. Stopped-flow kinetics studies indicate that the transition is second order and has a relatively low activation energy. The reaction rate increases with increasing ionic strength. These results are consistent with a three-step mechanism for the hairpin to duplex reaction: (i) fraying of the hairpin oligomers' terminal base pairs, (ii) a rate-determining bimolecular step involving formation of a cruciform-type intermediate from two hairpin oligomers with open terminal base pairs, and (iii) base-pair migration and formation in the intermediate to give the duplex.     

Influence of dangling thymidine residues on the stability and structure of two DNA duplexes

We have employed temperature-dependent UV spectroscopy, circular dichroism (CD), 400-MHz proton nuclear magnetic resonance (NMR), and computer modeling to characterize both structurally and thermodynamically the influence of unpaired, dangling thymidine residues (T) on the thermal stability and melting behavior of two DNA core duplexes. The specific DNA double helices that we have investigated in this work are core duplexes [d(GC)3]2 (I) and [d(CG)3]2 (IV), 3' dangling T derivatives [d[(GC)3TT]]2 (II) and [d[(CG)3TT]]2 (V), and 5' dangling T derivatives [d[TT(GC)3]]2 (III) and [d[TT(CG)3]]2 (VI). Our experimental data allow us to reach the following conclusions: (1) For both core duplexes (I and IV), the addition of dangling T residues on either the 5' or 3' end causes an increase in the optical melting temperature tm. (2) For both core duplexes, 5' dangling T residues induce a greater increase in the optical tm's than 3' dangling T residues. (3) For both cores duplexes, the increase in tm induced by the addition of dangling T residues is enthalpic in origin, with 5' dangling T residues inducing a greater increase in the van't Hoff transition enthalpy than 3' dangling T's. (4) Dangling T residues cause downfield shifts in all of the nonexchangeable aromatic protons of the [d(GC)3]2 core duplex (I), with the 5' T residues inducing the largest shifts. For the most part, this trend does not hold with the [d(CG)3]2 core duplex (IV). (5) For both core duplexes, the addition of dangling T residues causes an increase in the NMR tm's of almost all the nonexchangeable aromatic protons of the core duplex.(ABSTRACT TRUNCATED AT 250 WORDS)     

Hairpin Formation of d(CGCG-TA-CGCG), d(CGCG-TG-CGCG) and Their Cytosine Methylated Analogs

Abstract

Hairpin formations of decamers d(CGCG-TA-CGCG), d(CGCG-TG-CGCG), and their m⁵dC analogs are evidenced by the existence of biphasic absorbance melting profiles in which the lower transition temperature increases with increasing oligomer concentration, whereas the higher melting temperature is concentration independent. The corresponding temperature dependent CD intensity at 285 nm exhibits a maximum around 55°C. These observations are consistent with the interpretation that the lower temperature transition corresponds to the duplex to hairpin transformation while the melting of hairpins into single strands constitutes the higher temperature transition. The CD spectrum of the hairpin conformation appears to be characterized by a couplet with nearly equal positive and negative intensities at 285 and 255 nm, respectively, while a significantly smaller intensity at 285 nm is apparent for the duplex form. The hairpin conformation is suspected to contain a two-nucleotide loop. Titrations with NaCl further suggest that, in contrast to the TA sequence, the TG sequence with wobble base pairing favors Z formation under high salt conditions.     

Energetics of the hairpin to mismatched duplex transition of d(GCCGCAGC) on NaCl solution

We report Potential of Mean Force studies to describe the relative thermodynamic stabilities of d(GCCGCAGC) in a mismatched duplex and a hairpin monomer conformation in NaCl solution. The PMF calculations are combined with previous molecular mechanics and normal mode analysis in order to estimate the role of different components of the free energy in determining the relative stability of the duplex and hairpin structures. The high entropy associated with the loop region and the lack of minor groove phosphate-phosphate interactions in the hairpin compete against the gain in enthalpic contribution to the free energy due to base pairing in the mismatched duplex. The combined free energy calculations show that the hairpin is the most stable conformation at low salt and that a hairpin to duplex transition takes place at approximately 0.47 M NaCl. In addition, we studied the hairpin to partially stacked single helical conformation equilibrium at low salt. We found a small variation in transition temperature in salt concentration, delta Tm/delta log10(cs) approximately 2-3 degrees K/decade, in contrast to the duplex to hairpin or duplex to partially stacked single helix transition where the transition temperature exhibited marked dependence on salt concentration. This is in qualitative agreement with experimental data. Based on the Potential of Mean Force free energy calculation, the order of relative stability of the three-conformations studied varies with salt concentration. We observed the following orders of stability: stacked single helix greater than hairpin greater than duplex for cs less than 0.77 M NaCl; single helix greater than duplex greater than hairpin for 0.77 less than Cs less than 2.1 M; and duplex greater than hairpin greater than single strand for cs greater than 2.1 M. From the calculated PMF free energy curves in the NaCl concentration range, 0.012 less than cs less than 5.0 M, we can assign upper and lower bounds for the non-ionic differences in free energy between the duplex, hairpin, and stacked single helical states (at standard conditions: cs = 1.0 M, T = 25 degrees C, and 1 M oligomer concentration). We found that for delta G duplex single helix = G duplex - 2 x G single helix less than -7.38 Kcal/mol, the single helix is the least stable state. For the duplex-to-hairpin free energy difference in the range, -1.87 less than delta G duplex-hairpin less than 0.03 Kcal/mol, there will always be a salt-induced hairpin-to-duplex transition for 0.01 less than cs less than 1.6 M NaCl. If delta G duplex-hairpin less than -1.87, the duplex is always more stable than the hairpin; and for delta G duplex-hairpin greater than Kcal/mol, the hairpin state is always more stable than the duplex, for all salt concentrations.     

Melting behavior of a covalently closed, single-stranded, circular DNA

We synthesized the 26-residue deoxynucleotide sequence d(TTCCT5GGAATTCCT5GGAA) which folds intramolecularly to form a dumbbell-shaped, double-hairpin structure with a gap between the 3' and the 5' ends. We used T4 polynucleotide kinase to phosphorylate the 5' end followed by T4 DNA ligase to close the 3' and 5' ends. Melting of the dumbbell structure formed by this ligated sequence produces a covalently closed, single-stranded, circular final state. We employed calorimetric and spectroscopic techniques to characterize thermodynamically the melting behavior of the ligated molecule and compared it with the corresponding melting behavior of its unligated precursor. This comparison allowed us to characterize uniquely the influence of single-stranded ring closure on intramolecular duplex melting. The data reveal that ring closure produces a thermally more stable structure which exhibits significantly altered melting thermodynamics. We rationalize these thermodynamic differences in terms of differential solvation and differential counterion association between the ligated and unligated molecules. We also note the importance of such constrained dumbbell structures as models for hairpins, cruciforms, and locally melted domains within naturally occurring DNA polymers.     

A comparison of the hairpin stability of the palindromic d(CGCG(A/T)4CGCG) oligonucleotides.

The palindromic deoxyribonucleotides 5'-CGCGA-TATCGCG-3' and 5'-CGCGTTAACGCG-3' have been characterized by 1H NMR spectroscopy. The NMR data identified both B-DNA duplex conformations and hairpin conformations, the latter with loop regions consisting of the four central nucleotides. The resonances of the various conformations were assigned by use of two-dimensional NMR methods. The relative stability of the various conformations was investigated as a function of temperature, ionic strength and nucleotide concentration. The duplexes were found to be stabilized at high ionic strength and at low temperature, while the hairpins were stabilized at low ionic strength and at medium temperature. The thermodynamics of the duplex-hairpin and the hairpin-random coil transitions were examined, and compared to the other two oligonucleotide in the palindromic d(CGCG(A/T)4CGCG) oligonucleotide family. The relative stabilities of the duplex conformations with respect to the random coil conformations are similar for the d(CGCGAATTCGCG), d(CGCGATATCGCG) and d(CGCGTATACGCG) oligonucleotides. The duplex conformation of d(CGCGTTAACGCG) is less stable. The hairpin of d(CGCGTTAACGCG) seems also to be less stable relative to the random coil conformation than in the case of the other oligonucleotides at an equal oligonucleotide concentration. A cruciform intermediate between the duplex and hairpin conformations is suggested to explain some discrepancies observed in this work in case of the d(CGCGTTAACGCG) oligonucleotide. This is similar to what has been reported for the d(CGCGTATACGCG) oligonucleotide.     

Melting of a self-complementary DNA minicircle. Comparison of optical melting theory with exchange broadening of the nuclear magnetic resonance spectrum

Melting curves are calculated for the 16-base-pair duplex DNA sequence 5' GTATCCGTACGGATAC 3' linked on the ends by TTTT single-strand loops. The equilibrium statistical thermodynamic theory of DNA melting is modified to include effects of end-loops on the melting transition. An excellent fit of the experimental melting curve in 0.2 M-NaCl is obtained using two adjustable parameters, one for end-loop formation and the other for formation of the complete 40-base single-strand loop. The best-fit calculated melting curve permits evaluation of these parameters. The free energy to close a TTTT end-loop is 2.12 kcal/mol (1 cal = 4.184 J). A TTTT end-loop or hairpin loop is significantly more stable than an internal loop of comparable size sandwiched between two helical regions, even after allowing for the different stacking contributions. Reasons for this increased stability are presented. The loop free energy of the 40-base single-strand open minicircle is evaluated to be +1.27 kcal/mol, thus favoring the melting of two end-loops into the large open minicircle. The present results are compared with those of others for d(T-A) oligomers. The sequence TTTT forms a more stable end-loop, or hairpin, than TATA by about 2.0 kcal/mol. Theoretical rate constants for the proton-transfer step in the standard hydrogen-exchange model are calculated by extending the theory of diffusion-controlled reactions to take account of the electrostatic potential of the DNA. The predicted ratios of rate constants for different pairs of catalysts exchanging an A.T proton agree satisfactorily with the available experimental data for a 14-base-pair linear duplex, which confirms the diffusion-control of the proton-transfer step. Data presented here for the 16 base-pair duplex of the minicircle are consistent with catalysis-limited exchange in which the proton-transfer step is likewise diffusion-controlled. Under catalysis-limited conditions, the imino proton exchange rates are predicted from the catalytic rate constants, prevailing buffer catalyst concentrations, and the equilibrium constants to form the unstacked open state of optical melting theory. The observed exchange rates of the A.T base-pairs show no sign of the strong predicted end-melting trend, and exceed the predicted values by factors of 10 to 400. Moreover, the succession of "melting" in the nuclear magnetic resonance line-broadening deviates from that predicted by optical melting theory.(ABSTRACT TRUNCATED AT 400 WORDS)