Electrostatic stabilization of a thermophilic cold shock protein.

The cold shock protein Bc-Csp from the thermophile Bacillus caldolyticus differs from its mesophilic homolog Bs-CspB from Bacillus subtilis by 15.8 kJ mol(-1) in the Gibbs free energy of denaturation (DeltaG(D)). The two proteins vary in sequence at 12 positions but only two of them, Arg3 and Leu66 of Bc-Csp, which replace Glu3 and Glu66 of Bs-CspB, are responsible for the additional stability of Bc-Csp. These two positions are near the ends of the protein chain, but close to each other in the three-dimensional structure. The Glu3Arg exchange alone changed the stability by more than 11 kJ mol(-1). Here, we elucidated the molecular origins of the stability difference between the two proteins by a mutational analysis. Electrostatic contributions to stability were characterized by measuring the thermodynamic stabilities of many variants as a function of salt concentration. Double and triple mutant analyses indicate that the stabilization by the Glu3Arg exchange originates from three sources. Improved hydrophobic interactions of the aliphatic moiety of Arg3 contribute about 4 kJ mol(-1). Another 4 kJ mol(-1) is gained from the relief of a pairwise electrostatic repulsion between Glu3 and Glu66, as in the mesophilic protein, and 3 kJ mol(-1) originate from a general electrostatic stabilization by the positive charge of Arg3, which is not caused by a pairwise interaction. Mutations of all potential partners for an ion pair within a radius of 10 A around Arg3 had only marginal effects on stability. The Glu3-->Arg3 charge reversal thus optimizes ionic interactions at the protein surface by both local and global effects. However, it cannot convert the coulombic repulsion with another Glu residue into a corresponding attraction. Avoidance of unfavorable coulombic repulsions is probably a much simpler route to thermostability than the creation of stabilizing surface ion pairs, which can form only at the expense of conformational entropy.     

Stabilization of the cold shock protein CspB from Bacillus subtilis by evolutionary optimization of Coulombic interactions

The bacterial cold shock proteins (Csp) are used by both experimentalists and theoreticians as model systems for analyzing the Coulombic contributions to protein stability. We employ Proside, a method of directed evolution, to identify stabilized variants of Bs-CspB from Bacillus subtilis. Proside links the increased protease resistance of stabilized protein variants to the infectivity of a filamentous phage. Here, three cspB libraries were used for in vitro selections to explore the stabilizing potential of charged amino acids in Bs-CspB. In the first library codons for nine selected surface residues were partially randomized, in the second one random mutations were introduced non-specifically by error-prone PCR, and in the third one the spontaneous mutation rate of the phage in Escherichia coli was used. Stabilizing mutations were found at the surface positions 1, 3, 46, 48, 65, and 66. The contributions of these mutations to stability were characterized by analyzing them individually and in combination. The best combination (M1R, E3K, K65I, and E66L) increased the midpoint of thermal unfolding of Bs-CspB from 53.8 to 85.0 degrees C. The effects of most mutations are strongly context dependent. A good example is provided by the E3R mutation. It is strongly stabilizing (DeltaDeltaGD=11.1kJ mol(-1)) in the wild-type protein, but destabilizing (DeltaDeltaGD=-4.0kJ mol(-1)) in the A46K/S48R/E66L variant. The stabilizations by charge mutations did not correlate well with the corresponding changes in the protein net charge, and they could not be ascribed to the formation of ion pairs. Previous theoretical analyses did not identify the stabilization caused by the mutations at positions 1, 46, and 48. Also, electrostatics calculations based on protein net charge or charge asymmetry did not predict well the stability changes that occur when charged residues in Bs-CspB are mutated. It remains a challenge to model the Coulombic interactions of charged residues in a protein and to determine their contributions to the Gibbs free energy of protein folding.     

Optimized variants of the cold shock protein from in vitro selection: structural basis of their high thermostability

