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.     

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.     

Thermodynamic stability of the 5' dangling-ended DNA hairpins formed from sequences 5'-(XY)2GGATAC(T)4GTATCC-3', where X, Y = A, T, G, C

Expressions for the partition function Q (T) of DNA hairpins are presented. Calculations of Q (T), in conjunction with our previously reported numerically exact algorithm [T. M. Paner, M. Amaratunga, M. J. Doktycz, and A. S. Benight (1990) Biopolymers, 29, 1715-1734], yield a numerical method to evaluate the temperature dependence of the transition enthalpy, entropy, and free energy of a DNA hairpin directly from its optical melting curve. No prior assumptions that the short hairpins melt in a two-state manner are required. This method is then applied in a systematic manner to investigate the stability of the six basepair duplex stem 5'-GGATAC-3' having four-base dangling single-strand ends with the sequences (XY)2, where X, Y = A, T, G, C, on the 5' end and a T4 loop on the 3' end. Results show that all dangling ends of the sample set stabilize the hairpin against melting. Increases in transition temperatures as great as 4.0 degrees C above the blunt-ended control hairpin were observed. The hierarchy of the hairpin transition temperatures is dictated by the identity of the first base of the dangling end adjoining the duplex in the order: purine greater than T greater than C. Calculated melting curves of every hairpin were fit to experimental curves by adjustment of a single parameter in the numerically exact theoretical algorithm. Exact fits were obtained in all cases. Experimental melting curves were also calculated assuming a two-state melting process. Equally accurate fits of all dangling-ended hairpin melting curves were obtained with the two-state model calculation. This was not the case for the melting curve of the blunt-ended hairpin, indicating the presence of a four-base dangling-end drives hairpin melting to a two-state process. Q (T) was calculated as a function of temperature for each hairpin using the theoretical parameters that provided calculated curves in exact agreement with the experimentally obtained optical melting curves. From Q (T), the temperature dependence of the transition enthalpy delta H, entropy delta S, and free energy delta G were calculated for every hairpin providing a quantitative assessment of the effects of dangling ends on hairpin thermodynamics. Comparisons of our results are made with those of the Breslauer group [M. Senior, R. A. Jones, and K. J. Breslauer (1988) Biochemistry 27, 3879-3885] on the T2 5' dangling-ended d(GC)3 duplexes.(ABSTRACT TRUNCATED AT 400 WORDS)     

Studies of DNA dumbbells VII: evaluation of the next-nearest-neighbor sequence-dependent interactions in duplex DNA

Melting experiments were conducted on 22 DNA dumbbells as a function of solvent ionic strength from 25-115 mM Na(+). The dumbbell molecules have short duplex regions comprised of 16-20 base pairs linked on both ends by T(4) single-strand loops. Only the 4-8 central base pairs of the dumbbell stems differ for different molecules, and the six base pairs on both sides of the central sequence and adjoining loops on both ends are the same in every molecule. Results of melting analysis on the 22 new DNA dumbbells are combined with our previous results on 17 other DNA dumbbells, with stem lengths containing from 14-18 base pairs, reported in the first article of this series (Doktycz, Goldstein, Paner, Gallo, and Benight, Biopoly 32, 1992, 849-864). The combination of results comprises a database of optical melting parameters for 39 DNA dumbbells in ionic strengths from 25-115 mM Na(+). This database is employed to evaluate the thermodynamics of singlet, doublet, and triplet sequence-dependent interactions in duplex DNA. Analysis of the 25 mM Na(+) data reveals the existence of significant sequence-dependent triplet or next-nearest-neighbor interactions. The enthalpy of these interactions is evaluated for all possible triplets. Some of the triplet enthalpy values are less than the uncertainty in their evaluation, indicating no measurable interaction for that particular sequence. This finding suggests that the thermodynamic stability of duplex DNA depends on solvent ionic strength in a sequence-dependent manner. As a part of the analysis, the nearest-neighbor (base pair doublet) interactions in 55, 85, and 115 mM Na(+) are also reevaluated from the larger database.     

