|
|
||||||||
1 Departments of Biochemistry, 2 Molecular and Medical Genetics, and 3 Chemistry, University of Toronto, Toronto, Ontario, M5S-1A8, Canada
Reprint requests to: Alan R. Davidson, Department of Molecular and Medical Genetics, University of Toronto, Toronto, Ontario M5S-1A8, Canada; e-mail alan.davidson{at}utoronto.ca; fax (416) 978-6885.
(RECEIVED September 1, 2005; FINAL REVISION September 1, 2005; ACCEPTED October 11, 2005)
| Abstract |
|---|
|
|
|---|
Keywords: SH3 domain; protein folding kinetics; guanidinium-induced protein stabilization; peptide binding; arginineprotein interaction; specific guanidinium binding
Article and publication are at http://www.proteinscience.org/cgi/doi/10.1110/ps.051829106.
| Introduction |
|---|
|
|
|---|
The effects of GuHCl on protein stability are particularly complicated due to its ionic character. Ions can either bind to the folded or unfolded states of proteins or influence stability through "screening" coulombic interactions, the Hofmeister effect, or changing the structure of the solvent (for a recent review, see Collins 2004). Thus, GuHCl can modulate protein stability through mechanisms related to both its denaturant and ionic properties. Although GuHCl is thought of primarily as a denaturant, there have been many cases reported where low concentrations of GuHCl actually stabilize proteins (Mayr and Schmid 1993; Monera et al. 1994; Makhatadze et al. 1998; Bhuyan 2002). In the majority of these cases, it has been shown that the ionic nature of GuHCl is the main factor in its stabilizing effect (Monera et al. 1994; Makhatadze et al. 1998). However, some studies have suggested that GuHCl may cause protein stabilization in a more specific manner (Mayr and Schmid 1993). Solution of a crystal structure of Ribonuclease A in GuHCl demonstrated direct binding of this compound to the folded state of the enzyme, resulting in a decrease in temperature factors (Dunbar et al. 1997). However, it was observed that guanidine bound to sites comprised of various amino acid residues with no obvious sequence or structural similarity among these sites detected. Guanidinium has been recently shown to be able to mimic the interactions mediated by the guanidino (H2N-CNH-NH2)+ moiety of an arginine side chain and restabilize a T4 Lysozyme variant with an arginine to alanine mutation (Yousef et al. 2004).
In a previous study performed in our laboratory examining the thermodynamic stabilities of several site-directed mutants of the SH3 domain of the Fyn tyrosine kinase, we found that the use of GuHCl as a denaturant in unfolding experiments resulted in measured unfolding free energy (
Gu) values that were higher by ~1.1 kcal/mol compared with those obtained using urea-induced unfolding experiments or NMR spin relaxation dispersion spectroscopy in the absence of denaturants (Maxwell and Davidson 1998; Northey et al. 2002; Di Nardo et al. 2004). Our goal in the studies presented here is to elucidate the mechanism by which GuHCl stabilizes the Fyn SH3 domain. We have used folding kinetics studies, functional assays, mutagenesis, and NMR spectroscopy to systematically investigate the effects of GuHCl on the structure and stability of this domain, and we demonstrate the existence of a specific binding site for guanidinium on the surface of this protein.
| Results and Discussion |
|---|
|
|
|---|