The bacterial cold shock proteins (Csp) are widely used as models for the experimental and computational analysis of protein stability. In a previous study, in vitro evolution was employed to identify strongly stabilizing mutations in Bs-CspB from Bacillus subtilis. The best variant found by this approach contained the mutations M1R, E3K and K65I, which raised the midpoint of thermal unfolding of Bs-CspB from 53.8 degrees C to 83.7 degrees C, and increased the Gibbs free energy of stabilization by 20.9 kJ mol(-1). Another selected variant with the two mutations A46K and S48R was stabilized by 11.1 kJ mol(-1). To elucidate the molecular basis of these stabilizations, we determined the crystal structures of these two Bs-CspB variants. The mutated residues are generally well ordered and provide additional stabilizing interactions, such as charge interactions, additional hydrogen bonds and improved side-chain packing. Several mutations improve the electrostatic interactions, either by the removal of unfavorable charges (E3K) or by compensating their destabilizing interactions (A46K, S48R). The stabilizing mutations are clustered at a contiguous surface area of Bs-CspB, which apparently is critically important for the stability of the beta-barrel structure but not well optimized in the wild-type protein.     

In-vitro selection of highly stabilized protein variants with optimized surface

Thermostable proteins are of prime importance in protein science, but it has remained difficult to develop general strategies for stabilizing a protein. Site-directed mutagenesis based on comparisons with thermophilic homologs is rarely successful because the sequence differences are too numerous and dominated by neutral mutations. Here we used a method of directed evolution to increase the stability of a mesophilic protein, the cold shock protein Bs-CspB from Bacillus subtilis. It differs from its thermophilic counterpart Bc-Csp from Bacillus caldolyticus at 12 surface-exposed positions. To elucidate the stabilizing potential of exposed amino acid residues, six of these variant positions were randomized by saturation mutagenesis, the corresponding library of sequences was inserted into the gene-3-protein of the filamentous phage fd, and stabilized variants were selected by the Proside technique. Proside links the increased protease resistance of stabilized protein variants with the infectivity of the phage. Many strongly stabilized variants of Bs-CspB were identified in two selections, one in the presence of a denaturant and the other at elevated temperature. Several of them are significantly more stable than the naturally thermostable homolog Bc-Csp, and the best variant reaches Tm-Csp (the homolog from the hyperthermophile Thermotoga maritima) in stability. Remarkably, this variant differs from Tm-Csp at five and from Bc-Csp at all six randomized positions. This indicates that proteins can be strongly stabilized by many different sets of surface mutations, and Proside selects them efficiently from large libraries. The course of the selection could be directed by the conditions. In an ionic denaturant non-polar surface interactions were optimized, whereas at elevated temperature variants with improved electrostatics were selected, pointing to two different strategies for stabilization at protein surfaces.     

Electrostatic contributions to the stability of a thermophilic cold shock protein

The thermophilic Bacillus caldolyticus cold shock protein (Bc-Csp) differs from the mesophilic Bacillus subtilis cold shock protein B (Bs-CspB) in 11 of the 66 residues. Stability measurements of Schmid and co-workers have implicated contributions of electrostatic interactions to the thermostability. To further elucidate the physical basis of the difference in stability, previously developed theoretical methods that treat electrostatic effects in both the folded and the unfolded states were used in this paper to study the effects of mutations, ionic strength, and temperature. For 27 mutations that narrow the difference in sequence between Bc-Csp and Bs-CspB, calculated changes in unfolding free energy (Delta G) and experimental results have a correlation coefficient of 0.98. Bc-Csp appears to use destabilization of the unfolded state by unfavorable charge-charge interactions as a mechanism for increasing stability. Accounting for the effects of ionic strength and temperature on the electrostatic free energies in both the folded and the unfolded states, explanations for two important experimental observations are presented. The disparate ionic strength dependences of Delta G for Bc-Csp and Bs-CspB were attributed to the difference in the total charges (-2e and -6e, respectively). A main contribution to the much higher unfolding entropy of Bs-CspB was found to come from the less favorable electrostatic interactions in the folded state. These results should provide insight for understanding the thermostability of other thermophilic proteins.     

Role of the chain termini for the folding transition state of the cold shock protein.