Studies of DNA dumbbells. II. Construction and characterization of DNA dumbbells with a 16 base-pair duplex stem and Tn end loops (n = 2, 3, 4, 6, 8, 10, 14).

The preparation and characterization of DNA dumbbells that contain the 16 base-pair duplex sequences 5'G-C-A-T-A-G-A-T-G-A-G-A-A-T-G-C3' (set 1) and 5'G-C-A-T-C-A-T-C-G-A-T-G-A-T-G-C3' (set 2) are reported. The dumbbells of set 1 have the duplex stem nucleated on both ends by Tn (n = 2, 3, 4, 6, 8, 10, and 14) loops. The dumbbells of set 2 have Tn (n = 2, 4, 8, 10) end loops. For the molecules of set 1, effects of end loop size on the electrophoretic mobility, CD and UV absorbance spectra, and cleavage by restriction enzymes, were investigated. Effects of loop size on the CD spectra and restriction enzyme cleavage of the molecules of set 2 were also examined. Optical melting curves of the molecules of set 1 were collected as a function of sodium ion concentration from 30 to 120 mM. These investigations revealed that as loop size decreases, the electrophoretic mobilities, rates of enzyme cleavage, and optical melting temperatures increase. For end loops with at least three T's the observed increases are inversely proportional to loop size. The behavior of the dumbbell with T2 end loops departs from this linear dependence and is anomalous in every experimental context. For molecules with end loops comprised of at least four T's CD spectra were virtually indistinguishable. However, these spectra differed considerably from the CD spectrum of the T2-looped molecule. The CD spectrum of the dumbbell with T3 end loops displayed features common to the dumbbells with larger loops and T2 end loops. Thermodynamic evidence that the terminal G.C base pairs (bps) nucleating the T2 end loops were intact was obtained from a comparison of the melting temperature of this molecule with that of a DNA dumbbell containing the 14 central bps of the set 1 duplex sequence linked instead by end loops comprised of the four base sequence, C-T-T-C. The tm of this latter molecule was determined to be 9 degrees C less than that of the former dumbbell assumed to contain a 16-bp stem and T2 end loops.(ABSTRACT TRUNCATED AT 400 WORDS)     

Studies of DNA dumbbells VIII. Melting analysis of DNA dumbbells with dinucleotide repeat stem sequences

Melting curves and circular dichroism spectra were measured for a number of DNA dumbbell and linear molecules containing dinucleotide repeat sequences of different lengths. To study effects of different sequences on the melting and spectroscopic properties, six DNA dumbbells whose stems contain the central sequences (AA)(10), (AC)(10), (AG)(10), (AT)(10), (GC)(10), and (GG)(10) were prepared. These represent the minimal set of 10 possible dinucleotide repeats. To study effects of dinucleotide repeat length, dumbbells with the central sequences (AG)(n), n = 5 and 20, were prepared. Control molecules, dumbbells with a random central sequence, (RN)(n), n = 5, 10, and 20, were also prepared. The central sequence of each dumbbell was flanked on both sides by the same 12 base pairs and T(4) end-loops. Melting curves were measured by optical absorbance and differential scanning calorimetry in solvents containing 25, 55, 85, and 115 mM Na(+). CD spectra were collected from 20 to 45 degrees C and [Na(+)] from 25 to 115 mM. The spectral database did not reveal any apparent temperature dependence in the pretransition region. Analysis of the melting thermodynamics evaluated as a function of Na(+) provided a means for quantitatively estimating the counterion release with melting for the different sequences. Results show a very definite sequence dependence, indicating the salt-dependent properties of duplex DNA are also sequence dependent. Linear DNA molecules containing the (AG)(n) and (RN)(n), sequences, n = 5, 10, 20, and 30, were also prepared and studied. The linear DNA molecules had the exact sequences of the dumbbell stems. That is, the central repeat sequence in each linear duplex was flanked on both sides by the same 12-bp sequence. Melting and CD studies were also performed on the linear DNA molecules. Comparison of results obtained for the same sequences in dumbbell and linear molecular environments reveals several interesting features of the interplay between sequence-dependent structural variability, sequence length, and the unconstrained (linear) or constrained (dumbbell) molecular environments.     

