Binding of the recombinant proteinase inhibitor eglin c from leech Hirudo medicinalis to human leukocyte elastase, bovine alpha-chymotrypsin and subtilisin Carlsberg: thermodynamic study

The effect of pH and temperature on the apparent association equilibrium constant (Ka) for the binding of the recombinant proteinase inhibitor eglin c from leech Hirudo medicinalis to human leukocyte elastase (EC 3.4.21.37), bovine alpha-chymotrypsin (EC 3.4.21.1) and subtilisin Carlsberg (EC 3.4.21.14) has been investigated. On lowering the pH from 9.5 to 4.5, values of Ka for eglin c binding to the serine proteinases considered decrease thus reflecting the acid-pK shift of the invariant histidyl catalytic residue (His57 in human leukocyte elastase and bovine alpha-chymotrypsin, and His64 in subtilisin Carlsberg) from congruent to 6.9, in the free enzymes, to congruent to 5.1, in the enzyme:inhibitor adducts. At pH 8.0, values of the apparent thermodynamic parameters for eglin c binding are: human leukocyte elastase - Ka = 1.0 x 10(10) M-1, delta G phi = -13.4 kcal/mol, delta H phi = +1.8 kcal/mol, and delta S phi = +52 entropy units; bovine alpha-chymotrypsin -Ka = 5.0 x 10(9) M-1, delta G phi = -13.0 kcal/mol, delta H phi = +2.0 kcal/mol, and delta S phi = +51 entropy units; and subtilisin Carlsberg - Ka = 6.6 x 10(9) M-1, delta G phi = -13.1 kcal/mol, delta H phi = +2.0 kcal/mol, and delta S phi = +51 entropy units (values of Ka, delta G phi and delta S phi were obtained at 21 degrees C; values of delta H phi were temperature independent over the range explored, i.e. between 10 degrees C and 40 degrees C; 1 kcal = 4184J).(ABSTRACT TRUNCATED AT 250 WORDS)     

Thermodynamics of ribonuclease T1 denaturation.

Differential scanning calorimetry has been used to investigate the thermodynamics of denaturation of ribonuclease T1 as a function of pH over the pH range 2-10, and as a function of NaCl and MgCl2 concentration. At pH 7 in 30 mM PIPES buffer, the thermodynamic parameters are as follows: melting temperature, T1/2 = 48.9 +/- 0.1 degrees C; enthalpy change, delta H = 95.5 +/- 0.9 kcal mol-1; heat capacity change, delta Cp = 1.59 kcal mol-1 K-1; free energy change at 25 degrees C, delta G degrees (25 degrees C) = 5.6 kcal mol-1. Both T1/2 = 56.5 degrees C and delta H = 106.1 kcal mol-1 are maximal near pH 5. The conformational stability of ribonuclease T1 is increased by 3.0 kcal/mol in the presence of 0.6 M NaCl or 0.3 M MgCl2. This stabilization results mainly from the preferential binding of cations to the folded conformation of the protein. The estimates of the conformational stability of ribonuclease T1 from differential scanning calorimetry are shown to be in remarkably good agreement with estimates derived from an analysis of urea denaturation curves.     

Thermodynamic analysis of stacking hybridization of oligonucleotides with DNA template.