Residues Arg3 and Leu66 are crucially important for the enhanced stability of the cold shock protein Bc-Csp from the thermophile Bacillus caldolyticus relative to its homologue Bs-CspB from the mesophile Bacillus subtilis. Arg3, which replaces Glu3 of Bs-CspB, accounts for two-thirds of the stability difference and for the entire difference in Coulombic interactions between the two proteins. Leu66, which replaces Glu66 of Bs-CspB, contributes additional hydrophobic interactions. To elucidate the role of these two residues near the chain termini for the rapid folding of the cold shock proteins, we performed an extensive mutational analysis of the folding kinetics to characterize interactions between residues 3, 46, and 66 in the transition state of folding. We employed a pressure-jump apparatus which allows folding to be followed over a broad range of temperatures and urea concentrations in the time range of microseconds to minutes. The N-terminal region folds early, and the interactions that originate from residue 3 are present to a large extent in the transition state already. They include a hydrophobic contribution, a general electrostatic stabilization by the positive charge of Arg3 in Bc-Csp, and a pairwise Coulombic repulsion with Glu46 in the Arg3Glu variant. The C-terminus appears to be largely unfolded in the transition state. The interactions of Leu66, including those with the already structured N-terminal region, are established only after passage through the transition state. The N- and C-termini of the cold shock proteins thus contribute differently to the folding kinetics, although they are very close in space in the folded protein.     

Key role of coulombic interactions for the folding transition state of the cold shock protein

The cold shock protein CspB shows a five-stranded beta-sheet structure, and it folds rapidly via a native-like transition state. A previous Phi value analysis showed that most of the residues with Phi values close to one reside in strand beta1, and two of them, Lys5 and Lys7 are partially exposed charged residues. To elucidate how coulombic interactions of these two residues contribute to the energetic organisation of the folding transition state we performed comparative folding experiments in the presence of an ionic denaturant (guanidinium chloride) and a non-ionic denaturant (urea) and a double-mutant analysis. Lys5 contributes 6.6 kJ mol(-1) to the stability of the transition state, and half of it originates from screenable coulombic interactions. Lys7 contributes 5.3 kJ mol(-1), and 3.4 kJ mol(-1) of it are screened by salt. In the folded protein Lys7 interacts with Asp25, and the screenable coulombic interaction between these two residues is fully formed in the transition state. This suggests that long-range coulombic interactions such as those originating from Lys5 and Lys7 of CspB can be important for organizing and stabilizing native-like structure early in protein folding.     

In vitro evolution of a hyperstable Gbeta1 variant

An in-vitro selection strategy was used to obtain strongly stabilized variants of the beta1 domain of protein G (Gbeta1). In a two-step approach, first candidate positions with a high potential for stabilization were identified in Gbeta1 libraries that were created by error-prone PCR, and then, after randomization of these positions by saturation mutagenesis, strongly stabilized variants were selected. For both steps the in-vitro selection method Proside was employed. Proside links the stability of a protein with the infectivity of a filamentous phage. Ultimately, residues from the two best selected variants were combined in a single Gbeta1 molecule. This variant with the four mutations E15V, T16L, T18I, and N37L showed an increase of 35.1 degrees C in the transition midpoint and of 28.5 kJ mol(-1) (at 70 degrees C) in the Gibbs free energy of stabilization. It was considerably more stable than the best variant from a previous Proside selection, in which positions were randomized that had originally been identified by computational design. Only a single substitution (T18I) was found in both selections. The best variants from the present selection showed a higher cooperativity of thermal unfolding, as indicated by an increase in the enthalpy of unfolding by about 60 kJ mol(-1). This increase is apparently correlated with the presence of Leu residues that were selected at the positions 16 and 37.     

Thermal stability and atomic-resolution crystal structure of the Bacillus caldolyticus cold shock protein