Analysis of the two-state behavior of the thermal unfolding serum retinol binding protein containing a single retinol ligand.

Through the use of CD and DSC, the thermal unfolding of holo serum retinol binding protein containing a single, tightly bound retinol ligand was studied at pH 7.4. The DSC endotherm of the holoprotein ([retinol]/[protein] = 1) was asymmetric about the transition temperature of 78 degrees C. Using changes in ellipticity at 230 nm, the thermal unfolding curve was also asymmetric about the inflection point centered near 78 degrees C. van't Hoff enthalpies were determined by three means and compared to the calorimetric enthalpy (delta Hcal) of 200 kcal/mol. A van't Hoff enthalpy of 190 kcal/mol was determined from the dependence of transition temperature on the concentration of the ligand-bound protein. This value agreed well with the van't Hoff enthalpies found from fits of the DSC (delta HvH = 184 kcal/mol) and spectroscopic (delta HvH = 181 kcal/mol) curves to a two-state thermodynamic model that included ligand dissociation (NR in equilibrium with U+R, where NR is the native holoprotein, U is the unfolded apoprotein, and R is retinol). Poor agreement was obtained with a two-state model that ignored ligand dissociation (N in equilibrium with U). Furthermore, the NR in equilibrium with U+R model accounted for the asymmetry in both CD and DSC transitions and yielded a much improved fit of the data over the N in equilibrium with U model. From these considerations and simulations on other equilibrium models, it is suggested that the NR in equilibrium with U+R model is the simplest model that describes the thermal unfolding of this ligand-bound protein.(ABSTRACT TRUNCATED AT 250 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)     

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.     

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.     

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⁺.     

Studies of DNA dumbbells. I. Melting curves of 17 DNA dumbbells with different duplex stem sequences linked by T4 endloops: evaluation of the nearest-neighbor stacking interactions in DNA.

Seventeen DNA dumbbells were constructed that have duplex sequences ranging in length from 14 to 18 base pairs linked on the ends by T4 single-strand loops. Fifteen of the molecules have the core duplexes with the sequences 5'G-T-A-T-C-C-(W-X-Y-Z)-G-G-A-T-A-C3', where (W-X-Y-Z) represents a unique combination of A.T, T.A, G.C, and C.G base pairs. The remaining two molecules have the central sequence (W-X-Y-Z) = A-C and A-C-A-C-A-C. These duplex sequences were designed such that the central sequences include different combinations of the 10 possible nearest-neighbor (n-n) stacks in DNA. In this sense the set of molecules is complete and serves as a model system for evaluating sequence-dependent local stability of DNA. Optical melting curves of the samples were collected in 25, 55, 85, and 115 mM [Na+], and showed, regardless of solvent ionic strength, that the transition temperatures of the dumbbells vary by as much as 14 degrees for different molecules of the set. Results of melting experiments analyzed in terms of a n-n sequence-dependent model allowed evaluation of nine independent linear combinations of the n-n stacking interactions in DNA as a function of solvent ionic strength. Although there are in principle 10 possible different n-n interactions in DNA, these 10 are not linearly independent and therefore can not be uniquely determined. For molecules with ends, there are 9 linearly independent combinations, as opposed to circular or semiinfinite repeating copolymers where only 8 linear combinations of the 10 possible n-n interactions are linearly independent. The n-n interactions are presented as combinations of the deviations from average stacking for the 5'-3' base-pair doublets, delta Gi, and reveal several interesting features: (1) Titratable changes in the values of delta Gi with changing salt environment are observed. In all salts the most stable unique combination is delta G4 = (delta GGpC+delta GCpG)/2, and the least stable is the GpG/CpC stack, delta G2 = delta GGpG/CpC. (2) The chi 2 values of the fits of the evaluated delta Gi's to experimental data increased with decreasing [Na+], suggesting that significant interactions beyond nearest neighbors become more pronounced, particularly at 25 nM Na+.(ABSTRACT TRUNCATED AT 400 WORDS)     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in six 22-base pair deoxyoligonucleotides.

