A thermodynamic comparison of mesophilic and thermophilic ribonucleases H

The mechanisms by which thermophilic proteins attain their increased thermostability remain unclear, as usually the sequence and structure of these proteins are very similar to those of their mesophilic homologues. To gain insight into the basis of thermostability, we have determined protein stability curves describing the temperature dependence of the free energy of unfolding for two ribonucleases H, one from the mesophile Escherichia coli and one from the thermophile Thermus thermophilus. The circular dichroism signal was monitored as a function of temperature and guanidinium chloride concentration, and the resulting free energies of unfolding were fit to the Gibbs-Helmholtz equation to obtain a set of thermodynamic parameters for these proteins. Although the maximal stabilities for these proteins occur at similar temperatures, the heat capacity of unfolding for T. thermophilus RNase H is lower, resulting in a smaller temperature dependence of the free energy of unfolding and therefore a higher thermal melting temperature. In addition, the stabilities of these proteins are similar at the optimal growth temperatures for their respective organisms, suggesting that a balance of thermodynamic stability and flexibility is important for function.     

Conformational stability of HPr: the histidine-containing phosphocarrier protein from Bacillus subtilis.

The conformational stability of the histidine-containing phosphocarrier protein (HPr) from Bacillus subtilis has been determined using a combination of thermal unfolding and solvent denaturation experiments. The urea-induced denaturation of HPr was monitored spectroscopically at fixed temperatures and thermal unfolding was performed in the presence of fixed concentrations of urea. These data were analyzed in several different ways to afford a measure of the cardinal parameters (delta Hg, Tg, delta Sg, and delta Cp) that describe the thermodynamics of folding for HPr. The method of Pace and Laurents (Pace CN, Laurents DV, 1989, Biochemistry 28:2520-2525) was used to estimate delta Cp as was a global analysis of the thermal- and urea-induced unfolding data. Each method used to analyze the data gives a similar value for delta Cp (1,170 +/- 50 cal mol-1K-1). Despite the high melting temperature for HPr (Tg = 73.5 degrees C), the maximum stability of the protein, which occurs at 26 degrees C, is quite modest (delta Gs = 4.2 kcal mol-1). In the presence of moderate concentrations of urea, HPr exhibits cold denaturation, and thus a complete stability curve for HPr, including a measure of delta Cp, can be achieved using the method of Chen and Schellman (Chen B, Schellman JA, 1989, Biochemistry 28:685-691). A comparison of the different methods for the analysis of solvent denaturation curves is provided and the effects of urea on the thermal stability of this small globular protein are discussed. The methods presented will be of general utility in the characterization of the stability curve for many small proteins.     

The sarcosine effect on protein stability: A case of nonadditivity?

We have used differential scanning calorimetry to determine the effect of low concentrations (C = 0-2 M) of the osmolyte sarcosine on the Gibbs energy changes (deltaG) for the unfolding of hen-egg-white lysozyme, ribonuclease A, and ubiquitin, under the same buffer and pH conditions. We have also computed this effect on the basis of the additivity assumption and using published values of the transfer Gibbs energies for the amino acid side chains and the peptide backbone unit. The values thus predicted for the slope delta deltaG/deltaC agree with the experimental ones, but only if the unfolded state is assumed to be compact (that is, if the accessibility to solvent of the unfolded state is modeled using segments excised from native structures). The additivity-based calculations predict similar delta deltaG/deltaC values for the three proteins studied. We point out that, to the extent that this approximate constancy of delta deltaG/deltaC holds, osmolyte-induced increases in denaturation temperature will be larger for proteins with low unfolding enthalpy (small proteins that bury a large proportion of apolar surface). The experimental results reported here are consistent with this hypothesis.     

The HPr proteins from the thermophile Bacillus stearothermophilus can form domain-swapped dimers

The study of proteins from extremophilic organisms continues to generate interest in the field of protein folding because paradigms explaining the enhanced stability of these proteins still elude us and such studies have the potential to further our knowledge of the forces stabilizing proteins. We have undertaken such a study with our model protein HPr from a mesophile, Bacillus subtilis, and a thermophile, Bacillus stearothermophilus. We report here the high-resolution structures of the wild-type HPr protein from the thermophile and a variant, F29W. The variant proved to crystallize in two forms: a monomeric form with a structure very similar to the wild-type protein as well as a domain-swapped dimer. Interestingly, the structure of the domain-swapped dimer for HPr is very different from that observed for a homologous protein, Crh, from B.subtilis. The existence of a domain-swapped dimer has implications for amyloid formation and is consistent with recent results showing that the HPr proteins can form amyloid fibrils. We also characterized the conformational stability of the thermophilic HPr proteins using thermal and solvent denaturation methods and have used the high-resolution structures in an attempt to explain the differences in stability between the different HPr proteins. Finally, we present a detailed analysis of the solution properties of the HPr proteins using a variety of biochemical and biophysical methods.     

Interaction of epitope-related and -unrelated peptides with anti-p24 (HIV-1) monoclonal antibody CB4-1 and its Fab fragment

The binding of four epitope-related peptides and three library-derived, epitope-unrelated peptides of different lengths (10-14 amino acids) and sequence by anti-p24 (HIV-1) monoclonal antibody CB4-1 and its Fab fragment was studied by isothermal titration calorimetry. The binding constants K(A) at 25 degrees C vary between 5.1 x 10(7) M (-1) for the strongest and 1.4 x 10(5) M (-1) for the weakest binder. For each of the peptides complex formation is enthalpically driven and connected with unfavorable entropic contributions; however, the ratio of enthalpy and entropy contributions to deltaG(0) differs markedly for the individual peptides. A plot of -deltaH(0) vs -TdeltaS(0) shows a linear correlation of the data for a wide variety of experimental conditions as expected for a process with deltaC(p) much larger than deltaS(0). The dissimilarity of deltaC(p) and deltaS(0) also explains why deltaH(0) and TdeltaS(0) show similar temperature dependences resulting in relatively small changes of deltaG(0) with temperature. The heat capacity changes deltaC(p) upon antibody-peptide complex formation determined for three selected peptides vary only in a small range, indicating basic thermodynamic similarity despite different key residues interacting in the complexes. Furthermore, the comparison of van't Hoff and calorimetric enthalpies point to a non-two-state binding mechanism. Protonation effects were excluded by measurements in buffers of different ionization enthalpies. Differences in the solution conformation of the peptides as demonstrated by circular dichroic measurements do not explain different binding affinities of the peptides; specifically a high helix content in solution is not essential for high binding affinity despite the helical epitope conformation in the crystal structure of p24.     

Thermal and urea-induced unfolding of the marginally stable lac repressor DNA-binding domain: a model system for analysis of solute effects on protein processes

Thermodynamic and structural evidence indicates that the DNA binding domains of lac repressor (lacI) exhibit significant conformational adaptability in operator binding, and that the marginally stable helix-turn-helix (HTH) recognition element is greatly stabilized by operator binding. Here we use circular dichroism at 222 nm to quantify the thermodynamics of the urea- and thermally induced unfolding of the marginally stable lacI HTH. Van't Hoff analysis of the two-state unfolding data, highly accurate because of the large transition breadth and experimental access to the temperature of maximum stability (T(S); 6-10 degrees C), yields standard-state thermodynamic functions (deltaG(o)(obs), deltaH(o)(obs), deltaS(o)(obs), deltaC(o)(P,obs)) over the temperature range 4-40 degrees C and urea concentration range 0 相似文献   

Stabilization of proteins by enhancement of inter-residue hydrophobic contacts: lessons of T4 lysozyme and barnase

Although the hydrophobic interactions are considered as the main contributors to the protein stability, not much examples of protein stabilization by rational increasing of this type of interactions still can be found in literature. This is partly due to the lack of proper theoretical "measure" of hydrophobic interactions and their changes upon mutations. In the present paper the molecular hydrophobicity potential approach is used to assess how the changes in type and the strength of inter-residue contacts upon single amino acid mutations are correlated with the changes in thermodynamic stability of T4 lysozyme and barnase mutants, and which factors affect these correlations. Mutations changing unfavorable hydrophilic-to-hydrophobic contacts into favorable hydrophobic were found to enhance the thermodynamic stability in more than 81 % of cases, if these mutations do not create steric bumps and do not involve proline residues and hydrogen-bonded side-chains. Mutations increasing hydrophobic contributions (according to molecular hydrophobicity potential formalism) lead to increase of thermodynamic stability in more than 94% of cases for certain type of mutations (i.e., mutations not involving charged residues, Pro and residues with side-chain hydrogen bonds, when these mutations do not introduce steric bumps and do not involve strongly exposed residues and residues situated at helix N- and C-cap positions). For this type of mutations the correlation was found between the change in hydrophobic contributions of mutated residues deltaCphob and thermodynamic parameters deltaTm (change in melting temperature) and deltadeltaG (change in free energy of unfolding). Although the correlation coefficients were larger if the experimental structures of mutants were used for the calculations (correlation coefficients r(exp) deltaC,deltaT = .85 and r(exp) deltaC,deltadeltaG = 0.87) than if the modeled structures were used instead (r(mod) deltaC,deltaT = 0.74 and r(mod)deltaC,deltadeltaG = 0.76), the modelled structures of mutants in the vast majority of cases can be used for qualitative predictition of the protein stabilization. Basing on the analysis of mutations increasing hydrophobic contributions in T4 lysozyme the substitution matrix was derived, which can be used to decide which new residue should be put instead the old one to increase the stability of protein. The estimation shows that the number of potential mutation sites for enhancement of hydrophobic interactions in T4 lysozyme is quite large, and only approximately 10 per cent of them were studied thus far. Basing on the current analysis of T4 lysozyme and barnase mutations the algorithm for increasing of protein stability via increasing of hydrophobic interactions for the proteins with known spatial structure is proposed.     

Thermodynamic molecular switch in macromolecular interactions