Contiguous stacking hybridization of oligodeoxyribonucleotides with DNA as template was investigated using three types of complexes: oligonucleotide contiguously stacked with the stem of the preformed minihairpin (complexes I), oligonucleotide tandems containing two (complexes II) or three (complexes III) short oligomers with a common DNA template. Enthalpy Delta H degrees and entropy Delta S degrees of the coaxial stacking of adjacent duplexes were determined for GC/G*pC, GT/A*pC, AC/G*pT, AT/A*pT, CT/A*pG, AG/C*pT, AA/T*pT and TT/A*pA nicked (*) dinucleotide base pairs. The maximal efficiency of co-operative interaction was found for the GC/G*pC interface (Delta G degrees(NN/N*pN)=-2.7 kcal/mol) and the minimal one for the AA/T*pT interface (Delta G degrees(NN/N*pN)=-1.2 kcal/mol) at 37 degrees C. As a whole, the efficiency of the base pairs interaction Delta G degrees(NN/N*pN) in the nick is not lower than that within the intact DNA helix (Delta G degrees(NN/NN)).These observed Delta G degrees(NN/N*pN) values are proposed may include the effect of the partial removal of fraying at the adjacent helix ends additionally to the effect of the direct stacking of the terminal base pairs in the duplex junction (Delta G degrees(NN/NN). The thermodynamic parameters have been found to describe adequately the formation of all tandem complexes of the II and III types with oligonucleotides of various length and hybridization properties. The performed thermodynamic analysis reveals features of stacking oligonucleotide hybridization which allow one to predict the temperature dependence of association of oligonucleotides and the DNA template within tandem complexes as well as to determine optimal concentration for formation of these complexes characterized by high co-operativity level.     

Comparison of the base pairing properties of a series of nitroazole nucleobase analogs in the oligodeoxyribonucleotide sequence 5'-d(CGCXAATTYGCG)-3'.

The nucleoside analogs 1-(2'-deoxy-beta-D-ribofuranosyl)- 3-nitropyrrole (9), 1-(2'-deoxy-beta-D-ribofuranosyl)-4-nitropyrazole (10), 1-(2'-deoxy-beta-D-ribofuranosyl)-4-nitroimidazole (11) and 1-(2'-deoxy-beta-D-ribofuranosyl)-5-nitroindole (21) were incorporated into the oligonucleotide 5'-d(CGCXAATTYGCG)-3'in the fourth position from the 5'-end. Procedures for synthesis of two of the nitroazole nucleosides, 10 and 11, were developed for this study. Each of the nitroazoles was converted into a 3'-phosphoramidite for oligonucleotide synthesis by conventional automated protocols. Four oligonucleotides were synthesized for each modified nucleoside in order to obtain duplexes in which each of the four natural bases was placed opposite (position 9) the nitroazole. In order to assess the role of the nitro group on base stacking interaction, sequences were also synthesized in which the fourth base was 1-(2'-deoxy-beta-D-ribofuranosyl)pyrazole. Corresponding sequences containing an abasic site, as well as sequences containing inosine, were synthesized for comparison. Thermal melting studies yielded T m values and thermodynamic parameters. Each nucleoside analog displayed a unique pattern of base pairing preferences. The least discriminating analog was 3-nitropyrrole, for which T m values differed by 5 degrees C and Delta G 25 degrees C ranged from -6.1 to -6.5 kcal/mol. 5-Nitroindole gave duplexes with significantly higher thermal stability, with Tm values varying from 35.0 to 46.5 degrees C and -Delta G 25 degrees C ranging from 7.7 to 8.5 kcal/mol. Deoxyinosine (22), a natural analog which has found extensive use as a universal nucleoside, is far less non-discriminating than any of the nitroazole derivatives. Tm values ranged from 35.4 degrees C when paired with G to 62.3 degrees C when paired with C. The significance of the nitro substituent was determined by comparison of the base pairing properties of a simple azole nucleoside, 1-(2'-deoxy-beta-D-ribofuranosyl)pyrazole (12). The pyrazole-containing sequences melt at 10-20 degrees C lower than the corresponding nitropyrazole-containing sequences. On average, the pyrazole-containing sequences were equivalent in stability (average Delta G = -4.8 kcal/mol) to the sequences containing an abasic site (average Delta G = -4.7 kcal/mol).     

Physicochemical studies of processes involving potential photodynamic drugs on route to their targets: self-aggregation and membrane-binding of Zn-hematoporphyrin

The thermodynamics of zinc hematoporphyrin (ZnHP) dimerization and ZnHP-membrane binding were studied. The dimerization equilibrium was determined over the temperature range 19-40 degrees C, using fluorometric techniques. The dimerization constant obtained at 37 degrees C (neutral pH in phosphate-buffered saline) is 4.6 (+/- 0.6) X 10(4) M-1. The dimerization was found to decrease with temperature over the range 19-36 degrees C, the data allowing the extraction of the following thermodynamic parameters for the temperature range 19-31 degrees C: delta G0 = -9.3 kcal/mol, delta H0 = -7.4 kcal/mol, delta S0 = -6.4 eu. For temperatures above 36 degrees C the dimerization was found to be temperature independent, giving the following parameters: delta G0 = -6.6 kcal/mol, delta H0 = 0 kcal/mol, delta S0 = 21.2 eu. On the basis of the data the case is made for the existence of two types of ZnHP dimers, differing in the location of the fifth Zn2+ ligand and in the nature of the contribution of the solvent to the dimerization. For the membrane binding, large unilamellar liposomes served to model biological membranes. The binding of ZnHP to the liposomes was found to be similar, quantitatively, to the corresponding metal-free molecule, namely, fitting a case of one type of site and giving a binding constant of 1600 +/- 160 M (neutral pH and 37 degrees C) which is independent of the length of the porphyrin-liposome.     

Entropy drives integrin alphaIIbbeta3:echistatin binding--evidence from surface plasmon resonance spectroscopy.

This investigation examined the molecular mechanisms that enable the alphaIIbbeta3 integrin to bind efficiently, tightly, and selectively to echistatin, an RGD disintegrin. We used surface plasmon resonance spectroscopy to measure the rate, extent, and stability of complexes formed between micellar alphaIIbbeta3 and recombinant echistatin (rEch) mutants, immobilized on the surface of a biosensor chip. alphaIIbbeta3 bound readily and tightly to wild-type RGD-rEch and RGDF-rEch but not to RGA-rEch or AGD-rEch, demonstrating that both of those charged moieties contribute to integrin recognition. van't Hoff analysis of the temperature dependence of the RGD-rEch K d data yielded an unfavorable enthalpy change, Delta H degrees = 14 +/- 3 kcal/mol, offset by a favorable entropy term, TDelta S degrees = 23 +/- 3 kcal/mol. Eyring analysis of the temperature dependence of the kinetic parameters yielded Delta H a degrees (++) = 9 +/- 2 kcal/mol and TDelta S a degrees (++) = -4 +/- 2 kcal/mol, indicating that a substantial activation enthalpy barrier and a smaller activation entropy hinder assembly of the encounter complex. Thus, equilibrium thermodynamic data demonstrate that entropy is the dominant factor stabilizing integrin:echistatin binding, while transition-state thermodynamic parameters indicate that the rate of complex formation is enthalpy-limited. When electrostatic contacts are the major source of receptor:ligand stability, theory and experiment indicate that entropy-favorable ion-pair desolvation often provides the driving force for molecular recognition.     

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.     

Inhibition of lactate and malate dehydrogenase by dyes containing vinylsulfone and chlorotriazine groups

The difference spectra of lactate and malate dehydrogenase complexes with four native dyes containing vinylsulfonic and triazinic groups (light-resistant yellow 2KT, red-violet 2KT, etc.) were monitored in 0.1 M phosphate buffer pH 8.2 at 20 degrees C. The dissociation constants were calculated from the spectral data. The most stable complexes were lactate dehydrogenase--light-resistant yellow 2KT and malate dehydrogenase--light-resistant yellow 2KT ones. The values of delta H degree = 5.75 kcal/mole and standard thermodynamic parameters, delta G degree = -6.5 kcal/mole and delta S degree = 41.2 e. u., were calculated from the values of association constants for temperature dependence. The thermodynamic characteristics confirmed the key role of hydrophobic interactions in lactate dehydrogenase--reactive dye complex formation. All the dyes under study competitively inhibit lactate and malate oxidation by the corresponding dehydrogenases. The inhibition constants of both enzymes by the four dyes were determined at 20 degrees C in 0.1 M phosphate buffer pH 8.2. Light-resistant yellow 2KT appeared to be the most effective inhibitor of the enzymes.     

Thermodynamics of phospholipid bilayer assembly

Thermodynamic properties of bilayer assembly have been obtained from measurements of the solubility of the sodium salt of dimyristoylphosphatidylglycerol (DMPG) in water. The standard free energy of bilayer assembly delta G degree a is shown to be RT 1n Xs + zF psi 0 where Xs is the mole fraction of dissolved lipid, F is the Faraday constant, z is the valence of the counterion (Na+), and psi 0 is the electrical double-layer potential of the ionized bilayer. The function d 1n Xs/dT was found to be discontinuous at 24 degrees C, the gel-liquid-crystal transition temperature (Tm) for DMPG. This function was unaffected when solubilities were measured in 0.001 M NaCl solutions; thus, psi 0 is constant in the experimental temperature interval (4-40 degrees C). Using a value of psi 0 = -180 mV [Eisenberg et al. (1979) Biochemistry 18, 5213-5223], and the temperature dependence of delta G degrees a, values for delta H degrees a and delta S degree a at 24 degrees C were calculated for the gel and liquid-crystal states of DMPG. For the gel, delta H degrees a and T delta S a are -26.2 and 12.7 kcal/mol, respectively; for the liquid-crystal, delta H degrees a and T delta S degrees a are -19.2 and -5.7 kcal/mol, respectively. The calculated value for the latent heat of the gel-liquid-crystal transition is 7 kcal/mol, in agreement with calorimetric measurements.     

Thermal stability of DNA.

Tij and Delta Hij for stacking of pair i upon j in DNA have been obtained over the range 0.034-0.114 M Na+from high-resolution melting curves of well-behaved synthetic tandemly repeating inserts in recombinant pN/MCS plasmids. Results are consistent with neighbor-pair thermodynamic additivity, where the stability constant, sij , for different domains of length N depend quantitatively on the product of stability constants for each individual pair in domains, sijN . Unit transition enthalpies with average errors less than +/-5%, were determined by analysis of two-state equilibria associated with the melting of internal domains and verified from variations of Tij with [Na+]. Enthalpies increase with Tij , in close agreement with the empirical function: Delta Hij = 52.78@ Tij - 9489, and in parallel with a smaller increase in Delta Sij . Delta Hij and Delta Sij are in good agreement with the results of an extensive compilation of published Delta Hcal and Delta Scal for synthetic and natural DNAs. Neighbor-pair additivity was also observed for (dA@dT)-tracts at melting temperatures; no evidence could be detected of the familiar and unusual structural features that characterize tracts at lower temperatures. The energetic effects of loops were determined from the melting behavior of repeating inserts installed between (G+C)-rich barrier domains in the pN/MCS plasmids. A unique set of values for the cooperativity, loop exponent and stiffness parameters were found applicable to internal domains of all sizes and sequences. Statistical mechanical curves calculated with values of Tij([Na+]) , Delta Hij and these loop parameters are in good agreement with observation.     

1H NMR assignment and melting temperature study of cis-syn and trans-syn thymine dimer containing duplexes of d(CGTATTATGC).d(GCATAATACG)

The preparation and spectroscopic characterization of duplex decamers containing site-specific cis-syn and trans-syn thymine dimers are described. Three duplex decamers, d(CGTATTATGC).d(GCATAATACG), d(CGTAT[c,s]TATGC).d(GCATAATACG), and d(CGTAT[t,s]TATGC).d(GCATAATACG), were prepared by solid-phase phosphoramidite synthesis utilizing cis-syn and trans-syn cyclobutane thymine dimer building blocks (Taylor et al., 1987; Taylor & Brockie, 1988). NMR spectra (500 MHz 2D 1H and 202 MHz 1D 31P) were obtained in "100%" D2O at 10 degrees C, and 1D exchangeable 1H spectra were obtained in 10% D2O at 10 degrees C. 1H NMR assignments for H5, H6, H8, CH3, H1', H2', and H2" were made on the basis of standard sequential NOE assignment strategies and verified in part by DQF COSY data. Comparison of the chemical shift data suggests that the helix structure is perturbed more to the 3'-side of the cis-syn dimer and more to the 5'-side of the trans-syn dimer. Thermodynamic parameters for the helix in equilibrium coil equilibrium were obtained by two-state, all or none, analysis of the melting behavior of the duplexes. Analysis of the temperature dependence of the T5CH3 1H NMR signal gave delta H = 44 +/- 4 kcal and delta S = 132 +/- 13 eu for the trans-syn duplex. Analysis of the concentration and temperature dependence of UV spectra gave delta H = 64 +/- 6 kcal and delta S = 178 +/- 18 eu for the parent duplex and delta H = 66 +/- 7 kcal and delta S = 189 +/- 19 eu for cis-syn duplex. It was concluded that photodimerization of the dTpdT unit to give the cis-syn product causes little perturbation of the DNA whereas dimerization to give the trans-syn product causes much greater perturbation, possibly in the form of a kink or dislocation at the 5'-side of the dimer.     

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)     

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.     

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

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