Effects of different end sequences on stability, circular dichroism spectra (CD), and enzyme binding properties were investigated for six 22-base pair, non-self-complementary duplex DNA oligomers. The center sequences of these deoxyoligonucleotides have 8-14 base pairs in common and are flanked on both sides by sequences differing in context and A-T content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming a two-state melting transition. Melting free energies (20 degrees C) of the six DNAs evaluated from DSC experiments ranged from -18.7 to -32.7 kcal/mol. van't Hoff estimates of the free energies ranged from -18.5 to -48.0 kcal/mol. With either method, the trends in free energy as a function of sequence were identical. Equilibrium binding by BamHI restriction endonuclease to the 22-base pair DNAs was also investigated. The central eight base pairs of all six molecules, 5'-A-GGATCC-A-3', contained a BamHI recognition sequence bounded by A-T base pairs. Magnesium free binding assays were performed by titering BamHI against a constant concentration of each of the deoxyoligonucleotide substrates and analyzing reaction products by gel retardation. Binding isotherms of the total amount of bound DNA versus protein concentration were constructed which provided semiquantitative estimates of the equilibrium dissociation constants for dissociation of BamHI from the six DNA oligomers. Dissociation constants ranged from 0.5 x 10(-)(9) to 12.0 x 10(-)(9) M with corresponding binding free energies of -12.5 to -10.6 (+/-0. 1) kcal/mol. An inverse relationship is found when binding and stability are compared.     

Studies of DNA dumbbells. V. A DNA triplex formed between a 28 base-pair DNA dumbbell substrate and a 16 base linear single strand

CD spectra and melting curves were collected for a 28 base-pair DNA fragment in the form of a DNA dumbbell (linked on both ends by T₄ single-strand loops) and the same DNA sequence in the linear form (without end loops). The central 16 base pairs (bp) of the 28-bp duplex region is the poly(pu) sequence: 5′-AGGAAGGAGGAAAGAG-3′. Mixtures of the dumbbell and linear DNAs with the 16-base single-strand sequence 5′-TCCTTCCTCCTTTCTC-3′ were also prepared and studied. At 22°C, CD measurements of the mixtures in 950 mM NaCl, 10 mM sodium acetate, 1 mM EDTA, pH 5.5, at a duplex concentration of 1.8 μM, provided evidence for triplex formation. Spectroscopic features of the triplexes formed with either a dumbbell or linear substrate were quite similar. Melting curves of the duplex molecules alone and in mixtures with the third strand were collected as a function of duplex concentration from 0.16 to 2.15 μM. Melting curves of the dumbbell alone and mixtures with the third strand were entirely independent of DNA concentration. In contrast, melting curves of the linear duplex alone or mixed with the third strand were concentration dependent. At identical duplex concentrations, the dumbbell alone melts ～20°C higher than the linear duplex. The curve of the linear duplex displayed a significant pretransition probably due to end fraying. On melting curves of mixtures of the dumbbell or linear duplex with the third strand, a low temperature transition with much lower relative hyperchromicity change (～ 5%) was observed. This transition was attributed to the melting of a new molecular species, e.g., the triplex formed between the duplex and single-strand DNA molecules. In the case of the dumbbell/single-strand mixture, these melting transitions of the triplex and the dumbbell were entirely resolvable. In contrast, the melting transitions of the linear duplex and the triplex overlapped, thereby preventing their clear distinction. To analyze the data, a three-state equilibrium model is presented. The analysis utilizes differences in relative absorbance vs temperature curves of dumbbells (or linear molecules) alone and in mixtures with the third strand. From the model analysis a straightforward derivation of f_T(T), the fraction of triplex as a function of temperature, was obtained. Analysis of f_T vs temperature curves, in effect melting curves of the triplexes, provided evaluation of thermodynamic parameters of the melting transition. For the triplex formed with the dumbbell substrate, the total transition enthalpy is ΔH_T = 118.4 ± 12.8 kcal/mol (7.4 ± 0.8 kcal/mol per triplet unit) and the total transition entropy is ΔS_T = 344 ± 36.8 cal/K · mol (eu) (21.5 ± 2.3 eu per triple unit). The transition curves of the triplex formed with the linear duplex substrate displayed two distinct regions. A broad pretransition region from f_T = 0 to 0.55 and a higher, sharper transition above f_T = 0.55. The transition parameters derived from the lower temperature region of the curve are ΔH′_T = 44.8 ± 9.6 kcal/mol and ΔS′_T = 112 ± 33.6 eu (or ΔH′ = 2.8 ± 0.6 kcal/mol and ΔS′ = 7.0 ± 2.1 eu per triplet). These values are probably too small to correspond to actual melting of the triplex but instead likely reveal effects of end fraying of the duplex substrate on triplex stability. Transition parameters of the upper transition are ΔH′_T = 128.0 ± 2.3 kcal/mol and ΔS′_T = 379.2 ± 6.4 eu (ΔH′ = 8.0 ± 0.2 kcal/mol and ΔS′ = 23.7 ± 0.4 eu per triplet) in good agreement (within experimental error) with the transition parameters of the triplex formed with the dumbbell substrate. Supposing this upper transition reflects actual dissociation of the third strand from the linear duplex substrate this triplex is comparable in thermodynamic stability to the triplex formed with a dumbbell substrate. Even so, the biphasic melting character of the linear triplex obscures the whole analysis, casting doubt on its absolute reliability. Apparently triplexes formed with a dumbbell substrate offer technical advantages over triplexes formed from linear or hairpin duplex substrates for studies of DNA triplex stability. © 1993 John Wiley & Sons, Inc.     

