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

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

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.     

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.     

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.     

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.     

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.     

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.     

Impacts of the charged residues mutation S48E/N62H on the thermostability and unfolding behavior of cold shock protein: insights from molecular dynamics simulation with Gō model

The cold shock protein from the hyperthermophile Thermotoga maritima (Tm-Csp) exhibits significantly higher thermostability than its homologue from the thermophile Bacillus caldolyticus (Bc-Csp). Experimental studies have shown that the electrostatic interactions unique to Tm-Csp are responsible for improving its thermostability. In the present work, the favorable charged residues in Tm-Csp were grafted into Bc-Csp by a double point mutation of S48E/N62H, and the impacts of the mutation on the thermostability and unfolding/folding behavior of Bc-Csp were then investigated by using a modified Gō model, in which the electrostatic interactions between charged residues were considered in the model. Our simulation results show that this Tm-Csp-like charged residue mutation can effectively improve the thermostability of Bc-Csp without changing its two-state folding mechanism. Besides that, we also studied the unfolding kinetics and unfolding/folding pathway of the wild-type Bc-Csp and its mutant. It is found that this charged residue mutation obviously enhanced the stability of the C-terminal region of Bc-Csp, which decreases the unfolding rate and changes the unfolding/folding pathway of the protein. Our studies indicate that the thermostability, unfolding kinetics and unfolding/folding pathway of Bc-Csp can be artificially changed by introducing Tm-Csp-like favorable electrostatic interactions into Bc-Csp.

Open image in new window

Graphical abstract Tertiary structure of wild-type cold shock protein from the thermophile Bacillus caldolyticus

     

Origins of the high stability of an in vitro-selected cold-shock protein

In previous work, we had identified stabilized forms of the cold-shock protein Bs-CspB from Bacillus subtilis in a combinatorial library by an in vitro selection procedure. In this library, the sequence positions 2, 3, 46, 64, 66, and 67 had been randomized, because Bs-CspB differs from the naturally thermostable homolog Bc-Csp from Bacillus caldolyticus, among others, at these six positions. For the most stable selected variant, the midpoint of thermal unfolding (tM) increased by 28.2 deg. C and the Gibbs free energy of unfolding (deltaG(D)) by 19 kJ/mol. Here, we analyzed by site-directed mutagenesis how the selected residues contribute individually to this strong stabilization. Val3 and Val66, which replace Glu3 and Glu66 of wild-type Bs-CspB, each contribute about 7 kJ/mol to stability, the Thr64Arg substitution contributes 4.5 kJ/mol, and 3.2 kJ/mol originate from the Ala46Leu replacement. Gly67 at the carboxy terminus is unimportant for stability, the Arg selected at position 2 is overall slightly destabilizing but improves the coulombic interactions. The best variant differs from Bc-Csp at all six positions; nevertheless, natural and in vitro selection followed similar principles. In both cases, negatively charged residues at the adjacent positions 3 and 66 are avoided, and a positively charged residue is introduced into this area of the protein surface. Its exact location is unimportant. It can be at position 3, as in the thermophilic Bc-Csp, or at positions 2 or 64, as in the most stable selected variant. These positively charged residues contribute to stability not by engaging in pairwise coulombic interactions with a specific carboxyl group, but by generally improving the charge distribution in this particular region of the protein surface. These coulombic effects contribute significantly to the thermostability of the cold-shock proteins. They are only weakly interdependent and best explained by the presence of a flexible ion network at the protein surface. Our results emphasize that surface positions are very good candidates for optimizing protein stability.     

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.     

Thermophilicity of wild type and mutant cold shock proteins by molecular dynamics simulation

In the attempt to clarify possible mechanisms underlying thermal stability of proteins, we study through molecular dynamics thermophile Bc-Csp, mesophile Bs-CspB, and selected mutants. These proteins have been extensively characterized experimentally; researchers showed that differential thermostability among the wild type proteins is fundamentally linked to one or two mutated amino acids, and that the nature of the effect is electrostatic. They also inferred an atomistic mechanism related to removal of unfavorable interactions, rather than to the formation of salt bridges. Molecular dynamics allows us to confirm and support both hypotheses. Several other collective parameters have also been monitored in relation to thermophilicity, such as global and local rigidity, permanence and number of hydrogen bonds, or of salt links. None of these clearly correlates with the thermal stability of the presently studied proteins.     

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.     

The primary structure of Bacillus cereus neutral proteinase and comparison with thermolysin and Bacillus subtilis neutral proteinase

The complete amino-acid sequence of a neutral proteinase, produced by Bacillus cereus, was determined by protein sequencing. The neutral proteinase consists of 317 amino-acid residues. The primary structure is 70% homologous to thermolysin, a thermostable neutral proteinase and 45% homologous to Bacillus subtilis neutral proteinase. The zinc-binding site and the hydrophobic pocket of the active site are highly similar in all three proteinases. B. cereus neutral proteinase which is 20 degrees C less thermostable (60 degrees C) than thermolysin (80 degrees C) shows only minor differences in calcium binding sites and salt bridges compared to thermolysin (known from its X-ray diffraction analysis), whereas B. subtilis neutral proteinase (50 degrees C) differs considerably. Therefore it was assumed that the difference in thermostability between B. cereus neutral proteinase and thermolysin is not caused by different metal binding properties, or differences in the active site, but by changes within the rest of the molecule. Calculation of secondary structure potentials according to Chou & Fasman, hydrophobicity and bulkiness of the different structural elements and preferred cold----hot amino-acid residue exchanges indicated, that the thermostability of thermolysin compared to B. cereus neutral proteinase is caused by small effects contributed by numerous amino-acid exchanges distributed over the whole molecule, resulting in increased hydrophobicity of beta-pleated sheet and higher bulkiness of alpha-helical regions.     