It is known that most living systems can live and operate optimally only at a sharply defined temperature, or over a limited temperature range, at best, which implies that many basic biochemical interactions exhibit a well-defined Gibbs free energy minimum as a function of temperature. The Gibbs free energy change, deltaG(o) (T), for biological systems shows a complicated behavior, in which deltaG(o)(T) changes from positive to negative, then reaches a negative value of maximum magnitude (favorable), and finally becomes positive as temperature increases. The critical factor in this complicated thermodynamic behavior is a temperature-dependent heat capacity change (deltaCp(o)(T) of reaction, which 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, Keq, and/or spontaneity. The subsequent, mathematically predictable changes in deltaH(o)(T), deltaS(o)(T), deltaW(o)(T), and deltaG(o)(T) give rise to the classically observed behavior patterns in biological reactivity, as demonstrated in three interacting protein systems: the acid dimerization reaction of alpha-chymotrypsin at low pH, interaction of chromogranin A with the intraluminal loop peptide of the inositol 1,4,5-triphosphate receptor at pH 5.5, and the binding of L-arabinose and D-galactose to the L-arabinose binding protein of Escherichia coli. In cases of protein unfolding of four mutants of phage T4 lysozyme, no thermodynamic molecular switch is observed.     

Thermodynamic basis for the increased thermostability of CheY from the hyperthermophile Thermotoga maritima.

The CheY protein isolated from the hyperthermophile Thermotoga maritima is much more resistant to thermally induced unfolding than is its counterpart from the mesophile Bacillus subtilis. To determine the basis of this increased thermostability, the temperature dependence of the free energy of unfolding was determined for these CheY homologues using denaturant-induced unfolding experiments. This allowed comparison of T. maritima CheY with B. subtilis CheY and determination of the thermodynamic qualities responsible for the enhanced thermostability of T. maritima CheY. The stability of the thermophilic CheY protein is a direct result of the increased enthalpy contribution at the temperature of zero entropy, T(s), and the decreased heat capacity change upon unfolding, resulting in a decreased dependence of the free energy of unfolding on temperature. It was found that neither purely entropic nor purely enthalpic contributions alone (as reflected by T(s)) were sufficient to account for the increase in stability.     

The denaturation of alpha,beta and psi bovine trypsin at pH 3.0: evidence of intermediates

The conformational stability and the folding process of alpha, beta and Psi bovine trypsin at pH 3.0 followed by circular dichroism (CD) and size exclusion in HPLC have been analyzed as a function of urea concentration. The thermodynamic stability for a and b are deltaG = 15.91 +/- 0.28 kcal/mol, deltaG = 15.54 +/- 2.39 kcal/mol. respectively, and y trypsin is deltaG = 16.10 +/- 2.51 kcal/mol. The transition curves for alpha, beta and Psi forms suggest a molten globule state.     

Unfolding energetics and conformational stability of DLC8 monomer

To understand the rules governing the protein folding process it is essential to study the stability and unfolding of small monomeric proteins. Here, I present the pH dependent thermal unfolding energetics and conformational stability analysis of monomeric Dynein light chain protein (DLC8) in the pH range 3.5-2.0. DLC8 is the smallest and the most conserved light chain among the light chains of the dynein motor assembly. Thermal unfolding of DLC8 monomer is much complex with the presence of transient intermediates, which is in contrast to the notion that small proteins unfold via simple two-state process. The unfolding seems to be more cooperative at lower pH and the temperature of highest conformational stability (T(s)) is found to be maximum (295.7 K) at pH 2.76. Stability curves have been simulated to understand the thermodynamic parameters that govern the shapes of the experimentally obtained curves. Further, an effort has been made to correlate the observed differences in the denaturation energetics with the protein sequence in order to throw light on the structure-folding paradigm of the DLC8 monomer.     

Mechanism of thermostabilization in a designed cold shock protein with optimized surface electrostatic interactions

Using computational and sequence analysis of bacterial cold shock proteins, we designed a protein (CspB-TB) that has the core residues of mesophilic protein from Bacillus subtilis(CspB-Bs) and altered distribution of surface charged residues. This designed protein was characterized by circular dichroism spectroscopy, and found to have secondary and tertiary structure similar to that of CspB-Bs. The activity of the CspB-TB protein as measured by the affinity to a single-stranded DNA (ssDNA) template at 25 degrees C is somewhat higher than that of CspB-Bs. Furthermore, the decrease in the apparent binding constant to ssDNA upon increase in temperature is much more pronounced for CspB-Bs than for CspB-TB. Temperature-induced unfolding (as monitored by differential scanning calorimetry and circular dichroism spectroscopy) and urea-induced unfolding experiments were used to compare the stabilities of CspB-Bs and CspB-TB. It was found that CspB-TB is approximately 20 degrees C more thermostable than CspB-Bs. The thermostabilization of CspB-TB relative to CspB-Bs is achieved by decrease in the enthalpy and entropy of unfolding without affecting their temperature dependencies, i.e. these proteins have similar heat capacity changes upon unfolding. These changes in the thermodynamic parameters result in the global stability function, i.e. Gibbs energy, deltaG(T), that is shifted to higher temperatures with only small changes in the maximum stability. Such a mechanism of thermostabilization, although predicted from the basic thermodynamic considerations, has never been identified experimentally.     

Temperature, stability, and the hydrophobic interaction

Changes in free energy are normally used to track the effect of temperature on the stability of proteins and hydrophobic interactions. Use of this procedure on the aqueous solubility of hydrocarbons, a standard representation of the hydrophobic effect, leads to the conclusion that the hydrophobic effect increases in strength as the temperature is raised to approximately 140 degrees C. Acceptance of this interpretation leads to a number of far-reaching conclusions that are at variance with the original conception of the hydrophobic effect and add considerably to the complexity of interpretation. There are two legitimate thermodynamic functions that can be used to look at stability as a function of temperature: the standard Gibbs free energy change, deltaG degrees, and deltaG degrees/T. The latter is proportional to the log of the equilibrium constant and is sometimes called the Massieu-Planck function. Arguments are presented for using deltaG degrees/T rather than deltaG degrees for variations in stability with temperature. This makes a considerable difference in the interpretation of the hydrophobic interaction, but makes little change in the stability profile of proteins. Protein unfolding and the aqueous solubility of benzene are given as examples. The contrast between protein unfolding and the hydration of nonpolar molecules provides a rough estimate of the contribution of other factors that stabilize and destabilize protein structure.     

Thermodynamic differences among homologous thermophilic and mesophilic proteins.

Here, we analyze the thermodynamic parameters and their correlations in families containing homologous thermophilic and mesophilic proteins which show reversible two-state folding <--> unfolding transitions between the native and the denatured states. For the proteins in these families, the melting temperatures correlate with the maximal protein stability change (between the native and the denatured states) as well as with the enthalpic and entropic changes at the melting temperature. In contrast, the heat capacity change is uncorrelated with the melting temperature. These and additional results illustrate that higher melting temperatures are largely obtained via an upshift and broadening of the protein stability curves. Both thermophilic and mesophilic proteins are maximally stable around room temperature. However, the maximal stabilities of thermophilic proteins are considerably greater than those of their mesophilic homologues. At the living temperatures of their respective source organisms, homologous thermophilic and mesophilic proteins have similar stabilities. The protein stability at the living temperature of the source organism does not correlate with the living temperature of the protein. We tie thermodynamic observations to microscopics via the hydrophobic effect and a two-state model of the water structure. We conclude that, to achieve higher stability and greater resistance to high and low temperatures, specific interactions, particularly electrostatic, should be engineered into the protein. The effect of these specific interactions is largely reflected in an increased enthalpy change at the melting temperature.     

Temperature range of thermodynamic stability for the native state of reversible two-state proteins

The difference between the heat (T(G)) and the cold (T(G)') denaturation temperatures defines the temperature range (T(Range)) over which the native state of a reversible two-state protein is thermodynamically stable. We have performed a correlation analysis for thermodynamic parameters in a selected data set of structurally nonhomologous single-domain reversible two-state proteins. We find that the temperature range is negatively correlated with the protein size and with the heat capacity change (DeltaC(p)) but is positively correlated with the maximal protein stability [DeltaG(T(S))]. The correlation between the temperature range and maximal protein stability becomes highly significant upon normalization of the maximal protein stability with protein size. The melting temperature (T(G)) also shows a negative correlation with protein size. Consistently, T(G) and T(G)' show opposite correlations with DeltaC(p), indicating a dependence of the T(Range) on the curvature of the protein stability curve. Substitution of proteins in our data set with their homologues and arbitrary addition or removal of a protein in the data set do not affect the outcome of our analysis. Simulations of the thermodynamic data further indicate that T(Range) is more sensitive to variations in curvature than to the slope of the protein stability curve. The hydrophobic effect in single domains is the principal reason for these observations. Our results imply that larger proteins may be stable over narrower temperature ranges and that smaller proteins may have higher melting temperatures, suggesting why protein structures often differentiate into multiple substructures with different hydrophobic cores. Our results have interesting implications for protein thermostability.     

The histidine-phosphocarrier protein of Streptomyces coelicolor folds by a partially folded species at low pH.

The folding of a 93-residue protein, the histidine-phosphocarrier protein of Streptomyces coelicolor, HPr, has been studied using several biophysical techniques, namely fluorescence, 8-anilinonaphthalene-1-sulfate binding, circular dichroism, Fourier transform infrared spectroscopy, gel filtration chromatography and differential scanning calorimetry. The chemical-denaturation behaviour of HPr, followed by fluorescence, CD and gel filtration, at pH 7.5 and 25 degrees C, is described as a two-state process, which does not involve the accumulation of thermodynamically stable intermediates. Its conformational stability under those conditions is deltaG = 4.0 +/- 0.2 kcal x mol-1 (1 kcal = 4.18 kJ), which makes the HPr from S. coelicolor the most unstable member of the HPr family described so far. The stability of the protein does not change significantly from pH 7-9, as concluded from the differential scanning calorimetry and thermal CD experiments. Conformational studies at low pH (pH 2.5-4) suggest that, in the absence of cosmotropic agents, HPr does not unfold completely; rather, it accumulates partially folded species. The transition from those species to other states with native-like secondary and tertiary structure, occurs with a pKa = 3.3 +/- 0.3, as measured by the averaged measurements obtained by CD and fluorescence. However, this transition does not agree either with: (a) that measured by burial of hydrophobic patches (8-anilinonaphthalene-1-sulfate binding experiments); or (b) that measured by acquisition of native-like compactness (gel-filtration studies). It seems that acquisition of native-like features occurs in a wide pH range and it cannot be ascribed to a unique side-chain titration. These series of intermediates have not been reported previously in any member of the HPr family.     

Kinetic analysis of enhanced thermal stability of an alkaline protease with engineered twin disulfide bridges and calcium-dependent stability

The thermal stability of a cysteine-free alkaline protease (Alp) secreted by the eukaryote Aspergillus oryzae was improved both by the introduction of engineered twin disulfide bridges (Cys-69/Cys-101 and Cys-169/Cys-200), newly constructed as part of this study, and by the addition of calcium ions. We performed an extensive kinetic analysis of the increased thermal stability of the mutants as well as the role of calcium dependence. The thermodynamic activation parameters for irreversible thermal inactivation, the activation free energy (deltaG), the activation enthalpy (deltaH), and the activation entropy (deltaS) were determined from absolute reaction rate theory. The values of deltaH and deltaS were significantly and concomitantly increased as a result of introducing the twin disulfide bridges, for which the increase in the value of deltaH outweighed that of deltaS, resulting in significant increases in the value of deltaG. The enhancement of the thermal stability obtained by introducing the twin disulfide bridges is an example of the so-called low-temperature stabilization of enzymes. The stabilizing effect of calcium ions on wild-type Alp is similar to the results we obtained by introducing the engineered twin disulfide bridges.     

Thermostable variants of the recombinant xylanase A from Bacillus subtilis produced by directed evolution show reduced heat capacity changes

Directed evolution techniques have been used to improve the thermal stability of the xylanase A from Bacillus subtilis (XylA). Two generations of random mutant libraries generated by error prone PCR coupled with a single generation of DNA shuffling produced a series of mutant proteins with increasing thermostability. The most Thermostable XylA variant from the third generation contained four mutations Q7H, G13R, S22P, and S179C that showed an increase in melting temperature of 20 degrees C. The thermodynamic properties of a representative subset of nine XylA variants showing a range of thermostabilities were measured by thermal denaturation as monitored by the change in the far ultraviolet circular dichroism signal. Analysis of the data from these thermostable variants demonstrated a correlation between the decrease in the heat capacity change (deltaC(p)) with an increase in the midpoint of the transition temperature (T(m)) on transition from the native to the unfolded state. This result could not be interpreted within the context of the changes in accessible surface area of the protein on transition from the native to unfolded states. Since all the mutations are located at the surface of the protein, these results suggest that an explanation of the decrease in deltaC(p) should include effects arising from the protein/solvent interface.