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.     

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.     

The effects of ionic strength on protein stability: the cold shock protein family

Continuum electrostatic models are used to examine in detail the mechanism of protein stabilization and destabilization due to salt near physiological concentrations. Three wild-type cold shock proteins taken from mesophilic, thermophilic, and hyperthermophilic bacteria are studied using these methods. The model is validated by comparison with experimental data collected for these proteins. In addition, a number of single point mutants and three designed sequences are examined. The results from this study demonstrate that the sensitivity of protein stability toward salt is correlated with thermostability in the cold shock protein family. The calculations indicate that the mesophile is stabilized by the presence of salt while the thermophile and hyperthermophile are destabilized. A decomposition of the salt influence at a residue level permits identification of regions of the protein sequences that contribute toward the observed salt-dependent stability. This model is used to rationalize the effect of various point mutations with regard to sensitivity toward salt. Finally, it is demonstrated that designed cold shock protein variants exhibit electrostatic properties similar to the natural thermophilic and hyperthermophilic proteins.     

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.     

Role of the charge-charge interactions in defining stability and halophilicity of the CspB proteins

Charge-charge interactions on the surface of native proteins are important for protein stability and can be computationally redesigned in a rational way to modulate protein stability. Such computational effort led to an engineered protein, CspB-TB that has the same core as the mesophilic cold shock protein CspB-Bs from Bacillus subtilis, but optimized distribution of charge-charge interactions on the surface. The CspB-TB protein shows an increase in the transition temperature by 20 degrees C relative to the unfolding temperature of CspB-Bs. The CspB-TB and CspB-Bs protein pair offers a unique opportunity to further explore the energetics of charge-charge interactions as the substitutions at the same sequence positions are done in largely similar structural but different electrostatic environments. In particular we addressed two questions. What is the contribution of charge-charge interactions in the unfolded state to the protein stability and how amino acid substitutions modulate the effect of increase in ionic strength on protein stability (i.e. protein halophilicity). To this end, we experimentally measured the stabilities of over 100 variants of CspB-TB and CspB-Bs proteins with substitutions at charged residues. We also performed computational modeling of these protein variants. Analysis of the experimental and computational data allowed us to conclude that the charge-charge interactions in the unfolded state of two model proteins CspB-Bs and CspB-TB are not very significant and computational models that are based only on the native state structure can adequately, i.e. qualitatively (stabilizing versus destabilizing) and semi-quantitatively (relative rank order), predict the effects of surface charge neutralization or reversal on protein stability. We also show that the effect of ionic strength on protein stability (protein halophilicity) appears to be mainly due to the screening of the long-range charge-charge interactions.     

Complementation analysis of the cold-sensitive phenotype of the Escherichia coli csdA deletion strain

The cold shock response of Escherichia coli is elicited by downshift of temperature from 37 degrees C to 15 degrees C and is characterized by induction of several cold shock proteins, including CsdA, during the acclimation phase. CsdA, a DEAD-box protein, has been proposed to participate in a variety of processes, such as ribosome biogenesis, mRNA decay, translation initiation, and gene regulation. It is not clear which of the functions of CsdA play a role in its essential cold shock function or whether all do, and so far no protein has been shown to complement its function in vivo. Our screening of an E. coli genomic library for an in vivo counterpart of CsdA that can compensate for its absence at low temperature revealed only one protein, RhlE, another DEAD-box RNA helicase. We also observed that although not detected in our genetic screening, two cold shock-inducible proteins, namely, CspA, an RNA chaperone, and RNase R, an exonuclease, can also complement the cold shock function of CsdA. Interestingly, the absence of CsdA and RNase R leads to increased sensitivity of the cells to even moderate temperature downshifts. The correlation between the helicase activity of CsdA and the stability of mRNAs of cold-inducible genes was shown using cspA mRNA, which was significantly stabilized in the DeltacsdA cells, an effect counteracted by overexpression of wild-type CsdA or RNase R but not by that of the helicase-deficient mutant of CsdA. These results suggest that the primary role of CsdA in cold acclimation of cells is in mRNA decay and that its helicase activity is pivotal for promoting degradation of mRNAs stabilized at low temperature.     

The thermostability of DNA-binding protein HU from mesophilic, thermophilic, and extreme thermophilic bacteria

