Statistical thermodynamic analysis of differential scanning calorimetry data: structural deconvolution of heat capacity function of proteins

Energetics of β-lactoglobulin-flavor compounds interactions

Interaction of bovine β-lactoglobulin (BLG) with several flavor compounds (FC) (2-methylpyrazine, vanillin, 2-acetylpyridine, 2- and 3-acetylthiophene, methyl isoamyl ketone, heptanone, octanone, and nonanone) was studied by high-sensitivity differential scanning calorimetry. The denaturation temperature, enthalpy, and heat capacity increment were determined at different FC concentrations. It was found that the denaturation temperature and heat capacity increment do not depend on the FC concentration, while the denaturation enthalpy decreases linearly with the FC concentration. These thermodynamic effects disclose the preferential FC binding to the unfolded form of BLG. By the obtained calorimetric data, the free energies of FC binding vs. the FC concentrations were calculated. These dependences were shown to be linear. Their slope relates closely to the overall FC affinity for the unfolded BLG in terms of the Langmuir binding model. The overall BLG affinity for FC varies from 20 M-1 (2-methylpyrazine) up to 360 M-1(nonanone). The maximal stoichiometry of the BLG-FC complexes was roughly estimated as a ratio of the length of the unfolded BLG to the molecular length of FC. Using these estimates, the apparent BLG-FC binding constants were determined. They are in the range of 0.3-8.0 M-1 and correlated strictly with the FC lipophilicity descriptor (logP).

The Stability Landscape of de novo TIM Barrels Explored by a Modular Design Approach

•

The TIM barrel is a versatile fold to understand structure-stability relationships.

•

A collection of de novo TIM barrels with improved hydrophobic cores was designed.

•

DeNovoTIMs are reversible in chemical and thermal unfolding, which is uncommon in TIM barrels.

•

Epistatic effects play a central role in DeNovoTIMs stabilization.

•

DeNovoTIMs navigate a previously uncharted region of the stability landscape.

The ability to design stable proteins with custom-made functions is a major goal in biochemistry with practical relevance for our environment and society. Understanding and manipulating protein stability provide crucial information on the molecular determinants that modulate structure and stability, and expand the applications of de novo proteins. Since the (β/⍺)₈-barrel or TIM-barrel fold is one of the most common functional scaffolds, in this work we designed a collection of stable de novo TIM barrels (DeNovoTIMs), using a computational fixed-backbone and modular approach based on improved hydrophobic packing of sTIM11, the first validated de novo TIM barrel, and subjected them to a thorough folding analysis. DeNovoTIMs navigate a region of the stability landscape previously uncharted by natural TIM barrels, with variations spanning 60 degrees in melting temperature and 22 kcal per mol in conformational stability throughout the designs. Significant non-additive or epistatic effects were observed when stabilizing mutations from different regions of the barrel were combined. The molecular basis of epistasis in DeNovoTIMs appears to be related to the extension of the hydrophobic cores. This study is an important step towards the fine-tuned modulation of protein stability by design.

Aggregation-phase diagrams of β2-microglobulin reveal temperature and salt effects on competitive formation of amyloids versus amorphous aggregates

Several serious diseases are associated with crystal-like amyloid fibrils or glass-like amorphous aggregates of denatured proteins. However, protein aggregation involving both types of aggregates has not yet been elucidated in much detail. Using a protein associated with dialysis-related amyloidosis, β₂-microglobulin (β2m), we previously demonstrated that amyloid fibrils and amorphous aggregates form competitively depending on salt (NaCl) concentration. To examine the generality of the underlying competitive mechanisms, we herein investigated the effects of heat on acid-denatured β2m at pH 2. Using thioflavin fluorescence, CD, and light scattering analysis along with atomic force microscopy imaging, we found that the temperature-dependent aggregation of β2m markedly depends on NaCl concentration. Stepwise transitions from monomers to amyloids and then back to monomers were observed at low NaCl concentrations. Amorphous aggregates formed rapidly at ambient temperatures at high NaCl concentrations, but the transition from amorphous aggregates to amyloids occurred only as the temperature increased. Combining the data from the temperature- and NaCl-dependent transitions, we constructed a unified phase diagram of conformational states, indicating a parabolic solubility curve with a minimum NaCl concentration at ambient temperatures. Although amyloid fibrils formed above this solubility boundary, amorphous aggregates dominated in regions distant from this boundary. Kinetic competition between supersaturation-limited slow amyloid fibrillation and supersaturation-unlimited fast amorphous aggregation deformed the phase diagram, with amyloid regions disappearing with fast titration rates. We conclude that phase diagrams combining thermodynamics and kinetics data provide a comprehensive view of β2m aggregation exhibiting severe hysteresis depending on the heat- or salt-titration rates.

Differential scanning calorimetry: An invaluable tool for a detailed thermodynamic characterization of macromolecules and their interactions