Structural basis for the remarkable stability of Bacillus subtilis lipase (Lip A) at low pH

Understanding the structural basis of altered properties of proteins due to changes in temperature or pH provides useful insights in designing proteins with improved stability. Here we report the basis for the pH-dependent thermostability of the Bacillus subtilis lipase (Lip A) using spectroscopic and X-ray crystallographic studies. At pH values above 7, lipase denatures and aggregates when heated at temperatures above 45 degrees C. However, at pH below 6 lipase denatures upon heating but the activity and its native structure is completely recovered upon cooling. In order to obtain the structural basis of this unusual stability of lipase, we determined high-resolution crystal structures of the lipase in two different crystal forms at pH 4.5 and 5. These structures show linear oligomerization of lipase using only two types of dimeric associations and these inter-molecular interactions are completely absent in several crystal forms of wild-type and mutant proteins obtained at basic pH. In accordance with the crystallographic studies, spectroscopic investigations reveal an invariant secondary structure in the pH range of 4-10. Quaternary organization of lipase at low pH resulted in changes in the tryptophan environment and binding of 1-anilino-8-naphthalene sulfate (ANS) at low pH. Low pH stability of the lipase is not observed in the presence of sodium chloride (>0.2 M) indicating the importance of ionic interactions at low pH. Inter- and intra-molecular ionic interactions that occur at pH below 6.0 are proposed to trap the molecule in a conformation that allows its complete refolding upon cooling.     

The unique pentagonal structure of an archaeal Rubisco is essential for its high thermostability

We have previously determined the crystal structure of a novel pentagonal ribulose 1,5-bisphosphate carboxylase/oxygenase (Rubisco) from the hyperthermophilic archaeon, Thermococcus kodakaraensis KOD1. Here we have carried out biochemical studies to identify the necessities and/or advantages of this intriguing pentagonal structure. The structure indicated the presence of three neighboring residues (Glu-63, Arg-66, and Asp-69), participating in ionic interactions within unique dimer-dimer interfaces. We constructed three single mutant proteins (E63S, R66S, and D69S) and one triple mutant protein (E63S/R66S/D69S) by replacing the charged residues with serine. The wild type (WT) and all mutant proteins were purified and subjected to gel permeation chromatography at various temperatures. WT and D69S proteins were decameric at all temperatures examined between 30 and 90 degrees C. The majority of E63S and R66S were decamers at 30 degrees C but were found to gradually disassemble with the elevation in temperature. E63S/R66S/D69S was found in a dimeric form even at 30 degrees C. An interesting correlation was found between the subunit assembly and thermostability of the proteins. Circular dichroism and differential scanning calorimetry analyses indicated that the denaturation temperatures of dimeric enzymes (E63S, R66S, and E63S/R66S/D69S) were approximately 95 degrees C, whereas those of the enzymes retaining a decameric structure (WT and D69S) were approximately 110 degrees C. Disassembly into tetramer or dimer units did not alter the slopes of the Arrhenius plots, indicating that the decameric structure had no effect on catalytic performance per se. The results indicate that the decameric assembly of Tk-Rubisco contributes to enhance the thermostability of the enzyme. Taking into account the growth temperature of strain KOD1 (65-100 degrees C), the decameric structure of Tk-Rubisco can be considered essential for the stable presence of the enzyme in the host cells. This study provides an interesting example in which the thermostability of a protein can be enhanced by formation of a unique quaternary structure not found in mesophilic enzymes.     

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.     

Directed evolution of Thermus maltogenic amylase toward enhanced thermal resistance

The thermostability of maltogenic amylase from Thermus sp. strain IM6501 (ThMA) was improved greatly by random mutagenesis using DNA shuffling. Four rounds of DNA shuffling and subsequent recombination of the mutations produced the highly thermostable mutant enzyme ThMA-DM, which had a total of seven individual mutations. The seven amino acid substitutions in ThMA-DM were identified as R26Q, S169N, I333V, M375T, A398V, Q411L, and P453L. The optimal reaction temperature of the recombinant enzyme was 75 degrees C, which was 15 degrees C higher than that of wild-type ThMA, and the melting temperature, as determined by differential scanning calorimetry, was increased by 10.9 degrees C. The half-life of ThMA-DM was 172 min at 80 degrees C, a temperature at which wild-type ThMA was completely inactivated in less than 1 min. Six mutations that were generated during the evolutionary process did not significantly affect the specific activity of the enzyme, while the M375T mutation decreased activity to 23% of the wild-type level. The molecular interactions of the seven mutant residues that contributed to the increased thermostability of the mutant enzyme with other adjacent residues were examined by comparing the modeled tertiary structure of ThMA-DM with those of wild-type ThMA and related enzymes. The A398V and Q411L substitutions appeared to stabilize the enzyme by enhancing the interdomain hydrophobic interactions. The R26Q and P453L substitutions led potentially to the formation of genuine hydrogen bonds. M375T, which was located near the active site of ThMA, probably caused a conformational or dynamic change that enhanced thermostability but reduced the specific activity of the enzyme.