The bacterial cold shock proteins are small compact beta-barrel proteins without disulfide bonds, cis-proline residues or tightly bound cofactors. Bc-Csp, the cold shock protein from the thermophile Bacillus caldolyticus shows a twofold increase in the free energy of stabilization relative to its homolog Bs-CspB from the mesophile Bacillus subtilis, although the two proteins differ by only 12 out of 67 amino acid residues. This pair of cold shock proteins thus represents a good system to study the atomic determinants of protein thermostability. Bs-CspB and Bc-Csp both unfold reversibly in cooperative transitions with T(M) values of 49.0 degrees C and 77.3 degrees C, respectively, at pH 7.0. Addition of 0.5 M salt stabilizes Bs-CspB but destabilizes Bc-Csp. To understand these differences at the structural level, the crystal structure of Bc-Csp was determined at 1.17 A resolution and refined to R=12.5% (R(free)=17.9%). The molecular structures of Bc-Csp and Bs-CspB are virtually identical in the central beta-sheet and in the binding region for nucleic acids. Significant differences are found in the distribution of surface charges including a sodium ion binding site present in Bc-Csp, which was not observed in the crystal structure of the Bs-CspB. Electrostatic interactions are overall favorable for Bc-Csp, but unfavorable for Bs-CspB. They provide the major source for the increased thermostability of Bc-Csp. This can be explained based on the atomic-resolution crystal structure of Bc-Csp. It identifies a number of potentially stabilizing ionic interactions including a cation-binding site and reveals significant changes in the electrostatic surface potential.     

残基突变对冷休克蛋白热稳定性影响的计算分析研究(英文)

残基突变是提高蛋白质热稳定性最直接有效的方式。在本文中,我们选取一对冷休克蛋白质作为研究对象,其中一个来自嗜温的Bacillus subtilis(Bs-CspB),另一个来自嗜热的Bacillus caldolyticus(Bc-Csp),这两个蛋白质在序列和结构上具有高度的相似性,但两者的耐热能力却相差很大。我们利用全原子模型计算残基突变前后蛋白质的自由能和氨基酸之间相互作用能的变化,分析残基突变对冷休克蛋白热稳定性的影响。通过对比两个蛋白质对应位置上残基的能量,我们成功鉴别出对Bc-Csp的高热稳定性有突出贡献的残基。我们计算了这些残基突变前后,该残基的静电相互作用和范德华相互作用的变化,以分析该残基对Bc-Csp高热稳定性的主要贡献。同时,我们分析了离子键对蛋白质热稳定性的贡献。我们的计算结果和实验结果吻合得很好,关键在于利用该方法可以详细地说明残基突变影响蛋白质热稳定性的根本原因。本文为研究残基突变对蛋白质热稳定性的影响提供了一种计算思路和方法,并有助于设计具有高耐热能力的蛋白质。     

Thermodynamics of a diffusional protein folding reaction

The folding reactions of several proteins are well described as diffusional barrier crossing processes, which suggests that they should be analyzed by Kramers' rate theory rather than by transition state theory. For the cold shock protein Bc-Csp from Bacillus caldolyticus, we measured stability and folding kinetics, as well as solvent viscosity as a function of temperature and denaturant concentration. Our analysis indicates that diffusional folding reactions can be treated by transition state theory, provided that the temperature and denaturant dependence of the solvent viscosity is properly accounted for, either at the level of the measured rate constants or of the calculated activation parameters. After viscosity correction the activation barriers for folding become less enthalpic and more entropic. The transition from an enthalpic to an entropic folding barrier with increasing temperature is, however, apparent in the data before and after this correction. It is a consequence of the negative activation heat capacity of refolding, which is independent of solvent viscosity. Bc-Csp and its mesophilic homolog Bs-CspB from Bacillus subtilis differ strongly in stability but show identical enthalpic and entropic barriers to refolding. The increased stability of Bc-Csp originates from additional enthalpic interactions that are established after passage through the activated state. As a consequence, the activation enthalpy of unfolding is increased relative to Bs-CspB.     

Role of surface hydrophobic residues in the conformational stability of human lysozyme at three different positions

