Identification of pitfalls in the analysis of heat capacity changes of beta-lactoglobulin A

Information on changes in heat capacity (DeltaCp) of proteins upon unfolding is used frequently in literature to understand possible follow-up reactions of protein denaturation, like their aggregation propensity. This thermodynamic property is intrinsic to the protein's architecture and unfolding and should be independent of the approach used to evaluate it. However, for many proteins, the reported values for DeltaCp vary considerably. To identify whether the origin of these discrepancies lies within the experimental approach chosen and/or in the too simplified unfolding models used in the analysis of the data, we choose beta-lactoglobulin A, a relatively small protein, but disputed for its two-state unfolding, and established its DeltaCp from tryptophan fluorescence, near-UV circular dichroism and differential scanning calorimetric measurements. In view of the large variation for the obtained DeltaCp (between 3.2 and 10.1+/-0.8 kJ/(mol K)), it is evident that: (1) the sensitivity of different approaches to the structural changes; (2) irreversibility of unfolding; (3) non-ideal two-state unfolding behaviour need to be considered prior to interpretation. While the first two points can be addressed by using multiple approaches, the applicability of the selected unfolding behaviour for the analysis is often less easy to establish. In this work, we illustrate that by checking the wavelength-dependence used to detect protein conformational changes a tool is provided that gives a direct insight in the validity of the interpretation in these studies. An experimentally validated determination of DeltaCp allows a more proper use for the mechanistic understanding of protein denaturation and its follow-up reactions, avoiding pitfalls in the interpretation.     

Temperature dependence and thermodynamics of Klenow polymerase binding to primed-template DNA

DNA binding of Klenow polymerase has been characterized with respect to temperature to delineate the thermodynamic driving forces involved in the interaction of this polymerase with primed-template DNA. The temperature dependence of the binding affinity exhibits distinct curvature, with tightest binding at 25-30 degrees C. Nonlinear temperature dependence indicates Klenow binds different primed-template constructs with large heat capacity (DeltaCp) values (-870 to -1220 cal/mole K) and thus exhibits large temperature dependent changes in enthalpy and entropy. Binding is entropy driven at lower temperatures and enthalpy driven at physiological temperatures. Large negative DeltaCp values have been proposed to be a 'signature' of site-specific DNA binding, but type I DNA polymerases do not exhibit significant DNA sequence specificity. We suggest that the binding of Klenow to a specific DNA structure, the primed-template junction, results in a correlated thermodynamic profile that mirrors what is commonly seen for DNA sequence-specific binding proteins. Klenow joins a small number of other DNA-sequence independent DNA binding proteins which exhibit unexpectedly large negative DeltaCp values. Spectroscopic measurements show small conformational rearrangements of both the DNA and Klenow upon binding, and small angle x-ray scattering shows a global induced fit conformational compaction of the protein upon binding. Calculations from both crystal structure and solution structural data indicate that Klenow DNA binding is an exception to the often observed correlation between DeltaCp and changes in accessible surface area. In the case of Klenow, surface area burial can account for only about half of the DeltaCp of binding.     

Structural energetics of MgADP binding to the isolated beta subunit of F1-ATPase from thermophilic Bacillus PS3

The energetics of binding of MgADP to the isolated beta subunit of F(1)-ATPase from thermophilic Bacillus (Tbeta) was characterized by high-precision isothermal titration calorimetry. The reaction was enthalpically driven, with a DeltaCp of -36cal(molK)(-1). To gain insight into the molecular basis of this small DeltaCp, we analyzed the changes in accessible surface areas (DeltaASA) between the structures of empty and MgADP-filled beta subunits, extracted from the crystal structure of bovine heart F(1). Consistent with the experimental DeltaCp, the DeltaASA was small (-775A(2)). We used a reported surface area model developed for protein reactions to calculate DeltaCp and DeltaH from DeltaASA, obtaining good agreement with the experimental values. Conversely, using the same model, a DeltaASA of -770A(2) was estimated from experimental DeltaCp and DeltaH for the Tbeta-MgADP complex. Our structural-energetic study indicates that on MgADP binding the isolated Tbeta subunit exhibits intrinsic structural changes similar to those observed in F(1).     