Based on primary structure comparison between four highly homologous DNA-binding proteins (HUs) displaying differential thermostability, we have employed in vitro site-directed mutagenesis to decipher their thermostability mechanism at the molecular level. The contribution of the 11 amino acids that differ between the thermophilic HUBst from Bacillus stearothermophilus (Tm = 61.6 degrees C) and the mesophilic HUBsu from Bacillus subtilis (Tm = 39.7 degrees C) was evaluated by replacing these amino acids in HUBst with their mesophilic counterparts. Among 11 amino acids, three residues, Gly-15, Glu-34, and Val-42, which are highly conserved in the thermophilic HUs, have been found to be responsible for the thermostability of HUBst. These amino acids in combination (HUBst-G15E/E34D/V42I) reduce the thermostability of the protein (Tm = 45.1 degrees C) at the level of its mesophilic homologue HUBsu. By replacing these amino acids in HUBsu with their thermophilic counterparts, the HUBsu-E15G/D34E/142V mutant was generated with thermostability (Tm = 57.8 degrees C) at the level of thermophilic HUBst. Employing the same strategy, we generated several mutants in the extremely thermophilic HUTmar from Thermotoga maritima (Tm = 80.5 degrees C), and obtained data consistent with the previous results. The triplet mutant HUTmar-G15E/E34D/V421 (Tm = 35.9 degrees C) converted the extremely thermophilic protein HUTmar to mesophilic. The various forms of HU proteins were overproduced in Escherichia coli, highly purified, and the thermostability of the mutants confirmed by circular dichroism spectroscopy. The results presented here were elucidated on the basis of the X-ray structure of HUBst and HUTmar (our unpublished results), and their mechanism was proposed at the molecular level. The results clearly show that three individual local interactions located at the helix-turn-helix part of the protein are responsible for the stability of HU proteins by acting cooperatively in a common mechanism for thermostability.     

Stabilization of Escherichia coli ribonuclease HI by strategic replacement of amino acid residues with those from the thermophilic counterpart.

Thermus thermophilus ribonuclease H is exceptionally stable against thermal and guanidine hydrochloride denaturations as compared to Escherichia coli ribonuclease HI (Kanaya, S., and Itaya, M. (1992) J. Biol. Chem. 267, 10184-10192). The identity in the amino acid sequences of these enzymes is 52%. As an initial step to elucidate the stabilization mechanism of the thermophilic RNase H, we examined whether certain regions in its amino acid sequence confer the thermostability. A variety of mutant proteins of E. coli RNase HI were constructed and analyzed for protein stability. In these mutant proteins, amino acid sequences in loops or terminal regions were systematically replaced with the corresponding sequences from T. thermophilus RNase H. Of the nine regions examined, replacement of the amino acid sequence in each of four regions (R4-R7) resulted in an increase in protein stability. Simultaneous replacements of these amino acid sequences revealed that the effect of each replacement on protein stability is independent of each other and cumulative. Replacement of all four regions (R4-R7) gave the most stable mutant protein. The temperature of the midpoint of the transition in the thermal unfolding curve and the free energy change of unfolding in the absence of denaturant of this mutant protein were increased by 16.7 degrees C and 3.66 kcal/mol, respectively, as compared to those of E. coli RNase HI. These results suggest that individual local interactions contribute to the stability of thermophilic proteins in an independent manner, rather than in a cooperative manner.     

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.     

Distance-dependent statistical potentials for discriminating thermophilic and mesophilic proteins

Identification of the characteristic structural patterns responsible for protein thermostability is theoretically important and practically useful but largely remains an open problem. These patterns may be revealed through comparative study on thermophilic and mesophilic proteins that have distinct thermostability. In this study, we constructed several distance-dependant potentials from thermophilic and mesophilic proteins. These potentials were then used to evaluate the structural difference between thermophilic and mesophilic proteins. We found that using the subtraction or division of the potentials derived from thermophilic and mesophilic proteins can dramatically increase the discriminatory ability. This approach revealed that the ability to distinct the subtle structural features responsible for protein thermostability may be effectively enhanced through rationally designed comparative study.     

蛋白质空间结构属性与全基因组微生物耐热性的关系

 溶剂接触表面积、空腔个数和体积、紧密度、疏水性以及温度因子是影响蛋白质耐热的主要三维结构参数.挑选NCBI COG数据库中具有全基因组的单细胞微生物,选择其中三维结构已知的蛋白质作为研究对象,分析这些因素对细菌类和古细菌类微生物耐热性的影响.结果表明：（1）古细菌类蛋白质的空腔个数和体积与耐热性无关,极性面积和表面残基个数随耐热性增加而降低;（2）超高温细菌类蛋白质的分子量比较小,空腔个数和体积都小于常温蛋白质,而且空腔个数对稳定性的贡献大于空腔体积;（3）无论是古细菌还是细菌类蛋白质,疏水性和紧密度都不随耐热性变化,但暴露残基个数越多,蛋白质的耐热性越差;（4）两类蛋白质侧链的温度因子都高于主链,这与侧链的运动性（柔性）一般比主链高的实验结果一致;另外,超高温细菌类蛋白质的温度因子明显高于常温蛋白质.     