To evaluate the contribution of the amino acid residues on the surface of a protein to its stability, a series of hydrophobic mutant human lysozymes (Val to Gly, Ala, Leu, Ile, Met, and Phe) modified at three different positions on the surface, which are located in the alpha-helix (Val 110), the beta-sheet (Val 2), and the loop (Val 74), were constructed. Their thermodynamic parameters of denaturation and crystal structures were examined by calorimetry and by X-ray crystallography at 100 K, respectively. Differences in the denaturation Gibbs energy change between the wild-type and the hydrophobic mutant proteins ranged from 4.6 to -9.6 kJ/mol, 2.7 to -1.5 kJ/mol, and 3.6 to -0.2 kJ/mol at positions 2, 74, and 110, respectively. The identical substitution at different positions and different substitutions at the same position resulted in different degrees of stabilization. Changes in the stability of the mutant proteins could be evaluated by a unique equation considering the conformational changes due to the substitutions [Funahashi et al. (1999) Protein Eng. 12, 841-850]. For this calculation, secondary structural propensities were newly considered. However, some mutant proteins were not adapted to the equation. The hydration structures around the mutation sites of the exceptional mutant proteins were affected due to the substitutions. The stability changes in the exceptional mutant proteins could be explained by the formation or destruction of the hydration structures. These results suggest that the hydration structure mediated via hydrogen bonds covering the protein surface plays an important role in the conformational stability of the protein.     

Electrostatic interactions in leucine zippers: thermodynamic analysis of the contributions of Glu and His residues and the effect of mutating salt bridges

Electrostatic interactions play a complex role in stabilizing proteins. Here, we present a rigorous thermodynamic analysis of the contribution of individual Glu and His residues to the relative pH-dependent stability of the designed disulfide-linked leucine zipper AB(SS). The contribution of an ionized side-chain to the pH-dependent stability is related to the shift of the pK(a) induced by folding of the coiled coil structure. pK(a)(F) values of ten Glu and two His side-chains in folded AB(SS) and the corresponding pK(a)(U) values in unfolded peptides with partial sequences of AB(SS) were determined by 1H NMR spectroscopy: of four Glu residues not involved in ion pairing, two are destabilizing (-5.6 kJ mol(-1)) and two are interacting with the positive alpha-helix dipoles and are thus stabilizing (+3.8 kJ mol(-1)) in charged form. The two His residues positioned in the C-terminal moiety of AB(SS) interact with the negative alpha-helix dipoles resulting in net stabilization of the coiled coil conformation carrying charged His (-2.6 kJ mol(-1)). Of the six Glu residues involved in inter-helical salt bridges, three are destabilizing and three are stabilizing in charged form, the net contribution of salt-bridged Glu side-chains being destabilizing (-1.1 kJ mol(-1)). The sum of the individual contributions of protonated Glu and His to the higher stability of AB(SS) at acidic pH (-5.4 kJ mol(-1)) agrees with the difference in stability determined by thermal unfolding at pH 8 and pH 2 (-5.3 kJ mol(-1)). To confirm salt bridge formation, the positive charge of the basic partner residue of one stabilizing and one destabilizing Glu was removed by isosteric mutations (Lys-->norleucine, Arg-->norvaline). Both mutations destabilize the coiled coil conformation at neutral pH and increase the pK(a) of the formerly ion-paired Glu side-chain, verifying the formation of a salt bridge even in the case where a charged side-chain is destabilizing. Because removing charges by a double mutation cycle mainly discloses the immediate charge-charge effect, mutational analysis tends to overestimate the overall energetic contribution of salt bridges to protein stability.     

Contribution of surface salt bridges to protein stability: guidelines for protein engineering

The small globular protein, ubiquitin, contains a pair of oppositely charged residues, K11 and E34, that according to the three-dimensional structure are located on the surface of this protein with a spatial orientation characteristic of a salt bridge. We investigated the strength of this salt bridge and its contribution to the global stability of the ubiquitin molecule. Using the "double mutant cycle" analysis, the strength of the pairwise interactions between K11 and E34 was estimated to be favorable by 3.6kJ/mol. Further, the salt bridge of the reverse orientation, i.e. E11/K34, can be formed and is found to have a strength (3.8kJ/mol) similar to that of the K11/E34 pair. However, the global stability of the K11/E34 variant of ubiquitin is 2.2kJ/mol higher than that of the E11/K34 variant. The difference in the contribution of the opposing salt bridge orientations to the overall stability of the ubiquitin molecule is attributed to the difference in the charge-charge interactions between residues forming the salt bridge and the rest of the ionizable groups in this protein. On the basis of these results, we concluded that surface salt bridges are stabilizing, but their contribution to the overall protein stability is strongly context-dependent, with charge-charge interactions being the largest determinant. Analysis of 16 salt bridges from six different proteins, for which detailed experimental data on energetics have been reported, support the conclusions made from the analysis of the salt bridge in ubiquitin. Implications of these findings for engineering proteins with enhanced thermostability are discussed.     