Stability of recombinant Lys25-ribonuclease T1

The conformational stability of recombinant Lys25-ribonuclease T1 has been determined by differential scanning microcalorimetry (DSC), UV-monitored thermal denaturation measurements, and isothermal Gdn.HCl unfolding studies. Although rather different extrapolation procedures are involved in calculating the Gibbs free energy of stabilization, there is fair agreement between the delta G degrees values derived from the three different experimental techniques at pH 5, theta = 25 degrees C: DSC, 46.6 +/- 2.1 kJ/mol; UV melting curves, 48.7 +/- 5 kJ/mol; Gdn.HCl transition curves, 40.8 +/- 1.5 kJ/mol. Thermal unfolding of the enzyme is a reversible process, and the ratio of the van't Hoff and calorimetric enthalpy, delta HvH/delta Hcal, is 0.97 +/- 0.06. This result strongly suggests that the unfolding equilibrium of Lys25-ribonuclease T1 is adequately described by a simple two-state model. Upon unfolding the heat capacity increases by delta Cp degrees = 5.1 +/- 0.5 kJ/(mol.K). Similar values have been found for the unfolding of other small proteins. Surprisingly, this denaturational heat capacity change practically vanishes in the presence of moderate NaCl concentrations. The molecular origin of this effect is not clear; it is not observed to the same extent in the unfolding of bovine pancreatic ribonuclease A, which was employed in control experiments. NaCl stabilizes Lys25-ribonuclease T1. The transition temperature varies with NaCl activity in a manner that suggests two limiting binding equilibria to be operative. Below approximately 0.2 M NaCl activity unfolding is associated with dissociation of about one ion, whereas above that concentration about four ions are released in the unfolding reaction.(ABSTRACT TRUNCATED AT 250 WORDS)     

Equilibrium thermal transitions of collagen model peptides

The folding of collagen in vitro is very slow and presents difficulties in reaching equilibrium, a feature that may have implications for in vivo collagen function. Peptides serve as good model systems for examining equilibrium thermal transitions in the collagen triple helix. Investigations were carried out to ascertain whether a range of synthetic triple-helical peptides of varying sequences can reach equilibrium, and whether the triple helix to unfolded monomer transition approximates a two-state model. The thermal transitions for all peptides studied are fully reversible given sufficient time. Isothermal experiments were carried out to obtain relaxation times at different temperatures. The slowest relaxation times, on the order of 10-15 h, were observed at the beginning of transitions, and were shown to result from self-association limited by the low concentration of free monomers, rather than cis-trans isomerization. Although the fit of the CD equilibrium transition curves and the concentration dependence of T(m) values support a two-state model, the more rigorous comparison of the calorimetric enthalpy to the van't Hoff enthalpy indicates the two-state approximation is not ideal. Previous reports of melting curves of triple-helical host-guest peptides are shown to be a two-state kinetic transition, rather than an equilibrium transition.     