Effective factors in thermostability of thermophilic proteins

Thermostability of proteins in general and especially thermophilic proteins has been subject of a wide variety of studies based on theoretical and experimental investigation. Thermostability seems to be a property obtained through many minor structural modifications rather than certain amino acids substitution. In comparison with its mesophile homologue in a thermostable protein, usually a number of amino acids are exchanged. A wide variety of theoretical studies are based on comparative investigation of thermophilic proteins characteristics with their mesophilic counterparts in order to reveal their sequences, structural differences and consequently, to relate these observed differences to the thermostability properties. In this work we have compared a dataset of thermophilic proteins with their mesophilic homologues and furthermore, a mesophilic proteins dataset was also compared with its mesophilic homologue. This strategy enabled us first, to eliminate noise or background differences from signals and moreover, the important factors which were related to the thermostability were recognized too. Our results reveal that thermophilic and mesophilic proteins have both similar polar and nonpolar contribution to the surface area and compactness. On the other hand, salt bridges and main chain hydrogen bonds show an increase in the majority of thermophilic proteins in comparison to their mesophilic homologues. In addition, in thermophilic proteins hydrophobic residues are significantly more frequent, while polar residues are less. These findings indicate that thermostable proteins through evolution adopt several different strategies to withstand high temperature environments.     

Computational design of the Fyn SH3 domain with increased stability through optimization of surface charge charge interactions

Computational design of surface charge-charge interactions has been demonstrated to be an effective way to increase both the thermostability and the stability of proteins. To test the robustness of this approach for proteins with predominantly beta-sheet secondary structure, the chicken isoform of the Fyn SH3 domain was used as a model system. Computational analysis of the optimal distribution of surface charges showed that the increase in favorable energy per substitution begins to level off at five substitutions; hence, the designed Fyn sequence contained four charge reversals at existing charged positions and one introduction of a new charge. Three additional variants were also constructed to explore stepwise contributions of these substitutions to Fyn stability. The thermodynamic stabilities of the variants were experimentally characterized using differential scanning calorimetry and far-UV circular dichroism spectroscopy and are in very good agreement with theoretical predictions from the model. The designed sequence was found to have increased the melting temperature, DeltaT (m) = 12.3 +/- 0.2 degrees C, and stability, DeltaDeltaG(25 degrees C) = 7.1 +/- 2.2 kJ/mol, relative to the wild-type protein. The experimental data suggest that a significant increase in stability can be achieved through a very small number of amino acid substitutions. Consistent with a number of recent studies, the presented results clearly argue for a seminal role of surface charge-charge interactions in determining protein stability and suggest that the optimization of surface interactions can be an attractive strategy to complement algorithms optimizing interactions in the protein core to further enhance protein stability.     

Confirmation of the hierarchical folding of RNase H: a protein engineering study.

The kinetic intermediate of RNase H is structured in a core region of the protein. To probe the role of this intermediate in the folding of RNase H, the folding kinetics of mutant proteins with altered native state stabilities were investigated. Mutations within the folding core destabilize the kinetic intermediate and slow refolding in a manner consistent with an obligatory intermediate model. Mutations outside of the folding core, however, do not affect the stability of the kinetic intermediate but do perturb the native state and transition state. These results indicate that interactions formed in the intermediate persist in the transition and native states and that RNase H folds through a hierarchical mechanism.     

Different packing of external residues can explain differences in the thermostability of proteins from thermophilic and mesophilic organisms

MOTIVATION: Understanding the basis of protein stability in thermophilic organisms raises a general question: what structural properties of proteins are responsible for the higher thermostability of proteins from thermophilic organisms compared to proteins from mesophilic organisms? RESULTS: A unique database of 373 structurally well-aligned protein pairs from thermophilic and mesophilic organisms is constructed. Comparison of proteins from thermophilic and mesophilic organisms has shown that the external, water-accessible residues of the first group are more closely packed than those of the second. Packing of interior parts of proteins (residues inaccessible to water molecules) is the same in both cases. The analysis of amino acid composition of external residues of proteins from thermophilic organisms revealed an increased fraction of such amino acids as Lys, Arg and Glu, and a decreased fraction of Ala, Asp, Asn, Gln, Thr, Ser and His. Our theoretical investigation of folding/unfolding behavior confirms the experimental observations that the interactions that differ in thermophilic and mesophilic proteins form only after the passing of the transition state during folding. Thus, different packing of external residues can explain differences in thermostability of proteins from thermophilic and mesophilic organisms. AVAILABILITY: The database of 373 structurally well-aligned protein pairs is available at http://phys.protres.ru/resources/termo_meso_base.html. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.     

