The energetics of HMG box interactions with DNA: thermodynamic description of the target DNA duplexes

The thermal properties and energetics of formation of 10, 12 and 16 bp DNA duplexes, specifically interacting with the HMG box of Sox-5, have been studied by isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC). DSC studies show that the partial heat capacity of these short duplexes increases considerably prior to the cooperative process of strand separation. Direct extrapolation of the pre and post-transition heat capacity functions into the cooperative transition zone suggests that unfolding/dissociation of strands results in no apparent heat capacity increment. In contrast, ITC measurements show that the negative enthalpy of complementary strand association increases in magnitude with temperature rise, implying that strand association proceeds with significant decrease of heat capacity. Furthermore, the ITC-measured enthalpy of strand association is significantly smaller in magnitude than the enthalpy of cooperative unfolding measured by DSC. To resolve this paradox, the heat effects upon heating and cooling of the separate DNA strands have been measured by DSC. This showed that cooling of the strands from 100 degrees C to -10 degrees C proceeds with significant heat release associated with the formation of intra and inter-molecular interactions. When the enthalpy of residual structure in the strands and the temperature dependence of the heat capacity of the duplexes and of their unfolded strands have been taken into account, the ITC and DSC results are brought into agreement. The analysis shows that the considerable increase in heat capacity of the duplexes with temperature rise is due to increasing fluctuations of their structure (e.g. end fraying and twisting) and this effect obscures the heat capacity increment resulting from the cooperative separation of strands, which in fact amounts to 200(+/-40) JK(-1) (mol bp)(-1). Using this heat capacity increment, the averaged standard enthalpy, entropy and Gibbs energy of formation of fully folded duplexes from fully unfolded strands have been determined at 25 degrees C as -33(+/-2) kJ (mol bp)(-1), -93(+/-4) J K(-1) (mol bp)(-1) and -5.0(+/-0.5) kJ (mol bp)(-1), respectively.     

DNA binding of a non-sequence-specific HMG-D protein is entropy driven with a substantial non-electrostatic contribution

The thermal properties of two forms of the Drosophila melanogaster HMG-D protein, with and without its highly basic 26 residue C-terminal tail (D100 and D74) and the thermodynamics of their non-sequence-specific interaction with linear DNA duplexes were studied using scanning and titration microcalorimetry, spectropolarimetry, fluorescence anisotropy and FRET techniques at different temperatures and salt concentrations. It was shown that the C-terminal tail of D100 is unfolded at all temperatures, whilst the state of the globular part depends on temperature in a rather complex way, being completely folded only at temperatures close to 0 degrees C and unfolding with significant heat absorption at temperatures below those of the gross denaturational changes. The association constant and thus Gibbs energy of binding for D100 is much greater than for D74 but the enthalpies of their association are similar and are large and positive, i.e. DNA binding is a completely entropy-driven process. The positive entropy of association is due to release of counterions and dehydration upon forming the protein/DNA complex. Ionic strength variation showed that electrostatic interactions play an important but not exclusive role in the DNA binding of the globular part of this non-sequence-specific protein, whilst binding of the positively charged C-terminal tail of D100 is almost completely electrostatic in origin. This interaction with the negative charges of the DNA phosphate groups significantly enhances the DNA bending. An important feature of the non-sequence-specific association of these HMG boxes with DNA is that the binding enthalpy is significantly more positive than for the sequence-specific association of the HMG box from Sox-5, despite the fact that these proteins bend the DNA duplex to a similar extent. This difference shows that the enthalpy of dehydration of apolar groups at the HMG-D/DNA interface is not fully compensated by the energy of van der Waals interactions between these groups, i.e. the packing density at the interface must be lower than for the sequence-specific Sox-5 HMG box.     

Thermodynamic comparison of PNA/DNA and DNA/DNA hybridization reactions at ambient temperature

The thermodynamics of 13 hybridization reactions between 10 base DNA sequences of design 5'-ATGCXYATGC-3' with X, Y = A, C, G, T and their complementary PNA and DNA sequences were determined from isothermal titration calorimetry (ITC) measurements at ambient temperature. For the PNA/DNA hybridization reactions, the binding constants range from 1.8 x 10(6)M(-1)for PNA(TT)/DNA to 4.15 x 10(7)M(-1)for PNA(GA)/DNA and the binding enthalpies range from -194 kJ mol(-1)for PNA(CG)/DNA to -77 kJ mol(-1)for PNA(GT)/DNA. For the corresponding DNA/DNA binding reactions, the binding constants range from 2.9 x 10(5)M(-1)for DNA(GT)/DNA to 1.9 x 10(7)M(-1)for DNA(CC)/DNA and the binding enthalpies range from -223 kJ mol(-1)for DNA(CG)/DNA to -124 kJ mol(-1)for DNA(TT)/DNA. Most of the PNA sequences exhibited tighter binding affinities than their corresponding DNA sequences resulting from smaller entropy changes in the PNA/DNA hybridization reactions. van't Hoff enthalpies and extrapolated Delta G values determined from UV melting studies on the duplexes exhibited closer agreement with the ITC binding enthalpies and Delta G values for the DNA/DNA duplexes than for the PNA/DNA duplexes.     

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.     