Increasing protein stability by altering long-range coulombic interactions.

It is difficult to increase protein stability by adding hydrogen bonds or burying nonpolar surface. The results described here show that reversing the charge on a side chain on the surface of a protein is a useful way of increasing stability. Ribonuclease T1 is an acidic protein with a pI approximately 3.5 and a net charge of approximately -6 at pH 7. The side chain of Asp49 is hyperexposed, not hydrogen bonded, and 8 A from the nearest charged group. The stability of Asp49Ala is 0.5 kcal/mol greater than wild-type at pH 7 and 0.4 kcal/mol less at pH 2.5. The stability of Asp49His is 1.1 kcal/mol greater than wild-type at pH 6, where the histidine 49 side chain (pKa = 7.2) is positively charged. Similar results were obtained with ribonuclease Sa where Asp25Lys is 0.9 kcal/mol and Glu74Lys is 1.1 kcal/mol more stable than the wild-type enzyme. These results suggest that protein stability can be increased by improving the coulombic interactions among charged groups on the protein surface. In addition, the stability of RNase T1 decreases as more hydrophobic aromatic residues are substituted for Ala49, indicating a reverse hydrophobic effect.     

Contribution of proton linkage to the thermodynamic stability of the major cold-shock protein of Escherichia coli CspA

The stability of protein is defined not only by the hydrogen bonding, hydrophobic effect, van der Waals interactions, and salt bridges. Additional, much more subtle contributions to protein stability can arise from surface residues that change their properties upon unfolding. The recombinant major cold shock protein of Escherichia coli CspA an all-beta protein unfolds reversible in a two-state manner, and behaves in all other respects as typical globular protein. However, the enthalpy of CspA unfolding strongly depends on the pH and buffer composition. Detailed analysis of the unfolding enthalpies as a function of pH and buffers with different heats of ionization shows that CspA unfolding in the pH range 5.5-9.0 is linked to protonation of an amino group. This amino group appears to be the N-terminal alpha-amino group of the CspA molecule. It undergoes a 1.6 U shift in pKa values between native and unfolded states. Although this shift in pKa is expected to contribute approximately 5 kJ/mol to CspA stabilization energy the experimentally observed stabilization is only approximately 1 kJ/mol. This discrepancy is related to a strong enthalpy-entropy compensation due, most likely, to the differences in hydration of the protonated and deprotonated forms of the alpha-amino group.     

Thermodynamic analysis of the unfolding and stability of the dimeric DNA-binding protein HU from the hyperthermophilic eubacterium Thermotoga maritima and its E34D mutant.

We have studied the stability of the histone-like, DNA-binding protein HU from the hyperthermophilic eubacterium Thermotoga maritima and its E34D mutant by differential scanning microcalorimetry and CD under acidic conditions at various concentrations within the range of 2-225 micro m of monomer. The thermal unfolding of both proteins is highly reversible and clearly follows a two-state dissociation/unfolding model from the folded, dimeric state to the unfolded, monomeric one. The unfolding enthalpy is very low even when taking into account that the two disordered DNA-binding arms probably do not contribute to the cooperative unfolding, whereas the quite small value for the unfolding heat capacity change (3.7 kJ.K(-1).mol(-1)) stabilizes the protein within a broad temperature range, as shown by the stability curves (Gibbs energy functions vs. temperature), even though the Gibbs energy of unfolding is not very high either. The protein is stable at pH 4.00 and 3.75, but becomes considerably less so at pH 3.50 and below, to the point that a simple decrease in concentration will lead to unfolding of both the wild-type and the mutant protein at pH 3.50 and low temperatures. This indicates that various acid residues lose their charges leaving uncompensated positively charged clusters. The wild-type protein is more stable than its E34D mutant, particularly at pH 4.00 and 3.75 although less so at 3.50 (1.8, 1.6 and 0.6 kJ.mol(-1) at 25 degrees C for DeltaDeltaG at pH 4.00, 3.75 and 3.50, respectively), which seems to be related to the effect of a salt bridge between E34 and K13.     