Adenine base unstacking dominates the observed enthalpy and heat capacity changes for the Escherichia coli SSB tetramer binding to single-stranded oligoadenylates.

Isothermal titration calorimetry (ITC) was used to test the hypothesis that the relatively small enthalpy change (DeltaHobs) and large negative heat capacity change (DeltaCp,obs) observed for the binding of the Escherichia coli SSB protein to single-stranded (ss) oligodeoxyadenylates result from the temperature-dependent adenine base unstacking equilibrium that is thermodynamically coupled to binding. We have determined DeltaH1,obs for the binding of 1 mole of each of dT(pT)34, dC(pC)34, and dA(pA)34 to the SSB tetramer (20 mM NaCl at pH 8.1). For dT(pT)34 and dC(pC)34, we found large, negative values for DeltaH1,obs of -75 +/- 1 and -85 +/- 2 kcal/mol at 25 degrees C, with DeltaCp,obs values of -540 +/- 20 and -570 +/- 30 cal mol-1 K-1 (7-50 degrees C), respectively. However, for SSB-dA(pA)34 binding, DeltaH1,obs is considerably less negative (-14 +/- 1 kcal/mol at 25 degrees C), even becoming positive at temperatures below 13 degrees C, and DeltaCp,obs is nearly twice as large in magnitude (-1180 +/- 40 cal mol-1 K-1). These very different thermodynamic properties for SSB-dA(pA)34 binding appear to result from the fact that the bases in dA(pA)34 are more stacked at any temperature than are the bases in dC(pC)34 or dT(pT)34 and that the bases become unstacked within the SSB-ssDNA complexes. Therefore, the DeltaCp,obs for SSB-ssDNA binding has multiple contributions, a major one being the coupling to binding of a temperature-dependent conformational change in the ssDNA, although SSB binding to unstacked ssDNA still has an "intrinsic" negative DeltaCp,0. In general, such temperature-dependent changes in the conformational "end states" of interacting macromolecules can contribute significantly to both DeltaCp,obs and DeltaHobs.     

Structural insight into the kinetics and DeltaCp of interactions between TEM-1 beta-lactamase and beta-lactamase inhibitory protein (BLIP)

In a previous study, we examined thermodynamic parameters for 20 alanine mutants in beta-lactamase inhibitory protein (BLIP) for binding to TEM-1 beta-lactamase. Here we have determined the structures of two thermodynamically distinctive complexes of BLIP mutants with TEM-1 beta-lactamase. The complex BLIP Y51A-TEM-1 is a tight binding complex with the most negative binding heat capacity change (DeltaG = approximately -13 kcal mol(-1) and DeltaCp = approximately -0.8 kcal mol(-1) K(-1)) among all of the mutants, whereas BLIP W150A-TEM-1 is a weak complex with one of the least negative binding heat capacity changes (DeltaG = approximately -8.5 kcal mol(-1) and DeltaCp = approximately -0.27 kcal mol(-1) K(-1)). We previously determined that BLIP Tyr51 is a canonical and Trp150 an anti-canonical TEM-1-contact residue, where canonical refers to the alanine substitution resulting in a matched change in the hydrophobicity of binding free energy. Structure determination indicates a rearrangement of the interactions between Asp49 of the W150A BLIP mutant and the catalytic pocket of TEM-1. The Asp49 of W150A moves more than 4 angstroms to form two new hydrogen bonds while losing four original hydrogen bonds. This explains the anti-canonical nature of the Trp150 to alanine substitution, and also reveals a strong long distance coupling between Trp150 and Asp49 of BLIP, because these two residues are more than 25 angstroms apart. Kinetic measurements indicate that the mutations influence the dissociation rate but not the association rate. Further analysis of the structures indicates that an increased number of interface-trapped water molecules correlate with poor interface packing in a mutant. It appears that the increase of interface-trapped water molecules is inversely correlated with negative binding heat capacity changes.     

Break in the heat capacity change at 303 K for complex binding of netropsin to AATT containing hairpin DNA constructs