DNA binding and bending by HMG boxes: energetic determinants of specificity

To clarify the physical basis of DNA binding specificity, the thermodynamic properties and DNA binding and bending abilities of the DNA binding domains (DBDs) of sequence-specific (SS) and non-sequence-specific (NSS) HMG box proteins were studied with various DNA recognition sequences using micro-calorimetric and optical methods. Temperature-induced unfolding of the free DBDs showed that their structure does not represent a single cooperative unit but is subdivided into two (in the case of NSS DBDs) or three (in the case of SS DBDs) sub-domains, which differ in stability. Both types of HMG box, most particularly SS, are partially unfolded even at room temperature but association with DNA results in stabilization and cooperation of all the sub-domains. Binding and bending measurements using fluorescence spectroscopy over a range of ionic strengths, combined with calorimetric data, allowed separation of the electrostatic and non-electrostatic components of the Gibbs energies of DNA binding, yielding their enthalpic and entropic terms and an estimate of their contributions to DNA binding and bending. In all cases electrostatic interactions dominate non-electrostatic in the association of a DBD with DNA. The main difference between SS and NSS complexes is that SS are formed with an enthalpy close to zero and a negative heat capacity effect, while NSS are formed with a very positive enthalpy and a positive heat capacity effect. This indicates that formation of SS HMG box-DNA complexes is specified by extensive van der Waals contacts between apolar groups, i.e. a more tightly packed interface forms than in NSS complexes. The other principal difference is that DNA bending by the NSS DBDs is driven almost entirely by the electrostatic component of the binding energy, while DNA bending by SS DBDs is driven mainly by the non-electrostatic component. The basic extensions of both categories of HMG box play a similar role in DNA binding and bending, making solely electrostatic interactions with the DNA.     

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.     

The DNA bend angle and binding affinity of an HMG box increased by the presence of short terminal arms.

The HMG box of human LEF-1 (hLEF-1, formerly TCF1alpha) has been expressed in four forms: a parent box of 81 amino acids and constructs having either a 10 amino acid C-terminal extension, a 9 amino acid N-terminal extension, or both. These four species have been compared for DNA binding and bending ability using a 28 bp recognition sequence from the TCR alpha-chain enhancer. In the bending assay, whereas the parent box and that with the N-terminal extension bent the DNA by 57/58 degrees, the box extended at the C-terminus bent the DNA by 77/78 degrees, irrespective of the presence or absence of the N-terminal extension. A 6- fold increase in DNA affinity also resulted from addition of both terminal extensions. These observations redefine the functional boundaries of the HMG box. The structure of a mouse LEF-1/DNA complex recently published [Love et al. (1995) Nature 376, 791-795] implies that the higher DNA affinity and in particular the increased bend angle observed are consequences, at least in part, of the C-terminal extension spanning the major groove on the inside of the DNA bend.     

Thermodynamics of the interaction of aluminum ions with DNA: implications for the biological function of aluminum

Aluminum is a known neurotoxic agent and its neurotoxic effects may be due to its binding to DNA. However, the mechanism for the interaction of aluminum ions with DNA is not well understood. Here, we report the application of isothermal titration calorimetry (ITC), fluorescence spectroscopy, and UV spectroscopy to investigate the thermodynamics of the binding of aluminum ions to calf thymus DNA (CT DNA) under various pH and temperature conditions. The binding reaction is driven entirely by a large favorable entropy increase but with an unfavorable enthalpy increase in the pH range of 3.5-5.5 and at all temperatures examined. Aluminum ions show a strong and pH-dependent binding affinity to CT DNA, and a large positive molar heat capacity change for the binding, 1.57 kcal mol(-1) K(-1), demonstrates the burial of the polar surface of CT DNA upon groove binding. The fluorescence of ethidium bromide bound to CT DNA is quenched by aluminum ions in a dynamic way. Both Stern-Volmer quenching constant and the binding constant increase with the increase of the pH values, reaching a maximum at pH 4.5, and decline with further increasing the pH to 5.5. At pH 6.0 and 7.0, aluminum ions precipitate CT DNA completely and no binding of aluminum ions to CT DNA is observed by ITC. Combining the results from these three methods, we conclude that aluminum ions bind to CT DNA with high affinity through groove binding under aluminum toxicity pH conditions and precipitate CT DNA under physiological conditions.     

The DNA sequence specificity of HMG boxes lies in the minor wing of the structure.

