Enthalpy and heat capacity changes for formation of an oligomeric DNA duplex: interpretation in terms of coupled processes of formation and association of single-stranded helices.

The thermodynamics of self-assembly of a 14 base pair DNA double helix from complementary strands have been investigated by titration (ITC) and differential scanning (DSC) calorimetry, in conjunction with van't Hoff analysis of UV thermal scans of individual strands. These studies demonstrate that thermodynamic characterization of the temperature-dependent contributions of coupled conformational equilibria in the individual "denatured" strands and in the duplex is essential to understand the origins of duplex stability and to derive stability prediction schemes of general applicability. ITC studies of strand association at 293 K and 120 mM Na+ yield an enthalpy change of -73 +/- 2 kcal (mol of duplex)-1. ITC studies between 282 and 312 K at 20, 50, and 120 mM Na+ show that the enthalpy of duplex formation is only weakly salt concentration-dependent but is very strongly temperature-dependent, decreasing approximately linearly with increasing temperature with a heat capacity change (282-312 K) of -1.3 +/- 0.1 kcal K-1 (mol of duplex)-1. From DSC denaturation studies in 120 mM Na+, we obtain an enthalpy of duplex formation of -120 +/- 5 kcal (mol of duplex)-1 and an estimate of the corresponding heat capacity change of -0.8 +/- 0.4 kcal K-1 (mol of duplex)-1 at the Tm of 339 K. van't Hoff analysis of UV thermal scans on the individual strands indicates that single helix formation is noncooperative with a temperature-independent enthalpy change of -5.5 +/- 0.5 kcal at 120 mM Na+. From these observed enthalpy and heat capacity changes, we obtain the corresponding thermodynamic quantities for two fundamental processes: (i) formation of single helices from disordered strands, involving only intrastrand (vertical) interactions between neighboring bases; and (ii) formation of double helices by association (docking) of single helical strands, involving interstrand (horizontal and vertical) interactions. At 293 K and 120 mM Na+, we calculate that the enthalpy change for association of single helical strands is approximately -64 kcal (mol of duplex)-1 as compared to -210 kcal (mol of duplex)-1 calculated for duplex formation from completely unstructured single strands and to the experimental ITC value of -73 kcal (mol of duplex)-1. The intrinsic heat capacity change for association of single helical strands to form the duplex is found to be small and positive [ approximately 0.1 kcal K-1 (mol of duplex)-1], in agreement with the result of a surface area analysis, which also predicts an undetectably small heat capacity change for single helix formation.     

Thermodynamics of RNA duplexes with tandem mismatches containing a uracil-uracil pair flanked by C.G/G.C or G.C/A.U closing base pairs

The thermodynamics governing the denaturation of RNA duplexes containing 8 bp and a central tandem mismatch or 10 bp were evaluated using UV absorbance melting curves. Each of the eight tandem mismatches that were examined had one U-U pair adjacent to another noncanonical base pair. They were examined in two different RNA duplex environments, one with the tandem mismatch closed by G.C base pairs and the other with G.C and A.U closing base pairs. The free energy increments (Delta Gdegrees(loop)) of the 2 x 2 loops were positive, and showed relatively small differences between the two closing base pair environments. Assuming temperature-independent enthalpy changes for the transitions, (Delta Gdegrees(loop)) for the 2 x 2 loops varied from 0.9 to 1.9 kcal/mol in 1 M Na(+) at 37 degrees C. Most values were within 0.8 kcal/mol of previously estimated values; however, a few sequences differed by 1.2-2.0 kcal/mol. Single strands employed to form the RNA duplexes exhibited small noncooperative absorbance increases with temperature or transitions indicative of partial self-complementary duplexes. One strand formed a partial self-complementary duplex that was more stable than the tandem mismatch duplexes it formed. Transitions of the RNA duplexes were analyzed using equations that included the coupled equilibrium of self-complementary duplex and non-self-complementary duplex denaturation. The average heat capacity change (DeltaC(p)) associated with the transitions of two RNA duplexes was estimated by plotting DeltaH degrees and DeltaS degrees evaluated at different strand concentrations as a function of T(m) and ln T(m), respectively. The average DeltaC(p) was 70 +/- 5 cal K(-)(1) (mol of base pairs)(-)(1). Consideration of this heat capacity change reduced the free energy of formation at 37 degrees C of the 10 bp control RNA duplexes by 0.3-0.6 kcal/mol, which may increase Delta Gdegrees(loop) values by similar amounts.     