Differential Scanning Calorimetry (DSC) is a highly sensitive technique to study the thermotropic properties of many different biological macromolecules and extracts. Since its early development, DSC has been applied to the pharmaceutical field with excipient studies and DNA drugs. In recent times, more attention has been applied to lipid-based drug delivery systems and drug interactions with biomimetic membranes. Highly reproducible phase transitions have been used to determine values, such as, the type of binding interaction, purity, stability, and release from a drug delivery mechanism. This review focuses on the use of DSC for biochemical and pharmaceutical applications.

Changes in thermodynamic stability of von Willebrand factor differentially affect the force-dependent binding to platelet GPIbalpha

In circulation, plasma glycoprotein von Willebrand Factor plays an important role in hemostasis and in pathological thrombosis under hydrodynamic forces. Mutations in the A1 domain of von Willebrand factor cause the hereditary types 2B and 2M von Willebrand disease that either enhance (2B) or inhibit (2M) the interaction of von Willebrand factor with the platelet receptor glycoprotein Ibalpha. To understand how type 2B and 2M mutations cause clinically opposite phenotypes, we use a combination of protein unfolding thermodynamics and atomic force microscopy to assess the effects of two type 2B mutations (R1306Q and I1309V) and a type 2M mutation (G1324S) on the conformational stability of the A1 domain and the single bond dissociation kinetics of the A1-GPIbalpha interaction. At physiological temperature, the type 2B mutations destabilize the structure of the A1 domain and shift the A1-GPIbalpha catch to slip bonding to lower forces. Conversely, the type 2M mutation stabilizes the structure of the A1 domain and shifts the A1-GPIbalpha catch to slip bonding to higher forces. As a function of increasing A1 domain stability, the bond lifetime at low force decreases and the critical force required for maximal bond lifetime increases. Our results are able to distinguish the clinical phenotypes of these naturally occurring mutations from a thermodynamic and biophysical perspective that provides a quantitative description of the allosteric coupling of A1 conformational stability with the force dependent catch to slip bonding between A1 and GPIbalpha.

Amyloidogenic potential of alpha-chymotrypsin in different conformational states

Amyloid fibril formation is widely believed to be a generic property of polypeptide chains. In the present study, alpha-chymotrypsin, a well-known serine protease has been driven toward these structures by the use of two different conditions involving (I) high temperature, pH 2.5, and (II) low concentration of trifluoroethanol (TFE), pH 2.5. A variety of experimental methods, including fluorescence emission, dynamic quenching, steady-state fluorescence anisotropy, far-UV circular dichroism, nuclear magnetic resonance spectroscopy, and dynamic light scattering were employed to characterize the conformational states of alpha-chymotrypsin that precede formation of amyloid fibrils. The structure formed under Condition I was an unfolded monomer, whereas an alpha-helical rich oligomer was induced in Condition II. Both the amyloid aggregation-prone species manifested a higher solvent exposure of hydrophobic and aromatic residues compared with the native state. Upon incubation of the protein in these conditions for 48 h, amyloid-like fibrils were formed with diameters of about 10-12 nm. In contrast, at neutral pH and low concentration of TFE, a significant degree of amorphous aggregation was observed, suggesting that charge neutralization of acidic residues in the amyloid core region has a positive influence on amyloid fibril formation. In summary, results presented in this communication suggest that amyloid fibrils of alpha-chymotrypsin may be obtained from a variety of structurally distinct conformational ensembles highlighting the critical importance of protein evolution mechanisms related to prevention of protein misfolding.

Microcalorimetry of proteins and their complexes

Ultrasensitive microcalorimetric techniques for measuring the heat capacities of proteins in dilute solutions over a broad temperature range (DSC) and the heats of protein reactions at fixed temperatures (ITC) are described and the methods of working with these instruments are considered. Particular attention is paid to analyzing the thermal properties of individual proteins, their stability, the energetics of their folding, and their association with specific macromolecular partners. Use of these calorimetric methods is illustrated with examples of small compact globular proteins, small proteins having loose noncompact structure, multidomain proteins, and protein complexes, particularly with DNA.

Differential scanning calorimetry

Differential scanning calorimetry (DSC) has emerged as a powerful experimental technique for determining thermodynamic properties of biomacromolecules. The ability to monitor unfolding or phase transitions in proteins, polynucleotides, and lipid assemblies has not only provided data on thermodynamic stability for these important molecules, but also made it possible to examine the details of unfolding processes and to analyze the characteristics of intermediate states involved in the melting of biopolymers. The recent improvements in DSC instrumentation and software have generated new opportunities for the study of the effects of structure and changes in environment on the behavior of proteins, nucleic acids, and lipids. This review presents some of the details of application of DSC to the examination of the unfolding of biomolecules. After a brief introduction to DSC instrumentation used for the study of thermal transitions, the methods for obtaining basic thermodynamic information from the DSC curve are presented. Then, using DNA unfolding as an example, methods for the analysis of the melting transition are presented that allow deconvolution of the DSC curves to determine more subtle characteristics of the intermediate states involved in unfolding. Two types of transitions are presented for analysis, the first example being the unfolding of two large synthetic polynucleotides, which display high cooperativity in the melting process. The second example shows the application of DSC for the study of the unfolding of a simple hairpin oligonucleotide. Details of the data analysis are presented in a simple spreadsheet format.