Studies performed in our laboratory demonstrated the formation of two thermodynamically distinct complexes on binding of netropsin to a number of hairpin-forming DNA sequences containing AATT-binding regions. These two complexes were proposed to differ only by a bridging water molecule between the drug and the DNA in the lower affinity complex. A temperature-dependent isothermal titration calorimetry (ITC)-binding study was performed using one of these constructs (a 20-mer hairpin of sequence 5'-CGAATTCGTCTCCGAATTCG) and netropsin. This study demonstrated a break in the heat capacity change for the formation of the complex containing the bridging water molecule at approximately 303 K. In the plot of the binding enthalpy change versus temperature, the slope (DeltaCp) was -0.67 kcal mol-1 K-1 steeper after the break at 303 K. Because of the relatively low melting temperature of the 20-mer hairpin (341 K (68 degrees C)), the enthalpy change for complex formation might have included some energy of refolding of the partially denatured hairpin, giving the suggestion of a larger DeltaCp. Studies done on the binding of netropsin to similar constructs, a 24-mer and a 28-mer, with added GC basepairs in the hairpin stem to increase thermal stability, exhibit the same nonlinearity in DeltaCp over the temperature range of from 275 to 333 K. The slopes (DeltaCp) were -0.69 and -0.64 kcal mol-1 K-1 steeper after 303 K for the 24-mer and 28-mer, respectively. This observation strengthens the argument regarding the presence of a bridging water molecule in the lower affinity netropsin/DNA complex. The DeltaCp data seem to infer that because the break in the heat capacity change function for the lower affinity binding occurs at the isoequilibrium temperature for water, water may be included or trapped in the complex. The fact that this break does not occur in the heat capacity change function for formation of the higher affinity complex can similarly be taken as evidence that water is not included in the higher affinity complex.     

Structural reorganization and the cooperative binding of single-stranded telomere DNA in Sterkiella nova

In Sterkiella nova, alpha and beta telomere proteins bind cooperatively with single-stranded DNA to form a ternary alpha.beta.DNA complex. Association of telomere protein subunits is DNA-dependent, and alpha-beta association enhances DNA affinity. To further understand the molecular basis for binding cooperativity, we characterized several possible stepwise assembly pathways using isothermal titration calorimetry. In one path, alpha and DNA first form a stable alpha.DNA complex followed by the addition of beta in a second step. Binding energy accumulates with nearly equal free energy of association for each of these steps. Heat capacity is nonetheless dramatically different, with DeltaCp = -305 +/- 3 cal mol(-1) K(-1) for alpha binding with DNA and DeltaCp = -2010 +/- 20 cal mol(-1) K(-1) for the addition of beta to complete the alpha.beta.DNA complex. By examining alternate routes including titration of single-stranded DNA with a preformed alpha.beta complex, a significant portion of binding energy and heat capacity could be assigned to structural reorganization involving protein-protein interactions and repositioning of the DNA. Structural reorganization probably affords a mechanism to regulate high affinity binding of telomere single-stranded DNA with important implications for telomere biology. Regulation of telomere complex dissociation is thought to involve post-translational modifications in the lysine-rich C-terminal portion of beta. We observed no difference in binding energetics or crystal structure when comparing complexes prepared with full-length beta or a C-terminally truncated form, supporting interesting parallels between the intrinsically disordered regions of histones and this portion of beta.     

Compulsory order of substrate binding to herpes simplex virus type 1 thymidine kinase. A calorimetric study

Isothermal titration calorimetry has been used to investigate the thermodynamic parameters of the binding of thymidine (dT) and ATP to herpes simplex virus type 1 thymidine kinase (HSV1 TK). Binding follows a sequential pathway in which dT binds first and ATP second. The free enzyme does not bind ATP, whose binding site becomes only accessible in the HSV1 TK.dT complex. At pH 7.5 and 25 degrees C, the binding constants are 1.9 x 10(5) m(-1) for dT and 3.9 x 10(6) m(-1) for ATP binding to the binary HSV1 TK.dT complex. Binding of both substrates is enthalpy-driven and opposed by a large negative entropy change. The heat capacity change (DeltaCp) obtained from DeltaH in the range of 10-25 degrees C is -360 cal K(-1) mol(-1) for dT binding and -140 cal K(-1) mol(-1) for ATP binding. These large DeltaCp values are incompatible with a rigid body binding model in which the dT and ATP binding sites pre-exist in the free enzyme. Values of DeltaCp and TDeltaS strongly indicate large scale conformational adaptation of the active site in sequential substrate binding. The conformational changes seem to be more pronounced in dT binding than in the subsequent ATP binding. Considering the crystal structure of the ternary HSV1 TK.dT.ATP complex, a large movement in the dT binding domain and a smaller but substantial movement in the LID domain are proposed to take place when the enzyme changes from the substrate-free, presumably more open and less ordered conformation to the closed and compact conformation of the ternary enzyme-substrate complex.     

Evidence that structural rearrangements and/or flexibility during TCR binding can contribute to T cell activation

While in many cases the half-life of T cell receptor (TCR) binding to a particular ligand is a good predictor of activation potential, numerous exceptions suggest that other physical parameter(s) must also play a role. Accordingly, we analyzed the thermodynamics of TCR binding to a series of peptide-MHC ligands, three of which are more stimulatory than their stability of binding would predict. Strikingly, we find that during TCR binding these outliers show anomalously large changes in heat capacity, an indicator of conformational change or flexibility in a binding interaction. By combining the values for heat capacity (DeltaCp) and the half-life of TCR binding (t(1/2)), we find that we can accurately predict the degree of T cell stimulation. Structural analysis shows significant changes in the central TCR contact residue of the peptide-MHC, indicating that structural rearrangements within the TCR-peptide-MHC interface can contribute to T cell activation.     

Hydration heat capacity of nucleic acid constituents determined from the random network model

The heat capacities of hydration (dCp) of the five nucleic acid bases A, G, C, T, and U, the sugars ribose and deoxyribose, and the phosphate backbone were determined using Monte Carlo simulations and the random network model. Solute-induced changes in the mean length and root mean square angle of hydrogen bonds between hydration shell waters were used to compute dCp for these solutes. For all solutes the dCp is significantly more positive than predicted from accessible surface area (ASA) models of heat capacity. In ASA models, nitrogen, oxygen, and phosphorus atoms are considered as uniformly polar, therefore making a negative contribution to dCp. However, the simulations show that many of these polar atoms are hydrated by water whose hydrogen bonds are less distorted than in bulk, leading to a positive dCp. This is in contrast to the effect of polar groups seen previously in small molecules and amino acids, which increase the water H-bond distortion, giving negative dCp contributions. Our results imply that dCp accompanying DNA dehydration in DNA-ligand and DNA-protein binding reactions may be significantly more negative than previously believed and that dehydration is a significant contributor to the large decrease in heat capacity seen in experiments.     

Energetics of carbohydrate binding to Momordica charantia (bitter gourd) lectin: an isothermal titration calorimetric study

Physico-chemical and carbohydrate binding studies have been carried out on the Momordica charantia (bitter gourd) seed lectin (MCL). The lectin activity is maximal in the pH range 7.4-11.0, but decreases steeply below pH 7.0. The lectin activity is mostly unaffected in the temperature range 4-50 degrees C, but a sharp decrease is seen between 50 and 60 degrees C, which could be correlated to changes in the structure of the protein as seen by circular dichroism and fluorescence spectroscopy. Isothermal titration calorimetric studies show that the tetrameric MCL binds two sugar molecules and the binding constants (Kb), determined at 288.15 K, for various saccharides were found to vary between 7.3 x 10(3) and 1.52 x 10(4)M(-1). The binding reactions for all the saccharides investigated were essentially enthalpy driven, with the binding enthalpies (DeltaHb) at 288.15 K being in the range of -50.99 and -43.39 kJ mol(-1), whereas the contribution to the binding reaction from the entropy of binding was negative, with values of binding entropy (DeltaSb) ranging between -99.2 and -72.0 J mol(-1)K(-1) at 288.15 K. Changes in heat capacity (DeltaCp) for the binding of disaccharides, lactose and lactulose, were significantly larger in magnitude than those obtained for the monosaccharides, methyl-beta-D-galactopyranoside, and methyl-alpha-D-galactopyranoside, and could be correlated reasonably well with the surface areas of these ligands. Enthalpy-entropy compensation was observed for all the sugars studied, suggesting that water structure plays an important role in the overall binding reaction. CD spectroscopy indicates that carbohydrate binding does not lead to significant changes in the secondary and tertiary structures of MCL, suggesting that the carbohydrate binding sites on this lectin are mostly preformed.     

G148-GA3: a streptococcal virulence module with atypical thermodynamics of folding optimally binds human serum albumin at physiological temperatures

The third albumin binding domain of streptococcal protein G strain 148 (G148-GA3) belongs to a novel class of prokaryotic albumin binding modules that is thought to support virulence in several bacterial species. Here, we characterize G148-GA3 folding and albumin binding by using differential scanning calorimetry and isothermal titration calorimetry to obtain the most complete set of thermodynamic state functions for any member of this medically significant module. When buffered at pH 7.0 the 46-amino acid alpha-helical domain melts at 72 degrees C and exhibits marginal stability (15 kJ/mol) at 37 degrees C. G148-GA3 unfolding is characterized by small contributions to entropy from non-hydrophobic forces and a low DeltaCp (1.1 kJ/(deg mol)). Isothermal titration calorimetry reveals that the domain has evolved to optimally bind human serum albumin near 37 degrees C with a binding constant of 1.4 x 10 7 M(-1). Analysis of G148-GA3 thermodynamics suggests that the domain experiences atypically small per residue changes in structural dynamics and heat capacity while transiting between folded and unfolded states.     

Unfolding thermodynamics of Trp-cage, a 20 residue miniprotein, studied by differential scanning calorimetry and circular dichroism spectroscopy

Small proteins provide convenient models for computational studies of protein folding and stability, which are usually compared with experimental data. Until recently, the unfolding of Trp-cage was considered to be a two-state process. However, no direct experimental evidence for this has been presented, and in some cases, the contrary has been suggested. To elucidate a detailed unfolding mechanism, we studied the thermodynamics of unfolding of Trp-cage by differential scanning calorimetry (DSC) and circular dichroism (CD) spectroscopy. The observation that at low temperatures only approximately 90-95% of Trp-cage exists in the native conformation presented an analytical challenge. Nevertheless, it was found that the DSC and CD data can be fitted simultaneously to the same set of thermodynamic parameters. The major uncertainty in such a global fit is the heat capacity change upon unfolding, DeltaCp. This can be circumvented by obtaining DeltaCp directly from the difference between heat capacity functions of the native and unfolded states. Using such an analysis it is shown that Trp-cage unfolding can be represented by a two-state model with the following thermodynamic parameters: Tm = 43.9 +/- 0.8 degrees C, DeltaH(Tm) = 56 +/- 2 kJ/mol, DeltaCp = 0.3 +/- 0.1 kJ/(mol.K). Using these thermodynamic parameters it is estimated that Trp-cage is marginally stable at 25 degrees C, DeltaG(25 degrees C) = 3.2 +/- 0.2 kJ/mol, which is only 30% more than the thermal fluctuation energy at this temperature.     

Linear free-energy model description of the conformational stability of uracil-DNA glycosylase inhibitor A thermodynamic characterization of interaction with denaturant and cold denaturation.

The equilibrium unfolding of uracil DNA glycosylase inhibitor (Ugi), a small acidic protein of molecular mass 9474 Da, has been studied by a combination of thermal-induced and guanidine hydrochloride (GdnCl)-induced denaturation. The analysis of the denaturation data provides a measure of the changes in conformational free energy, enthalpy, entropy and heat capacity DeltaCp that accompany the equilibrium unfolding of Ugi over a wide range of temperature and GdnCl concentration. The unfolding of Ugi is a simple two-state, reversible process. The protein undergoes both low-temperature and high-temperature unfolding even in the absence of GdnCl but more so in the presence of denaturant. The data are consistent with the linear free-energy model and with a temperature independent DeltaCp over the large temperature range of unfolding. The small DeltaCp (6.52 kJ.mol-1.K-1) for the unfolding of Ugi, is perhaps a reflection of a relatively small, buried hydrophobic core in the folded form of this small monomeric protein. Despite a relatively low value of DeltaG(H2O), 7.40 kJ.mol-1 at pH 8.3, Ugi displays considerable stability with the temperature of maximum stability being 301.6 K.     

Interaction of human apolipoprotein A-I with model membranes exhibiting lipid domains

Several mechanisms for cell cholesterol efflux have been proposed, including membrane microsolubilization, suggesting that the existence of specific domains could enhance the transfer of lipids to apolipoproteins. In this work isothermal titration calorimetry, circular dichroism spectroscopy, and two-photon microscopy are used to study the interaction of lipid-free apolipoprotein A-I (apoA-I) with small unilamellar vesicles (SUVs) of 1-palmitoyl, 2-oleoyl phosphatidylcholine (POPC) and sphingomyelin (SM), with and without cholesterol. Below 30 degrees C the calorimetric results show that apoA-I interaction with POPC/SM SUVs produces an exothermic reaction, characterized as nonclassical hydrophobic binding. The heat capacity change (DeltaCp degrees ) is small and positive, whereas it was larger and negative for pure POPC bilayers, in the absence of SM. Inclusion of cholesterol in the membranes induces changes in the observed thermodynamic pattern of binding and counteracts the formation of alpha-helices in the protein. Above 30 degrees C the reactions are endothermic. Giant unilamellar vesicles (GUVs) of identical composition to the SUVs, and two-photon fluorescence microscopy techniques, were utilized to further characterize the interaction. Fluorescence imaging of the GUVs indicates coexistence of lipid domains under 30 degrees C. Binding experiments and Laurdan generalized-polarization measurements suggest that there is no preferential binding of the labeled apoA-I to any particular domain. Changes in the content of alpha-helix, binding, and fluidity data are discussed in the framework of the thermodynamic parameters.     

Thermodynamics of nucleotide binding to the chaperonin GroEL studied by isothermal titration calorimetry: evidence for noncooperative nucleotide binding.

We characterized the thermodynamics of binding reactions of nucleotides ADP and ATPgammaS (a nonhydrolyzable analog of ATP) to GroEL in a temperature range of 5 degrees C to 35 degrees C by isothermal titration calorimetry. Analysis with a noncooperative binding model has shown that the bindings of nucleotides are driven enthalpically with binding constants of 7x103 M-1 and 4x104 M-1 for ADP and ATPgammaS, respectively, at 26 degrees C and that the heat capacity change DeltaCp is about 100 cal/mol.K for both the nucleotides. The stoichiometries of binding were about 8 and 9 molecules for ADP and ATPgammaS, respectively, per GroEL tetradecamer at 5 degrees C, and both increased with temperature to reach about 14 (ADP) and 12 (ATPgammaS) for both nucleotides at 35 degrees C. The absence of initial increase of binding heat as well as Hill coefficient less than 1.2, which were obtained from the fitting to the model curve by assuming positive cooperativity, showed that there was virtually no positive cooperativity in the nucleotide bindings. Incorporating a difference in affinity for the nucleotide (ADP and ATPgammaS) between the two rings of GroEL into the noncooperative binding model improved the goodness of fitting and the difference in the affinity increased with decreasing temperature.     

Heat capacity effects in protein folding and ligand binding: a re-evaluation of the role of water in biomolecular thermodynamics

Large "anomalous" heat capacity (DeltaC(p)) effects are a common feature of the thermodynamics of biomolecular interactions in aqueous solution and, as a result of the improved facility for direct calorimetric measurements, there is a growing body of experimental data for such effects in protein folding, protein-protein and protein-ligand interactions. Conventionally such heat capacity effects have been ascribed to hydrophobic interactions, and there are some remarkably convincing demonstrations of the usefulness of this concept. Nonetheless, there is also increasing evidence that hydrophobic interactions are not the only possible source of such effects. Here we re-evaluate the possible contributions of other interactions to the heat capacity changes to be expected for cooperative biomolecular folding and binding processes, with particular reference to the role of hydrogen bonding and solvent water interactions. Simple models based on the hydrogen-bonding propensity of water as a function of temperature give quantitative estimates of DeltaC(p) that compare well with experimental observations for both protein folding and ligand binding. The thermodynamic contribution of bound waters in protein complexes is also estimated. The prediction from simple lattice models is that trapping of water in a complex should give more exothermic binding (DeltaDeltaH-6 to -12 kJ mol(-1)) with lower entropy (DeltaDeltaS(0) approximately -11 J mol(-1) K(-1)) and more negative DeltaC(p) (by about -75 J mol(-1) K(-1)) per water molecule. More generally, it is clear that significant DeltaC(p) effects are to be expected for any macromolecular process involving a multiplicity of cooperative weak interactions of whatever kind.     

Experiment-guided thermodynamic simulations on reversible two-state proteins: implications for protein thermostability

Here, we perform protein thermodynamic simulations within a set of boundary conditions, effectively blanketing the experimental data. The thermodynamic parameters, melting temperature (TG), enthalpy change at the melting temperature (DeltaHG) and heat capacity change (DeltaCp) were systematically varied over the experimentally observed ranges for small single domain reversible two-state proteins. Parameter sets that satisfy the Gibbs-Helmholtz equation and yield a temperature of maximal stability (TS) around room temperature were selected. The results were divided into three categories by arbitrarily chosen TG ranges. The TG ranges in these categories correspond to typical values of the melting temperatures observed for the majority of the proteins from mesophilic, thermophilic and hyperthermophilic organisms. As expected, DeltaCp values tend to be high in mesophiles and low in hyperthermophiles. An increase in TG is accompanied by an up-shift and broadening of the protein stability curves, however, with a large scatter. Furthermore, the simulations reveal that the average DeltaHG increases with TG up to approximately 360 K and becomes constant thereafter. DeltaCp decreases with TG with different rates before and after approximately 360 K. This provides further justification for the separate grouping of proteins into thermophiles and hyperthermophiles to assess their thermodynamic differences. This analysis of the Gibbs-Helmholtz equation has allowed us to study the interdependence of the thermodynamic parameters TG, DeltaHG and DeltaCp and their derivatives in a more rigorous way than possible by the limited experimental protein thermodynamics data available in the literature. The results provide new insights into protein thermostability and suggest potential strategies for its manipulation.     

A thermodynamic molecular switch in biological systems: ribonuclease S' fragment complementation reactions

It is well known that essentially all biological systems function over a very narrow temperature range. Most typical macromolecular interactions show DeltaH degrees (T) positive (unfavorable) and a positive DeltaS degrees (T) (favorable) at low temperature, because of a positive (DeltaCp degrees /T). Because DeltaG degrees (T) for biological systems shows a complicated behavior, wherein DeltaG degrees (T) changes from positive to negative, then reaches a negative value of maximum magnitude (favorable), and finally becomes positive as temperature increases, it is clear that a deeper-lying thermodynamic explanation is required. This communication demonstrates that the critical factor is a temperature-dependent DeltaCp degrees (T) (heat capacity change) of reaction that is positive at low temperature but switches to a negative value at a temperature well below the ambient range. Thus the thermodynamic molecular switch determines the behavior patterns of the Gibbs free energy change and hence a change in the equilibrium constant, K(eq), and/or spontaneity. The subsequent, mathematically predictable changes in DeltaH degrees (T), DeltaS degrees (T), DeltaW degrees (T), and DeltaG degrees (T) give rise to the classically observed behavior patterns in biological reactivity, as may be seen in ribonuclease S' fragment complementation reactions.     

Thermodynamics of glutathione binding to the tyrosine 7 to phenylalanine mutant of glutathione S-transferase from Schistosoma japonicum

The binding of glutathione (GSH) to the tyrosine 7 to phenylalanine mutant of Schistosoma japonicum glutathione S-transferase (SjGST-Y7F) has been studied by isothermal titration calorimetry (ITC). At pH 6.5 and 25 °C this mutant shows a higher affinity for glutathione than wild type enzyme despite an almost complete loss of activity in the presence of 1-chloro-2,4-dinitrobenzene (CDNB) as second substrate. The enthalpy change upon binding of GSH is more negative for the mutant than for the wild type GST (SjGST). Changes in accessible solvent areas (ASA) have been calculated based on enthalpy and heat capacity changes. ASA values indicated the burial of apolar surfaces of protein and ligand upon binding. A more negative ΔC_p value has been obtained for the mutant enzyme, suggesting a more hydrophobic interaction, as may be expected from the change of a tyrosine residue to phenylalanine.