Differences in thermal stability between reduced and oxidized cytochrome b562 from Escherichia coli

The thermal stabilities of ferri- and ferrocytochrome b562 were examined. Thermally induced spectral changes, monitored by absorption and second-derivative spectroscopies, followed the dissociation of the heme moiety and the increased solvation of tyrosine residue(s) located in close proximity to the heme binding site. All observed thermal transitions were independent of the rate of temperature increase (0.5-2 degrees C/min), and the denatured protein exhibited partial to near-complete reversibility upon return to ambient temperature. The extent of renaturation of cytochrome b562 is dependent on the amount of time the unfolded conformer is exposed to temperatures above the transition temperature, Tm. All thermally induced spectra changes fit a simple two-state model, and the thermal transition was assumed to be reversible. The thermal transition for ferrocytochrome b562 yielded Tm and van't Hoff enthalpy (delta HvH) values of 81.0 degrees C and 137 kcal/mol, respectively. In contrast, Tm and delta HvH values obtained for the ferricytochrome were 66.7 degrees C and 110 kcal/mol, respectively. The estimated increase in the stabilization free energy at the Tm of ferricytochrome b562 following the one-electron reduction to the ferrous form, where delta delta G = delta Tm delta Sm [delta Sm = 324 cal/(K.mol), delta Tm = 14.3 degrees C] [Becktel, W. J., & Schellman, J. A. (1987) Biopolymers 26, 1859-1877], is 4.6 kcal/mol.     

Thermodynamics of DNA hairpins: contribution of loop size to hairpin stability and ethidium binding.

A combination of calorimetric and spectroscopic techniques was used to evaluate the thermodynamic behavior of a set of DNA hairpins with the sequence d(GCGCTnGCGC), where n = 3, 5 and 7, and the interaction of each hairpin with ethidium. All three hairpins melt in two-state monomolecular transitions, with tm's ranging from 79.1 degrees C (T3) to 57.5 degrees C (T7), and transition enthalpies of approximately 38.5 kcal mol-1. Standard thermodynamic profiles at 20 degrees C reveal that the lower stability of the T5 and T7 hairpins corresponds to a delta G degree term of +0.5 kcal mol-1 per thymine residue, due to the entropic ordering of the thymine loops and uptake of counterions. Deconvolution of the ethidium-hairpin calorimetric titration curves indicate two sets of binding sites that correspond to one ligand in the stem with binding affinity, Kb, of approximately 1.8 x 10(6) M-1, and two ligands in the loops with Kb of approximately 4.3 x 10(4) M-1. However, the binding enthalpy, delta Hb, ranges from -8.6 (T3) to -11.6 kcal mol-1 (T7) for the stem site, and -6.6 (T3) to -12.7 kcal mol-1 (T7) for the loop site. Relative to the T3 hairpin, we obtained an overall thermodynamic contribution (per dT residue) of delta delta Hb = delta(T delta Sb) = -0.7(5) kcal mol-1 for the stem sites and delta delta Hb = delta(T delta Sb) = -1.5 kcal mol-1 for the loop sites. Therefore, the induced structural perturbations of ethidium binding results in a differential compensation of favorable stacking interactions with the unfavorable ordering of the ligands.     

Thermodynamic, spectroscopic, and equilibrium binding studies of DNA sequence context effects in four 40 base pair deoxyoligonucleotides

Effects of different end sequences on melting, circular dichroism spectra (CD), and enzyme binding properties were investigated for four 40 base pair, non-self-complementary duplex DNA oligomers. The center sequences of these oligoduplexes have either of two 22 base pair modules flanked on both sides by sequences differing in AT content. Temperature-induced melting transitions monitored by differential scanning calorimetry (DSC) and ultraviolet absorbance were measured for the six duplexes in buffered 115 mM Na(+) solutions. Values of the melting transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal), were obtained directly from DSC experiments. Melting transition parameters, DeltaH(vH) and DeltaS(vH), were also estimated from a van't Hoff analysis of optical melting curves collected as a function of DNA concentration, assuming that the melting transition is two-state. Melting free energies (20 degrees C) evaluated from DSC melting experiments on the four duplex DNAs ranged from -52.2 to -77.5 kcal/mol. Free energies based on the van't Hoff analysis were -37.9 to -58.8 kcal/mol. Although the values are different, trends in the melting free energies of the four duplex DNAs as a function of sequence were identical in both DSC and optical analyses. Subject to several assumptions, values for the initiation free energy were estimated for each duplex, defined as DeltaG(int) = DeltaG(cal) - DeltaG(pred), where DeltaG(cal) is the experimental free energy at 20 degrees C determined from the experimentially measured values of the transition enthalpy, DeltaH(cal), and entropy, DeltaS(cal). The predicted free energy of the sequence, DeltaG(pred)(20 degrees C), is obtained using published nearest-neighbor sequence stability values. For three of the four duplexes, values of DeltaG(int) are essentially nil. In contrast, the duplex with 81.8% GC has a considerably higher estimate of DeltaG(int) = 7.1 kcal/mol. The CD spectra for the six duplexes collected over the wavelength range from 200 to 320 nm are also sequence-dependent. Factor analysis of the CD spectra by singular value decomposition revealed that the experimental CD spectra could be reconstructed from linear combinations of two minor and one major subspectra. Changes in the coefficients of the major subspectrum for different sequences reflect incremental sequence-dependent variations of the CD spectra. Equilibrium binding by BamHI restriction endonuclease to the 40 base pair DNAs whose central eight base pairs contain the recognition sequence for BamHI restriction enzyme bounded by A.T base pairs, 5'-A-GGATCC-A-3' was investigated. Binding assays were performed by titering BamHI against a constant concentration of each of the duplex DNA substrates, in the absence of Mg(2+), followed by analysis by gel retardation. Under the conditions employed, the enzyme binds but does not cleave the DNAs. Results of the assays revealed two binding modes with retarded gel mobilities. Binding isotherms for the fraction of bound DNA species versus enzyme concentration for each binding mode were constructed and analyzed with a simple two-step equilibrium binding model. This analysis provided semiquantitative estimates on the equilibrium binding constants for BamHI to the four DNAs. Values obtained for the binding constants varied only 7-fold and ranged from 6 x 10(-)(8) to 42 x 10(-)(8) M, with binding free energies from -8.6 to -9.7 (+/- 0.2) kcal/mol depending on the sequence that flanks the enzyme binding site. Unlike what was found earlier in binding studies of the 22 base pair duplexes that constitute the core modules of the present 40-mers [Riccelli, P. V., Vallone, P. M., Kashin, I., Faldasz, B. D., Lane, M. J., and Benight, A. S. (1999) Biochemistry 38, 11197-11208], no obvious relationship between binding and stability was found for these longer DNAs. Apparently, effects of flanking sequence stability on restriction enzyme binding may only be measurable in very short duplex deoxyoligonucl     

Is it always possible to distinguish two- and three-state systems by evaluating the van't Hoff enthalpy?

Many small globular proteins are traditionally classified as thermodynamical two-state systems, i.e., the protein is either in the native, active state (folded) or in the denatured state (unfolded). We challenge this view and show that there may exist (protein) systems for which a van't Hoff analysis of experimental data cannot determine whether the system corresponds to two or three thermodynamical states when only temperatures in a narrow temperature region around the transition are considered. We generalize a widely employed two-state protein folding model to include a third, transition state. For this three-state system we systematically study the deviation of the calorimetric enthalpy (heat of transition) from the van't Hoff enthalpy, a measure of the two-stateness of a transition. We show that under certain conditions the heat capacity of the three-state system can be almost indistinguishable from the heat capacity for the two-state system over a broad temperature interval. The consequence may be that some three-state (or even more than three-states) systems have been misinterpreted as two-state systems when the conclusion is drawn solely upon the van't Hoff enthalpy. These findings are important not only for proteins, but also for the interpretation of thermodynamical systems in general.