Microcalorimetry of biological macromolecules

The capabilities of contemporary differential scanning and isothermal titration microcalorimetry for studying the thermodynamics of protein unfolding/refolding and their association with partners, particularly target DNA duplexes, are considered. It is shown that the predenaturational changes of proteins must not be ignored in studying the thermodynamics of formation of their native structure and their complexes with partners, particularly their cognate DNA duplexes.

Role of the N-terminal region for the conformational stability of esterase 2 from Alicyclobacillus acidocaldarius

In order to clarify the role played by the N-terminal region for the conformational stability of the thermophilic esterase 2 (EST2) from Alicyclobacillus acidocaldarius, two mutant forms have been investigated: a variant obtained by deleting the first 35 residues at the N-terminus (EST2-36del), and a variant obtained by mutating Lys102 to Gln (K102Q) to perturb the N-terminus by destroying the salt bridge E43-K102. The temperature- and denaturant-induced unfolding of EST2 and the two mutant forms have been studied by means of circular dichroism (CD), differential scanning calorimetry (DSC) and fluorescence measurements. In line with its thermophilic origin, the denaturation temperature of EST2 is high: T(d)=91 degrees C and 86 degrees C if detected by recording the CD signal at 222 nm and 290 nm, respectively. This difference suggests that the thermal denaturation process, even though reversible, is more complex than a two-state Nright arrow over left arrowD transition. The non-two-state behaviour is more pronounced in the case of the two mutant forms. The complex DSC profiles of EST2 and both mutant forms have been analysed by means of a deconvolution procedure. The thermodynamic parameters characterizing the two transitions obtained in the case of EST2 are: T(d,1)=81 degrees C, Delta(d)H(1)=440 kJ mol(-1), Delta(d)C(p,1)=7 kJ K(-1)mol(-1), T(d,2)=86 degrees C, Delta(d)H(2)=710 kJ mol(-1), and Delta(d)C(p,2)=9 kJ K(-1)mol(-1). The first transition occurs at lower temperatures in the two mutant forms, whereas the second transition is always centred at 86 degrees C. The results indicate that EST2 possesses two structural domains whose coupling is tight in the wild-type protein, but markedly weakens in the two mutant forms as a consequence of the perturbations in the N-terminal region.

Calorimetric and crystallographic analysis of the oligomeric structure of Escherichia coli GMP kinase

Guanosine monophosphate kinases (GMPKs), which catalyze the phosphorylation of GMP and dGMP to their diphosphate form, have been characterized as monomeric enzymes in eukaryotes and prokaryotes. Here, we report that GMPK from Escherichia coli (ecGMPK) assembles in solution and in the crystal as several different oligomers. Thermodynamic analysis of ecGMPK using differential scanning calorimetry shows that the enzyme is in equilibrium between a dimer and higher order oligomers, whose relative amounts depend on protein concentration, ionic strength, and the presence of ATP. Crystallographic structures of ecGMPK in the apo, GMP and GDP-bound forms were solved at 3.2A, 2.9A and 2.4A resolution, respectively. ecGMPK forms a hexamer with D3 symmetry in all crystal forms, in which the two nucleotide-binding domains are able to undergo closure comparable to that of monomeric GMPKs. The 2-fold and 3-fold interfaces involve a 20-residue C-terminal extension and a sequence signature, respectively, that are missing from monomeric eukaryotic GMPKs, explaining why ecGMPK forms oligomers. These signatures are found in GMPKs from proteobacteria, some of which are human pathogens. GMPKs from these bacteria are thus likely to form the same quaternary structures. The shift of the thermodynamic equilibrium towards the dimer at low ecGMPK concentration together with the observation that inter-subunit interactions partially occlude the ATP-binding site in the hexameric structure suggest that the dimer may be the active species at physiological enzyme concentration.

HEAT CAPACITY IN PROTEINS

Heat capacity (Cp) is one of several major thermodynamic quantities commonly measured in proteins. With more than half a dozen definitions, it is the hardest of these quantities to understand in physical terms, but the richest in insight. There are many ramifications of observed Cp changes: The sign distinguishes apolar from polar solvation. It imparts a temperature (T) dependence to entropy and enthalpy that may change their signs and which of them dominate. Protein unfolding usually has a positive deltaCp, producing a maximum in stability and sometimes cold denaturation. There are two heat capacity contributions, from hydration and protein-protein interactions; which dominates in folding and binding is an open question. Theoretical work to date has dealt mostly with the hydration term and can account, at least semiquantitatively, for the major Cp-related features: the positive and negative Cp of hydration for apolar and polar groups, respectively; the convergence of apolar group hydration entropy at T approximately 112 degrees C; the decrease in apolar hydration Cp with increasing T; and the T-maximum in protein stability and cold denaturation.

Deciphering the function of an ORF: Salmonella enterica DeoM protein is a new mutarotase specific for deoxyribose

We identified in Salmonella enterica serovar Typhi a cluster of four genes encoding a deoxyribokinase (DeoK), a putative permease (DeoP), a repressor (DeoQ), and an open reading frame encoding a 337 amino acid residues protein of unknown function. We show that the latter protein, called DeoM, is a hexamer whose synthesis is increased by a factor over 5 after induction with deoxyribose. The CD spectrum of the purified recombinant protein indicated a dominant contribution of betatype secondary structure and a small content of alpha-helix. Temperature and guanidinium hydrochloride induced denaturation of DeoM indicated that the hexamer dissociation and monomer unfolding are coupled processes. DeoM exhibits 12.5% and 15% sequence identity with galactose mutarotase from Lactococcus lactis and respectively Escherichia coli, which suggested that these three proteins share similar functions. Polarimetric experiments demonstrated that DeoM is a mutarotase with high specificity for deoxyribose. Site-directed mutagenesis of His183 in DeoM, corresponding to a catalytically active residue in GalM, yielded an almost inactive deoxyribose mutarotase. DeoM was crystallized and diffraction data collected for two crystal systems, confirmed its hexameric state. The possible role of the protein and of the entire gene cluster is discussed in connection with the energy metabolism of S. enterica under particular growth conditions.

A general criterion based on the implicit function theorem to model aggregation and adsorption in colloidal dispersions

This paper, which may interest not only colloid scientists and physical chemists but also applied mathematicians, completes some previous results on aqueous silicon nitride dispersions. Experimental data on adsorption from liquid solution were first obtained by a titration method and then used to derive the number of solid particles from an equilibrium constraint. To discuss the complex mechanisms affecting simultaneous solid particle aggregation and small ion adsorption at the solid/liquid interface, the Dini implicit function theorem (DT) has been applied to the equilibrium condition for a former suspension Gibbs free energy. It was able to relate the average particle number to the ion concentration adsorbed, but not to unequivocally specify their dependence on the liquid phase pH. We attempt here to model aggregation both through bulk and interfacial quantities. The generalized DT-based criterion has first been formulated in all generality, and then adopted according to a wider investigation. The results obtained confirm the original guess, i.e., to regard solid aggregation as dominated by interfacial mechanisms.

Temperature dependence of the thermodynamics of helix-coil transition

The temperature-induced helix to coil transition in a series of host peptides was monitored using circular dichroism spectroscopy (CD) and differential scanning calorimetry (DSC). Combination of these two techniques allowed direct determination of the enthalpy of helix-coil transition for the studied peptides. It was found that the enthalpy of the helix-coil transition differs for different peptides and this difference is related to the difference in the temperature for the midpoint of helix-coil transition. The enthalpy of the helix-coil transition decreases with the increase in temperature, thus providing the first experimental estimate for the heat capacity changes upon helix-coil transition, DeltaC(p). The values for DeltaC(p) of helix-coil transition are found to be negative, which is in contrast to the positive DeltaC(p) for protein unfolding. Analysis suggests that this negative DeltaC(p) of helix-coil transition is due to the exposure of the polar peptide backbone to solvent upon helix unfolding.

Thermodynamic characterization of variants of mesophilic cytochrome c and its thermophilic counterpart

Thermal stability was measured for variants of cytochrome c-551 (PA c-551) from a mesophile, Pseudomonas aeruginosa, and a thermophilic counterpart, Hydrogenobacter thermophilus cytochrome c-552 (HT c-552), by differential scanning calorimetry (DSC) at pH 3.6. The mutated residues in PA c-551, selected with reference to the corresponding residues in HT c-552, were located in three spatially separated regions: region I, Phe7 to Ala/Val13 to Met; region II, Glu34 to Tyr/Phe43 to Tyr; and region III, Val78 to Ile. The thermodynamic parameters determined indicated that the mutations in regions I and III caused enhanced stability through not only enthalpic but also entropic contributions, which reflected improved packing of the side chains. Meanwhile, the mutated region II made enthalpic contributions to the stability through electrostatic interactions. The obtained differences in the Gibbs free energy changes of unfolding [Delta(DeltaG)] showed that the three regions contributed to the overall stability in an additive manner. HT c-552 had the smallest heat capacity change (DeltaC(P)), resulting in higher DeltaG values over a wide temperature range (0-100 degrees C), compared to the PA c-551 variants; this contributed to the highest stability of HT c-552. Our DSC measurement results, in conjunction with mutagenesis and structural studies on the homologous mesophilic and thermophilic cytochromes c, provided an extended thermodynamic view of protein stabilization.

Understanding thermostability in cytochrome P450 by combinatorial mutagenesis

The cytochromes P450 are an important class of mono-oxygenases involved in xenobiotic metabolism and steroid biosynthesis in a diverse set of life forms. Discovery of CYP-119, a P450 from the archea Sulfolobus solfataricus has provided a means for understanding nature's method of stabilizing this important protein superfamily. To identify classes of stabilizing interactions used by CYP-119, we have generated a randomized library of point mutants and screened for mutants that are less thermostable than the wild type by monitoring the characteristic Soret band in the visible region of the cell lysis. The selected mutants were characterized by differential scanning calorimetry to compare the temperatures of the melting transitions of the various mutants. The identified mutations suggested that electrostatic interactions involving salt links and charge-charge interactions, as well as contributions from other interactions such as aromatic stacking, and side chain volume of hydrophobic residues contribute to enhanced thermostability in this cytochrome P450.

Selected mutations in a mesophilic cytochrome c confer the stability of a thermophilic counterpart

Mesophilic cytochrome c(551) of Pseudomonas aeruginosa (PA c(551)) became as stable as its thermophilic counterpart, Hydrogenobacter thermophilus cytochrome c(552) (HT c(552)), through only five amino acid substitutions. The five residues, distributed in three spatially separated regions, were selected and mutated with reference to the corresponding residues in HT c(552) through careful structure comparison. Thermodynamic analysis indicated that the stability of the quintuple mutant of PA c(551) could be partly attained through an enthalpic factor. The solution structure of the mutant showed that, as in HT c(552), there were tighter side chain packings in the mutated regions. Furthermore, the mutant had an increased total accessible surface area, resulting in great negative hydration free energy. Our results provide a novel example of protein stabilization in that limited amino acid substitutions can confer the overall stability of a natural highly thermophilic protein upon a mesophilic molecule.

Isothermal titration calorimetry and differential scanning calorimetry as complementary tools to investigate the energetics of biomolecular recognition

The principles of isothermal titration calorimetry (ITC) and differential scanning calorimetry (DSC) are reviewed together with the basic thermodynamic formalism on which the two techniques are based. Although ITC is particularly suitable to follow the energetics of an association reaction between biomolecules, the combination of ITC and DSC provides a more comprehensive description of the thermodynamics of an associating system. The reason is that the parameters DeltaG, DeltaH, DeltaS, and DeltaCp obtained from ITC are global properties of the system under study. They may be composed to varying degrees of contributions from the binding reaction proper, from conformational changes of the component molecules during association, and from changes in molecule/solvent interactions and in the state of protonation.

Ethanol-induced conformational transitions in holo-alpha-lactalbumin: spectral and calorimetric studies

Conformational transitions of holo-alpha-lactalbumin in a hydro-ethanolic cosolvent system was studied by spectrofluorescence, CD in near- and far-uv regions, and high-sensitivity differential scanning calorimetry. Experimental results allow us to propose that in isothermal conditions alpha-lactalbumin undergoes a number of conformational transitions with increasing ethanol concentration: N<=>I<=>D<=>H. The existence of I-state was deduced from spectrofluorometric and near-uv CD data. In this state the aromatic chromophores of the amino acid side chains are more accessible to the solvent displaying higher local mobility. The H-state was detected from far-uv CD spectra as a state corresponding to the content of alpha-helices higher than originally found in native protein. However, calorimetric measurements provide data revealing only the two-state mechanism of alpha-lactalbumin unfolding in both water and in aqueous ethanol solutions. This indicates that the energy levels of N- and I-states as well as of D- and H-states are similar. Thermodynamics of the unfolding of alpha-lactalbumin in hydroethanolic solutions was analyzed with the help of the linear model of solvent denaturation. Unfolding increments of enthalpy, entropy, and Gibbs energy of transfer of the protein from a reference aqueous solution to hydro-ethanolic solutions of different concentrations were determined from the calorimetric data. They are linear functions of molar ethanol fraction. The slope of the unfolding increment of Gibbs energy of transfer was calculated from data on transfer of amino acid residues taking into account the average solvent accessibility of amino acid residues in the native structure of small globular proteins, using the additive group contribution method.

The energetics of HMG box interactions with DNA. Thermodynamic description of the box from mouse Sox-5

The structural energetics of the HMG box from the DNA-binding protein mouse Sox-5 were examined calorimetrically. It was found that this box, notwithstanding its small size (molecular mass about 10 kDa), does not behave as a single cooperative unit and, on heating, the box reversibly unfolds in two separate stages. The first transition (tt approximately 34 degrees C) involves about 40% of the total enthalpy and the second (tt approximately 46 degrees C) the remainder. Both transitions proceed with significant heat capacity increment, showing that they are associated with the unfolding of two sub-domains having non-polar cores. According to heat capacity, ellipticity, fluorescence and NMR criteria, this HMG box is in a fully compact native state only below 5 degrees C. HMG boxes consist of two approximately orthogonal wings: the minor wing comprises helix 3 and its associated antiparallel N-terminal strand, whilst the major wing is composed of helices I and II. Analysis of the fluorescence and NMR spectra for this box obtained at different temperatures shows that the lower melting transition can be assigned to the minor wing and the upper transition to the major wing. Under physiological conditions (37 degrees C), the minor wing is considerably unfolded, whilst the major wing is essentially fully folded. DNA binding in vivo therefore involves refolding of the minor wing.

Solvent isotope effect on thermodynamics of hydration

Partial molar heat capacities of five linear alcohols (methanol, ethanol, n-propanol, n-butanol, n-pentanol) and five N-substituted amides (n-propionamide, N-methylformamide, N-methylacetamide, N-methylpropionamide, N-ethylacetamide) in aqueous D(2)O solution have been measured at 25 degrees C. The heat capacities of transfer of these compounds from H(2)O to D(2)O were calculated using previously reported (Makhatadze et al., Biophys. Chem. 64 (1997) 93) values of partial heat capacities of alcohols and amides in aqueous H(2)O solutions. It is shown that the sign and magnitude of the heat capacity change upon transfer from H(2)O to D(2)O depends on the relative amount of polar and non-polar solvent accessible surface areas of solute. Analysis shows that transfer of non-polar surface from H(2)O to D(2)O is accompanied by a positive heat capacity change. In contrast, transfer of polar surface from H(2)O to D(2)O occurs with negative heat capacity change. Estimates show that the solvent isotope effect on the heat capacity changes upon protein unfolding can be predicted using the changes of the polar and non-polar surface area changes upon protein unfolding and the transfer data of model compounds. Analysis of the thermodynamic functions of transfer of non-polar compounds from H(2)O to D(2)O shows puzzling behavior which contradicts current definitions of the hydrophobic effect.

Thermodynamic characterization of an intermediate state of human growth hormone

The thermal denaturation of recombinant human growth hormone (rhGH) was studied by differential scanning calorimetry and circular dichroism spectroscopy (CD). The thermal unfolding is reversible only below pH 3.5, and under these conditions a single two-state transition was observed between 0 and 100 degrees C. The magnitudes of the deltaH and deltaCp of this transition indicate that it corresponds to a partial unfolding of rhGH. This is also supported by CD data, which show that significant secondary structure remains after the unfolding. Above pH 3.5 the thermal denaturation is irreversible due to the aggregation of rhGH upon unfolding. This aggregation is prevented in aqueous solutions of alcohols such as n-propanol, 2-propanol, or 1,2-propanediol (propylene glycol), which suggests that the self-association of rhGH is caused by hydrophobic interactions. In addition, it was found that the native state of rhGH is stable in relatively high concentrations of propylene glycol (up to 45% v/v at pH 7-8 or 30% at pH 3) and that under these conditions the thermal unfolding is cooperative and corresponds to a transition from the native state to a partially folded state, as observed at acidic pH in the absence of alcohols. In higher concentrations of propylene glycol, the tertiary structure of rhGH is disrupted and the cooperativity of the unfolding decreases. Moreover, the CD and DSC data indicate that a partially folded intermediate with essentially native secondary structure and disordered tertiary structure becomes significantly populated in 70-80% propylene glycol.

Structure-based statistical thermodynamic analysis of T4 lysozyme mutants: structural mapping of cooperative interactions

The recent development of a structural parameterization of the energetics of protein folding has permitted the incorporation of the functions that describe the enthalpy, entropy and heat capacity changes, i.e. the individual components of the Gibbs energy, into a statistical thermodynamic formalism that describes the distribution of conformational states under equilibrium conditions. The goal of this approach is to construct with the computer a large ensemble of conformational states, and then to derive the most probable population distribution, i.e. the distribution of states that best accounts for a wide array of experimental observables. This analysis has been applied to four different mutants of T4 lysozyme (S44A, S44G, V131A, V131G). It is shown that the structural parameterization predicts well the stability of the protein and the effects of the mutations. The entire set of folding constants per residue has been calculated for the four mutants. In all cases, the effect of the mutations propagates beyond the mutation site itself through sequence and three-dimensional space. This phenomenon occurs despite the fact that the mutations are at solvent-exposed locations and do not directly affect other interactions in the protein. These results suggest that single amino acid mutations at solvent-exposed locations, or other locations that cause a minimal perturbation, can be used to identify the extent of cooperative interactions. The magnitude and extent of these effects and the accuracy of the algorithm can be tested by means of NMR-detected hydrogen exchange.

Conformational stability of pGEX-expressed Schistosoma japonicum glutathione S-transferase: a detoxification enzyme and fusion-protein affinity tag

A glutathione S-transferase (Sj26GST) from Schistosoma japonicum, which functions in the parasite's Phase II detoxification pathway, is expressed by the Pharmacia pGEX-2T plasmid and is used widely as a fusion-protein affinity tag. It contains all 217 residues of Sj26GST and an additional 9-residue peptide linker with a thrombin cleavage site at its C-terminus. Size-exclusion HPLC (SEC-HPLC) and SDS-PAGE studies indicate that purification of the homodimeric protein under nonreducing conditions results in the reversible formation of significant amounts of 160-kDa and larger aggregates without a loss in catalytic activity. The basis for oxidative aggregation can be ascribed to the high degree of exposure of the four cysteine residues per subunit. The conformational stability of the dimeric protein was studied by urea- and temperature-induced unfolding techniques. Fluorescence-spectroscopy, SEC-HPLC, urea- and temperature-gradient gel electrophoresis, differential scanning microcalorimetry, and enzyme activity were employed to monitor structural and functional changes. The unfolding data indicate the absence of thermodynamically stable intermediates and that the unfolding/refolding transition is a two-state process involving folded native dimer and unfolded monomer. The stability of the protein was found to be dependent on its concentration, with a delta G degree (H2O) = 26.0 +/- 1.7 kcal/mol. The strong relationship observed between the m-value and the size of the protein indicates that the amount of protein surface area exposed to solvent upon unfolding is the major structural determinant for the dependence of the protein's free energy of unfolding on urea concentration. Thermograms obtained by differential scanning microcalorimetry also fitted a two-state unfolding transition model with values of delta Cp = 7,440 J/mol per K, delta H = 950.4 kJ/mol, and delta S = 1,484 J/mol.

Temperature-induced denaturation of ribonuclease S: a thermodynamic study

In this paper the thermal denaturation of ribonuclease S, the product of mild digestion of ribonuclease A by subtilisin, is deeply investigated by means of DSC and CD measurements. It results that at whatever pH in the range 4-7.5 the process if fully reversible but not well represented by the simple two-state N<-->D transition. Actually, a two-state model that considers both unfolding and dissociation, NL<-->D + L*, well accounts for the main features of the process: the tail present in the low-temperature side of DSC peaks and the marked dependence of denaturation temperature on protein concentration. This mechanism is strictly linked to the exact stoichiometry of RNase S. An excess of the protein component of RNase S, the so-called S-protein, shifts the system toward a more complex behavior, that deserves a separate treatment in the accompanying paper [Graziano, G., Catanzano, F., Giancola, C., & Barone, G. (1996) Biochemistry 35, 13386-13392]. The thermodynamic analysis leads to the conclusion that the difference in thermal stability between RNase S and RNase A is due to entropic effects, i.e., a greater conformational flexibility of both backbone and side chains in RNase S. The process becomes irreversible at pH 8.0-8.5, probably due to side-reactions occurring at high temperature. Finally, the influence of phosphate ion on the stability of RNase A and RNase S at pH 7.0 is studied and explained in terms of its binding on the active site of ribonuclease. The analysis enables us to obtain an estimate of the apparent association constant and binding enthalpy also.

Thermodynamic effects of mutations on the denaturation of T4 lysozyme

We investigated the folding of substantially destabilized mutant forms of T4 lysozyme using differential scanning calorimetry and circular dichroism measurements. Three mutations in an alpha-helix in the protein's N-terminal region, the alanine insertion mutations S44[A] and K48[A], and the substitution A42K had previously been observed to result in unexpectedly low apparent enthalpy changes of melting, compared to a pseudo-wild-type reference protein. The pseudo-wild-type reference protein thermally unfolds in an essentially two-state manner. However, we found that the unfolding of the three mutant proteins has reduced cooperativity, which partially explains their lower apparent enthalpy changes. A three-state unfolding model including a discrete intermediate is necessary to describe the melting of the mutant proteins. The reduction in cooperativity must be considered for accurate calculation of the energy changes of folding. Unfolding in two stages reflects the underlying two-subdomain structure of the lysozyme protein family.

Two partially unfolded states of Torpedo californica acetylcholinesterase

Chemical modification with sulfhydryl reagents of the single, nonconserved cysteine residue Cys231 in each subunit of a disulfide-linked dimer of Torpedo californica acetylcholinesterase produces a partially unfolded inactive state. Another partially unfolded state can be obtained by exposure of the enzyme to 1-2 M guanidine hydrochloride. Both these states display several important features of a molten globule, but differ in their spectroscopic (CD, intrinsic fluorescence) and hydrodynamic (Stokes radii) characteristics. With reversal of chemical modification of the former state or removal of denaturant from the latter, both states retain their physiochemical characteristics. Thus, acetylcholinesterase can exist in two molten globule states, both of which are long-lived under physiologic conditions without aggregating, and without either intraconverting or reverting to the native state. Both states undergo spontaneous intramolecular thioldisulfide exchange, implying that they are flexible. As revealed by differential scanning calorimetry, the state produced by chemical modification lacks any heat capacity peak, presumably due to aggregation during scanning, whereas the state produced by guanidine hydrochloride unfolds as a single cooperative unit, thermal transition being completely reversible. Sucrose gradient centrifugation reveals that reduction of the interchain disulfide of the native acetylcholinesterase dimer converts it to monomers, whereas, after such reduction, the two subunits remain completely associated in the partially unfolded state generated by guanidine hydrochloride, and partially associated in that produced by chemical modification. It is suggested that a novel hydrophobic core, generated across the subunit interfaces, is responsible for this noncovalent association. Transition from the unfolded state generated by chemical modification to that produced by guanidine hydrochloride is observed only in the presence of the denaturant, yielding, on extrapolation to zero guanidine hydrochloride, a high free energy barrier (ca. 23.8 kcal/mol) separating these two flexible, partially unfolded states.

The N-terminal domain of Escherichia coli enzyme I of the phosphoenolpyruvate/glycose phosphotransferase system: molecular cloning and characterization

The bacterial phosphoenolpyruvate/glycose phosphotransferase system (PTS) comprises a group of proteins that catalyze the transfer of the phosphoryl group from phosphoenolpyruvate (PEP) to sugars concomitant with their translocation. The first two steps of the phosphotransfer sequence are PEP <--> Enzyme I (EI) <--> HPr (the histidine-containing phosphocarrier protein). We have proposed that many functions of the PTS are regulated by EI, which undergoes a monomer/dimer transition. EI monomer (63.5 kDa) comprises two major domains: a flexible C-terminal domain (EI-C) and a protease-resistant, structurally stable N-terminal domain (EI-N) containing the active site His. Trypsin treatment of Salmonella typhimurium EI yielded EI-N, designated EI-N(t). Homogeneous recombinant Escherichia coli EI-N [i.e., EI-N(r)], has now been prepared in quantity, shows the expected thermodynamic unfolding properties and, similarly to EI-N(t), is phosphorylated by phospho-HPr, but not by PEP. In addition, binding of EI-N(r) to HPr was studied by isothermal titration calorimetry: K/a = 1.4 x 10(5) M(-1) and delta H = +8.8 kcal x mol(-1). Both values are comparable to those for HPr binding to intact EI. Fluorescence anisotropy [dansyl-EI-N(r)] and gel filtration of EI-N(r) show that it does not dimerize. These results emphasize the role of EI-C in dimerization and the regulation of intact EI.

Thermodynamic studies of the core histones: pH and ionic strength effects on the stability of the (H3-H4)/(H3-H4)2 system

The self-associative behavior and the thermal stability of the H3/H4 histone complex was studied in low-ionic strength conditions by several physicochemical techniques, including differential scanning calorimetry and circular dichroism spectroscopy. At neutrality, the major molecular species present in solution is the (H3-H4)2 tetramer. Its thermodynamic properties cannot be studied directly though, since its thermal denaturation is completely irreversible even at the lowest salt concentrations. However, a complete thermodynamic analysis can be performed at low ionic strength and pH 4.5, where the (H3-H4)2 tetramer is quantitatively dissociated into two H3-H4 dimers and where almost complete reversibility of the thermal transitions is attained. The unfolding transition temperature of the 26.5 kDa H3-H4 dimer increases as a function of both the ionic strength of the solvent and the total protein concentration. The thermal denaturation of the H3-H4 dimer is characterized by the presence of a single calorimetric peak, centered at 58 degrees C, with a corresponding enthalpy change of 25 kcal/mol of a 13 kDa monomer unit and a change in heat capacity upon unfolding of about 0.6 kcal/(K mol of 13 kDa monomer unit). The complex between histones H3 and H4 (tetramer or dimer) is stable between pH 9.5 and 3.0. At pH 1.5, the system is almost completely unfolded at all temperatures. At low ionic strengths and pH values between 5.0 and 2.5, the H3-H4 dimer behaves as a highly cooperative system, melting as a single unit; i.e. individual H3 and H4 folded monomers are not detectable during the treatment. The two-state mechanism accounting for the unfolding of the H3-H4 dimer at pH 4.5 is the same as that described for the H2A-H2B dimer at neutrality. Just like for the H2A and H2B histones, the H3 and H4 polypeptides are properly folded only when assembled as H3-H4 dimers or in higher-order histone assemblies. Therefore, coupling along the interfaces of the two chains within the heterodimer is the major factor contributing to the stabilization of the secondary and tertiary structures of the chains as well as of the histone dimers.

Calorimetry of proteins and nucleic acids

The availability of sensitive calorimetric instrumentation has led to a considerable increase in thermodynamic studies of proteins, nucleic acids, and their interactions. This article reviews some of the recent contributions of calorimetry to characterizing the thermodynamic origins of protein and nucleic acid stability and conformational preferences, as well as the interactions of proteins with each other, with small molecules, and with nucleic acids.