Optimized and automated protocols for high-throughput screening of amylosucrase libraries

This article describes the design and validation of a general procedure for the high-throughput isolation of amylosucrase variants displaying higher thermostability or increased resistance to organic solvents. This procedure consists of 2 successive steps: an in vivo selection that eliminates inactive variants followed by automated screening of active variants to isolate mutants displaying enhanced features. The authors chose an Escherichia coli expression vector, allowing a high production rate of the recombinant enzyme in miniaturized culture conditions. The screening assay was validated by minimizing variability for various parameters of the protocol, especially bacterial growth and protein production in cultures in 96-well microplates. Recombinant amylosucrase production was normalized by decreasing the coefficient of variance from 27% to 12.5%. Selective screening conditions were defined to select variants displaying higher thermostability or increased resistance to organic solvents. A first-generation amylosucrase variant library, constructed by random mutagenesis, was subjected to this procedure, yielding a mutant displaying a 25-fold increased stability at 50 degrees C compared to the parental wild-type enzyme.     

Solution structure and backbone dynamics of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin mutant: a molecular analysis of its reduced thermal stability

No general strategy for thermostability has been yet established, because the extra stability of thermophiles appears to be the sum of different cumulative stabilizing interactions. In addition, the increase of conformational rigidity observed in many thermophilic proteins, which in some cases disappears when mesophilic and thermophilic proteins are compared at their respective physiological temperatures, suggests that evolutionary adaptation tends to maintain corresponding states with respect to conformational flexibility. In this study, we accomplished a structural analysis of the K18G/R82E Alicyclobacillus acidocaldarius thioredoxin (BacTrx) mutant, which has reduced heat resistance with respect to the thermostable wild-type. Furthermore, we have also achieved a detailed study, carried out at 25, 45, and 65 degrees C, of the backbone dynamics of both the BacTrx and its K18G/R82E mutant. Our findings clearly indicate that the insertion of the two mutations causes a loss of energetically favorable long-range interactions and renders the secondary structure elements of the double mutants more similar to those of the mesophilic Escherichia coli thioredoxin. Moreover, protein dynamics analysis shows that at room temperature the BacTrx, as well as the double mutant, are globally as rigid as the mesophilic thioredoxins; differently, at 65 degrees C, which is in the optimal growth temperature range of A. acidocaldarius, the wild-type retains its rigidity while the double mutant is characterized by a large increase of the amplitude of the internal motions. Finally, our research interestingly shows that fast motions on the pico- to nanosecond time scale are not detrimental to protein stability and provide an entropic stabilization of the native state. This study further confirms that protein thermostability is reached through diverse stabilizing interactions, which have the key role to maintain the structural folding stable and functional at the working temperature.     

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.     

Directed evolution of the thermostable xylanase from Thermomyces lanuginosus

The thermostability of the endo-beta-1,4-xylanase from Thermomyces lanuginosus (xynA) was improved by directed evolution using error-prone PCR. Transformants expressing the variant xylanases were first selected on 0.4% Remazol Brilliant Blue-xylan and then exposed to 80 degrees C. Whereas the wild type XynA lost 90% activity after 10 min at 80 degrees C, five mutants displayed both higher stabilities and activities than XynA. Four mutants were subjected to further mutagenesis to improve the stability and activity of the xylanase. Subsequent screening revealed three mutants with enhanced thermostability. Mutant 2B7-10 retained 71% of its activity after treatment at 80 degrees C for 60 min and had a half-life of 215 min at 70 degrees C, which is higher than that attained by XynA. Sequence analysis of second generation mutants revealed that mutations were not concentrated in any particular region of the protein and exhibited much variation. The best mutant obtained from this study was variant 2B7-10, which had a single substitution (Y58F) in beta-sheet A of the protein, which is the hydrophilic, solvent-accessible outer surface of the enzyme. Most of the mutants obtained in this study displayed a compromise between stability and activity, the only exception being mutant 2B7-10. This variant showed increased activity and thermostability.