The vnd/NK-2 homeodomain: thermodynamics of reversible unfolding and DNA binding for wild-type and with residue replacements H52R and H52R/T56W in helix III

The conformational stabilities of the vnd (ventral nervous system defective)/NK-2 homeodomain [HD(wt); residues 1-80 that encompass the 60-residue homeodomain] and those harboring mutations in helix III of the DNA recognition site [HD(H52R) and HD(H52R/T56W)] have been investigated by differential scanning calorimetry (DSC) and ellipticity changes at 222 nm. Thermal unfolding reactions at pH 7.4 are reversible and repeatable in the presence of 50-500 mM NaCl with DeltaC(p) = 0.52 +/- 0.04 kcal K(-1) mol(-1). A substantial stabilization of HD(wt) is produced by 50 mM phosphate or by the addition of 100-500 mM NaCl to 50 mM Hepes, pH 7.4, buffer (from T(m) = 35.5 degrees C to T(m) 43-51 degrees C; DeltaH(vH) congruent with 47 +/- 5 kcal mol(-1)). The order of stability is HD(H52R/T56W) > HD(H52R) > HD(wt), irrespective of the anions present. Progress curves for ellipticity changes at 222 nm as a function of increasing temperature are fitted well by a two-state unfolding model, and the cooperativity of secondary structure changes is greater for mutant homeodomains than for HD(wt) and also is increased by adding 100 mM NaCl to Hepes buffer. A 33% quench of the intrinsic tryptophanyl residue fluorescence of HD(wt) by phosphate binding (K(D)' = 2.6 +/- 0.3 mM phosphate) is reversed approximately 60% by DNA binding. Thermodynamic parameters for vnd/NK-2 homeodomain proteins binding sequence-specific 18 bp DNA have been determined by isothermal titration calorimetry (10-30 degrees C). Values of DeltaC(p) are +0.25, -0.17, and -0.10 +/- 0.04 kcal K(-1) mol(-1) for HD(wt), HD(H52R), and HD(H52R/T56W) binding duplex DNA, respectively. Interactions of homeodomains with DNA are enthalpically controlled at 298 K and pH 7.4 with corresponding DeltaH values of -6.6 +/- 0.5, -10.8 +/- 0.1, and -9.0 +/- 0.6 kcal mol(-1) and DeltaG' values of -11.0 +/- 0.1, -11.0 +/- 0.1, and -11.3 +/- 0.3 kcal mol(-1) with a binding stoichiometry of 1.0 +/- 0.1. Thermodynamic parameters for DNA binding are not predicted from homeodomain structural changes that occur upon complexing to DNA and must reflect also solvent and possibly DNA rearrangements.     

Disulfide bond effects on protein stability: designed variants of Cucurbita maxima trypsin inhibitor-V

Attempts to increase protein stability by insertion of novel disulfide bonds have not always been successful. According to the two current models, cross-links enhance stability mainly through denatured state effects. We have investigated the effects of removal and addition of disulfide cross-links, protein flexibility in the vicinity of a cross-link, and disulfide loop size on the stability of Cucurbita maxima trypsin inhibitor-V (CMTI-V; 7 kD) by differential scanning calorimetry. CMTI-V offers the advantage of a large, flexible, and solvent-exposed loop not involved in extensive intra-molecular interactions. We have uncovered a negative correlation between retention time in hydrophobic column chromatography, a measure of protein hydrophobicity, and melting temperature (T(m)), an indicator of native state stabilization, for CMTI-V and its variants. In conjunction with the complete set of thermodynamic parameters of denaturation, this has led to the following deductions: (1) In the less stable, disulfide-removed C3S/C48S (Delta Delta G(d)(50 degrees C) = -4 kcal/mole; Delta T(m) = -22 degrees C), the native state is destabilized more than the denatured state; this also applies to the less-stable CMTI-V* (Delta Delta G(d)(50 degrees C) = -3 kcal/mole; Delta T(m) = -11 degrees C), in which the disulfide-containing loop is opened by specific hydrolysis of the Lys(44)-Asp(45) peptide bond; (2) In the less stable, disulfide-inserted E38C/W54C (Delta Delta G(d)(50 degrees C) = -1 kcal/mole; Delta T(m) = +2 degrees C), the denatured state is more stabilized than the native state; and (3) In the more stable, disulfide-engineered V42C/R52C (Delta Delta G(d)(50 degrees C) = +1 kcal/mole; Delta T(m) = +17 degrees C), the native state is more stabilized than the denatured state. These results show that a cross-link stabilizes both native and denatured states, and differential stabilization of the two states causes either loss or gain in protein stability. Removal of hydrogen bonds in the same flexible region of CMTI-V resulted in less destabilization despite larger changes in the enthalpy and entropy of denaturation. The effect of a cross-link on the denatured state of CMTI-V was estimated directly by means of a four-state thermodynamic cycle consisting of native and denatured states of CMTI-V and CMTI-V*. Overall, the results show that an enthalpy-entropy compensation accompanies disulfide bond effects and protein stabilization is profoundly modulated by altered hydrophobicity of both native and denatured states, altered flexibility near the cross-link, and residual structure in the denatured state.     

Thermodynamic stability of DNA tandem mismatches

The thermodynamics of nine hairpin DNAs were evaluated using UV-monitored melting curves and differential scanning calorimetry (DSC). Each DNA has the same five-base loop and a stem with 8-10 base pairs. Five of the DNAs have a tandem mismatch in the stem, while four have all base pairs. The tandem mismatches examined (ga/ga, aa/gc, ca/gc, ta/ac, and tc/tc) spanned the range of stability observed for this motif in a previous study of 28 tandem mismatches. UV-monitored melting curves were obtained in 1.0 M Na(+), 0.1 M Na(+), and 0.1 M Na(+) with 5 mM Mg(2+). DSC studies were conducted in 0.1 M Na(+). Transition T(m) values were unchanged over a 50-fold range of strand concentration. Model-independent enthalpy changes (DeltaH degrees ) evaluated by DSC were in good agreement (+/-8%) with enthalpy values determined by van't Hoff analyses of the melting curves in 0.1 M Na(+). The average heat capacity change (DeltaC(p)) associated with the hairpin to single strands transitions was estimated from plots of DeltaH degrees and DeltaS degrees with T(m) and ln T(m), respectively, and from profiles of DSC curves. The average DeltaC(p) values (113 +/- 9 and 42 +/- 27 cal x K(-1) x mol(-1) of bp), were in the range of values reported in previous studies. Consideration of DeltaC(p) produced large changes in DeltaH degrees and DeltaS degrees extrapolated from the transition region to 37 degrees C and smaller but significant changes to free energies. The loop free energy of the five tandem mismatches at 37 degrees C varied over a range of approximately 4 kcal x mol(-1) for each solvent.     

Thermodynamics of glucagon aggregation.

Heats of dilution of concentrated glucagon solutions have been measured calorimetrically at 10 and 25 degrees C in 0.2 M potassium phosphate buffer of pH 10.6. Analysis of the data in terms of a monomer-trimer equilibrium gives the following thermodynamic parameters for the association reaction at 25 degrees C: delta G degrees = 7.34 kcal/mol of trimer, delta H degrees = -31.2 kcal/mol, deltaS degrees = -80 cal/(K mol), deltaCp = 430 cal/(K mol). The sensitivity of heat of dilution data to the association constant and stoichiometry of the reaction is discussed.