To establish the basis of sequence-specific DNA recognition by HMG boxes we separately transferred the minor and major wings from the sequence-specific HMG box of TCF1 alpha into their equivalent position in the non-sequence-specific box 2 of HMG1. Thus chimera THT1 contains the minor wing (of 11 N-terminal and 25 C-terminal residues) from the HMG box of TCF1 alpha and the major wing (the 45 residue central section) from HMG1 box 2, whilst the situation is reversed in chimera HTH1. The structural integrity of the two chimeric proteins was established by CD, NMR and their binding to four-way junction DNA. Gel retardation and circular permutation assays showed that only chimera THT1, containing the TCF1 alpha minor wing, formed a sequence-specific complex and bent the DNA. The bend angle was estimated to be 59 degrees for chimera THT1 and 52 degrees for the HMG box of TCF1 alpha. Our results, in combination with mutagenesis and other data, suggests a model for the DNA binding of HMG boxes in which the N-terminal residues and part of helix 1 contact the minor groove on the outside of a bent DNA duplex.     

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

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

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.     

Forces driving the binding of homeodomains to DNA

Homeodomains are helix-turn-helix type DNA-binding domains that exhibit sequence-specific DNA binding by insertion of their "recognition" alpha helices into the major groove and a short N-terminal arm into the adjacent minor groove without inducing substantial distortion of the DNA. The stability and DNA binding of four representatives of this family, MATalpha2, engrailed, Antennapedia, and NK-2, and truncated forms of the last two lacking their N-terminal arms have been studied by a combination of optical and microcalorimetric methods at different temperatures and salt concentrations. It was found that the stability of the free homeodomains in solution is rather low and, surprisingly, is reduced by the presence of the N-terminal arm for the Antennapedia and NK-2 domains. Their stabilities depend significantly upon the presence of salt: strongly for NaCl but less so for NaF, demonstrating specific interactions with chloride ions. The enthalpies of association of the homeodomains with their cognate DNAs are negative, at 20 degrees C varying only between -12 and -26 kJ/mol for the intact homeodomains, and the entropies of association are positive; i.e., DNA binding is both enthalpy- and entropy-driven. Analysis of the salt dependence of the association constants showed that the electrostatic component of the Gibbs energy of association resulting from the entropy of mixing of released ions dominates the binding, being about twice the magnitude of the nonelectrostatic component that results from dehydration of the protein/DNA interface, van der Waals interactions, and hydrogen bonding. A comparison of the effects of NaCl/KCl with NaF showed that homeodomain binding results in a release not only of cations from the DNA phosphates but also of chloride ions specifically associated with the proteins. The binding of the basic N-terminal arms in the minor groove is entirely enthalpic with a negative heat capacity effect, i.e., is due to sequence-specific formation of hydrogen bonds and hydrophobic interactions rather than electrostatic contacts with the DNA phosphates.     

Thermodynamics of HMGB1 interaction with duplex DNA

The high mobility group protein HMGB1 is a small, highly abundant protein that binds to DNA in a non-sequence-specific manner. HMGB1 consists of 2 DNA binding domains, the HMG boxes A and B, followed by a short basic region and a continuous stretch of 30 glutamate or aspartate residues. Isothermal titration calorimetry was used to characterize the binding of HMGB1 to the double-stranded model DNAs poly(dAdT).(dTdA) and poly(dGdC).(dCdG). To elucidate the contribution of the different structural motifs to DNA binding, calorimetric measurements were performed comparing the single boxes A and B, the two boxes plus or minus the basic sequence stretch (AB(bt) and AB), and the full-length HMGB1 protein. Thermodynamically, binding of HMGB1 and all truncated constructs to duplex DNA was characterized by a positive enthalpy change at 15 degrees C. From the slopes of the temperature dependence of the binding enthalpies, heat capacity changes of -0.129 +/- 0.02 and -0.105 +/- 0.05 kcal mol(-1) K(-1) were determined for box A and full-length HMGB1, respectively. Significant differences in the binding characteristics were observed using full-length HMGB1, suggesting an important role for the acid tail in modulating DNA binding. Moreover, full-length HMGB1 binds differently these two DNA templates: binding to poly(dAdT).(dTdA) was cooperative, had a larger apparent binding site size, and proceeded with a much larger unfavorable binding enthalpy than binding to poly(dGdC).(dCdG).     

Dual binding modes for an HMG domain from human HMGB2 on DNA

High mobility group B (HMGB) proteins contain two HMG box domains known to bind without sequence specificity into the DNA minor groove, slightly intercalating between basepairs and producing a strong bend in the DNA backbone. We use optical tweezers to measure the forces required to stretch single DNA molecules. Parameters describing DNA flexibility, including contour length and persistence length, are revealed. In the presence of nanomolar concentrations of isolated HMG box A from HMGB2, DNA shows a decrease in its persistence length, where the protein induces an average DNA bend angle of 114 +/- 21 degrees for 50 mM Na+, and 87 +/- 9 degrees for 100 mM Na+. The DNA contour length increases from 0.341 +/- 0.003 to 0.397 +/- 0.012 nm per basepair, independent of salt concentration. In 50 mM Na+, the protein does not unbind even at high DNA extension, whereas in 100 mM Na+, the protein appears to unbind only below concentrations of 2 nM. These observations support a flexible hinge model for noncooperative HMG binding at low protein concentrations. However, at higher protein concentrations, a cooperative filament mode is observed instead of the hinge binding. This mode may be uniquely characterized by this high-force optical tweezers experiment.