Translational Entropy and DNA Duplex Stability

Investigation of folding/unfolding DNA duplexes of various size and composition by superprecise calorimetry has revised several long-held beliefs concerning the forces responsible for the formation of the double helix. It was established that: 1) the enthalpy and the entropy of duplex unfolding are temperature dependent, increasing with temperature rise and having the same heat capacity increment for CG and AT pairs; 2) the enthalpy of AT melting is greater than that of the CG pair, so the stabilizing effect of the CG pair in comparison with AT results not from its larger enthalpic contribution (as expected from its extra hydrogen bond), but from the larger entropic contribution of the AT pair that results from its ability to fix ordered water in the minor groove and release it upon duplex unfolding; 3) the translation entropy, resulting from the appearance of a new kinetic unit on duplex dissociation, determines the dependence of duplex stability on its length and its concentration (it is an order-of-magnitude smaller than predicted from the statistical mechanics of gases and is fully expressed by the stoichiometric correction term); 4) changes in duplex stability on reshuffling the sequence (the “nearest-neighbor effect”) result from the immobilized water molecules fixed by AT pairs in the minor groove; and 5) the evaluated thermodynamic components permit a quantitative expression of DNA duplex stability.     

The energetic basis of the DNA double helix: a combined microcalorimetric approach

Microcalorimetric studies of DNA duplexes and their component single strands showed that association enthalpies of unfolded complementary strands into completely folded duplexes increase linearly with temperature and do not depend on salt concentration, i.e. duplex formation results in a constant heat capacity decrement, identical for CG and AT pairs. Although duplex thermostability increases with CG content, the enthalpic and entropic contributions of an AT pair to duplex formation exceed that of a CG pair when compared at the same temperature. The reduced contribution of AT pairs to duplex stabilization comes not from their lower enthalpy, as previously supposed, but from their larger entropy contribution. This larger enthalpy and particularly the greater entropy results from water fixed by the AT pair in the minor groove. As the increased entropy of an AT pair exceeds that of melting ice, the water molecule fixed by this pair must affect those of its neighbors. Water in the minor groove is, thus, orchestrated by the arrangement of AT groups, i.e. is context dependent. In contrast, water hydrating exposed nonpolar surfaces of bases is responsible for the heat capacity increment on dissociation and, therefore, for the temperature dependence of all thermodynamic characteristics of the double helix.     

Effect of mutation of the Sac7d intercalating residues on the temperature dependence of DNA distortion and binding thermodynamics

Sac7d is a small chromatin protein from the hyperthermophile Sulfolobus acidocaldarius which kinks duplex DNA by approximately 66 degrees at a single base pair step with intercalation of V26 and M29 side chains. Site-directed mutagenesis coupled with calorimetric and spectroscopic data has been used to characterize the influence of the intercalating side chains on the structure and thermodynamics of the DNA complex from 5 to 85 degrees C. Two single-alanine substitutions (V26A and M29A) and five double-glycine, -alanine, -leucine, -phenylalanine, and -tryptophan substitutions of the surface residues have been created. NMR and fluorescence titrations indicated that the substitutions had little effect on the structure of the protein or DNA binding site size. Each of the mutant proteins demonstrated a temperature-dependent binding enthalpy which was correlated with a similar temperature dependence in the structure of the complex reflected by changes in fluorescence and circular dichroism. A positive heat capacity change (DeltaC(p)) for DNA binding was observed for only those mutants which also demonstrated a thermotropic structural transition in the complex, and the temperature range for the positive DeltaC(p) coincided with that observed for the structural transition. The thermodynamic data are interpreted using a model in which binding is linked to an endothermic distortion of the DNA in the complex. The results support the proposal that the unfavorable enthalpy of binding of Sac7d at 25 degrees C is due in part to the distortion of DNA.     