Gu value of 2.61 kcal/mol with respect to wild type. The transition midpoint of urea denaturation of this mutant was shifted to a lower concentration of urea so that more data points could be collected along the unfolding arm of the chevron plot. For both the wild type and R40N mutant, it can be seen that GuHCl is indeed highly stabilizing at low concentrations with free energy of unfolding (
Gu) of the wild-type domain increased by 1 kcal/mol in 0.4 M GuHCl and that of the R40N mutant increased by 0.6 kcal/mol at 0.2 M GuHCl (Table 1
Gu) of the R40N mutant assumes a strong linear relationship (r = 0.99) with GuHCl concentration (plot not shown). The slope of the best-fit line (1.71 kcal/mol.M) obtained using the R40N data is almost identical to the previously reported "meq value" (1.73 kcal/mol.M) for the Fyn SH3 domain (Maxwell and Davidson 1998).
|
|
The deceleration of the unfolding rate in the presence of GuHCl was accompanied by no drastic change in the folding rate, which suggests that guanidinium exclusively interacts with the native state of this protein. Consistent with this notion, the mkf values remain unaltered upon increasing the concentrations of GuHCl (Table 1
), whereas the mku values increase with increasing overall stability (
Gu). As illustrated in Figure 2
, the strong correlation (r = 0.99 for wild type and r = 0.92 for R40N) between mku and
Gu argues for exclusive binding of guanidinium to the native state of the Fyn SH3 domain. Since the mku value is believed to be proportional to the difference in the exposed surface area between the native state and the folding transition state (Sanchez and Kiefhaber 2003), these results suggest that guanidinium binding buries part of the surface of the native state. Thus, a relatively larger surface area would become solvent exposed upon unfolding to the transition state of the guanidinium-bound Fyn SH3 domain. In addition, the lack of change of the mkf value suggests that GuHCl does not preferentially interact with either the unfolded or the folding transition state of this protein.
|
|
|
strand. Interestingly, these residues are part of a conserved functional surface on SH3 domains that is often involved in binding to the side chain of a conserved arginine residue found in many SH3 target peptides. For example, the structure of a complex between the c-Src tyrosine kinase SH3 domain and the phage display-derived Vsl12 peptide ligand (Feng et al. 1995) reveals that the guanidino moiety of the target peptide arginine at the P3 position (see Fig. 5C
|
|
|
Does GuHCl stabilize other SH3 domains?
To address the question of whether stabilization by GuHCl might be a general feature of proteins that interact with arginine, we examined the ability of GuHCl to stabilize two other SH3 domains, one from the yeast Actin Binding Protein 1 (Abp1) protein and the other from the yeast Sho1 protein. As a simple means to test the stabilizing effect of GuHCl on these domains, we performed temperature-induced unfolding experiments in varying amounts of GuHCl. As expected, the wild-type Fyn SH3 domain exhibited an increase in the midpoint of its temperature-induced unfolding transition (Tm) in the presence of low concentrations of GuHCl, which was maximal at 0.2 M GuHCl, where the Tm was increased by 3.6°C (Fig. 8
). The Sho1 SH3 domain exhibited an even greater degree of stabilization in the presence of GuHCl with a Tm increase of more than 7°C at 0.1 M GuHCl. In contrast, the Abp1 SH3 domain did not exhibit an increase in Tm at any concentration of GuHCl. Since the Sho1 SH3 domain has been shown to bind strongly to target peptides possessing arginine at the P3 position (Marles et al. 2004), the stabilization of this domain by GuHCl is not surprising. On the other hand, the Abp1 SH3 domain binds tightly to a different class of peptide ligands containing lysine at the P3 position (Rath and Davidson 2000; Zarrinpar et al. 2004), and this domain has never been directly shown to tightly bind peptides containing arginine at the P3 position4. These data suggest that the ability of SH3 domains to be stabilized by GuHCl correlates with their ability to bind to peptides possessing arginine at the P3 position.
|
| Conclusions |
|---|
|
|
|---|
| Materials and methods |
|---|
|
|
|---|
Mutagenesis and protein purification
All of the proteins and protein variants used in this study were expressed as C-terminal hexahistidine fusions in Escherichia coli strain Bl21* (DE3). Recombinant proteins were purified using nickel affinity chromatography under denaturing conditions as described previously (Maxwell and Davidson 1998). Proteins were folded through dialysis and used as such without cleaving the hexahistidine tag. The dialysis buffer consisted of a 50 mM sodium phosphate solution (pH 7.0). Protein folding and functional assays were performed in this buffer.
Folding kinetics
Folding and unfolding rates were measured at 298 K on a Bio-Logic SFM-4 stopped-flow device equipped with a Photo Multiplier Tube (PMT) monitoring the recovery of intrinsic Trp fluorescence upon folding of urea denatured protein. Excitation was carried out at 295 nm and all of the fluorescence above 309 nm was collected. Traces were fit to appropriate single exponential functions using BioKine. At each concentration of urea, at least five separate shots of 35 µM protein were averaged. Assuming a linear dependence of ln kobs on urea concentration, kinetic chevron plots were fit to chevron equation:
![]() |
where kf and ku are the folding and unfolding rates at 0Murea and mkf and mku denote the dependence of kf and ku on the concentration of denaturant. Fitting was performed by a least squares method using Kaleidagraph (Synergy software).
NMR spectroscopy and quantitative peak shift studies
NMR experiments were performed on 0.25 mM uniformly 15N isotopically enriched protein samples containing 50 mM sodium phosphate (pH 7.0), 100 mM NaCl, 0.2 mM EDTA, and 4 mM DSS at 25°C using a Varian Inova spectrometer operating at 500 MHz 1H frequency. GuHCl was added to final concentrations of 0, 5, 10, 15, 20, 29, 38, 57, 83, 115, 153, 198, 260, 350, 436, 538, 646, 750, and 851 mM. Spectra were processed and peak positions quantitated using the NMRPipe/NMRDraw suite of programs (Delaglio et al. 1995), referencing 1H and 15N chemical shifts relative to the 1H methyl resonances of DSS (Wishart et al. 1995). Assignment of NMR spectra was accomplished with a 15N, 13C-labeled protein sample under identical buffer conditions using HNCACB (Wittekind and Mueller 1993) and TOCSY-HSQC (Montelione et al. 1992; Logan et al. 1993) pulse sequences.
Amide 15N and 1H chemical shifts showing changes greater than 0.4 or 0.08 ppm, respectively, were fit separately to the equation
![]() | ([1]) |
where
(C) is the chemical shift at GuHCl concentration C,
0 =
(0),
F =
(
), and m is a baseline correction factor discussed below. fB is the fraction of protein bound to GuHCl, which is given by
![]() | ([2]) |
The parameters
0,
F, m, and Kd were allowed to vary in least-squares fits of the data. The weighted mean of the dissociation constant was calculated according to
![]() | ([3]) |
where
j is the experimental uncertainty in the dissociation constant, Kj, calculated from a single set of 1H or 15N chemical shifts, and the sum extends over all individual Kd estimates. The uncertainty in the weighted mean was computed using a boostrap simulation (Efron and Tibshirani 1986).
Linear baseline corrections were found to be necessary for adequate fits of experimental data. Extracted values range from m = 0.30 ppm/M for A12 to 0.06 ppm/M for T14 1H chemical shifts; in the case of 15N chemical shifts, the range extends from 1.29 ppm/M for T14 to 0.39 ppm/M for D35. This linear dependence on GuHCl concentration is likely due to either slight noncooperative changes in protein conformation or very weak guanidinium binding in addition to the stronger (Kd
90 mM) interaction.
Binding studies
Dissociation constants were measured using the changes in the intrinsic fluorescence of the Trp residue located in the target binding pocket of the Fyn SH3 domain as a probe for binding. A
CI fusion construct of a previously characterized Vsl12 target peptide (Maxwell and Davidson 1998) was used in these studies. Fluorescence signal was collected at equilibrium on an Aviv Spectrofluorometer Model ATF 105 (Aviv Associates). After subtraction of the fluorescence contribution from target alone from that of the complex, fluorescence data were fit as described previously (Maxwell and Davidson 1998).
Temperature melts
The changes in the ellipticity (mdeg) at 220 nm were recorded at equilibrium at different temperatures on an Aviv Circular Dichroism spectrometer Model 62A DS (Aviv Associates). Sigmoidal unfolding curves were fit by nonlinear least-squares methods using the program Igor Pro to the standard equations to obtain the Tm as described previously (Maxwell and Davidson 1998).
Error analysis
The errors reported for kobs, mkf, mku, Kd, and Tm are uncertainties associated with the least-squares fitting of the data to appropriate equations. Folding kinetics experiments, temperature melts, and binding assays were repeated at least twice. In all cases, the variations observed in kobs, mkf, mku, Kd, and Tm parameters were <5%. Uncertainties reported for
G,
Tm are by error propagation according to the general equation:
![]() |
where u = f(x,y) and dx and dy are the residual errors of x and y, respectively.
| Footnotes |
|---|
| Acknowledgments |
|---|
| References |
|---|
|
|
|---|
Bogan, A.A. and Thorn, K.S. 1998. Anatomy of hot spots in protein interfaces. J. Mol. Biol. 280: 19.[CrossRef][Medline]
Collins, K.D. 2004. Ions from the Hofmeister series and osmolytes: Effects on proteins in solution and in the crystallization process. Methods 34: 300311.[CrossRef][Medline]
de Los Rios, M.A. and Plaxco, K.W. 2005. Apparent Debye-Huckel electrostatic effects in the folding of a simple, single domain protein. Biochemistry 44: 12431250.[CrossRef][Medline]
Delaglio, F., Grzesiek, S., Vuister, G.W., Zhu, G., Pfeifer, J., and Bax, A. 1995. NMRPipe: A multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 6: 277293.[Medline]
Di Nardo, A.A., Korzhnev, D.M., Stogios, P.J., Zarrine-Afsar, A., Kay, L.E., and Davidson, A.R. 2004. Dramatic acceleration of protein folding by stabilization of a nonnative backbone conformation. Proc. Natl. Acad. Sci. 101: 79547959.
Dunbar, J., Yennawar, H.P., Banerjee, S., Luo, J., and Farber, G.K. 1997. The effect of denaturants on protein structure. Protein Sci. 6: 17271733.[Abstract]
Efron, B. and Tibshirani, R. 1986. Bootstrap methods for standard errors, confidence intervals and other measures of statistical accuracy. Stat. Sci. 1: 5477.
Fazi, B., Cope, M.J., Douangamath, A., Ferracuti, S., Schirwitz, K., Zucconi, A., Drubin, D.G., Wilmanns, M., Cesareni, G., and Castagnoli, L. 2002. Unusual binding properties of the SH3 domain of the yeast actin-binding protein Abp1: Structural and functional analysis. J. Biol. Chem. 277: 52905298.
Feng, S., Kasahara, C., Rickles, R.J., and Schreiber, S.L. 1995. Specific interactions outside the proline-rich core of two classes of Src homology 3 ligands. Proc. Natl. Acad. Sci. 92: 1240812415.
Jones, S. and Thornton, J.M. 1996. Principles of proteinprotein interactions. Proc. Natl. Acad. Sci. 93: 1320.
Kay, L.E., Keifer, P., and Saarinen, T. 1992. Pure absorption gradient enhanced heteronuclear single quantum correlation spectroscopy with improved sensitivity. J. Am. Chem. Soc. 114: 1066310665.[CrossRef]
Koradi, R., Billeter, M., and Wuthrich, K. 1996. MOLMOL: A program for display and analysis of macromolecular structures. J. Mol. Graph. 14: 5155.[CrossRef][Medline]
Logan, T.M., Olejniczak, E.T., Xu, R.X., and Fesik, S.W. 1993. A general method for assigning NMR spectra of denatured proteins using 3D HC(CO)NH-TOCSY triple resonance experiments. J. Biomol. NMR 3: 225231.[Medline]
Makhatadze, G.I. 1999. Thermodynamics of protein interactions with Urea and Guanidinium Hydrochloride. J. Phys. Chem. B 103: 47814785.
Makhatadze, G.I. and Privalov, P.L. 1992. Protein interactions with urea and guanidinium chloride. A calorimetric study. J. Mol. Biol. 226: 491505.[CrossRef][Medline]
Makhatadze, G.I., Lopez, M.M., Richardson 3rd, J.M., and Thomas, S.T. 1998. Anion binding to the ubiquitin molecule. Protein Sci. 7: 689697.[Abstract]
Marles, J.A., Dahesh, S., Haynes, J., Andrews, B.J., and Davidson, A.R. 2004. Proteinprotein interaction affinity plays a crucial role in controlling the Sho1p-mediated signal transduction pathway in yeast. Mol. Cell 14: 813823.[CrossRef][Medline]
Maxwell, K.L. and Davidson, A.R. 1998. Mutagenesis of a buried polar interaction in an SH3 domain: Sequence conservation provides the best prediction of stability effects. Biochemistry 37: 1617216182.[CrossRef][Medline]
Mayr, L.M. and Schmid, F.X. 1993. Stabilization of a protein by guanidinium chloride. Biochemistry 32: 79947998.[CrossRef][Medline]
Moglich, A., Krieger, F., and Kiefhaber, T. 2005. Molecular basis for the effect of urea and guanidinium chloride on the dynamics of unfolded polypeptide chains. J. Mol. Biol. 345: 153162.[CrossRef][Medline]
Monera, O.D., Kay, C.M., and Hodges, R.S. 1994. Protein denaturation with guanidine hydrochloride or urea provides a different estimate of stability depending on the contributions of electrostatic interactions. Protein Sci. 3: 19841991.[Abstract]
Montelione, G.T., Lyons, B.A., Emerson, S.D., and Tashiro, M. 1992. An efficient triple resonance experiment using carbon 13 isotropic mixing for determining sequence-specific resonance assignments of isotopically- enriched proteins. J. Am. Chem. Soc. 114: 1097410975.[CrossRef]
Noble, M.E., Musacchio, A., Saraste, M., Courtneidge, S.A., and Wierenga, R.K. 1993. Crystal structure of the SH3 domain in human Fyn; comparison of the three-dimensional structures of SH3 domains in tyrosine kinases and spectrin. EMBO J. 12: 26172624.[Medline]
Northey, J.G., Di Nardo, A.A., and Davidson, A.R. 2002. Hydrophobic core packing in the SH3 domain folding transition state. Nat. Struct. Biol. 9: 126130.[CrossRef][Medline]
Rath, A. and Davidson, A.R. 2000. The design of a hyperstable mutant of the Abp1p SH3 domain by sequence alignment analysis. Protein Sci. 9: 24572469.[Abstract]
Rickles, R.J., Botfield, M.C., Zhou, X.M., Henry, P.A., Brugge, J.S., and Zoller, M.J. 1995. Phage display selection of ligand residues important for Src homology 3 domain binding specificity. Proc. Natl. Acad. Sci. 92: 1090910913.
Sanchez, I.E. and Kiefhaber, T. 2003. Hammond behavior versus ground state effects in protein folding: Evidence for narrow free energy barriers and residual structure in unfolded states. J. Mol. Biol. 327: 867884.[CrossRef][Medline]
Schellman, J.A. 2002. Fifty years of solvent denaturation. Biophys. Chem. 96: 91101.[CrossRef][Medline]
Wishart, D.S., Bigam, C.G., Yao, J., Abildgaard, F., Dyson, H.J., Oldfield, E., Markley, J.L., and Sykes, B.D. 1995. 1H, 13C and 15N chemical shift referencing in biomolecular NMR. J. Biomol. NMR 6: 135140.[Medline]
Wittekind, M. and Mueller, L. 1993. HNCACB, a high sensitivity 3D NMR experiment to correlate amide proton and nitrogen resonances with the
- and
-carbon resonances in proteins. J. Magn. Reson. Series B 101: 201205.[CrossRef]
Yousef, M.S., Baase, W.A., and Matthews, B.W. 2004. Use of sequence duplication to engineer a ligand-triggered, long-distance molecular switch in T4 lysozyme. Proc. Natl. Acad. Sci. 101: 1158311586.
Zarrinpar, A., Bhattacharyya, R.P., Nittler, M.P., and Lim, W.A. 2004. Sho1 and Pbs2 act as coscaffolds linking components in the yeast high osmolarity MAP kinase pathway. Mol. Cell 14: 825832.[CrossRef][Medline]
![]()
CiteULike
Connotea
Del.icio.us
Digg
Reddit
Technorati What's this?
This article has been cited by other articles:
![]() |
R. F. Latypov, D. Liu, K. Gunasekaran, T. S. Harvey, V. I. Razinkov, and A. A. Raibekas Structural and thermodynamic effects of ANS binding to human interleukin-1 receptor antagonist Protein Sci., April 1, 2008; 17(4): 652 - 663. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. L. Schweiker, A. Zarrine-Afsar, A. R. Davidson, and G. I. Makhatadze Computational design of the Fyn SH3 domain with increased stability through optimization of surface charge charge interactions Protein Sci., December 1, 2007; 16(12): 2694 - 2702. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |