Evolution and thermodynamics of the slow unfolding of hyperstable monomeric proteins
© Okada et al; licensee BioMed Central Ltd. 2010
Received: 29 March 2010
Accepted: 9 July 2010
Published: 9 July 2010
The unfolding speed of some hyperthermophilic proteins is dramatically lower than that of their mesostable homologs. Ribonuclease HII from the hyperthermophilic archaeon Thermococcus kodakaraensis (Tk-RNase HII) is stabilized by its remarkably slow unfolding rate, whereas RNase HI from the thermophilic bacterium Thermus thermophilus (Tt-RNase HI) unfolds rapidly, comparable with to that of RNase HI from Escherichia coli (Ec-RNase HI).
To clarify whether the difference in the unfolding rate is due to differences in the types of RNase H or differences in proteins from archaea and bacteria, we examined the equilibrium stability and unfolding reaction of RNases HII from the hyperthermophilic bacteria Thermotoga maritima (Tm-RNase HII) and Aquifex aeolicus (Aa-RNase HII) and RNase HI from the hyperthermophilic archaeon Sulfolobus tokodaii (Sto-RNase HI). These proteins from hyperthermophiles are more stable than Ec-RNase HI over all the temperature ranges examined. The observed unfolding speeds of all hyperstable proteins at the different denaturant concentrations studied are much lower than those of Ec-RNase HI, which is in accordance with the familiar slow unfolding of hyperstable proteins. However, the unfolding rate constants of these RNases H in water are dispersed, and the unfolding rate constant of thermophilic archaeal proteins is lower than that of thermophilic bacterial proteins.
These results suggest that the nature of slow unfolding of thermophilic proteins is determined by the evolutionary history of the organisms involved. The unfolding rate constants in water are related to the amount of buried hydrophobic residues in the tertiary structure.
Proteins from thermophiles and hyperthermophiles generally exhibit higher stability than their mesostable counterparts [1–4]. Recent research in this field has indicated that some proteins from hyperthermophiles are stabilized by their remarkably slow unfolding rates [5–17]. At present, the molecular mechanism of slow unfolding is not completely clear. Recently, using hydrophobic mutant proteins, Dong et al.  demonstrated that the hydrophobic effect is one of the reasons for the slow unfolding of ribonuclease HII from the hyperthermophilic archaeon Thermococcus kodakaraensis (Tk-RNase HII). However, the attributes of proteins that exhibit unusually slow unfolding are unclear. For example, Tk-RNase HII is kinetically robust , whereas RNase HI from the thermophilic bacterium Thermus thermophilus (Tt-RNase HI), which has equivalent thermostability but no sequence similarity to Tk-RNase HII, unfolds rapidly. Moreover, the unfolding rate of Tt-RNase HI is comparable with to that of RNase HI from Escherichia coli (Ec-RNase HI) [20, 21]. The reason for the difference in the unfolding rate is unknown.
Both equilibrium stability and kinetic unfolding/refolding experiments can improve our understanding of the nature of hyperstable proteins. Luke et al.  listed the hyperthermophilic proteins for which kinetic folding/unfolding parameters have been reported. These are generally oligomers or monomers that show irreversible unfolding. Monomeric proteins with reversible unfolding are suitable for detailed analysis of the essentials of slow unfolding. Only RNase H, cold shock protein (Csp), and ribosomal protein S16 (S16) satisfy these prerequisites . Csp from the hyperthermophilic bacterium Thermotoga maritima (Tm-Csp) and ribosomal protein S16 from the hyperthermophilic bacterium Aquifex aeolicus (Aa-S16) do not exhibit extremely slow unfolding [24–27]. Data for the unfolding kinetics of RNase H, Csp, and S16 suggest that differences in the type of protein, i.e., whether it is RNase HI or RNase HII, or differences in the organism kingdom, i.e., whether the host organism to which the protein belongs originates from archaea or bacteria, may be the reason for the differences in the unfolding reactions of Tk-RNase HII and Tt-RNase HI.
Based on their sequence similarity, RNases H are classified into two major types--Type 1 RNase H and Type 2 RNase H [28, 29]. RNase HI is a member of Type 1 RNase H, whereas RNase HII is a member of Type 2 RNase H. Comparison of the crystal structures of Ec-RNase HI, Tt-RNase HI, and Tk-RNase HII indicates that these enzymes share a folding motif referred to as the RNase H-fold due to its characterization in Ec-RNase HI , although there is no sequence similarity between RNase HI and RNase HII. The stability and folding of Ec-RNase HI, Tt-RNase HI, and Tk-RNase HII have been well studied [18–21, 31–40].
Archaea and bacteria have divided at an early stage of evolution . Hyperthermophiles and thermophiles in archaea originated in a hot environment but those in bacteria recolonized at a later stage under extreme conditions. The differences in the unfolding characteristics of hyperstable proteins may be ascribed to such differences in the evolutionary history of thermophilic archaea and bacteria.
Types and origins of RNases H.
Sto; Sulfolobus tokodaii
Ec; Escherichia coli
Tt; Thermus thermophilus
Tk; Thermococcus kodakaraensis
Aa; Aquifex aeolicus
Tm; Thermotoga maritima
CD spectra and reversible unfolding in GdnHCl
The far-UV spectra of Tm-RNase HII, Aa-RNase HII, and Sto-RNase HI were recorded in the absence and presence of GdnHCl at 25°C (see additional file 2). The spectra represent the proteins in their native and GdnHCl denatured states in the absence and presence of GdnHCl. The spectra of the folded condition obtained from the denatured state by dilution of GdnHCl were also recorded. The results demonstrate that GdnHCl-induced unfolding of these proteins is almost completely reversible under each condition (see additional file 2). Experiments on GdnHCl-induced unfolding of these proteins were carried out under conditions of reversible unfolding (see below). It has been shown that Tk-RNase HII, Tt-RNase HI, and Ec-RNase HI undergo reversible unfolding in the presence of denaturants [19, 20, 31]. RNases H are considered to be good models for clarifying the mechanisms of conformational stability and folding of hyperstable proteins. RNase HI from the psychotrophic bacterium Shewanella onedensis MR-1 was also shown to unfold reversibly in the presence of chemical denaturants [46–48].
Equilibrium stability of GdnHCl-induced unfolding and stability profile
Comparison of the thermodynamic parameters of Tm-, Aa-, Tk-RNases HII and Sto-, Tt-, Ec-RNases HI.
(kJ mol-1 K-1)
7.3 ± 0.6
10.4 ± 0.6
14.5 ± 1.8
11.9 ± 1.2
7.5 ± 0.4
11.3 ± 0.8
471 ± 19
549 ± 21
745 ± 49
858 ± 45
548 ± 21
502 ± 17
RNases H from thermophiles and hyperthermophiles are more stable than Ec-RNase HI over most of the temperature ranges examined. In particular, Sto-RNase HI has much higher stability. However, the thermodynamic parameters and stabilization mechanism from the viewpoint of equilibrium stability for the hyperstable RNases H differ depending on the proteins (see Discussion below).
Thermal stability of heat-induced unfolding
Heat-induced unfolding of Sto-RNase HI has already been analyzed by differential scanning calorimetry (DSC) . The Tm of Sto-RNase HI is 102.0°C. The thermal denaturation curves of Tm-RNase HII and Aa-RNase HII were recorded by measuring the change in CD. The thermal denaturation of these proteins was almost completely reversible in the presence of various concentrations of GdnHCl (see Methods). Tm was calculated by curve fitting the resultant CD spectrum versus temperature based on a least-squares analysis using Eq. (4). The Tm values of Tm-RNase HII and Aa-RNase HII in the absence of GdnHCl were estimated to be 78.6°C and 84.0°C, respectively, and were comparable to those of Tt-RNase HII (86.0°C) and Tk-RNase HII (82.8°C) [19, 20].
Kinetics of GdnHCl-induced unfolding
Comparison of the kinetic parameters of Tm-, Aa-, Tk-RNases HII and Sto-, Tt-, Ec-RNases HI at 25°C.
6.5 × 10-7
3.7 × 10-5
6.0 × 10-10
5.7 × 10-11
(1.3 × 10-5)c
1.1 × 10-5
(1.7 × 10-5)d
Proteins from hyperthermophiles often exhibit slower unfolding than those from organisms that grow at moderate temperatures [5–17]. However, the extent of the decrease in the unfolding rate of hyperthermophilic proteins relative to that of their mesophilic counterparts varies, depending on the type of protein involved. The attributes of proteins that exhibit unusually slow unfolding have not been elucidated. In this study, we attempted to characterize super slow unfolding and elucidate its mechanism by examining several RNase H proteins from thermophiles and hyperthermophiles.
Equilibrium stability of hyperstable RNases H
The stability of proteins in solution is quantitatively evaluated by the Gibbs energy changes (ΔG) observed upon unfolding when the reaction is reversible under experimental conditions. The temperature dependence of ΔG (stability profile) provides some information on the thermodynamic stability of proteins and can be expressed by Eq. (3).
Three thermodynamic models have been proposed to explain the high stability of thermostable proteins . In the first model, the entire stability curve is raised to a higher ΔG; in the second, the stability curve is flattened; and in the third, the curve is shifted to a higher temperature. Adherence to all three models is observed in nature, and in many cases, proteins from hyperthermophiles utilize various combinations of these mechanisms [50–55]. The outcome depends on the model used to adapt the proteins to a high temperature.
RNases H from thermophiles and hyperthermophiles are more stable than Ec-RNase HI over most of the temperature ranges examined (Figure 1). In the case of hyperstable RNases HI, Tt-RNase HI adopts the first and second models, and Sto-RNase HI uses the first and third models (Figure 1 and Table 2). There is no data on mesophilic RNase HII. Among the hyperthermophilic RNases HII, Tm-RNase HII has a flatter curve than Aa-RNase HII and Tk-RNase HII as a result of its lower ΔCp. The Ts value of Tm-RNase HII is also lower than that of Aa-RNase HII and Tk-RNase HII. However, ΔG(Ts) of Tk-RNase HII is higher than that of Tm-RNase HII and Aa-RNase HII. These results suggest that the thermodynamic stabilization of RNases H from thermophiles and hyperthermophiles is achieved through different mechanisms.
Unfolding of hyperstable RNases H induced by GdnHCl concentration jump
Estimation of the unfolding speed of hyperstable RNases H at zero denaturant concentration
All the RNases H from the hyperthermophiles examined showed a remarkably slow rate of unfolding in the presence of GdnHCl, as described above. However, the dependencies of the apparent rate constants of the unfolding reactions on the GdnHCl concentration varied according to the type of protein (Figure 3). By extrapolating the unfolding rate constants to zero, the denaturant concentration yields the semilogarithmic unfolding rate constants in water (ln ku(H2O)). It has been reported that ku(H2O) at 25°C is 1.1 × 10-5 s-1 for GdnHCl-induced unfolding and 1.7 × 10-5 s-1 for urea-induced unfolding of Ec-RNase HI. Furthermore, the values are 1.3 × 10-5 s-1 for urea-induced unfolding of Tt-RNase HI and 6.0 × 10-10 s-1 for GdnHCl-induced unfolding of Tk-RNase HII [19, 21, 31, 33]. These values indicate that Tt-RNase HI unfolds in water on a time scale similar to that of Ec-RNase HI; however, the unfolding rate of Tk-RNase HII is dramatically lower than that of Ec-RNase HI. Therefore, the unfolding of Tk-RNase HII is much slower than that of Tt-RNase HI, although the equilibrium stabilities of Tt-RNase HI and Tk-RNase HII are comparable. Furthermore, although there are no unfolding kinetics data for Ec-RNase HII, it is possible that Ec-RNase HII may not exhibit slow unfolding because Ec-RNase HII is less stable than Ec-RNase HI .
The ku(H2O) value for GdnHCl-induced unfolding at 25°C is 6.5 × 10-7 s-1 for Tm-RNase HII, 3.7 × 10-5 s-1 for Aa-RNase HII, and 5.7 × 10-11 s-1 for Sto-RNase HI (Table 3). The results indicate that Aa-RNase HII unfolds in water on a time scale similar to that of Ec-RNase HI and Tt-RNase HI. In contrast, Sto-RNase HI has remarkably slow unfolding characteristics, as in the case of Tk-RNase HII. The unfolding of Tm-RNase HII is slower than that of Ec-RNase HI, but the difference is much less than that between Ec-RNase HI and Sto-RNase HI. Thus, Tm-RNase HII does not show the extremely slow unfolding reactions.
For Tk-RNase HII and Sto-RNase HI, which exhibit slow unfolding, the experimental conditions were pH 9.0  and pH 3.0, respectively, because the proteins unfold reversibly at each of these pH conditions. At neutral pH where these proteins do not show completely reversible unfolding, the ku(H2O) value due to GdnHCl-induced unfolding at 25°C was not very different from the value obtained under reversible conditions. For example, the ku(H2O) value was 9.3 × 10-10 s-1 for Tk-RNase HII at pH 7.5 and 3.3 × 10-11 s-1 for Sto-RNase HI at pH 5.5. Therefore, Tk-RNase HII and Sto-RNase HI unfold extremely slowly at neutral pH.
It should be noted here that the differences in the unfolding rates of hyperstable RNases H reflect the organism kingdom from which the host microorganisms of RNases H originate, i.e., archaea or bacteria. RNases H from thermophilic bacteria, such as Tt-RNase HI, Aa-RNase HII, and Tm-RNase HII unfold faster than similar proteins from thermophilic archaea, Sto-RNase HI, and Tk-RNase HII. This finding is supported by data from other monomeric proteins with reversible unfolding. Tm-Csp and Aa-S16 from hyperthermophilic bacteria do not exhibit extremely slow unfolding [24–27].
Origin of super slow unfolding
As described above, the nature of slow unfolding depends on the evolutionary history of the organisms. What is the difference between thermophilic archaeal and bacterial proteins? Berezovsky & Shakhnovich  reported that proteins from archaea that originated in a hot environment are more compact and hydrophobic than their mesophilic homologs. In contrast, proteins from some bacteria, such as Tm, that recolonized at a later stage under extreme conditions are stabilized by specific interactions such as salt bridges. Mizuguchi et al.  have shown differences in the amino acid composition of archaeal and bacterial proteins. They concluded that thermal adaptation is achieved in different ways in the archaeal and bacterial kingdoms.
We also calculated the buried area of the polar (O and N) and nonpolar (C and S) atoms in the interior of the proteins because the buried area contributes to protein stability . The buried areas of nonpolar atoms in Tt-RNase HI, Sto-RNase HI, and Tk-RNase HII were 62.6%, 67.3%, and 64.5%, respectively (see additional file 6). Surprisingly, there was no significant difference between bacterial and archaeal proteins, although the hydrophobicity of archaeal RNases H (Sto-RNase HI and Tk-RNase HII) was higher than that of bacterial RNase H (Tt-RNase HI). Therefore, we suggest that the overall hydrophobicity, calculated from the buried area of atoms, does not differ greatly between bacterial and archaeal proteins and contributes to their equilibrium stability. The equlibrium stability of Sto-RNase HI, which is the most hydrophobic RNase H, is considerably higher than that of Tt-RNase HI and Tk-RNase HII (Figure 1). It has been reported that 100 Å2 of buried nonpolar atoms contributes to an equilibrium stability of about 15 kJ mol-1 . In contrast, the estimate from the sequence/structure data (the fractions of amino acid residues in the interior of proteins) differs between bacterial and archaeal proteins [57, 58], as described above. This suggests that the hydrophobic interaction arising from the partially exposed hydrophobic side chains may affect the unfolding rate (see below).
Figures 4D and 4E represent the amino acid sequence alignments, buried residues (Tt, Sto and Tk), and hydrophobic residues of thermophilic RNases H. In RNases HI, Sto-RNase HI has 24 buried hydrophobic residues with nonburied or nonhydrophobic counterparts in Tt-RNase HI, whereas Tt-RNase HI has 11 such residues. Tk-RNase HII also possesses 10 buried hydrophobic residues with nonhydrophobic counterparts in Aa-RNase HII and Tm-RNase HII. These residues in Sto-RNase HI and Tk-RNase HII may be possible causes of the slow unfolding rate. These results will be confirmed in the near future by amino acid substitution experiments. Figures 4F, 4G, and 4H show the crystal structures of the buried hydrophobic residues of Tt-RNase HI, Sto-RNase HI, and Tk-RNase HII. Among RNases HI, Sto-RNase HI has more buried hydrophobic residues than Tt-RNase HI. Furthermore, it seems that the the 24 buried hydrophobic residues in Sto-RNase HI, which have nonburied or nonhydrophobic counterparts in Tt-RNase HI, are located near the surface of the protein. Tk-RNase HII also shows a similar tendency. The results suggest that buried hydrophobic side-chains near the protein surface may contribute to slow unfolding. Unfortunately, at present, the mechanism of slow unfolding by buried hydrophobic side-chains near the protein surface remains unclear.
Energy diagram of the stability and folding of hyperstable proteins
It has been reported that some hyperthermophilic proteins have a remarkably slow unfolding rate but a normal refolding rate in comparison with their homologous mesophilic proteins [6, 8]. Tk-RNase HII unfolds extremely slowly but folds rapidly . In contrast, Tt-RNase HI, which has thermostability equivalent to that of Tk-RNase HII, unfolds and folds quickly, similar to the behavior of Ec-RNase HI [19, 21]. We can consider the energy diagrams of the folding process (see additional file 7). The first is for a mesophilic protein and is a two-state model. The second is for hyperstable proteins with super slow unfolding, whereas the third is for hyperstable proteins without super slow unfolding. In the second model, the increase in the equilibrium stability (ΔG) between the native and denatured states is accompanied by an increase in the unfolding barrier (N to T). In contrast, the third model includes an intermediate state in the folding process, which results in an increase in ΔG without an increase in either the folding or unfolding barriers. It has been reported that although Tt-RNase HI has an intermediate folding state, Tk-RNase HII folds according to the two-state model [19, 36, 60]. Therefore, Tt-RNase HI and Tk-RNase HII are represented by the third and second models, respectively. For Aa-RNase HI and Tm-RNase HI, the refolding experiments show that there is an intermediate state in the folding process (data not shown), suggesting that these proteins can be represented by the third model. When Sto-RNase HI folds according to the two-state model, kf(H2O) = 9.3 × 10 s-1, mf = -2.98 M-1 s-1, and ΔGkinetics(H2O) = 73.9 kJ mol-1 (unpublished data). This indicates that Sto-RNase HI is represented by the second model. On the other hand, the tetrameric protein Pf-PCP both folds and unfolds quite slowly . Therefore, this protein is represented by the fourth model, in which the increase in equilibrium stability (ΔG) between the native and denatured states is accompanied by increases in both the unfolding (N to T) and folding (D to T) barriers. As a result, the unfolding barrier of the fourth model is higher than that of the second model. The unfolding rate of Pf-PCP (1.6 × 10-15 s-1) is lower than those of Tk-RNase HII and Sto-RNase HI (6.0 × 10-10 s-1 and 5.7 × 10-11 s-1) .
In summary, we have presented the stability and unfolding kinetics of hyperstable RNases H. This is a good protein family for investigating the stability and folding mechanisms from the viewpoints of thermal adaptation and evolution because the proteins are in the monomeric form and undergo reversible unfolding. Although different mechanisms are responsible for the equilibrium stabilization of hyperstable RNases H, unfolding depends on the evolutional history of the host organisms from which the proteins originate. This has led to the use of different strategies for thermal adaptation by thermophilic archaea and bacteria.
Construction of plasmids
Plasmids pET800TM and pET600AA, used for overproduction of Tm-RNase HII and Aa-RNase HII, were constructed by ligating 714bp and 588bp DNA fragments, respectively, into the NdeI-SalI sites of pET25b (Novagen). The DNA fragments were amplified by PCR using the genomic DNA of Tm and Aa as a template. The sequences of the PCR primers for Tm-RNase HII are 5'-GAAAGGAGACGGCATATGGGAATAGATGA-3' for the 5' primer and 5'-GCGGTCGACCTTTGTGTTCGGAGATAAAG-3' for the 3' primer, and the PCR primers for Aa-RNase HII are 5'-AACCATATGCTTAATTACGAATTAGAACTTTG-3' for the 5' primer and 5'-GCGGTCGACCTAAAAAAGTTTTTTCTTAAGGGG-3' for the 3' primer, where the bold bases indicate the positions of the NdeI and SalI sites. The DNA sequence of the PCR products was confirmed with a Prism 310 DNA sequencer (Perkin-Elmer).
Overproduction and purification
E. coli BL21(DE3) was transformed with pET800TM and pET600AA then grown in NZCYM medium (Novagen) containing 50 mgL-1 ampicillin 30 mgL-1 chloramphenicol at 30°C. When the optical density at 660 nm of the culture reached 0.5, 1 mM IPTG was added to the culture medium and cultivation continued at 30°C for four hours. Cells were transformed with pET800TM, harvested by centrifugation at 8000 rpm for 10 min, suspended in 20 mM Tris-HCl (pH 8) containing 10 mM EDTA, disrupted by sonication lysis, and centrifuged at 30,000 g for 30 min. The supernatant was subjected to heat treatment at 75°C for 20 min and centrifuged at 30,000 g for 30 min. The supernatant was applied to an anion-exchange column (5 mL) of a HiTrap Q HP column (Amersham Biosciences) equilibrated with 20 mM Tris-HCl (pH 8). The protein was eluted from the column with 20 mM Tris-HCl (pH 8) containing 10 mM EDTA. The flow-through containing the protein was pooled and applied to a Hitrap Heparin and SP column (Amersham Biosciences) equilibrated with the same buffer. After the column was washed, the protein was eluted from the column with a linear gradient of NaCl from 0 to 1 M in the same buffer. The fractions containing the protein were pooled and dialyzed against 20 mM Tris-HCl (pH 8) containing 200 mM NaCl. Cells were transformed with pET600AA, harvested by centrifugation at 8000 rpm for 10 min, suspended in 20 mM Acetate-NaOH (pH 5.5) containing 10 mM EDTA, disrupted by French press lysis, and centrifuged at 30,000 g for 30 min. The supernatant was applied to an anion-exchange column (5 mL) of a HiTrap Q HP column equilibrated with 20 mM Acetate-NaOH (pH 5.5). The protein was eluted from the column with 20 mM Acetate-NaOH (pH 5.5) containing 10 mM EDTA. The flow-through containing the protein was pooled and applied to a Hitrap SP column equilibrated with the same buffer. After the column was washed, the protein was eluted from the column with a linear gradient of NaCl from 0 to 1 M in the same buffer. The fractions containing the protein were pooled and dialyzed against 20 mM Acetate-NaOH (pH 5.0). The fractions containing the protein were pooled and applied to a Hitrap Heparin column equilibrated with the same buffer. After the column was washed, the protein was eluted from the column with a linear gradient of NaCl from 0 to 1 M in the same buffer. The fractions containing the protein were pooled and dialyzed against 20 mM Acetate-NaOH (pH 5.0).
Sto-RNase HI was overproduced and purified as previously described [44, 61]. The purity of the protein was analyzed by SDS-PAGE using a 15% polyacrylamide gel, followed with Coomassie Brilliant Blue staining.
The protein concentrations were estimated by assuming absorbances at 280 nm of 0.24 for 1 mgmL-1 Tm-RNase HII and 0.89 for 1 mgmL-1 Aa-RNase HII. These values were calculated using 1576 M-1cm-1 for Tyr and 5225 M-1cm-1 for Trp at 280 nm .
CD spectra measurements
The CD spectra of Tm-RNase HII, Aa-RNase HII, and Sto-RNase HI in the absence and presence of GdnHCl were measured on a J-725 automatic spectropolarimeter (Japan Spectroscopic). The optical path length was 2 mm and the protein concentration was 0.16 mgml-1. The buffers we used were 20 mM Tris-HCl, 200 mM NaCl at pH 7.5 for Tm-RNase HII, 20 mM Acetate-NaOH at pH 5.0 for Aa-RNase HII, and 20 mM Glycine-NaOH at pH 3.0 for Sto-RNase HI. The measurements were done at 25°C. For the spectra of a refolded protein, the protein that was completely unfolded at a 7.0 M GdnHCl concentration was diluted with buffer for refolding, and the diluted protein solution was incubated at the selected temperature until refolding reached equilibrium. The mean residue ellipticity, θ, in units of degrees square centimeters per decimole, was calculated using an average amino acid molecular weight of 110.
Equilibrium experiments on GdnHCl-induced unfolding
where y is the observed CD signal at a given concentration of GdnHCl, [D] is the concentration of GdnHCl, b0n is the CD signal for the native state, b0u is the CD signal for the unfolded states, an is the slope of the pre-transition of the baseline, and au is the slope of the posttransition of the baseline. ΔG(H2O) is the Gibbs energy change (ΔG) of the unfolding in the absence of GdnHCl, m is the slope of the linear correlation between ΔG and the GdnHCl concentration [D], and Cm is the GdnHCl concentration at the midpoint of the curve. The raw experiment data were directly fitted to Eq. (1) using SigmaPlot (Jandel Scientific). The equilibrium experiments were conducted in 20 mM Tris-HCl, 200 mM NaCl at pH 7.5 for Tm-RNase HII, 20 mM Acetate-NaOH at pH 5.0 for Aa-RNase HII, and 20 mM Glycine-NaOH at pH 3.0 for Sto-RNase HI. To use Tris-HCl buffer at different temperatures, the pH was adjusted at 25°C so that the final pH would be 7.5 at the desired temperature . The protein concentration when the CD measurements were taken was 0.16 mgmL-1.
where ΔH(To) and ΔS(To) are the enthalpy and entropy of unfolding at the reference temperature To, and ΔCp is the difference in heat capacity between the native and unfolded states.
Heat-induced unfolding experiments
where y is the observed CD signal at a given temperature [T], bn is the CD signal for the native state, bu is the CD signal for the unfolded states, an is the slope of the pretransition of the baseline, au is the slope of the posttransition of the baseline, ΔHm is the enthalpy of unfolding at the transition midpoint temperature (Tm), T is the temperature, and R is the gas constant. Curve fitting was performed using SigmaPlot. Because the Tm values of these proteins increased in proportion to a decrease in GdnHCl concentration, the Tm values in the absence of GdnHCl were estimated by linear extrapolation .
Kinetic experiments on GdnHCl-induced unfolding
where [D] is the concentration of GdnHCl and mu represents the slopes of the linear correlations of ln ku with the GdnHCl concentration.
cold shock protein
Gibbs energy change (ΔG) of the unfolding in water
pyrrolidone carboxyl peptidase
ribosomal protein S16
thermal denaturation temperature
This work was supported in part by an Industrial Technology Research Grant Program from the New Energy and Industrial Technology Development Organization (NEDO) of Japan, and by a Grant from the Ministry of Education, Culture, Sports, Science, and Technology of Japan.
- Jaenicke R, Bohm G: The stability of proteins in extreme environments. Curr Opin Struct Biol. 1998, 8: 738-748. 10.1016/S0959-440X(98)80094-8.View ArticlePubMedGoogle Scholar
- Vieille C, Zeikus GJ: Hyperthermophilic enzymes: sources, uses, and molecular mechanisms for thermostability. Microbiol Mol Biol Rev. 2001, 65: 1-43. 10.1128/MMBR.65.1.1-43.2001.PubMed CentralView ArticlePubMedGoogle Scholar
- Perl D, Schmid FX: Some like it hot: The molecular determinants of protein thermostability. ChemBioChem. 2002, 3: 39-44. 10.1002/1439-7633(20020104)3:1<39::AID-CBIC39>3.0.CO;2-D.View ArticlePubMedGoogle Scholar
- Razvi A, Scholtz JM: Lessons in stability from thermophilic proteins. Protein Sci. 2006, 15: 1569-1578. 10.1110/ps.062130306.PubMed CentralView ArticlePubMedGoogle Scholar
- Cavagnero S, Debe DA, Zhou ZH, Adams MWW, Chan SI: Kinetic role of electrostatic interactions in the unfolding of hyperthermophilic and mesophilic rubredoxins. Biochemsitry. 1998, 37: 3369-3376. 10.1021/bi9721795.View ArticleGoogle Scholar
- Ogasahara K, Nakamura M, Nakura S, Tsunasawa S, Kato I, Yoshimoto T, Yutani K: The unusually slow unfolding rate causes the high stability of pyroolidone carboxyl peptidase from a hyperthermophilile Pyrococcus furiosus: Equilibrium and kinetic studies of guanidine hydrochloride-induced unfolding and refolding. Biochemistry. 1998, 37: 17537-17544. 10.1021/bi9814585.View ArticlePubMedGoogle Scholar
- Dams T, Jaenicke R: Stability and folding of dihydrofolate reductase from the hyperthermophilic bacterium Thermotoga maritima. Biochemistry. 1999, 38: 9169-9178. 10.1021/bi990635e.View ArticlePubMedGoogle Scholar
- Kaushik JK, Ogasahara K, Yutani K: The unusually slow relaxation kinetics of the folding-unfolding of pyrrolidone carboxyl peptidase from a hyperthermophile. Pyrococcus furiosus. J Mol Biol. 2002, 316: 991-1003. 10.1006/jmbi.2001.5355.View ArticlePubMedGoogle Scholar
- Jaswal SS, Sohl JL, Dans JH, Agard DA: Energetic landscape of alpha-lytic protease optimizes longevity through kinetic stability. Nature. 2002, 415: 343-346. 10.1038/415343a.View ArticlePubMedGoogle Scholar
- Iimura S, Yagi H, Ogasahara K, Akutsu H, Noda Y, Segawa S, Yutani K: Unusually slow denaturation and refolding process of pyrrolidone carcoxyl peptidase from a hyper-thermophile are highly cooperative: Real-time NMR studies. Biochemistry. 2004, 43: 11906-11915. 10.1021/bi048762k.View ArticlePubMedGoogle Scholar
- Zeeb M, Lipps G, Lilie H, Balbach J: Folding and association of an extremely stable dimeric protein from Sulfolobus islandicus. J Mol Biol. 2004, 336: 227-240. 10.1016/j.jmb.2003.12.003.View ArticlePubMedGoogle Scholar
- Forrer P, Chang C, Ott D, Wlodawer A, Plückthun A: Kinetic stability and crystal structure of the viral capside protein SHP. J Mol Biol. 2004, 344: 179-193. 10.1016/j.jmb.2004.09.030.View ArticlePubMedGoogle Scholar
- Wittung-Stafshede P: Slow unfolding explains high stability of thermostable ferredoxins: common mechanism governing thermostability?. Biochim Biophys Acta. 2004, 1700: 1-4.View ArticlePubMedGoogle Scholar
- Duy C, Fitter J: Thermostability of irreversible unfolding alpha-amylases analyzed by unfolding kinetics. J Biol Chem. 2005, 280: 37360-37365. 10.1074/jbc.M507530200.View ArticlePubMedGoogle Scholar
- Kaushik JK, Iimura S, Ogasahara K, Yamagata Y, Segawa S, Yutani K: Completely buried, non-ion-paired glutamic acid contributes favorably to the conformational stability of pyrrolidone carboxyl peptidases from hyperthermophiles. Biochemistry. 2006, 45: 7100-7112. 10.1021/bi052610n.View ArticlePubMedGoogle Scholar
- Mishra R, Olofsson L, Karlsson M, Carlsson U, Nicholls IA, Hammarström P: A conformationally isoformic thermophilic protein with high kinetic unfolding barriers. Cell Mol Life Sci. 2008, 65: 827-839. 10.1007/s00018-008-7517-4.View ArticlePubMedGoogle Scholar
- Costas M, Rodríguez-Larrea D, De Maria L, Borchert TV, Gómez-Puyou A, Sanchez-Ruiz JM: Between-species variation in the kinetic stability of TIM proteins linked to solvation-barrier free energies. J Mol Biol. 2009, 385: 924-937. 10.1016/j.jmb.2008.10.056.View ArticlePubMedGoogle Scholar
- Dong H, Mukaiyama A, Tadokoro T, Koga Y, Takano K, Kanaya S: Hydrophobic effect on the stability and folding of a hyperthermophilic protein. J Mol Biol. 2008, 378: 264-272. 10.1016/j.jmb.2008.02.039.View ArticlePubMedGoogle Scholar
- Mukaiyama A, Takano K, Haruki M, Morikawa M, Kanaya S: Kinetically robust monomeric protein from a hyperthermophile. Biochemistry. 2004, 43: 13859-13866. 10.1021/bi0487645.View ArticlePubMedGoogle Scholar
- Hollien J, Marqusee S: A thermodynamic comparison of mesophilic and thermophilic ribonucleases H. Biochemistry. 1999, 38: 3831-3836. 10.1021/bi982684h.View ArticlePubMedGoogle Scholar
- Hollien J, Marqusee S: Comparison of the folding processes of T. thermophilus and E. coli ribonucleases H. J Mol Biol. 2002, 316: 327-340. 10.1006/jmbi.2001.5346.View ArticlePubMedGoogle Scholar
- Luke KA, Higgins CL, Wittung-Stafshede P: Thermodynamic stability and folding of proteins from hyperthermophilic organisms. FEBS J. 2007, 274: 4023-4033. 10.1111/j.1742-4658.2007.05955.x.View ArticlePubMedGoogle Scholar
- Mukaiyama A, Takano K: Slow unfolding of monomeric proteins from hyperthermophiles with reversible unfolding. Int J Mol Sci. 2009, 10: 1369-1385. 10.3390/ijms10031369.PubMed CentralView ArticlePubMedGoogle Scholar
- Perl D, Welker C, Schindler T, Schröder K, Marahiel MA, Janicke R, Schmid FX: Conservation of rapid two-state folding in mesophilic, thermophilic and hyperthermophilic cold shock proteins. Nature Struct Biol. 1998, 5: 229-235. 10.1038/nsb0398-229.View ArticlePubMedGoogle Scholar
- Dominy BN, Perl D, Schmid FX, Brooks CLIII: The effect of ionic strength on protein stability: The cold shock protein family. J Mol Biol. 2002, 319: 541-554. 10.1016/S0022-2836(02)00259-0.View ArticlePubMedGoogle Scholar
- Schuler B, Kremer W, Kalbitzer HR, Jaenicke R: Role of entropy in protein thermostability: folding kinetics of a hyperthermophilic cold shock protein at high temperatures using 19F NMR. Biochemistry. 2002, 41: 11670-11680. 10.1021/bi026293l.View ArticlePubMedGoogle Scholar
- Wallgren M, Adén J, Pylypenko O, Mikaelsson T, Johansson LB, Rak A, Wolf-Watz M: Extreme temperature tolerance of a hyperthermophilic protein coupled to residual structure in the unfolded state. J Mol Biol. 2008, 379: 845-858. 10.1016/j.jmb.2008.04.007.View ArticlePubMedGoogle Scholar
- Ohtani N, Haruki M, Morikawa M, Crouch RJ, Itaya M, Kanaya S: Identifiacation of the genes encoding Mn2+-dependent RNase HII and Mg2+-dependent RNase HIII from Bacillus subtilis: Classification of RNase H into these families. Biochemistry. 1999, 38: 605-618. 10.1021/bi982207z.View ArticlePubMedGoogle Scholar
- Ohtani N, Haruki M, Morikawa M, Kanaya S: Molecular diversities of RNase H. J Biosci Bioeng. 1999, 88: 12-19. 10.1016/S1389-1723(99)80168-6.View ArticlePubMedGoogle Scholar
- Yang W, Steitz TA: Recombining the structures of HIV integrase, RuvC and RNase H. Structure. 1995, 3: 131-134. 10.1016/S0969-2126(01)00142-3.View ArticlePubMedGoogle Scholar
- Yamasaki K, Ogasahara K, Yutani K, Oobatake M, Kanaya S: Folding pathway of Escherichia coli ribonucrease H: A circular dichroism, fluorescence and NMR study. Biochemistry. 1995, 34: 16552-16562. 10.1021/bi00051a003.View ArticlePubMedGoogle Scholar
- Haruki M, Tanaka M, Motegi T, Tadokoro T, Koga Y, Takano K, Kanaya S: Structural and thermodynamic analyses of Escherichia coli RNase HI variant with quintuple thermostabilizing mutations. FEBS J. 2007, 274: 5815-5825. 10.1111/j.1742-4658.2007.06104.x.View ArticlePubMedGoogle Scholar
- Raschke TM, Marqusee S: The kinetic folding intermediate of ribonuclease H resembles the acid molten globule and partially unfolded molecules detected under native conditions. Naure Struct Biol. 1997, 4: 298-304. 10.1038/nsb0497-298.View ArticleGoogle Scholar
- Raschke TM, Kho J, Marqusee S: Confirmation of the hierarchical folding of RNase H: A protein engineering study. Naure Struct Biol. 1999, 6: 825-831. 10.1038/12277.View ArticleGoogle Scholar
- Spudich GM, Miller EJ, Marquesee S: Destabilization of the Escherichia coli RNase H kinetic intermediate: Switching between a two-state and three-state folding mechanism. J Mol Biol. 2004, 335: 609-618. 10.1016/j.jmb.2003.10.052.View ArticlePubMedGoogle Scholar
- Zhou Z, Feng H, Ghirlando R, Bai Y: The high-resolution NMR structure of the early folding intermediate of the Thermus thermophilus ribonuclease H. J Mol Biol. 2008, 384: 531-539. 10.1016/j.jmb.2008.09.044.PubMed CentralView ArticlePubMedGoogle Scholar
- Mukaiyama A, Haruki M, Ota M, Koga Y, Takano K, Kanaya S: A hyperthermophilic protein acquires function at the cost of stability. Biochemistry. 2006, 45: 12673-12679. 10.1021/bi060907v.View ArticlePubMedGoogle Scholar
- Mukaiyama A, Koga Y, Takano K, Kanaya S: Osmolyte effect on the stability and folding of a hyperthermophilic protein. Proteins: Struct Funct Bioinf. 2008, 71: 110-118. 10.1002/prot.21660.View ArticleGoogle Scholar
- Takano K, Higashi R, Okada J, Mukaiyama A, Tadokoro T, Koga Y, Kanaya S: Proline effect on the thermostability and slow unfolding of a hyperthermophilic protein. J Biochem. 2008, 145: 79-85. 10.1093/jb/mvn144.View ArticlePubMedGoogle Scholar
- Rohman MS, Koga Y, Takano K, Chon H, Crouch RJ, Kanaya S: Effect of the disease-causing mutations identified in human ribonuclease (RNase) H2 on the activities and stabilities of yeast RNase H2 and archaeal RNase HII. FEBS J. 2008, 275: 4836-4849. 10.1111/j.1742-4658.2008.06622.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Woese CR, Kandler O, Wheelis ML: Towards a natural system of organisms: Proposal for the domains Archaea, Bacteria, and Eucarya. Proc Natl Acad Sci USA. 1990, 87: 4576-4579. 10.1073/pnas.87.12.4576.PubMed CentralView ArticlePubMedGoogle Scholar
- Katayanagi K, Miyagawa M, Matsushima M, Ishikawa M, Kanaya S, Nakamura H, Ikehara M, Matsuzaki T, Morikawa K: Structural details of ribonuclease H from Escherichia coli as refined to an atomic resolution. J Mol Biol. 1992, 223: 1029-1052. 10.1016/0022-2836(92)90260-Q.View ArticlePubMedGoogle Scholar
- Ishikawa K, Okumura M, Katayanagi K, Kimura S, Kanaya S, Nakamura H, Morikawa K: Crystal structure of ribonuclease H from Thermus thermophilus HB8 refined at 2.8 Å resolution. J Mol Biol. 1993, 230: 529-542. 10.1006/jmbi.1993.1169.View ArticlePubMedGoogle Scholar
- You DJ, Chon H, Koga Y, Takano K, Kanaya S: Crystal structure of type 1 ribonuclease H from hyperthermophilic archaeon Sulfolobus tokodaii: role of arginine 118 and C-terminal anchoring. Biochemistry. 2007, 461: 1494-1503.Google Scholar
- Muroya A, Tsuchiya D, Ishikawa M, Haruki M, Morikawa M, Kanaya S, Morikawa K: Catalytic center of an archaeal type2 ribonuclease H as revealed by X-ray crystallographic and mutational analyses. Protein Sci. 2001, 10: 707-714. 10.1110/ps.48001.PubMed CentralView ArticlePubMedGoogle Scholar
- Tadokoro T, You DJ, Abe Y, Chon H, Matsumura H, Koga Y, Takano K, Kanaya S: Structural, thermodynamic, and mutational analyses of a psychrotrophic RNase HI. Biochemistry. 2007, 46: 7460-7468. 10.1021/bi7001423.View ArticlePubMedGoogle Scholar
- Tadokoro T, Matsushita K, Abe Y, Rohman MS, Koga Y, Takano K, Kanaya S: Remarkable stabilization of a psychrotrophic RNase HI by a combination of thermostabilizing mutations identified by the suppressor mutation method. Biochemistry. 2008, 47: 8040-8047. 10.1021/bi800246e.View ArticlePubMedGoogle Scholar
- Rohman MS, Tadokoro T, Angkawidjaja C, Abe Y, Matsumura H, Koga Y, Takano K, Kanaya S: Destabilization of psychrotrophic RNase HI in a localized fashion as revealed by mutational and X-ray crystallographic analyses. FEBS J. 2009, 276: 603-613. 10.1111/j.1742-4658.2008.06811.x.View ArticlePubMedGoogle Scholar
- Nojima H, Ikai A, Oshima T, Noda H: Reversible thermal unfolding of thermostable phosphoglycerate kinase. Thermostability associated with mean zero enthalpy change. J Mol Biol. 1977, 116: 429-442. 10.1016/0022-2836(77)90078-X.View ArticlePubMedGoogle Scholar
- Razvi A, Scholtz JM: A thermodynamic comparison of HPr proteins from extremophilic organisms. Biochemistry. 2006, 45: 4084-4092. 10.1021/bi060038+.View ArticlePubMedGoogle Scholar
- McCrary BS, Edmondson SP, Shriver JW: Hyperthermophile protein folding thermodynamics: Differential scanning calorimetry and chemical denaturation of Sac7d. J Mol Biol. 1996, 264: 784-805. 10.1006/jmbi.1996.0677.View ArticlePubMedGoogle Scholar
- Li WT, Grayling RA, Sandman K, Edmondson S, Shriver JW, Reeve JN: Thermodynamic stability of archaeal histones. Biochemistry. 1998, 37: 10563-10572. 10.1021/bi973006i.View ArticlePubMedGoogle Scholar
- Grättinger M, Dankesreiter A, Schurig H, Jaenicke R: Recombinant phosphoglycerate kinase from the hyperthermophilic bacterium Thermotoga maritima: Catalytic, spectral and thermodynamic properties. J Mol Biol. 1998, 280: 525-533. 10.1006/jmbi.1998.1861.View ArticlePubMedGoogle Scholar
- Deutschman WA, Dahlquist FW: Thermodynamic basis for the increased thermostability of CheY from the hyperthermophile Thermotoga maritima. Biochemistry. 2001, 40: 13107-13113. 10.1021/bi010665t.View ArticlePubMedGoogle Scholar
- Lee CF, Allen MD, Bycroft M, Wong KD: Electrostatic interactions contribute reduced heat capacity change of unfolding in a thermophilic ribosomal protein L30e. J Mol Biol. 2005, 348: 419-431. 10.1016/j.jmb.2005.02.052.View ArticlePubMedGoogle Scholar
- Chon H, Tadokoro T, Ohtani N, Koga Y, Takano K, Kanaya S: Identification of RNase HII from psychrotrophic bacterium. Shewanella sp, SIB1 as a high-activity type RNase H. FEBS J. 2006, 273: 2264-2275. 10.1111/j.1742-4658.2006.05241.x.View ArticlePubMedGoogle Scholar
- Berezovsky IN, Shakhnovich EI: Physics and evolution of thermophilic adaptation. Proc Natl Acad Sci USA. 2005, 102: 12742-12747. 10.1073/pnas.0503890102.PubMed CentralView ArticlePubMedGoogle Scholar
- Mizuguchi K, Sele M, Cubellis MV: Environment specific substitution tables for thermophilic proteins. BMC Bioinformatics. 2007, 8 (Suppl 1): S15-10.1186/1471-2105-8-S1-S15.PubMed CentralView ArticlePubMedGoogle Scholar
- Funahashi J, Takano K, Yutani K: Are the parameters of various stabilization factors estimated from mutant human lysozymes compatible with other proteins?. Protein Eng. 2001, 14: 127-134. 10.1093/protein/14.2.127.View ArticlePubMedGoogle Scholar
- Ratcliff K, Corn J, Marqusee S: Structure, stability, and folding of ribonuclease H1 from the moderately thermophilic Chlorobium tepidum: comparison with thermophilic and mesophilic homologues. Biochemistry. 2009, 48: 5890-5898. 10.1021/bi900305p.PubMed CentralView ArticlePubMedGoogle Scholar
- You DJ, Chon H, Koga Y, Takano K, Kanaya S: Crystallization and preliminary crystallographic analysis of type 1 RNase H from the hyperthermophilic archaeon Sulfolobus tokodaii 7. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2006, 62: 781-784. 10.1107/S1744309106024420.PubMed CentralView ArticlePubMedGoogle Scholar
- Goodwin TW, Morton RA: The spectrophotometric determination of tyrosine and tryptophan in proteins. Biochem J. 1946, 40: 628-632.PubMed CentralView ArticlePubMedGoogle Scholar
- Pace CN: Measuring and increasing protein stability. Trends Biotechnol. 1990, 8: 93-98. 10.1016/0167-7799(90)90146-O.View ArticlePubMedGoogle Scholar
- Good NE, Winget GD, Winter W, Connolly TN, Izawa S, Singh RMM: Hydrogen ion buffers for biological research. Biochemistry. 1966, 5: 467-477. 10.1021/bi00866a011.View ArticlePubMedGoogle Scholar
- Takano K, Saito M, Morikawa M, Kanaya S: Mutational and structural-based analyses of the osmolyte effect on protein stability. J Biochem. 2004, 135: 701-708. 10.1093/jb/mvh085.View ArticlePubMedGoogle Scholar
- Santoro MM, Bolen BW: Unfolding free energy changes determined by the linear extrapolation method. 1. Unfolding of phenylmethanesufonyl alpha-chymotrypsin using different denaturants. Biochemistry. 1988, 27: 8063-8068. 10.1021/bi00421a014.View ArticlePubMedGoogle Scholar
- You DJ, Fukuchi S, Nishikawa K, Koga Y, Takano K, Kanaya S: Protein thermostabilization requires a fine-tuned placement of surface-charged residues. J Biochem. 2007, 142: 507-516. 10.1093/jb/mvm157.View ArticlePubMedGoogle Scholar
- DeLano WL: PyMOL User's Guide. 2004, DeLano Scientific San Carlos, CaliforniaGoogle Scholar