PEDS Advance Access published online on January 5, 2008
Protein Engineering Design and Selection, doi:10.1093/protein/gzm077
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
The folding pathway of an engineered circularly permuted PDZ domain
Istituto di Biologia e Patologia Molecolari del CNR, Dipartimento di Scienze Biochimiche A. Rossi Fanelli, Sapienza - Università di Roma La Sapienza, Piazzale A. Moro 5, 00185 Rome, Italy
1 To whom correspondence should be addressed. E-mail: stefano.gianni{at}uniroma1.it
| Abstract |
|---|
|
|
|---|
To understand the role of sequence connectivity in the folding pathway of a multi-state protein, we have analysed the folding kinetics of an engineered circularly permuted PDZ domain. This variant has been designed with the specific aim of posing two of the strands participating in the stabilisation of an early folding nucleus as contiguous elements in the primary structure. Folding of the circularly permuted PDZ2 has been explored by a variety of different experimental approaches including stopped-flow and continuous-flow kinetics, as well as ligand-induced folding experiments. Data reveal that although circular permutation introduces a significant destabilisation of the native state, a folding intermediate is stabilised and accumulated prior folding. Furthermore, quantitative analysis of the observed kinetics indicates an acceleration of the early folding events by more than two orders of magnitude. The results support the importance of sequence connectivity both in the mechanism and the speed of protein folding.
Keywords: circular permutation/intermediates/kinetics/protein folding/transition state
| Introduction |
|---|
|
|
|---|
Native state topology has been proposed to be one of the main factors in controlling both speed (Munoz and Eaton, 1999
PDZ domains are small (80–100 amino acids) globular protein–protein interacting modules, with important roles in molecular recognition (Harris and Lim, 2001
; Jemth and Gianni, 2007
). The typical fold of the canonical PDZ domains involves six β-strands and two
-helices (Fig. 1A). PDZ domains are generally part of multi-domain proteins and exert their function by recognising small target motifs, exploiting a specific binding groove. The folding of murine PDZ2 from PTP-BL has been thoroughly investigated both by experimental and theoretical approaches. PDZ2 folds through a three-state mechanism involving two different transition states (namely a denatured like TS1 and a native-like TS2) with an obligatory high-energy on-pathway folding intermediate (Gianni et al., 2005a
, 2005b
; Ivarsson et al., 2007
). Molecular dynamics simulations using
-values as structural restrains (Vendruscolo et al., 2001
) clarified that TS1 consists of a heterogeneous ensemble of quite expanded structures (Gianni et al., 2007
); nevertheless, in all structures a specific cluster of sufficiently stable interactions between residues in β-strands 1, 4 and 6 can be identified (Fig. 1B). The formation of a β-sheet involving contacts between residues from the N- and C-termini can hence be seen as a key-event in PDZ2 folding. Here, we investigate the folding mechanism of an engineered circularly permuted PDZ2 (cpPDZ2), where the canonical N-terminal strand is shifted to the C-terminus position (Fig. 1A); thus two of the strands participating in the stabilisation of the early folding nucleus are contiguous elements in the primary structure. We studied cpPDZ2 folding using a variety of different experimental approaches including stopped-flow and continuous-flow kinetics as well as ligand-induced folding experiments. The results on cpPDZ2 suggest that although this circular permutation significantly destabilises the native state, a folding intermediate is stabilised and folding turns from single exponential to double exponential. Furthermore, a comparison with PDZ2 folding suggests that such mutation speeds up the early events of protein folding, as revealed by a lower energy barrier for TS1.
|
| Materials and methods |
|---|
|
|
|---|
Design of circularly permuted PDZ2
An engineered circular permutant of PDZ2 (cpPDZ2) was designed by cleaving between D22 and G23 [numbering according the murine PDZ2 structure 1GM1 (Walma et al., 2002
)], and linking residues G99 and G11 together. Residues D22 and G23 occur in the structurally variable, solvent accessible loop between the structurally conserved regions including β-strands 2 and 3. The introduction of a cleavage site within this loop is therefore not expected to significantly affect the PDZ domain structure. In the structure of the PDZ-like domain of the stress sensor protease DegS (PDB ID: 1TE0), a natural circular permutant of mouse PDZ2, residues D338 and D339 (structurally equivalent to G99 and G11 of mouse PDZ2, respectively) are linked together and the main-chain conformation of all structurally equivalent secondary structure elements between the two domains is conserved. This further suggests that the linkage of residues G99 and G11 of PDZ2 should not affect significantly the rest of the structure.
As previously described for other PDZ domains (Chi et al., 2007
), a tryptophan residue was introduced to function as a fluorescent probe (Y43W). A synthetic gene encoding cpPDZ2 was purchased from GENEART (Germany).
Protein expression and purification
The gene was subcloned into the expression-vector pET28(c) and protein was expressed in BL21(DE3) cells (Invitrogen). After induction with 1 mM IPTG, the cells were grown for 24 h at 25°C. The resulting hexahis-tagged protein formed inclusion bodies and was resuspended at high [urea]. The protein was purified on an S-Sepahorose ion exchange column (Amersham Biosciences) and eluted with a gradient of 0–2 M sodium chloride in 50 mM sodium phosphate pH 6.3.
Ligand-binding kinetics of cpPDZ2 were monitored by Förster resonance energy transfer (FRET) between the fluorescence donor Y43W and a dansyl-group covalently attached to the N-terminal of the target peptide D-EQVSAV, as previously described (Gianni et al., 2005b
; Gianni et al., 2006
). The experiments were performed at pH 7.2, 0.4 M sodium sulphate and 10°C with a protein concentration of 1.5 µM and varying the peptide concentration between 10 and 80 µM. The excitation wavelength was 280 nm and the fluorescence emission was measured using a 330 ± 20 nm band-pass glass filter.
Thermal denaturation of cpPDZ2 in 50 mM sodium phosphate buffer pH 7.2, in absence or presence of stabilising 0.4 M sodium sulphate, was monitored by far-UV CD (202 nm) on a JASCO circular dichroism spectrophotometer using a 0.1 cm quartz cuvette (Hellma) and a protein concentration of 16 µM.
Urea-induced equilibrium denaturation of cpPDZ2 was measured by following the decrease in Trp emission at different wavelengths on increasing urea concentration at 25°C, pH 7.2 in 50 mM sodium phosphate buffer in presence of 0.4 M sodium sulphate at a constant protein concentration of 1.5 µM.
Kinetic (un)folding experiments were performed using a Pi-star stopped-flow apparatus (Applied Photopysics, Leatherhead, UK) in 50 mM sodium phosphate, pH 7.2, in absence (10°C) or presence (25°C) of 0.4 M sodium sulphate. The protein sample was excited at 280 nm and the folding reaction was followed by the change in fluorescence using 320 or 360 nm cut-off filters. An 11-fold dilution of denatured or native protein in appropriate buffer initiated refolding and unfolding.
Peptide-induced folding experiments were carried out at pH 7.2, 25°C and 0.4 M sodium sulfate by mixing unfolded cpPDZ2 (i.e. in the presence of 1.0–4.8 M urea) with different concentrations of the peptide D-EQVSAV. The excitation wavelength was 280 nm and the fluorescence emission was measured using a 330 ± 20 nm band-pass glass filter.
Ultra-fast kinetic refolding experiments were carried out using the in-house built continuous-flow apparatus equipped with an ultra-fast mixing device (dead-time 50 µs) as described elsewhere (Jemth et al., 2004
; Gianni et al., 2006
). Folding was initiated by 11-fold dilution in solutions containing 0.44 M sodium sulphate and different urea concentrations. Final protein concentration was typically 15 µM. The folding reaction was monitored by exciting the Trp at 280 nm and following the increase in fluorescence signal using a 360 nm cut-off filter.
| Results |
|---|
|
|
|---|
An engineered circularly permuted PDZ2 was designed posing two of the three β-strands participating in the early folding nucleus (β1 and β6) as contiguous elements in the primary structure (Fig. 1A). We tested the permissibility of such circular permutation in PDZ2 by performing classic ligand-binding experiments under pseudo-first-order conditions (Fig. 2). A dansylated target peptide D-EQVSAV has previously been shown to be a valid ligand for wild-type PDZ2 (kon = 5.0 ± 0.3 µM–1 s–1 and koff = 45 ± 3 s–1, Gianni et al., 2006
|
Thermal denaturation of cpPDZ2 at pH 7.2 in absence and in the presence of 0.4 M sodium sulphate was monitored by far-UV CD as reported in Fig. 3A. The midpoint of the thermal denaturation in absence of sodium sulphate is 303 ± 1 K (the Tm of PDZ2 under the same experimental conditions being 321 ± 1 K), with an on-set of the transition at
290 K. Hence, cpPDZ2 is partially unfolded at pH 7.2 and 25°C. By addition of 0.4 M sodium sulphate, the midpoint of the transition is increased to 326 ± 1 K.
|
The urea-induced equilibrium denaturation of cpPDZ2 in presence of 0.4 M sodium sulphate is reported in Fig. 3B. Assuming a standard two-state model (Jackson and Fersht, 1991
|
|
Gd represents the free energy of folding at a concentration D of denaturant, mD-N is the slope of the transition and D50 is the midpoint of the denaturation transition. The equation used takes into account the pre- and post-transition baselines (Santoro and Bolen, 1988
The kinetics of the cpPDZ2 (un)folding was investigated at pH 7.2 in absence and presence of 0.4 M sodium sulphate at 10°C and 25°C, respectively. The semi-logarithmic plots of the folding/unfolding rate constants of cpPDZ2 versus denaturant concentrations (chevron plots) measured by stopped-flow are shown in Fig. 4, together with the chevron plots of wild-type PDZ2 (also in absence and presence of 0.4 M sodium sulphate) at 25°C. The chevron plot of the circularly permuted protein indicates a significant destabilisation. In fact, when the (un)folding kinetics was investigated in the absence of stabilising salt, only the unfolding rate constants of cpPDZ2 were measurable. On the other hand, in presence of 0.4 M sodium sulphate, the chevron plot of cpPDZ2 is clearly V-shaped, indicating a simple two-state system without detectable intermediates. The data set was therefore fitted to an equation for apparent two-state kinetics. The observed rate constant kobs is the sum of the folding and unfolding rates, kobs = kf + ku and it is dependent on urea concentration according to Eq. (1):
|
|
50 µs (Shastry and Roder, 1998
|
|
Owing to the low stability of cpPDZ2, we decided to strengthen the quantitative analysis of (un)folding kinetics by performing ligand-induced refolding experiments. In these experiments, a large excess of the ligand peptide is added to the denatured protein, thereby trapping newly folded protein by binding, draining the equilibrium towards the folded state. It has been shown that when and if (i) the ligand interacts only with the native protein and (ii) the binding reaction is much faster than the folding process, a determination of folding rates can be obtained independently of the unfolding even under strongly denaturing conditions (Sanz and Fersht, 1994
600 µM to circa 2.5 mM at low and high [urea], respectively). All the observed binding traces followed single exponential behaviour. Data were analysed as previously described (Ivarsson et al., 2007
The urea dependence of the relaxation rate constants for the fast (
1) and slow (
2) phases together with the folding rate independently determined by ligand induced binding experiments (
3) are reported in Fig. 5. The combined data set was globally fitted to the roots of a quadratic equation describing a three-state folding mechanism.
|
|
|
|
|
|
|
|
The resulting rate constants of formation and breakage of the intermediate are k1 = 6000 ± 500 s–1 and k–1 = 300 ± 100 s–1 and the formation and breakage of the natively folded protein have the rate constants k2 = 130 ± 8 s–1 and k–2 = 5.3 ± 0.1 s–1, returning a total mD–N value of 1.15 kcal mol–1 M–1. However, since the two phases do not appear to be kinetically coupled (i.e. they take place on different time scales at all urea concentrations), it is at this point not possible to discriminate between the on- and off-pathway schemes, contrary to what previously observed for other systems (Bai, 1999
; Capaldi et al., 2001
; Travaglini-Allocatelli et al., 2003
).
| Discussion |
|---|
|
|
|---|
The docking of the N- and C-terminal segments appears to be a key event in the folding of PDZ2. Docking of terminal elements as a first step in folding has been suggested for several small proteins (Krishna and Englander, 2005
The strategy of circular permutation has been applied to a number of small protein domains aiming at the investigation of the stability and folding mechanisms upon change of sequence connectivity (Viguera et al., 1996
; Otzen and Fersht, 1998
; Li and Shakhnovich, 2001
; Lindberg et al., 2002
; Miller et al., 2002
). A specific emphasis has been put on the concept of contact order, according to which, in general terms, a protein with many local contacts should fold faster than a protein with more long-range contacts (Plaxco et al., 1998
). Here, we analysed the folding pathway of an engineered circular permutant of a PDZ domain, where two segments of the previously identified folding nucleus are posed as contiguous elements in the primary structure. We have previously shown that PDZ2 folds via a complex folding mechanism involving two energy barriers (a denatured-like TS1 and a native-like TS2) and a high-energy on-pathway intermediate. A comparative study on five different PDZ domains showed that members of this protein family fold via a consensus mechanism where the transition state TS1 is rate limiting in water for all PDZ domain tested. Surprisingly, although our circular permutation destabilises the native state (
Gwt-cp
2 kcal mol–1), the folding kinetics of cpPDZ2 reveal a remarkable stabilisation of the folding intermediate, which accumulates transiently during folding. Furthermore, examination of all four microscopic rate constants allowed us to conclude that cpPDZ2 displays a lower TS1 energy barrier, the microscopic rate constant kDI turning from 35 ± 2 s–1 (in the case of PDZ2) to 6000 ± 500 s–1 (in the case of cpPDZ2). Thus, in agreement with theoretical predictions (Baker, 2000
; Fersht, 2000
) and with direct measurements of loop formation in disordered polypeptides (Buscaglia et al., 2006
), reducing the sequence separation between elements constituting an early folding nucleus produces a remarkable increase of the folding rate constants.
The kinetic folding pathway of cpPDZ2 raises the intriguing question of the nature of the folding intermediate. In fact, by assuming that the cpPDZ2 folding intermediate represents a productive species on path to the native state, as is the case in the wild-type PDZ2 folding, it may be difficult to reconcile the observed stabilisation of the intermediate with a destabilisation of the native state. Lindberg and Oliveberg (2007)
recently proposed pathway malleability to be a consensus feature of protein folding. In simple terms, when and if the denatured chain may reach its native conformation by means of different independent units (foldons), circular permutation involving the cleavage of a folding nucleus present in a dominant foldon may result in a different dominant folding pathway. In line with this hypothesis, it is tempting to speculate that although circular permutation stabilises the early events of PDZ folding (as mirrored by stabilisation of both TS1 and the folding intermediate), the main folding trajectory may be re-routed to different folding pathways. Under such conditions, contrary to what has been described for PDZ2, the observed intermediate may act as an off-pathway kinetic trap, which competes with productive folding.
In summary, we have here outlined the effects of a circular permutation on the folding mechanism of PDZ2. The most intriguing effect is the stabilisation of a folding intermediate that may be either on- or off-pathway. Whether or not the observed folding intermediate represents a stabilised version of the on-pathway intermediate detected in the PDZ2 folding remains to be found by extensive
-value analysis.
| Funding |
|---|
|
|
|---|
Y.I. is supported by a grant from the Wenner-Gren Foundations (Sweden). Work partly supported by grants from the Italian Ministero dell'Università e della Ricerca (2005027330_005) to C.T.A.; University of Rome, Sapienza, project Ateneo 2006 to M. B.
| Footnotes |
|---|
Abbreviations: PDZ, post-synaptic density-95/discs large/zonula occludens-1; PDZ2, the second PDZ domain from murine protein tyrosine phosphatase.
| References |
|---|
|
|
|---|
Bai Y. Proc. Natl Acad. Sci. USA (1999) 96:477–480.
Baker D. Nature (2000) 405:39–42.[CrossRef][Medline]
Buscaglia M., Lapidus L.J., Eaton W.A., Hofrichter J. Biophys. J. (2006) 91:276–288.[CrossRef][Web of Science][Medline]
Capaldi A.P., Shastry M.C., Kleanthous C., Roder H., Radford S.E. Nat. Struct. Biol. (2001) 8:68–72.[CrossRef][Web of Science][Medline]
Chi C.N., Gianni S., Calosci N., Travaglini-Allocatelli C., Engstrom K., Jemth P. FEBS Lett. (2007) 581:1109–1113.[CrossRef][Web of Science][Medline]
Christopher J.A., Baldwin T.O. J. Mol. Biol. (1996) 257:175–187.[CrossRef][Web of Science][Medline]
Fersht A.R. Proc. Natl Acad. Sci. USA (2000) 97:1525–1529.
Friel C.T., Capaldi A.P., Radford S.E. J. Mol. Biol. (2003) 326:293–305.[CrossRef][Web of Science][Medline]
Gianni S., Calosci N., Aelen J.M., Vuister G.W., Brunori M., Travaglini-Allocatelli C. Prot. Eng. Des. Sel. (2005) a18:389–395.[CrossRef]
Gianni S., Engstrom A., Larsson M., Calosci N., Malatesta F., Eklund L., Ngang C.C., Travaglini-Allocatelli C., Jemth P. J. Biol. Chem. (2005) b280:34805–34812.
Gianni S., Walma T., Arcovito A., Calosci N., Bellelli A., Engstrom A., Travaglini-Allocatelli C., Brunori M., Jemth P., Vuister G.W. Structure (2006) 14:1801–1809.[Medline]
Gianni S., Geierhaas C.D., Calosci N., Jemth P., Vuister G.W., Travaglini-Allocatelli C., Vendruscolo M., Brunori M. Proc. Natl Acad. Sci. USA (2007) 104:128–133.
Harris B.Z., Lim W.A. J. Cell Sci. (2001) 114:3219–3231.[Web of Science][Medline]
Ivarsson Y., Travaglini-Allocatelli C., Jemth P., Malatesta F., Brunori M., Gianni S. J. Biol. Chem. (2007) 282:8568–8572.
Jackson S.E., Fersht A.R. Biochemistry (1991) 30:10436–10443.[CrossRef][Web of Science][Medline]
Jemth P., Gianni S. Biochemistry (2007) 46:8701–8708.[CrossRef][Web of Science][Medline]
Jemth P., Gianni S., Day R., Li B., Johnson C.M., Daggett V., Fersht A.R. Proc. Natl Acad. Sci. USA (2004) 101:6450–6455.
Krishna M.M., Englander S.W. Proc. Natl Acad. Sci. USA (2005) 102:1053–1058.
Li L., Shakhnovich E.I. J. Mol. Biol. (2001) 306:121–132.[CrossRef][Web of Science][Medline]
Lindberg M.O., Oliveberg M. Curr. Opin. Struct. Biol. (2007) 17:21–29.[CrossRef][Web of Science][Medline]
Lindberg M., Tangrot J., Oliveberg M. Nat. Struct. Biol. (2002) 9:818–822.[Web of Science][Medline]
Mayor U., Guydosh N.R., Johnson C.M., Grossmann J.G., Sato S., Jas G.S., Freund S.M.V., Alonso D.O.V., Daggett V., Fersht A.R. Nature (2003) 421:863–867.[CrossRef][Medline]
Miller E.J., Fischer K.F., Marqusee S. Proc. Natl Acad. Sci. USA (2002) 99:10359–10363.
Munoz V., Eaton W.A. Proc. Natl Acad. Sci. USA (1999) 96:11311–11316.
Nishimura C., Prytulla S., Dyson J.H., Wright P.E. Nat. Struct. Biol. (2000) 7:679–686.[CrossRef][Web of Science][Medline]
Otzen D.E., Fersht A.R. Biochemistry (1998) 37:8139–8146.[CrossRef][Web of Science][Medline]
Plaxco K.W., Simons K.T., Baker D. J. Mol. Biol. (1998) 277:985–994.[CrossRef][Web of Science][Medline]
Santoro M.M., Bolen D.W. Biochemistry (1988) 27:8063–8068.[CrossRef][Web of Science][Medline]
Sanz J.M., Fersht A.R. FEBS Lett. (1994) 344:216–220.[CrossRef][Web of Science][Medline]
Shastry M.C., Roder H. Nat. Struct. Biol. (1998) 5:385–392.[CrossRef][Web of Science][Medline]
Teilum K., Maki K., Kragelund B.B., Poulsen F.M., Roder H. Proc. Natl Acad. Sci. USA (2002) 99:9807–9812.
Travaglini-Allocatelli C., Gianni S., Morea V., Tramontano A., Soulimane T., Brunori M. J. Biol. Chem. (2003) 278:41136–41140.
Travaglini-Allocatelli C., Gianni S., Brunori M. Trends Biochem. Sci. (2004) 29:535–541.[CrossRef][Web of Science][Medline]
Vendruscolo M., Paci E., Dobson C.M., Karplus M. Nature (2001) 409:641–645.[CrossRef][Medline]
Viguera A.R., Serrano L., Wilmanns M. Nat. Struct. Biol. (1996) 3:874–880.[CrossRef][Web of Science][Medline]
Walma T., Spronk C.A., Tessari M., Aelen J., Schepens J., Hendriks W., Vuister G.W. J. Mol. Biol. (2002) 316:1101–1110.[CrossRef][Web of Science][Medline]
Received November 19, 2007; revised November 19, 2007; accepted November 19, 2007.
![]()
CiteULike
Connotea
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||




