Salt Bridges in the Hyperthermophilic Protein Ssh10b Are Resilient to Temperature Increases*
Meng Ge,Xia-Yu Xia,Xian-Ming Pan
DOI: https://doi.org/10.1074/jbc.m805750200
2008-01-01
Abstract:A double mutant cycle (DMC) approach was employed to estimate the effect of temperature on the contribution of two highly conserved salt bridges to protein stability in the hyperthermophilic protein Ssh10b. The coupling free energy were 2.4 ± 0.4 kJ/mol at 298 K and 2.2 ± 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 ± 0.2 kJ/mol at 298 K and 5.9 ± 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridges are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures. The reason for the difference in coupling free energy between salt bridges Glu-54/Arg-57 and Glu-36/Lys-68 is discussed. Comparing our results with published DMC data for the contribution of salt bridges to stability in other proteins, we found that the energy contribution of a salt bridge formed by two charged residues far apart in the primary sequence is higher than that of those formed between two very close ones. Implications of this finding are useful for engineering proteins with enhanced thermostability. A double mutant cycle (DMC) approach was employed to estimate the effect of temperature on the contribution of two highly conserved salt bridges to protein stability in the hyperthermophilic protein Ssh10b. The coupling free energy were 2.4 ± 0.4 kJ/mol at 298 K and 2.2 ± 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 ± 0.2 kJ/mol at 298 K and 5.9 ± 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridges are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures. The reason for the difference in coupling free energy between salt bridges Glu-54/Arg-57 and Glu-36/Lys-68 is discussed. Comparing our results with published DMC data for the contribution of salt bridges to stability in other proteins, we found that the energy contribution of a salt bridge formed by two charged residues far apart in the primary sequence is higher than that of those formed between two very close ones. Implications of this finding are useful for engineering proteins with enhanced thermostability. Hyperthermophiles are a group of microorganisms with an optimum growth temperature between 80 and 100 °C. Remarkably, hyperthermophilic proteins and their mesophilic homologues typically show 40–85% sequence similarity, and their three-dimensional structures are highly superimposable (1Ladenstein R. Antranikian G. Adv. Biochem. Eng. Biotechnol. 1998; 61: 37-85PubMed Google Scholar, 2Woese C.R. Kandler O. Wheelis M.L. Proc. Natl. Acad. Sci. U. S. A. 1990; 87: 4576-4579Crossref PubMed Scopus (4535) Google Scholar, 3Vieille C. Zeikus G.J. Microbiol. Mol. Biol. Rev. 2001; 65: 1-43Crossref PubMed Scopus (1671) Google Scholar, 4Adams M.W. Annu. Rev. Microbiol. 1993; 47: 627-658Crossref PubMed Scopus (303) Google Scholar). These facts indicate that the molecular basis of their extreme thermal tolerance is associated with the delicate balance of non-covalent interactions. Considerable efforts have been invested in recent years to understand how proteins from hyperthermophiles can maintain stability at high temperatures. Several factors, such as improvement in packing density, strengthening of hydrophobic cores, and decreased length of surface loops, are thought to correlate with the increased thermostability of hyperthermophilic proteins (5Britton K.L. Baker P.J. Borges K.M. Engel P.C. Pasquo A. Rice D.W. Robb F.T. Scandurra R. Stillman T.J. Yip K.S. Eur. J. Biochem. 1995; 229: 688-695Crossref PubMed Scopus (122) Google Scholar, 6Russell R.J. Hough D.W. Danson M.J. Taylor G.L. Structure. 1994; 2: 1157-1167Abstract Full Text Full Text PDF PubMed Scopus (182) Google Scholar, 7Spassov V.Z. Karshikoff A.D. Ladenstein R. Protein Sci. 1995; 4: 1516-1527Crossref PubMed Scopus (153) Google Scholar, 8Schumann J. Bohm G. Schumacher G. Rudolph R. Jaenicke R. Protein Sci. 1993; 10: 1612-1620Crossref Scopus (59) Google Scholar, 9Vieille C. Zeikus J.G. Trends Biotechnol. 1996; 14: 183-191Abstract Full Text PDF Scopus (298) Google Scholar). The most commonly cited feature, however, is an increased number of salt bridges and salt bridge networks distributed over the protein surface (10Yip K.S. Stillman T.J. Britton K.L. Artymiuk P.J. Baker P.J. Sedelnikova S.E. Engel P.C. Pasquo A. Chiaraluce R. Consalvi V. Structure. 1995; 3: 1147-1158Abstract Full Text Full Text PDF PubMed Scopus (436) Google Scholar, 11Hennig M. Darimont B. Sterner R. Kirschner K. Jansonius J.N. Structure. 1995; 3: 1295-1306Abstract Full Text Full Text PDF PubMed Scopus (227) Google Scholar, 12Jaenicke R. Schurig H. Beaucamp N. Ostendorp R. Protein Chem. 1996; 48: 181-269Crossref PubMed Google Scholar, 13Korndorfer I. Steipe B. Huber R. Tomschy A. Jaenicke R. J. Mol. Biol. 1995; 246: 511-521Crossref PubMed Scopus (225) Google Scholar, 14Kawamura S. Tanaka I. Yamasaki N. Kimura M. J. Biochem. 1997; 121: 448-455Crossref PubMed Scopus (32) Google Scholar, 15Vetriani C. Maeder D.L. Tolliday N. Yip K.S. Stillman T.J. Britton K.L. Rice D.W. Klump H.H. Robb F.T. Proc. Natl. Acad. Sci. U. S. A. 1998; 95: 12300-12305Crossref PubMed Scopus (243) Google Scholar). A comparison of 13 structural parameters derived from the tertiary structures of 64 proteins from mesophiles and 29 proteins from thermophiles and hyperthermophiles revealed that the only feature shared in common is an increase in the number of salt bridges with increasing growth temperatures (16Szilágyi A. Závodszky P. Structure. 2000; 8: 493-504Abstract Full Text Full Text PDF PubMed Scopus (570) Google Scholar). However, experimental estimates of the free energy contribution of salt bridges to protein stability have led to conflicting conclusions, ranging from them having a stabilizing, insignificant, or even a destabilizing effect (17Marqusee S. Sauer R.T. Protein Sci. 1994; 3: 2217-2225Crossref PubMed Scopus (131) Google Scholar, 18Lounnas V. Wade R.C. Biochemistry. 1997; 36: 5402-5417Crossref PubMed Scopus (93) Google Scholar, 19Lebbink J.H. Consalvi V. Chiaraluce R. Berndt K.D. Ladenstein R. Biochemistry. 2002; 41: 15524-15535Crossref PubMed Scopus (26) Google Scholar, 20Barril X. Aleman C. Orozco M. Luque F.J. Proteins. 1998; 32: 67-79Crossref PubMed Scopus (81) Google Scholar, 21Marti D.N. Bosshard H.R. J. Mol. Biol. 2003; 330: 621-637Crossref PubMed Scopus (47) Google Scholar, 22Sun D.P. Sauer U. Nicholson H. Matthews B.W. Biochemistry. 1991; 30: 7142-7153Crossref PubMed Scopus (264) Google Scholar, 23Waldburger C.D. Schildbach J.F. Sauer R.T. Nat. Struct. Biol. 1995; 2: 122-128Crossref PubMed Scopus (339) Google Scholar). The association of two charged residues to form a salt bridge is thought to incur a substantial desolvation penalty that is seldom completely compensated for by favorable interactions within the salt bridge and the rest of the protein (24Hendsch Z.S. Tidor B. Protein Sci. 1994; 3: 211-226Crossref PubMed Scopus (590) Google Scholar, 25Elcock A.H. J. Mol. Biol. 1998; 284: 489-502Crossref PubMed Scopus (277) Google Scholar, 26Thomas A.S. Elcock A.H. J. Am Chem. Soc. 2004; 126: 2208-2214Crossref PubMed Scopus (77) Google Scholar). A continuum solvation model was used to investigate how this same argument applies at the high temperatures hyperthermophiles experience and predicted that at such high temperatures the desolvation penalty for the formation of a salt bridge is markedly reduced (25Elcock A.H. J. Mol. Biol. 1998; 284: 489-502Crossref PubMed Scopus (277) Google Scholar). Furthermore, molecular simulations suggest salt bridges are extremely resilient to temperature increases and thus are specially suited to promoting protein stability at high temperatures (26Thomas A.S. Elcock A.H. J. Am Chem. Soc. 2004; 126: 2208-2214Crossref PubMed Scopus (77) Google Scholar). However, attempts to obtain thermodynamic information at high temperature have been plagued by irreversibility of unfolding and/or inaccessibility of the unfolding transition to physical measurements. So far there are no reports of experiments to measure the stability contribution of salt bridges at high temperatures and to provide evidence for these theoretical predictions. The DNA-binding protein Ssh10b from the archaeon Sulfolobus shibatae is a member of the Sac10b family that is thought to be involved in chromosomal organization or DNA repair/recombination. Ssh10b is a highly thermostable dimeric protein composed of two identical subunits, each monomer consisting of 97 amino acid residues with no disulfide bonds (27Xue H. Guo R. Wen Y. Liu D. Huang L. J. Bacteriol. 2000; 182: 3929-3933Crossref PubMed Scopus (78) Google Scholar, 28Xu S. Qin S. Pan X.M. Biochem. J. 2004; 382: 433-440Crossref PubMed Scopus (23) Google Scholar). Ssh10b constitutes about 4–5% of total cellular protein, and binds dsDNA without apparent sequence specificity. Ssh10b is also capable of constraining negative DNA supercoils in a temperature-dependent fashion (29Cui Q. Tong Y. Xue H. Huang L. Feng Y. Wang J. J. Biol. Chem. 2003; 278: 51015-51022Abstract Full Text Full Text PDF PubMed Scopus (25) Google Scholar). The crystal structure of Ssh10b reveals that the monomer is a mixed α/β structure comprised of four β-strands and two α-helices (Fig. 1), each monomer containing four isolated salt bridges: Glu-36/Lys-68, Glu-54/Arg-57, Asp-63/Lys-97, Glu-66/Arg-95, and an ion-pair network, Lys-40/Glu-91/Arg-71/Glu-69 (30Wardleworth B.N. Russell R.J. Bell S.D. Taylor G.L. White M.F. EMBO J. 2002; 21: 4654-4662Crossref PubMed Scopus (132) Google Scholar), evidently more than the average of about 5 salt bridges per 150 amino acid residues (31Barlow D.J. Thornton J.M. J. Mol. Biol. 1983; 168: 867-885Crossref PubMed Scopus (627) Google Scholar). However, only the Glu-36/Lys-68 and Glu-54/Arg-57 salt bridges are highly conserved across the Sac10b family (Fig. 2). All of these salt bridges and ion-pair networks are surface-exposed, located on the back of the Ssh10b dimer and thought to play an important role in bracing and stabilizing the structure (30Wardleworth B.N. Russell R.J. Bell S.D. Taylor G.L. White M.F. EMBO J. 2002; 21: 4654-4662Crossref PubMed Scopus (132) Google Scholar).FIGURE 2Sequence alignment of archeal Ssh10b homologues.Ssh (S. shibatae), Sso (S. solfataricus), Sac (S. acidocaldarius), Afu (Archaeoglobus fulgidus), Hbu (Hyperthermus butylicus), Mja (Methanocaldococcus jannaschii), Mka (Methanopyrus kandleri), Mst (Methanosphaera stadtmanae), Mth (Methanosaeta thermophila), Pab (Pyrococcus abyssi), Pfu (Pyrococcus furiosus), Pho (Pyrococcus horikoshii), Pko (Pyrococcus kodakaraensis), Uncmar (Uncultured methanogenic archaeon RC-I), Tac (Thermoplasma acidophilum), Pae (Pyrobaculum aerophilum). The four residues involved in the two conserved salt bridges are marked with plus signs. The figure was rendered using the program DNAMAN.View Large Image Figure ViewerDownload Hi-res image Download (PPT) Electrostatic interactions of salt bridges can be studied by several methods. A particularly useful approach is to employ a method called double mutant cycle (DMC) 2The abbreviation used is: DMC, double mutant cycle. analysis to obtain quantitative information about the stability contribution of a salt bridge by mutating the residues involved separately and concurrently (32Serrano L. Horovitz A. Avron B. Bycroft M. Fersht A.R. Biochemistry. 1990; 29: 9343-9352Crossref PubMed Scopus (358) Google Scholar, 33Krylov D. Barchi J. Vinson C. J. Mol. Biol. 1998; 279: 959-972Crossref PubMed Scopus (107) Google Scholar, 34Vaughan C.K. Harryson P. Buckle A.M. Fersht A.R. Acta Crystallogr. Sect. D Biol. Crystallogr. 2002; 58: 591-600Crossref PubMed Scopus (35) Google Scholar, 35Luisi D.L. Snow C.D. Lin J.J. Hendsch Z.S. Tidor B. Raleigh D.P. Biochemistry. 2003; 42: 7050-7060Crossref PubMed Scopus (65) Google Scholar, 36Makhatadze G.I. Loladze V.V. Ermolenko D.N. Chen X. Thomas S.T. J. Mol. Biol. 2003; 327: 1135-1148Crossref PubMed Scopus (193) Google Scholar), thus cancelling all interaction effects except the direct interaction between the two charged residues. A DMC gives the so-called coupling free energy of a salt bridge associated with the direct interaction between the two charged residues under ideal conditions. In the present work, we analyzed the net strength of the two highly conserved salt bridges, Glu-36/Lys-68 and Glu-54/Arg-57, by DMC analysis. The coupling free energy were 2.4 ± 0.4 kJ/mol at 298 K and 2.2 ± 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 ± 0.2 kJ/mol at 298 K and 5.9 ± 0.6 kJ/mol at 353 K for Glu-36/Lys-68. The stability free energy of Ssh10b decrease greatly with increasing temperature, while the direct contribution of these two salt bridges to protein stability remain almost constant, providing evidence supporting the theoretical prediction that salt bridge interactions are extremely resilient to temperature increases and thus are specially suited to improving protein stability at high temperatures (26Thomas A.S. Elcock A.H. J. Am Chem. Soc. 2004; 126: 2208-2214Crossref PubMed Scopus (77) Google Scholar). Mutagenesis, Expression, and Purification of Ssh10b and Its Variants—All expression plasmids of Ssh10b variants were constructed from the parental plasmid pET11a-ssh10b. All mutations were introduced by site-directed mutagenesis through overlap extension PCR. The construct for each variant was verified by DNA sequencing. Expression plasmids of Ssh10b and its variants were transformed into the Escherichia coli BL21 (DE3) host strain. For protein expression, a single colony was grown in 100 ml of LB media containing 100 μg/ml ampicillin by shaking (∼220 rpm) at 37 °C overnight. Cultures were diluted 1:50 in fresh LB antibiotic-containing media and shaken for about 3 h at 37 °C. Target protein expression was then induced at A600 = 0.8–1.0 by the addition of isopropyl-1-thio-β-d-galactopyranoside to a final concentration of 0.3 mm, and growth was allowed to continue for about 5 h at 37 °C with constant shaking. After harvesting by centrifugation (4,000 rpm for 30 min), cell pellets from 1-liter cultures were re-suspended in 25 ml of buffer A (20 mm Tris-HCl, pH 7.5) and then disrupted by ultrasonication on ice. The lysate was maintained at 60 °C for 20 min to precipitate the E. coli proteins, and then centrifuged at 16,000 rpm for 30 min at 4 °C. The supernatant was applied to a 20-ml Source 30S column, which had been equilibrated with buffer A, and proteins were eluted with a 0–50% gradient (120 ml) of buffer B (1.5 m NaCl, 20 mm Tris-HCl, pH7.5). Fractions containing the target protein were identified by 15% SDS-PAGE and dialyzed overnight against buffer A. After centrifugation at 16,000 rpm for 30 min, the supernatant was applied to a 6-ml Resource-S column, which had been equilibrated with buffer A, and proteins were eluted with a 0–50% gradient (60 ml) of buffer B. Fractions containing the target protein were identified by 15% SDS-PAGE, dialyzed overnight against deionized water, and then lyophilized. Unfolding Studies—Unfolding of Ssh10b and its variants was studied by taking circular dichroism (CD) measurements with a π*-pistar 180 spectrometer (Applied Photophysics Ltd, UK), performed in buffer H (10 mm HEPES, pH 7.0), at a protein monomer concentration of about 25 μm. The CD signal was monitored using a rectangular quartz cuvette with a path length of 1 mm. For urea-induced unfolding, the urea solution was freshly prepared on the day of use. The samples containing various concentrations of urea were equilibrated at 25 °C overnight and then measured by far-UV CD at 222 nm. For heat-induced unfolding, each sample was heated from 40 to 98 °C and then cooled from 98 to 40 °C using stepwise changes of 2 °C, and the CD signal was recorded after equilibration for 2 min at each temperature point. All unfolding experiments were repeated 3–4 times. Analysis of the Denaturation Data—Previous studies in our laboratory have shown that both denaturant- and heat-induced unfolding of Ssh10b are fully reversible and follow a two-state mechanism involving a native dimer and two denatured monomers (28Xu S. Qin S. Pan X.M. Biochem. J. 2004; 382: 433-440Crossref PubMed Scopus (23) Google Scholar). Therefore, in this work the thermodynamic properties of Ssh10b and its variants were calculated assuming a two-state denaturation process in Reaction 1. N2Kobs↔2 U REACTION 1 The observed equilibrium constant (Kobs) and the corresponding free energy change (ΔG) at temperature T or denaturant concentration [D] were calculated according to Equations 1 and 2 (28Xu S. Qin S. Pan X.M. Biochem. J. 2004; 382: 433-440Crossref PubMed Scopus (23) Google Scholar), Kobs=2×PtyN+mNT[D]−yU+mUT[D]×[yN+mNT[D]−y]2y−yU+mUT[D](Eq. 1) ΔG=−RT ln Kobs (Eq. 1) where Pt is the total protein concentration in monomer units; R is the gas constant; T is the absolute temperature; y is the experimentally measured signal value at a given temperature (T) or given denaturant concentration ([D]); yN and yU are the intercepts; and mN and mU are the slopes of the native and unfolded baselines, respectively. According to the linear free energy model (37Schellman J.A. Annu. Rev. Biophys. Chem. 1987; 16: 115-137Crossref PubMed Scopus (299) Google Scholar, 38Chen B.L. Schellman J.A. Biochemistry. 1989; 28: 685-691Crossref PubMed Scopus (152) Google Scholar, 39Agashe V.R. Udgaonkar J.B. Biochemistry. 1995; 34: 3286-3299Crossref PubMed Scopus (209) Google Scholar, 40Kamal J.K. Nazeerunnisa M. Behere D.V. J. Biol. Chem. 2002; 277: 40717-40721Abstract Full Text Full Text PDF PubMed Scopus (25) Google Scholar), changes in free energy (ΔG) that occur on unfolding are expected to vary linearly with denaturant concentration in Equation 3 ([D]), ΔG=ΔG(H2O)−mG[D] (Eq. 3) where ΔG(H2O) represents the free energy change of unfolding in the absence of denaturant and mG is the slope of the transition for the free energy. For thermal unfolding, assuming that the heat capacity change (ΔCp) between the native and unfolded states of the system is relatively independent of temperature, gives us Equation 4 (41Santoro M.M. Bolen D.W. Biochemistry. 1988; 27: 8063-8068Crossref PubMed Scopus (1609) Google Scholar, 42Freire E. Methods Enzymol. 1995; 259: 144-168Crossref PubMed Scopus (98) Google Scholar), ΔG(T)=−RT ln(Kobs)=ΔHm−TΔSm+ΔCp((T−Tm)+T ln(1−T−TmT)) (Eq. 4) where ΔHm and ΔSm are the enthalpy and entropy changes, respectively, of the protein at the transition midpoint, where T = Tm. Within the transition range, where T × ln(1 - (T-Tm)/T) ∼ Tm - T, Equation 4 can be simplified to the van't Hoff plot in Equation 5. ΔG(T)=−RT ln(Kobs)=ΔHm−TΔSm(Eq. 5) The temperature of the transition midpoint (Tm) can be calculated according to Equation 6. Tm=ΔHmΔSm−R×ln(Pt)(Eq. 6) The uncertainties of the data were represented as two-sided 95% confidence intervals, which are given by Equation 7, b±tv,0.975×se(b)(Eq. 7) where b is the data (ΔG, ΔH, and Tm, etc.), se(b) is the standard error, tv,0.975 is the value of the t distribution with the number of degrees of freedom v for the two-sided 95% confidence interval. The direct stability contribution of each salt bridge was represented by using the average values of the coupling free energy of the four different DMCs and uncertainties estimated by Equation 7. Mutagenesis, Expression, and Purification of Ssh10b and Its Variants—Eight mutants for each salt bridge (Ala-36/Lys-68, Gln-36/Lys-68, Glu-36/Ala-68, Glu-36/Leu-68, Ala-36/Ala-68, Ala-36/Leu-68, Gln-36/Ala-68, and Gln-36/Leu-68 for salt bridge Glu-36/Lys-68; Ala-54/Arg-57, Gln-54/Arg-57, Glu-54/Ala-57, Glu-54/Leu-57, Ala-54/Ala-57, Ala-54/Leu-57, Gln-54/Ala-57, and Gln-54/Leu-57 for salt bridge Glu-54/Arg-57) were constructed by site-directed mutagenesis. Wild-type Ssh10b and its variants were overexpressed in E. coli BL21 (DE3) and purified to homogeneity. Protein purity was higher than 95% as confirmed by 15% SDS-PAGE. Unfolding Studies—Ssh10b is resistant to urea-induced denaturation in phosphate buffer (28Xu S. Qin S. Pan X.M. Biochem. J. 2004; 382: 433-440Crossref PubMed Scopus (23) Google Scholar), but is more susceptible to urea-induced denaturation in monovalent ion buffers (43Mao Y.J. Sheng X.R. Pan X.M. BMC Biochemistry. 2007; 8: 28Crossref PubMed Scopus (22) Google Scholar). And for any one of above proteins, when the urea-denatured protein solution was diluted to a lower concentration of denaturant, the CD spectrum of the renatured protein was identical to that of the native protein, indicating that the urea-induced unfolding processes of Ssh10b and its variants are fully reversible (data not shown). Fig. 3 shows representative urea-induced unfolding profiles for Ssh10b and its partial variants in buffer H. The linear free energy model was used to analyze the urea-induced unfolding profiles. Results are presented in Table 1 and reveal that the difference in stability between different variants is marked: relative to the Ssh10b wild-type, some variants show an increase in stability free energy at 298 K (ΔG(298)) of up to 7 kJ/mol, while others show a decrease in ΔG(298) of more than 10 kJ/mol.TABLE 1Parameters obtained from the analysis of unfolding experiments ΔG(298) was calculated from the results of urea-induced unfolding, and other parameters were calculated from the results of heat-induced unfolding. The uncertainties of the data were estimated according to Eq. 7.ProteinsΔG(298)ΔHmΔSmTmΔG(353)kJ mol–1kJ mol–1J mol–1 deg–1KkJ mol–1Wild type51.1 ± 0.7379.8 ± 15.9985.1 ± 45.1353.9 ± 0.232.1 ± 0.3A36/K6836.0 ± 0.4304.5 ± 10.2795.1 ± 29.5344.8 ± 0.222.8 ± 0.2Q36/K6839.7 ± 0.5348.1 ± 12.8910.1 ± 36.7348.7 ± 0.226.8 ± 0.3E36/A6846.8 ± 0.6384.0 ± 13.9996.4 ± 39.3354.1 ± 0.232.2 ± 0.3E36/L6858.1 ± 0.7332.0 ± 14.2848.1 ± 40.0354.6 ± 0.232.6 ± 0.3A36/A6843.6 ± 0.6318.3 ± 14.4806.5 ± 40.8355.8 ± 0.333.6 ± 0.3A36/L6854.9 ± 0.7402.8 ± 17.81042.1 ± 50.0356.4 ± 0.235.0 ± 0.4Q36/A6847.7 ± 0.5486.0 ± 22.91264.5 ± 64.0359.3 ± 0.239.6 ± 0.4Q36/L6858.4 ± 0.8458.1 ± 17.91186.9 ± 50.0359.3 ± 0.239.1 ± 0.4A54/R5742.0 ± 0.7331.5 ± 11.5856.3 ± 32.9351.0 ± 0.229.2 ± 0.3Q54/R5742.1 ± 0.6291.7 ± 11.0748.9 ± 31.7348.5 ± 0.227.4 ± 0.3E54/A5745.8 ± 0.7317.5 ± 12.2817.7 ± 35.0350.5 ± 0.228.8 ± 0.3E54/L5756.9 ± 0.6356.3 ± 19.8907.2 ± 55.2358.0 ± 0.336.1 ± 0.4A54/A5740.9 ± 0.5366.0 ± 15.9952.1 ± 45.0351.9 ± 0.230.0 ± 0.3A54/L5753.2 ± 0.9395.8 ± 16.21012.5 ± 45.5359.6 ± 0.238.4 ± 0.4Q54/A5741.3 ± 0.5316.6 ± 14.0816.7 ± 39.9349.8 ± 0.228.2 ± 0.3Q54/L5752.9 ± 0.7400.7 ± 17.61034.4 ± 49.4357.0 ± 0.235.6 ± 0.4 Open table in a new tab Ssh10b is high resistant to heat-induced denaturation in phosphate buffer (28Xu S. Qin S. Pan X.M. Biochem. J. 2004; 382: 433-440Crossref PubMed Scopus (23) Google Scholar), but we found that its thermostability decrease greatly in a monovalent ion buffer resulting in accessibility of the heat-induced unfolding transition to physical measurement. The heat-induced unfolding and refolding curves of Ssh10b and its variants in buffer H monitored by far-UV CD at 222 nm are nearly superimposable, and for each one the CD spectrum of the refolded protein is nearly identical to that of the native protein (data not shown), indicating that the processes of heat-induced unfolding of Ssh10b and its variants are fully reversible. Fig. 4 shows representative heat-induced unfolding profiles for Ssh10b and its partial variants in buffer H. The heat-induced unfolding profiles of Ssh10b and its variants were analyzed by using the two-state denaturation model described under “Experimental Procedures”: the parameters ΔHm and ΔSm were obtained by fitting Kobs to the van't Hoff plot, the parameter Tm was then obtained by substituting ΔHm and ΔSm into Equation 6. The Tm values of these variants were all around 353 K, and so the stability free energy values of these proteins at 353 K (ΔG(353)) were calculated using the van't Hoff plot (Table 1). As with the results for urea-induced unfolding, the difference between different variants in stability free energy at 353 K was also marked: the difference in ΔG(353) between mutants and the wild-type ranged from -8.3 to 7.5 kJ/mol. Estimating the Contribution of Salt Bridges to Stability by Using DMC Analysis—DMC analysis was used to estimate the strength of the two conserved salt bridges (Glu-36/Lys-68 and Glu-54/Arg-57) in protein Ssh10b. Fig. 5 shows a general scheme for a DMC consisting of two single and one double mutant. A thermodynamic cycle is set up between the wild type, each of the single mutants, and the double mutant. The effects of single amino acid substitutions of each residue involved in the salt bridge and the effect of simultaneous replacement at both positions on the stability of the protein were measured. If the two charged residues do not interact with each other, the effect of substitution of either of the two residues will be independent of the replacement of the other, in other words, the effect on protein stability resulting from their simultaneous replacement will be equal to the sum of the effects of the two single mutations. By contrast, if the two residues do interact with each other, the effect of substituting either of the two residues will depend on the substitution of the other. By using the DMC method, the so-called coupling free energy of a salt bridge can be obtained. As shown in Fig. 5, the coupling free energy (ΔΔGcoup) is defined in Equations 8, 9, 10, 11 (44Bosshard H.R. Marti D.N. Jelesarov I. J. Mol. Recognit. 2004; 17: 1-16Crossref PubMed Scopus (225) Google Scholar). ΔΔGCoup=ΔΔGwt→dm−ΔΔGwt→m1−ΔΔGwt→m2 (Eq. 8) where, ΔΔGwt→dm=ΔGdm−ΔGwt(Eq. 9) ΔΔGwt→m1=ΔGm1−ΔGwt(Eq. 10) ΔΔGwt→m2=ΔGm2−ΔGwt(Eq. 11) So Equation 8 can be simplified to Equation 12. ΔΔGCoup=ΔGdm−ΔGm1−ΔGm2+ΔGwt(Eq. 12) Here, because Ssh10b is a homodimer, the coupling free energy for each salt bridge in the molecule can be represented as Equation 13. ΔΔGCoup=12(ΔGdm−ΔGm1−ΔGm2+ΔGwt)(Eq. 13) Fig. 6, a and b shows the results of different DMCs for the two conserved salt bridges (Glu-36/Lys-68 and Glu-54/Arg-57) in protein Ssh10b at 298 K. For each salt bridge, four different DMCs were constructed to evaluate the interaction strength between the two charged residues involved in the salt bridge. For Glu-36/Lys-68, the reference proteins (doubly substituted variants) for the four DMCs were Ala-36/Ala-68, Ala-36/Leu-68, Gln-36/Ala-68, and Gln-36/Leu-68. Though the reference protein for each DMC has a different stability, all four DMCs gave identical, within the experimental error, values for the coupling free energy of the salt bridge Glu-36/Lys-68 (5.9, 6.0, 6.2, and 5.8 kJ/mol, respectively) at 298 K. For Glu-54/Arg-57, the reference proteins for the four DMCs were Ala-54/Ala-57, Ala-54/Leu-57, Gln-54/Ala-57, and Gln-54/Leu-57. The four DMCs of Glu-54/Arg-57 also gave identical, within the experimental error, coupling free energy values (2.1, 2.7, 2.3, and 2.5 kJ/mol, respectively) at 298 K. Thus, we can assign favorable Gibbs free energy of 6.0 ± 0.2 kJ/mol and 2.4 ± 0.4 kJ/mol to the salt bridges Glu-36/Lys-68 and Glu-54/Arg-57, respectively, at 298 K. Fig. 6, c and d shows the results of different DMCs for the two conserved salt bridges in protein Ssh10b at 353 K. For Glu-36/Lys-68, the four coupling free energy values at 353 K were 5.4, 5.9, 6.3, and 5.9 kJ/mol, and for Glu-54/Arg-57, the coupling free energy values of the four DMCs at 353 K were 2.0, 2.6, 2.1, and 2.1 kJ/mol. We can assign favorable Gibbs free energy of 5.9 ± 0.6 kJ/mol and 2.2 ± 0.4 kJ/mol to the salt bridges Glu-36/Lys-68 and Glu-54/Arg-57, respectively at 353 K. Remarkably, these results were highly consistent with those at 298 K. Experimental estimates of the energy contribution of a salt bridge to protein stability have led to conflicting conclusions, ranging from them having stabilizing, insignificant, or even destabilizing effects (17Marqusee S. Sauer R.T. Protein Sci. 1994; 3: 2217-2225Crossref PubMed Scopus (131) Google Scholar, 18Lounnas V. Wade R.C. Biochemistry. 1997; 36: 5402-5417Crossref PubMed Scopus (93) Google Scholar, 19Lebbink J.H. Consalvi V. Chiaraluce R. Berndt K.D. Ladenstein R. Biochemistry. 2002; 41: 15524-15535Crossref PubMed Scopus (26) Google Scholar, 20Barril X. Aleman C. Orozco M. Luque F.J. Proteins. 1998; 32: 67-79Crossref PubMed Scopus (81) Google Scholar, 21Marti D.N. Bosshard H.R. J. Mol. Biol. 2003; 330: 621-637Crossref PubMed Scopus (47) Google Scholar, 22Sun D.P. Sauer U. Nicholson H. Matthews B.W. Biochemistry. 1991; 30: 7142-7153Crossref PubMed Scopus (264) Google Scholar, 23Waldburger C.D. Schildbach J.F. Sauer R.T. Nat. Struct. Biol. 1995; 2: 122-128Crossref PubMed Scopus (339) Google Scholar). In some reports, electrostatic interactions of salt bridges were measured by mutating a charged residue to a non-charged residue and measuring the impact of this change on protein stability. This type of approach is unsuitable as such mutations not only remove charge-charge interactions but also alter a number of other interactions, including the desolvation penalty and the background interactions of the residue within the protein. However, DMC analysis provides an elegant method to estimate the strength of a salt bridge. The contribution of a salt bridge (ΔΔGbri) to protein stability involving two charged residues can be represented as Equation 14, ΔΔGbri=ΔΔGdir+ΔΔGdesol+ΔΔGbackg(Eq. 14) where ΔΔGdir is the direct contribution of the salt bridge to protein stability. In most situations, it is the direct charge-charge interaction, but it sometimes also includes other direct interactions, e.g. van der Waals interactions between atoms of the two charged residues. The other two components are all indirect contributions of the salt bridges: ΔΔGdesol is the desolvation penalty of the two charged residues, and ΔΔGbackg is the stability contributions due to background interactions of the two charged residues with other components of the protein. Generally, the desolvation penalty is regarded as the dominant component of the indirect contributions (44Bosshard H.R. Marti D.N. Jelesarov I. J. Mol. Recognit. 2004; 17: 1-16Crossref PubMed Scopus (225) Google Scholar). DMC analysis is designed to cancel all interactions except the direct interactions between the two residues involved in the salt bridge. However, the validity of the method is subject to the following two assumptions. First, all indirect interactions are simply additive in nature from the single to the double mutation. Secondly, there is no interaction between the two residues in the doubly substituted mutant. Then, in an ideal DMC, according to Equations 12 and 13, indirect contributions are first removed for each residue (-ΔGm1 -ΔGm2) and then added back again by (+ΔGwt). Hence, in an ideal DMC, the coupling free energy is the direct contribution of the salt bridge. To test the above assumptions, some DMC analyses have been carried out alongside structural studies (34Vaughan C.K. Harryson P. Buckle A.M. Fersht A.R. Acta Crystallogr. Sect. D Biol. Crystallogr. 2002; 58: 591-600Crossref PubMed Scopus (35) Google Scholar, 35Luisi D.L. Snow C.D. Lin J.J. Hendsch Z.S. Tidor B. Raleigh D.P. Biochemistry. 2003; 42: 7050-7060Crossref PubMed Scopus (65) Google Scholar). Another effective method is to construct several independent DMCs for a specific salt bridge. If the coupling free energy is independent of different DMCs, the assumptions described above are likely to be correct (36Makhatadze G.I. Loladze V.V. Ermolenko D.N. Chen X. Thomas S.T. J. Mol. Biol. 2003; 327: 1135-1148Crossref PubMed Scopus (193) Google Scholar, 44Bosshard H.R. Marti D.N. Jelesarov I. J. Mol. Recognit. 2004; 17: 1-16Crossref PubMed Scopus (225) Google Scholar). Here, we used the DMC approach to estimate the net strength of the two highly conserved salt bridges, Glu-54/Arg-57 and Glu-36/Lys-68 in protein Ssh10b. Four independent DMCs were generated for each salt bridge to test the feasibility of this method. The coupling free energies were 2.4 ± 0.4 kJ/mol at 298 K and 2.2 ± 0.4 kJ/mol at 353 K for Glu-54/Arg-57, and 6.0 ± 0.2 kJ/mol at 298 K and 5.9 ± 0.6 kJ/mol at 353 K for Glu-36/Lys-68. Our results demonstrate that, though the stability free energy of Ssh10b decrease greatly with increasing temperature, the direct stability contribution of these two salt bridges remain almost constant, providing evidence supporting the theoretical prediction that salt bridge interactions are extremely resilient to temperature increases and thus are specially suited to promoting protein stability at high temperatures (26Thomas A.S. Elcock A.H. J. Am Chem. Soc. 2004; 126: 2208-2214Crossref PubMed Scopus (77) Google Scholar). This may be the reason why thermophilic and hyperthermophilic proteins have increasing numbers of salt bridges with increasing growth temperature. It has been shown that a higher transition temperature in a hyperthermophilic protein can be obtained in three theoretical ways: 1) by shifting the stability curve to increase the overall unfolding free energy at any temperature, 2) by decreasing the ΔCp between the folded and unfolded states to flatten the stability curve, or 3) by shifting the stability curve toward higher temperatures (45Nojima H. Ikai A. Oshima T. Noda H. J. Mol. Biol. 1977; 116: 429-442Crossref PubMed Scopus (130) Google Scholar, 46Nojima H. Ikai A. Noda H. Hon-nami K. Oshima T. Friedman S.M. Biochemistry of Thermophily. Academic Press, New York1978: 305-323Crossref Google Scholar). The Ssh10b structure is stabilized by eight isolated salt bridges and two ion-pair networks (each involves three salt bridges). Because salt bridge interactions are extremely resilient to temperature increases, increasing the number of salt bridges in Ssh10b should decrease ΔCp between the folded and unfolded states of the protein and thus flatten its stability curve. According to a previous report, Glu-54/Arg-57 and Glu-36/Lys-68 form two salt bridges, which are located on the surface of the Ssh10b molecule (30Wardleworth B.N. Russell R.J. Bell S.D. Taylor G.L. White M.F. EMBO J. 2002; 21: 4654-4662Crossref PubMed Scopus (132) Google Scholar). However according to solvent accessible surface area (ASA) calculations, the four residues involved in these salt bridges are only partially solvent exposed. The solvent ASA of side chains of these residues are: Glu-54: 5%, Arg-57: 28%, Glu-36: 40%, and Lys-68: 58%. The two residues of the Glu-54/Arg-57 salt bridge are both located in the second α-helix, while those of the Glu-36/Lys-68 salt bridge are far apart in primary sequence and located in two separate β-strands (Fig. 1). Because the dielectric constant inside a protein is lower than that of the protein surface, it is generally accepted that the direct stability contribution of a salt bridge buried in a protein should be higher than that of a salt bridge located on the surface (47Tissot A.C. Vuilleumier S. Fersht A.R. Biochemistry. 1996; 35: 6786-6794Crossref PubMed Scopus (126) Google Scholar). However, though the residues involved in the Glu-36/Lys-68 salt bridge are more exposed than those of the Glu-54/Arg-57 salt bridge, the coupling free energies of Glu-36/Lys-68 (6.0 kJ/mol at 298 K and 5.9 kJ/mol at 353 K) are much higher than those of Glu-54/Arg-57 (2.4 kJ/mol at 298 K and 2.2 kJ/mol at 353 K). The contribution of a salt bridge to protein stability is defined as the free energy of unfolding the protein containing the salt bridge minus that of unfolding the protein without the salt bridge (26Thomas A.S. Elcock A.H. J. Am Chem. Soc. 2004; 126: 2208-2214Crossref PubMed Scopus (77) Google Scholar). If a salt bridge interaction could exist in both a native and unfolded state, the stability contribution of this salt bridge would be lower than that if the salt bridge only existed in the native state. In the primary sequence of Ssh10b, the Glu-54 and Arg-57 residues are separated by two hydrophobic residues, Val-55 and Ile-56. In the unfolded state, the most energy favorable conformation should arise when Glu-54 and Arg-57 interact with water, while residues Val-55 and Ile-56 bind to hydrophobic areas of the protein. This conformation would increase the frequency of salt bridge contacts formed by Glu-54 and Arg-57. However, the situation for the Glu-36/Lys-68 salt bridge is different: the two residues are far apart in the primary sequence, so the frequency of salt bridge contacts formed by Glu-36 and Lys-68 in the unfolded state is much lower. To determine whether conclusions reached from our experiments on salt bridges in Ssh10b are applicable to other proteins, we analyzed published DMC experimental data for the contribution of salt bridges to protein stability in other proteins (Table 2). This confirmed that the average coupling free energy of intrahelical salt bridges (typically formed by two very close charged residues) is lower than that of salt bridges formed between residues far apart in sequence. Five of the proteins included in Table 2 have intrahelical salt bridges; coupling free energies (ΔΔGcoup) for these salt bridges determined by DMC experiments are: 0.3 and 0.9 kJ/mol for Asp-12/Arg-16 and Glu-28/Lys-32 respectively in Barnase, 3.4 kJ/mol for Asp-14/Arg-17 in λ Repressor, 0.4 kJ/mol for Asp-116/Arg-119 in T4 lysozyme, 1.2 and 2.1 kJ/mol for Lys-8/Glu-11 and Glu-11/Lys-15, respectively, in GCN4-p1, and 2.4 kJ/mol for Glu-54/Lys-57 in Ssh10b as reported in this work. The average coupling free energy of salt bridges located in an α-helix is only about 1.5 kJ/mol, much lower than that of salt bridges found in other structures (about 5.4 kJ/mol). This finding provides further insight into how to engineer proteins for improved stability at high temperatures; salt bridges introduced between two residues, which are close to each other in tertiary structure but far from each other in primary structure will make a greater contribution to protein stability. Further experimental evidence is required to substantiate these findings.TABLE 2The contribution of salt bridges to protein stability (DMC experimental data)Salt bridgeΔΔGcoup kJ mol–1Intrahelix (yes/no)Barnase (32,47–49)D8/R1103.5NoD12/R1101.4NoD12/R160.9YesE28/K320.3YesR69/D9314.3NoRubredoxin (50Strop P. Mayo S.L. Biochemistry. 2000; 39: 1251-1255Crossref PubMed Scopus (151) Google Scholar)K6/E491.2NoNt/E146.3NoZinc finger (51Blasie C.A. Berg J.M. Biochemistry. 1997; 36: 6218-6222Crossref PubMed Scopus (62) Google Scholar)K3/E140.8NoGCN4-p1 (52Ibarra-Molero B. Zitzewitz J.A. Matthews C.R. J. Mol. Biol. 2004; 336: 989-996Crossref PubMed Scopus (46) Google Scholar)K8/E111.2YesE11/K152.1YesNTL9 (35Luisi D.L. Snow C.D. Lin J.J. Hendsch Z.S. Tidor B. Raleigh D.P. Biochemistry. 2003; 42: 7050-7060Crossref PubMed Scopus (65) Google Scholar)Nt/D237.1NoSsh10bE36/K686.0NoE54/R572.4YesProtein G (53Merkel J.S. Sturtevant J.M. Regan L. Structure. 1999; 7: 1333-1343Abstract Full Text Full Text PDF PubMed Scopus (92) Google Scholar,54Lassila K.S. Datta D. Mayo S.L. Protein Sci. 2002; 11: 688-690Crossref PubMed Scopus (20) Google Scholar)R6/E53(A44)2.5NoR6/E53(R44)2.7NoE6/R533.1NoK6/E531.6NoR44/E53(16)2.2NoR44/E53(R6)2.4Noλ repressor (17)D14/R173.4YesT4 lysozyme (55Nicholson H. Anderson D.E. Sun D.P. Matthews B.W. Biochemistry. 1991; 30: 9816-9828Crossref PubMed Scopus (160) Google Scholar)D116/R1190.4YesUbiquitin (36)E34/K113.6NoOmpA (56Hong H. Szabo G. Tamm L.K. Nat. Chem. Biol. 2006; 2: 627-635Crossref PubMed Scopus (103) Google Scholar)E52/K8214.7NoE52/R13823.5NoK68/E1287.1NoE128/R1382.5No Open table in a new tab