Two exposed amino acid residues confer thermostability on a cold shock protein

Thermophilic organisms produce proteins of exceptional stability. To understand protein thermostability at the molecular level we studied a pair of cold shock proteins, one of mesophilic and one of thermophilic origin, by systematic mutagenesis. Although the two proteins differ in sequence at 12 positions, two surface-exposed residues are responsible for the increase in stability of the thermophilic protein (by 15.8 kJ mol-1 at 70 degrees C). 11.5 kJ mol-1 originate from a predominantly electrostatic contribution of Arg 3 and 5.2 kJ mol-1 from hydrophobic interactions of Leu 66 at the carboxy terminus. The mesophilic protein could be converted to a highly thermostable form by changing the Glu residues at positions 3 and 66 to Arg and Leu, respectively. The variation of surface residues may thus provide a simple and powerful approach for increasing the thermostability of a protein.     

Solution structure and thermal stability of ribosomal protein L30e from hyperthermophilic archaeon Thermococcus celer

To understand the structural basis of thermostability, we have determined the solution structure of a thermophilic ribosomal protein L30e from Thermococcus celer by NMR spectroscopy. The conformational stability of T. celer L30e was measured by guanidine and thermal-induced denaturation, and compared with that obtained for yeast L30e, a mesophilic homolog. The melting temperature of T. celer L30e was 94 degrees C, whereas the yeast protein denatured irreversibly at temperatures >45 degrees C. The two homologous proteins also differ greatly in their stability at 25 degrees C: the free energy of unfolding was 45 kJ/mole for T. celer L30e and 14 kJ/mole for the yeast homolog. The solution structure of T. celer L30e was compared with that of the yeast homolog. Although the two homologous proteins do not differ significantly in their number of hydrogen bonds and the amount of solvent accessible surface area buried with folding, the thermophilic T. celer L30e was found to have more long-range ion pairs, more proline residues in loops, and better helix capping residues in helix-1 and helix-4. A K9A variant of T. celer L30e was created by site-directed mutagenesis to examine the role of electrostatic interactions on protein stability. Although the melting temperatures of the K9A variant is approximately 8 degrees C lower than that of the wild-type L30e, their difference in T(m) is narrowed to approximately 4.2 degrees C at 0.5 M NaCl. This salt-dependency of melting temperatures strongly suggests that electrostatic interactions contribute to the thermostability of T. celer L30e.     

Thermodynamic and kinetic destabilization of triosephosphate isomerase resulting from the mutation of conserved and non-conserved cysteines

Several variants of Saccharomyces cerevisiae triosephosphate isomerase (yTIM) were studied to determine how mutations of conserved and non-conserved Cys residues affect the enzyme. Wild-type yTIM has two buried free cysteines: Cys 41 (non-conserved) and the invariant Cys 126. Single-site mutants, containing substitutions of these cysteines with Ala, Val, or Ser (the three most conservative changes for a buried Cys, according to substitution matrices), were examined for stability and enzymatic activity. Neither of the Cys residues was found to be essential for enzyme catalysis. Determination of the global stability of the mutants indicated that, regardless of which Cys was substituted, individual Cys→Ala and Cys→Val mutations, as well as the C41S substitution, all decrease the unfolding free energy of the dimeric protein by less than 23 kJ mol(-1) (at 37 °C, pH 7.4), as compared to the wild-type enzyme. In contrast, a substantially larger destabilization (37 kJ mol(-1)) was found in the C126S mutant. These results suggest that, with the exception of C126S, all of these mutations can be regarded as neutral (i.e., mutations that do not impair the reproductive success of the organism). Accordingly, Cys 126 has remained invariant across evolution because its neutral substitutions by Ala or Val would require a highly unlikely, concerted double mutation at any of the Cys codons. Furthermore, detrimental effects to a cell expressing the C126S TIM mutant more likely arise from the high unfolding rate of this enzyme.