Thermodynamics of DNA branching.

Branched DNA molecules arise transiently as intermediates in genetic recombination or on extrusion of cruciforms from covalent circular DNA duplexes that contain palindromic sequences. The free energy of these structures relative to normal DNA duplexes is of interest both physically and biologically. Oligonucleotide complexes that can form stable branched structures, DNA junctions, have made it possible to model normally unstable branched states of DNA such as Holliday recombinational intermediates. We present here an evaluation of the free energy of creating four-arm branch points in duplex DNA, using a system of two complementary junctions and four DNA duplexes formed from different combinations of the same set of eight 16-mer strands. The thermodynamics of formation of each branched structure from the matching pair of intact duplexes have been estimated in two experiments. In the first, labeled strands are allowed to partition between duplexes and junctions in a competition assay on polyacrylamide gels. In the second, the heats of forming branched or linear molecules from the component strands have been determined by titration microcalorimetry at several temperatures. Taken together these measurements allow us to determine the standard thermodynamic parameters for the process of creating a branch in an otherwise normal DNA duplex. The free energy for reacting two 16-mer duplexes to yield a four-arm junction in which the branch site is incapable of migrating is + 1.1 (+/- 0.4) kcal mol-1 (at 18 degrees C, 10 mM-Mg2+). Analysis of the distribution of duplex and tetramer products by electrophoresis confirms that the free energy difference between the four duplexes and two junctions is small at this temperature. The associated enthalpy change at 18 degrees C is +27.1 (+/- 1.3) kcal mol-1, while the entropy is +89 (+/- 30) cal K-1 mol-1. The free energy for branching is temperature dependent, with a large unfavorable enthalpy change compensated by a favorable entropy term. Since forming one four-stranded complex from two duplexes should be an entropically unfavorable process, branch formation is likely to be accompanied by significant changes in hydration and ion binding. A significant apparent delta Cp is also observed for the formation of one mole of junction, +0.97 (+/-0.05) kcal deg-1 mol-1.     

Energetics of sequence-specific protein-DNA association: conformational stability of the DNA binding domain of integrase Tn916 and its cognate DNA duplex

Sequence-specific DNA recognition by bacterial integrase Tn916 involves structural rearrangements of both the protein and the DNA duplex. Energetic contributions from changes of conformation, thermal motions and soft vibrational modi of the protein, the DNA, and the complex significantly influence the energetic profile of protein-DNA association. Understanding the energetics of such a complicated system requires not only a detailed calorimetric investigation of the association reaction but also of the components in isolation. Here we report on the conformational stability of the integrase Tn916 DNA binding domain and its cognate 13 base pair target DNA duplex. Using a combination of temperature and denaturant induced unfolding experiments, we find that the 74-residue DNA binding domain is compact and unfolds cooperatively with only small deviation from two-state behavior. Scanning calorimetry reveals an increase of the heat capacity of the native protein attributable to increased thermal fluctuations. From the combined calorimetric and spectroscopic experiments, the parameters of protein unfolding are T(m) = 43.8 +/- 0.3 degrees C, DeltaH(m) = 255 +/- 18 kJ mol(-1), DeltaS(m) = 0.80 +/- 0.06 kJ mol(-1), and DeltaC(p) = 5.0 +/- 0.8 kJ K(-1) mol(-1). The DNA target duplex displays a thermodynamic signature typical of short oligonucleotide duplexes: significant heat absorption due to end fraying and twisting precedes cooperative unfolding and dissociation. The parameters for DNA unfolding and dissociation are DeltaH(m) = 335 +/- 4 kJ mol(-1) and DeltaC(p) = 2.7 +/- 0.9 kJ K(-(1) mol(-1). The results reported here have been instrumental in interpreting the thermodynamic features of the association reaction of the integrase with its 13 base pair target DNA duplex reported in the accompanying paper [Milev et al. (2003) Biochemistry 42, 3481-3491].     

Stability of DNA duplexes containing GG, CC, AA, and TT mismatches

We employed salt-dependent differential scanning calorimetric measurements to characterize the stability of six oligomeric DNA duplexes (5'-GCCGGAXTGCCGG-3'/5'-CCGGCAYTCCGGC-3') that contain in the central XY position the GC, AT, GG, CC, AA, or TT base pair. The heat-induced helix-to-coil transitions of all the duplexes are associated with positive changes in heat capacity, DeltaC(p), ranging from 0.43 to 0.53 kcal/mol. Positive values of DeltaC(p) result in strong temperature dependences of changes in enthalpy, DeltaH degrees, and entropy, DeltaS degrees , accompanying duplex melting and cause melting free energies, DeltaG degrees, to exhibit characteristically curved shapes. These observations suggest that DeltaC(p) needs to be carefully taken into account when the parameters of duplex stability are extrapolated to temperatures distant from the transition temperature, T(M). Comparison of the calorimetric and van't Hoff enthalpies revealed that none of the duplexes studied in this work exhibits two-state melting. Within the context of the central AXT/TYA triplet, the thermal and thermodynamic stabilities of the duplexes in question change in the following order: GC > AT > GG > AA approximately TT > CC. Our estimates revealed that the thermodynamic impact of the GG, AA, and TT mismatches is confined within the central triplet. In contrast, the thermodynamic impact of the CC mismatch propagates into the adjacent helix domains and may involve 7-9 bp. We discuss implications of our results for understanding the origins of initial recognition of mismatched DNA sites by enzymes of the DNA repair machinery.     

Thermodynamic dependence of DNA/DNA and DNA/RNA hybridization reactions on temperature and ionic strength

The thermodynamics of 5'-ATGCTGATGC-3' binding to its complementary DNA and RNA strands was determined in sodium phosphate buffer under varying conditions of temperature and salt concentration from isothermal titration calorimetry (ITC). The Gibbs free energy change, DeltaG degrees of the DNA hybridization reactions increased by about 6 kJ mol(-1) from 20 degrees C to 37 degrees C and exhibited heat capacity changes of -1.42 +/- 0.09 kJ mol(-1) K(-1) for DNA/DNA and -0.87 +/- 0.05 kJ mol(-1) K(-1) for DNA/RNA. Values of DeltaG degrees decreased non-linearly by 3.5 kJ mol(-1) at 25 degrees C and 6.0 kJ mol(-1) at 37 degrees C with increase in the log of the sodium chloride concentration from 0.10 M to 1.0 M. A near-linear relationship was observed, however, between DeltaG degrees and the activity coefficient of the water component of the salt solutions. The thermodynamic parameters of the hybridization reaction along with the heat capacity changes were combined with thermodynamic contributions from the stacking to unstacking transitions of the single-stranded oligonucleotides from differential scanning calorimetry (DSC) measurements, resulting in good agreement with extrapolation of the free energy changes to 37 degrees C from the melting transition at 56 degrees C.     

Enthalpy distribution functions for the unwinding of a short DNA duplex

In this article we use the published heat capacity data of Dragan et al. (J Mol Biol 2003, 327, 293-411) for a short DNA duplex to calculate the enthalpy probability distribution for this species as a function of temperature. Our approach is based on a procedure that we developed (Poland, D. J Chem Phys 2000, 112, 6554) whereby one obtains moments of the enthalpy distribution from the temperature dependence of the heat capacity. One then uses the maximum-entropy method to construct the enthalpy probability distribution from the set of enthalpy moments. For the DNA duplex treated here the heat capacity goes through a maximum as a function of temperature reflecting the unwinding of the duplex structure. In the neighborhood of the heat capacity maximum, the enthalpy distribution functions show a clear bimodal structure, indicating the coexistence of two distinct states, the duplex and the single-strand state. The probabilities of theses two states can be estimated from the enthalpy distribution functions and can be used to calculate the temperature dependence of the equilibrium constant for the unwinding of the DNA duplex. This example illustrates that the temperature dependence of the heat capacity can be used to give a detailed picture of conformational transitions in biological macromolecules. In particular, the structure of the enthalpy distribution in this case allows one to see the temperature evolution of the two-state distribution in detail.     

Energetics of sequence-specific protein-DNA association: binding of integrase Tn916 to its target DNA

The DNA binding domain of the transposon Tn916 integrase (INT-DBD) binds to its DNA target site by positioning the face of a three-stranded antiparallel beta-sheet within the major groove. Binding of INT-DBD to a 13 base pair duplex DNA target site was studied by isothermal titration calorimetry, differential scanning calorimetry, thermal melting followed by circular dichroism spectroscopy, and fluorescence spectroscopy. The observed heat capacity change accompanying the association reaction (DeltaC(p)) is temperature-dependent, decreasing from -1.4 kJ K(-1) mol(-1) at 4 degrees C to -2.9 kJ K(-1) mol(-1) at 30 degrees C. The reason is that the partial molar heat capacities of the free protein, the free DNA duplex, and the protein-DNA complex are not changing in parallel when the temperature increases and that thermal motions of the protein and the DNA are restricted in the complex. After correction for this effect, DeltaC(p) is -1.8 kJ K(-1) mol(-1) and temperature-independent. However, this value is still higher than DeltaC(p) of -1.2 kJ K(-1) mol(-1) estimated by semiempirical methods from dehydration of surface area buried at the complex interface. We propose that the discrepancy between the measured and the structure-based prediction of binding energetics is caused by incomplete dehydration of polar groups in the complex. In support, we identify cavities at the interface that are large enough to accommodate approximately 10 water molecules. Our results highlight the difficulties of structure-based prediction of DeltaC(p) (and other thermodynamic parameters) and emphasize how important it is to consider changes of thermal motions and soft vibrational modi in protein-DNA association reactions. This requires not only a detailed investigation of the energetics of the complex but also of the folding thermodynamics of the protein and the DNA alone, which are described in the accompanying paper [Milev et al. (2003) Biochemistry 42, 3492-3502].     

Thermodynamic parameters for DNA sequences with dangling ends

The thermodynamic contributions to duplex formation of all 32 possible single-nucleotide dangling ends on a Watson-Crick pair are reported. In most instances, dangling ends are stabilizing with free energy contributions ranging from +0.48 (GT(A)) to-0.96 kcal/mol (). In comparison, Watson-Crick nearest-neighbor increments range from -0. 58 (TA/AT) to -2.24 (GC/CG) kcal/mol. Hence, in some cases, a dangling end contributes as much to duplex stability as a Watson-Crick A-T base pair. The implications of these results for DNA probe design are discussed. Analysis of the sequence dependence of dangling-end stabilities show that the nature of the closing base pair largely determines the stabilization. For a given closing base pair, however, adenine dangling ends are always more or equally as stable as the other dangling nucleotides. Moreover, 5' dangling ends are more or equally as stabilizing as their 3' counterparts. Comparison of DNA with RNA dangling-end motifs shows that DNA motifs with 5' dangling ends contribute to stability equally or more than their RNA counterparts. Conversely, RNA 3' dangling ends contribute to stability equally or more than their DNA counterparts. This data set has been incorporated into a DNA secondary structure prediction algorithm (DNA MFOLD) (http://mfold2.wustl.edu/mfold/dna/for m1.cgi) as well as a DNA hybridization prediction algorithm (HYTHERtrade mark) (http://jsl1.chem.wayne.edu/Hyther/hythermenu .html).     

Temperature dependence of the gel electrophoretic mobility of superhelical DNA.

We have determined the gel electrophoretic behavior of closed circular plasmid pSM1 DNA (5420 bp) as a function of both temperature and of linking number (Lk). At temperatures below 37 degrees, the electrophoretic mobility first increases, then becomes constant as Lk is decreased below that of the relaxed closed DNA. As the temperature is increased above 37 degrees the electrophoretic mobility first increases as Lk decreases and then varies in a cyclic manner with further decreases in Lk. As the temperature is increased over the range 37 degrees - 65 degrees the cyclic behavior is manifested at progressively smaller decreases in Lk and the amplitude of the cycles increases. We interpret the results in terms of the early melting of superhelical DNA, in which the free energy associated with superhelix formation is progressively transferred to local denaturation. Using a two state approximation, we estimate the free energy change in the first cyclic transition to be 35 Kcal/mole DNA at 37 degrees and to decrease linearly with temperature. The free energy becomes equal to zero at a temperature of 71.6 degrees, which lies within 3 degrees of the melting temperature for the corresponding nicked circular DNA. From the slope of this relationship we estimate the apparent entropy and enthalpy of the first mobility transition to be 6.0 Kcal/mole base pair and 17.3 cal/mole base pair/degree, values consistent with duplex melting.     

On the stability of nucleic acid structures in solution: enthalpy-entropy compensations, internal rotations and reversibility.

The thermodynamics of self-association (stacking) of free bases and nucleotides, intramolecular stacking in dinucleotides, nearest-neighbour base pair stacking interactions in duplex DNA and RNA, and the formation of hairpin loops illustrate enthalpy/entropy compensations. Large stacking exothermicities are associated with large negative entropy changes that ensure that delta G is small, permitting readily reversible associations in solution. We rationalise enthalpy/entropy compensations with reference to residual motions and torsional vibrations which make a larger entropic contribution to binding when - delta H approximately kT (thermal energy at room temperature), than when - delta H >> kT. We present a factorisation of experimental free energies for helix formation in terms of approximate contributions from the restriction of rotations, hydrophobic interactions, electrostatic interactions due to base stacking, and contributions from hydrogen bonding, and estimate the adverse free energy cost per rotor (mainly entropy) of ordering the phosphate backbone as between 1.9 and 5.4 kJ mol-1 [averaged over 12 rotors per base pair for A-U on A-U stacking (lower limit), and G-C on C-G stacking (upper limit)]. The largest cost is associated with the most exothermic stacking interactions, while the range of values is consistent with earlier conclusions from data on the fusion of hydrocarbon chains (lower value), and with entropy changes in covalent isomerisations of small molecules involving severe restrictions (upper value).     

Specific DNA binding by the homeodomain Nkx2.5(C56S): detection of impaired DNA or unfolded protein by isothermal titration calorimetry

Titrations of specific 18-bp duplex DNA with the cardiac-specific homeodomain Nkx2.5(C56S) have utilized an ultrasensitive isothermal titration calorimeter (ITC). As the free DNA nears depletion, we observe large apparent decreases in the binding enthalpy when the DNA is impaired or when the temperature is sufficiently high to produce some unfolding of the free protein. Either effect can be attributed to refolding of the biopolymer that occurs as a result of stabilization due to the large favorable change in free energy on the homeodomain binding to DNA (-49.4 kJ/mol at 298 K). In either case, thermodynamic parameters obtained in such ITC experiments are unreliable. By using a lower temperature (85 vs. 95 degrees C) during the annealing of complementary DNA strands, damage of the 18-bp duplex DNA (T(m) = 72 degrees C) is avoided, and titrations with the homeodomain are normal at temperatures from 10 to 40 degrees C when >95% of the protein is folded. Under the latter conditions, the heat capacity plot is linear with a DeltaC(p) value of -0.80 +/- 0.03 kJ K(-1) mol(-1), which is more negative than that calculated from the burial of solvent accessible surface areas (-0.64 +/- 0.05 kJ K(-1) mol(-1)), consistent with water structures being at the protein-DNA interfaces.     

Thermodynamic properties of the specific binding between Ag+ ions and C:C mismatched base pairs in duplex DNA

Metal-mediated base pairs formed by the interaction between metal ions and artificial bases in oligonucleotides have been developed for potential applications in nanotechnology. We recently found that a natural C:C mismatched base pair bound to an Ag(+) ion to generate a novel metal-mediated base pair in duplex DNA. Preparation of the novel C-Ag-C base pair involving natural bases is more convenient than that of metal-mediated base pairs involving artificial bases because time-consuming base synthesis is not required. Here, we examined the thermodynamic properties of the binding between the Ag(+) ion and each of single and double C:C mismatched base pair in duplex DNA by isothermal titration calorimetry. The Ag(+) ion specifically bound to the C:C mismatched base pair at a 1:1 molar ratio with 10(6) M(-1) binding constant, which was significantly larger than those for nonspecific metal ion-DNA interactions. The specific binding between the Ag(+) ion and the single C:C mismatched base pair was mainly driven by the positive dehydration entropy change and the negative binding enthalpy change. In the interaction between the Ag(+) ion and each of the consecutive and interrupted double C:C mismatched base pairs, stoichiometric binding at a 1:1 molar ratio was achieved in each step of the first and second Ag(+) binding. The binding affinity for the second Ag(+) binding was similar to that for the first Ag(+) binding. Stoichiometric binding without interference and negative cooperativity may be favorable for aligning multiple Ag(+) ions in duplex DNA for applications of the metal-mediated base pairs in nanotechnology.     

Thermodynamic parameters for loop formation in RNA and DNA hairpin tetraloops.

We determined the melting temperatures (Tm) and thermodynamic parameters of 15 RNA and 19 DNA hairpins at 1 M NaCl, 0.01 M sodium phosphate, 0.1 mM EDTA, at pH 7. All these hairpins have loops of four bases, the most common loop size in 16S and 23S ribosomal RNAs. The RNA hairpins varied in loop sequence, loop-closing base pair (A.U, C.G, or G.C), base sequence of the stem, and stem size (four or five base pairs). The DNA hairpins varied in loop sequence, loop-closing base pair (C.G, or G.C), and base sequence of the four base-pair stem. Thermodynamic properties of a hairpin may be represented by nearest-neighbor interactions of the stem plus contributions from the loop. Thus, we obtained thermodynamic parameters for the formation of RNA and DNA tetraloops. For the tetraloops we studied, a free energy of loop formation (at 37 degrees C) of about +3 kcal/mol is most common for either RNA or DNA. There are extra stable loops with delta G degrees 37 near +1 kcal/mol, but the sequences are not necessarily the same for RNA and DNA. The closing base pair is also important; changing from C.G to G.C lowered the stability of several tetraloops in both RNA and DNA. These values will be useful in predicting RNA and DNA secondary structures.     

Heats of thermally induced helix-coil transitions of DNA in aqueous solution

Thermal denaturation of calf thymus DNA at both alkaline and neutral pH values was studied by differential scanning calorimetry. It was shown that the dependence of the enthalpy of transition on pH and salt concentration could be accounted for on the basis of a heat capacity change of +40 cal deg^?1(base pair)^?1. In the pH range between 10.3 and 11.3, a release of 0.6 proton per base pair was calculated from the pH dependence of the melting temperature. The heat effect associated with the release of this proton was calculated to be 5 kcal mole^?1.     

Spectroscopic and thermodynamic studies on the binding of sanguinarine and berberine to triple and double helical DNA and RNA structures

A comparative study on the interaction of sanguinarine and berberine with DNA and RNA triplexes and their parent duplexes was performed, by using a combination of spectrophotometric, UV thermal melting, circular dichroic and thermodynamic techniques. Formation of the DNA and RNA triplexes was confirmed from UV-melting and circular dichroic measurements. The interaction process was characterized by increase of thermal melting temperature, perturbation in circular dichroic spectrum and the typical hypochromic and bathochromic effects in the absorption spectrum. Scatchard analysis indicated that both the alkaloids bound to the triplex and duplex structures in a non-cooperative manner and the binding was stronger to triplexes than to parent duplexes. Thermal melting studies further indicated that sanguinarine stabilized the Hoogsteen base paired third strand of both DNA and RNA triplexes more tightly compared to their Watson-Crick strands, while berberine stabilized the third strand only without affecting the Watson-Crick strand. However, sanguinarine stabilized the parent duplexes while no stabilization was observed with berberine under identical conditions. Circular dichroic studies were also consistent with the observation that perturbations of DNA and RNA triplexes were more compared to their parent duplexes in presence of the alkaloids. Thermodynamic data revealed that binding of sanguinarine and berberine to triplexes (T.AxT and U.AxU) and duplexes (A.T and A.U) showed negative enthalpy changes and positive entropy changes but that of sanguinarine to C.GxC(+) triplex and G.C duplex exhibited negative enthalpy and negative entropy changes. Taken together, these results suggest that both sanguinarine and berberine can bind and stabilize the DNA and RNA triplexes more strongly than their respective parent duplexes.