PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 45174.
Published online 2017 March 24. doi:  10.1038/srep45174
PMCID: PMC5364529

Entropic stabilization of a deubiquitinase provides conformational plasticity and slow unfolding kinetics beneficial for functioning on the proteasome

Abstract

Human ubiquitin C-terminal hydrolyase UCH-L5 is a topologically knotted deubiquitinase that is activated upon binding to the proteasome subunit Rpn13. The length of its intrinsically disordered cross-over loop is essential for substrate recognition. Here, we showed that the catalytic domain of UCH-L5 exhibits higher equilibrium folding stability with an unfolding rate on the scale of 10−8 s−1, over four orders of magnitudes slower than its paralogs, namely UCH-L1 and -L3, which have shorter cross-over loops. NMR relaxation dynamics analysis confirmed the intrinsic disorder of the cross-over loop. Hydrogen deuterium exchange analysis further revealed a positive correlation between the length of the cross-over loop and the degree of local fluctuations, despite UCH-L5 being thermodynamically and kinetically more stable than the shorter UCHs. Considering the role of UCH-L5 in removing K48-linked ubiquitin to prevent proteasomal degradation of ubiquitinated substrates, our findings offered mechanistic insights into the evolution of UCH-L5. Compared to its paralogs, it is entropically stabilized to withstand mechanical unfolding by the proteasome while maintaining structural plasticity. It can therefore accommodate a broad range of substrate geometries at the cost of unfavourable entropic loss.

Ubiquitin C-terminal hydrolases (UCHs) are papain-like cysteine proteases that deubiquitinate by hydrolysing the isopeptide bond between the C-terminus of ubiquitin and the lysine side-chain of its adduct, countering ubiquitination in protein homeostasis1,2. Part of the much larger superfamily of deubiqutinases (DUBs), there are four UCHs in the human genome, namely UCH-L1, -L3, -L5 (also known as UCH37) and BAP1 (BRCA1-associated protein 1)3. While UCH-L1 and -L3 are single-domain proteins, UCH-L5 and BAP1 have an additional UCH37-like domain, ULD, at their C-termini. The DUB activity of UCH37 can be activated upon binding to a subunit of 26S proteasome, Rpn134,5,6,7. As a proteasome-associated DUB, UCH-L5 acts at the last checkpoint of the ubiquitin-proteasome pathway (UPS) before the ubiquitinated substrates are degraded by the proteasome. Deubiquitination by UCH-L5 prevents ubiquitinated substrates from being degraded by the UPS and thus, this protein editing process can be likened to a tug of war between DUBs and UPS for maintaining protein homeostasis. Despite its importance, however, little is known about how and when UCH-L5 rescues specific substrates by exerting DUB activity.

While the UCH domain of UCH-L5 shares limited sequence homology with UCH-L1 (23%) and UCH-L3 (28%) compared to the much high sequence homology between UCH-L1 and -L3 (55%), their three-dimensional structures are essentially identical. The few differences in local dynamics result in some disordered secondary structures (Supplementary Results)1,8,9,10,11,12,13,14. In addition to the conserved catalytic triad, C90, H169, and D184 in UCH-L5, the residues lining the ubiquitin binding interface are also highly conserved, forming a deep crater to accommodate the ubiquitin C-terminus. Part of the ubiquitin-binding interface of UCHs is formed by the so-called cross-over loop12,13. The cross-over loop of UCH-L5 is six and three residues longer than that of UCH-L1 and -L3, respectively. The electron density of the cross-over loops and preceding α-helix 6 of UCH-L3 and UCH-L5 in their apo forms could not be resolved due to intrinsic disorder. The flexibility, however, is important for substrate recognition, as UCH-L5 is the only UCH amongst the three that is capable of hydrolysing K48-linked di-ubiquitin. Grafting the cross-over loop of UCH-L5 onto UCH-L1 renders the chimera capable of hydrolysing K48-linked di-ubiquitin, indicating that conformational plasticity of the ubiquitin binding pocket is critical for substrate recognition of UCHs15. Indeed, comparing the reported crystal structures of UCH-L5 in its apo forms and its co-factor bound forms reveals that the UCH-L5 catalytic domain folds into a well-defined structure with limited variation in terms of its overall structures under different states (Supplementary Results), though there exists a high degree of structural plasticity in the cross-over loop region as well as in α-helix 7 where ubiquitin makes direct contacts (Supplementary Results). In contrast, the impacts of ubiquitin binding to UCH-L1 and -L3 are mostly limited to the conformational rearrangements in the cross-over loops and flanking regions, and, in the case of UCH-L1, a cascade of side-chain movements to align catalytic residues into a productive configuration9,13.

Enzyme kinetics analyses of UCH-L1, -L3, and -L5 using a fluorogenic ubiquitin derivative, ubiquitin C-terminal 7-amido-4- methylcoumarin (UbAMC), as a model substrate indicate that UCH-L1 and -L3 form a productive Michelis-Menten complex with ubiquitin (KM value) at a low nM range16,17,18,19, while that of UCH-L5 (a truncated form of residues 1–240 that retains the catalytic domain, hereafter denoted as UCH-L5N240)18 is three orders of magnitude higher, which may be attributed to the unfavourable entropic loss upon ubiquitin binding. Nonetheless, UCH-L5 exhibits a robust turnover rate (kcat) for which UbAMC that is over three orders of magnitude higher than that of UCH-L1 and slightly higher than that of UCH-L31. Considering that UCH-L5 is associated with the proteasome, where the local concentration of ubiquitin is much higher than the average cellular concentration, UCH-L5 might effectively deubiquitinate its substrates for editing purposes before the substrates are degraded by the proteasome.

In addition to their potential role in protein homeostasis, UCHs contain a complex Gordian knot with five projected crossings in their folding topologies, rendering them an intriguing system for protein folding studies20,21. Knotted proteins can be found in all kingdoms of life and many knotted proteins are enzymes. The question of how knotted proteins attain their intricate knotted topologies have attracted the attention of experimental and computational biophysicists in recent years22. Remarkably, two bacterial RNA methyltransferases, YibK and YbeA, have been shown to knot spontaneously immediately after protein synthesis during de novo folding without the aid of chaperonins23, and they remain knotted in the presence of highly concentrated chemical denaturant24 while behaving like random coils25. The knotted structural element of bacterial RNA methyltransferases has been shown to be of functional importance in substrate recognition26. Indeed, the cross-over loop of UCHs is also essential in ubiquitin recognition despite the variations in loop length and sequence composition.

The folding dynamics and kinetics of UCH-L1 and -L3 were recently delineated27,28. Both UCHs exhibit two well-defined kinetic folding intermediates connected by two parallel folding pathways between the native and denatured states. We further demonstrated that UCH-L1 exhibits partially unfolded forms (PUFs) under native conditions that share structural characteristics with those of the chemically-induced intermediate states28. As part of our endeavour to systematically characterize the folding of knotted proteins28,29,30,31,32, particularly those within the same family, we applied the same experimental strategy for studying UCH-L1 and -L3 to investigate in detail the folding dynamics and kinetics of the UCH domain of UCH-L5. Our folding analyses indicate that UCH-L5 is thermodynamically more stable than UCH-L1 and -L3, with unfolding kinetics over four orders of magnitude slower than those of UCH-L1 and -L3. Despite the slow global unfolding kinetics, however, UCH-L5 displays abundant local fluctuations within its secondary structure elements manifested in rapid hydrogen-deuterium exchange (HDX). Our results therefore represent a rare example of an entropically stabilized biological system that requires a significant degree of structural plasticity for substrate recognition and accommodation in a hostile environment that is constantly subject to mechanical unfolding by the proteasome for protein degradation.

Results

Folding equilibrium of UCH-L5N240 monitored by intrinsic fluorescence and far-UV CD spectroscopy

To compare the folding stability of UCH-L5N240 with that of UCH-L133 and -L327, we analysed its chemical stability by urea-induced equilibrium unfolding analysis through intrinsic fluorescence and far-UV CD spectroscopy. UCH-L5N240 has four tryptophan residues that are strategically distributed in the backbone topology that serve as ideal structural probes (Fig. 1). W8 and W36 are in proximity to the N- and C-termini that are involved in knot formation. W58 is adjacent to the catalytic site, encompassing C88, H164, and D179. W196 is located at the turn between β5 and the longest helix, α7. SVD analyses of the intrinsic fluorescence spectra of UCH-L5N240 as a function of urea concentration indicate that UCH-L5N240 unfolds in a two-state fashion without significant contribution from folding intermediates as in the case for UCH-L1 (Fig. 2). This is further confirmed by far-UV CD spectroscopy that monitors changes in secondary structures. The intrinsic fluorescence and far-UV CD unfolding data were global-fit to a two-state model (Table 1). The free energy of unfolding of UCH-L5N240 (9.21 ± 0.39 kcal mol−1) is higher than that of UCH-L1 (8.83 ± 0.10 kcal mol−1)33 and UCH-L3 (7.11 ± 0.20 kcal mol−1)27, and the transition point of UCH-L5N240, [D]N-D,50%, is about 0.5 M higher than those of UCH-L1 and -L3 with comparable m-values (Table 1).

Figure 1
Three-dimensional structure of UCH-L5N240 and its knotted folding topology.
Figure 2
Urea-induced equilibrium unfolding of UCH-L5N240.
Table 1
Thermodynamic parameters of equilibrium unfolding of UCH-L1, -L3 and -L5N240 induced by urea.

Folding kinetics of UCH-L5N240

We next evaluated the folding kinetics of UCH-L5N240 in comparison with those of UCH-L1 and -L327,28. Single-jump (SJ) stopped-flow fluorescence measurements of UCH-L5N240 as a function of urea concentration yielded four linear refolding arms that we assigned as phases 1, 2, 3, and 4, corresponding to the slowest to the fastest kinetic phases, respectively. The slowest refolding phase 1 has the largest amplitude (A1), indicating its association with global folding events (Fig. 3). In the case of unfolding kinetics, we only observed three linear unfolding arms that correspond to the kinetic phases 2–4 because the kinetics of phase 1 is too slow to be reliably measured by the stopped-flow mixing device. We therefore used manual mixing followed by far-UV CD measurements (Supplementary Results) to extract two unfolding kinetic phases, with the faster one being consistent with the kinetic phase 2 up to 4.4 M urea; above 4.4 M urea, the kinetic CD traces were fit to a single exponential function due to the limited time resolution to resolve the contribution of the faster kinetic phase. The slower kinetic phase observed in the CD measurements was then associated with unfolding phase 1 (open black circles in Fig. 3a).

Figure 3
Chevron plot analysis of UCH-L5N240.

By fitting to a two-state model, we obtained four refolding rates and four unfolding rates associated with the four kinetic phases. The fastest and slowest refolding rates were separated by four orders of magnitude, and the fastest and slowest unfolding rates were separated by six orders of magnitude, reflecting a very large dynamic range for the different folding events (Table 2). Amongst the four kinetic phases, the kinetic m-value (mkin) and the associated free energy of unfolding of the phase 1 are essentially the same as those derived from equilibrium unfolding (Table 1), while the values for the other three kinetic phases are markedly smaller, suggesting that phase 1 is the dominant term that gives rise to the observed spectral changes in the equilibrium unfolding measurements. In other words, phase 1 reflects the unfolding of the native state (N) into one of the intermediates. Indeed, CD-derived A1 associated with unfolding phase 1 is much larger than those of the second slowest phase (Supplementary Results); likewise, SJ fluorescence-derived A2 is 4–10 times larger than A3 and A4 (Fig. 3b).

Table 2
Kinetic parameters of UCH-L5N240 derived from stopped-flow fluorescence measurements.

To compare the folding kinetics of UCH-L1, -L3, and -L5N240, we repeated the SJ stopped-flow fluorescence measurements for UCH-L1 and -L3 under the same experimental conditions and focused on the slowest and the most dominant kinetic phases to compare with that of UCH-L5N240 (Fig. 3c and Table 3). While the m-values associated with the refolding arms of UCHs are similar, the m-value associated with the unfolding arm of UCH-L5N240 is about twice that of UCH-L1 and -L3, resulting in a much slower unfolding rate in water (An external file that holds a picture, illustration, etc.
Object name is srep45174-m1.jpg = 8.1 * 10−9 s−1) that is five and four orders of magnitude slower than those of UCH-L1 and -L3. The folding rate of UCH-L5N240 in water (An external file that holds a picture, illustration, etc.
Object name is srep45174-m2.jpg = 0.14 s−1) is comparable with that of UCHL1 and is about 100 times slower than that of UCH-L3.

Table 3
Comparison of the folding kinetics of UCH-L1, -L3 and -L5N240.

In an effort to establish the relationships between the four kinetic phases in the context of folding pathways, we used double-jump (DJ) interrupted refolding fluorescence analyses to obtain additional refolding rates through global fitting of the kinetic traces obtained with different aging times (tage) in the presence of different urea concentrations (open blue circles for k3 and open green triangles for k4 in Fig. 3a). Unlike UCH-L1 and -L3, however, amplitude analyses of the kinetic phases 3 and 4 of UCH-L5N240 as a function of the aging time (tage) can be best fit to a linear three-state reaction with the microscopic forward and back reaction rates extracted through global fitting, although the results are noisy and a lag phase in the amplitude build up of the second fastest phase is barely visible. The results therefore suggest the presence of transiently populated intermediate corresponds to the fastest kinetic phase of UCH-L5N240 (Supplementary Results). Completion of the kinetic folding pathway of UCH-L5N240 for direct comparison with those of UCH-L1 and -L3 will require better instrument stability for kinetic measurements with a timescale that spans over five orders of magnitudes (Fig. 3a).

NMR HDX reveals abundant local backbone fluctuations in UCH-L5N240

Having established the slow global folding kinetics of UCH-L5N240, we sought to obtain further structural insights into the origin of stability by solution state NMR spectroscopy. We first obtained near complete assignments of the fingerprint backbone 15N-1H HSQC spectrum of UCH-L5N240, which display highly dispersed cross-peaks corresponding to a well-folded structure, while a group of poorly dispersed cross-peaks with very strong intensities are ascribed to the highly disordered cross-over loop (residues 150–160) and the C-terminal extension beyond the UCH domain (residues 230–240) for which we did not have complete assignments (Supplementary Results). The structural disorder in the cross-over loop, the flanking residues, and the C-terminal region was confirmed by their low order parameters (S2) reflecting fast local motions on the ps to ns timescale (Supplementary Results). This timescale however, is more than six orders of magnitudes faster than the fluorescence-based kinetics of UCH-L5N240. We therefore resorted to NMR hydrogen-deuterium exchange (HDX) analysis to probe molecular dynamics on a much slower timescale by determining the rate of HDX for a specific amide group and comparing it with the expected value for a random coil to derive the corresponding PF that reflects the stability of the corresponding hydrogen bond (see Material and Methods)34.

Contrary to expectations of a slower HDX of UCH-L5N240 compared to UCH-L1 and -L3 based on the fluorescence-based equilibrium and kinetic data (Tables 1 and and2),2), UCH-L5N240 exhibits significantly lower PF values, i.e., faster HDX, throughout its backbone amide groups compared to those of UCH-L1 and -L3 (Fig. 4). For instance, α-helix 3 of UCH-L1, which is packed against the central β-sheet as part of the hydrophobic core, is highly protected against HDX, and the C-terminal half of α-helix 3 of UCH-L3 exhibits partial protection against HDX; in contrast, α-helix 3 of UCH-L5N240 shows negligible protection against HDX. Likewise, the PFs of the central β-sheet are higher for UCH-L1 than UCH-L3, both of which are higher than the minimal HDX protections of UCH-L5N240. In line with the crystallographic findings that indicated highly disordered structures in α-helix 6 for UCH-L3 and α-helix 2 for UCH-L5N240, we did not observe appreciable PF for any of the residues within these regions (shown in dashed lines in the top panel of Fig. 4). Consistent with our recent finding of the PUFs of UCH-L1 under native conditions, our current data on UCH-L3 and UCH-L5N240 suggest the presence of highly populated PUFs as well: the degree of fluctuations for UCH-L5N240 is much higher than that of UCH-L3, and even more so than that of UCH-L1, indicating a high degree of conformational plasticity for UCH-L5N240 (Supplementary Results).

Figure 4
NMR HDX analyses of UCH-L1, UCH-L3 and UCH-L5N240.

Thermodynamics associated with ubiquitin binding to UCH variants

Earlier studies on the enzyme kinetics of UCH-L1, -L3 and -L5N240 have established that both UCHs bind to ubiquitin with sub μM KM values16,17,18,19,35,36. While UCH-L1 exhibits higher binding affinity for ubiquitin compared to that of UCH-L3, its hydrolysis turnover rate (kcat) is about two orders of magnitude slower than that of UCH-L3. We repeated and confirmed the previously reported parameters for UCH-L1, -L3, and -L5N240 (Supplementary Results). Although the reported values have a rather broad range of distributions and our kcat value is at the higher end, the overall kcat/KM ratios are similar across different reports.

Having established that our recombinant UCHs are enzymatically active, we subsequently carried out isothermal titration calorimetry at four different temperatures (20, 25, 30, and 37 °C) to gain further insights into the ubiquitin binding process (Supplementary Results). For UCH-L1 and -L3, their binding affinities are not sensitive to temperature changes but the enthalpic contributions are clearly more favourable upon increasing temperature. For UCH-L1, the entropic contributions are favourable up to 30 °C and become unfavourable at body temperature, 37 °C. In contrast, for UCH-L3, the entropic contributions become unfavourable at 25 °C. ITC analysis is only applicable to UCH-L1 and -L3 due to the poor ubiquitin binding affinity of UCH-L5; the binding affinity between UCH-L5 and ubiquitin could only be determined by ITC in the presence of Rpn136. Indeed, the addition of 1.2 molar equivalent of ubiquitin slightly increased the thermal stability of UCH-L1 and -L3, raising the melting temperature by 2.2 ± 0.1 °C and 0.9 ± 0.2 °C, respectively, while no significant change in thermal stability was observed for UCH-L5N240 (Supplementary Results). The lack of significant thermal stability shift in UCH-L5N240 upon the addition of its substrate is consistent with the expectation for weaker binding events37. As the CD melting analyses were carried out with only 5 μM protein concentration, the addition of equal molar concentration of ubiquitin may not be sufficient to saturate UCH-L5N240. Even after the addition of six-fold higher amount of ubiquitin (30 μM), we did not observe appreciable changes in melting temperature (data not shown). We next carried out NMR titration analyses with 0.1 mM UCH-L5N240 and up to 0.15 mM ubiquitin to evaluate the binding processes under conditions that should significantly increase the population of ubiquitin-bound UCH-L5N240. The end NMR titration points nonetheless showed severe line-broadening of UCH-L5N240 backbone amide and side-chain methyl resonances as opposed to the well-defined backbone and side-chain resonances of ubiquitin-bound UCH-L138, suggesting fast to intermediate exchange processes and therefore unfavourable line-broadening (Supplementary Results).

Despite our inability to obtain a stable binary complex of UCH-L5N240 and ubiquitin in solution, we focused on the global analysis of the variable temperature ITC data to extract the heat capacity ΔCp in addition to the enthalpy, entropy, and thus free energy associated with ubiquitin binding for the UCH-L1 and -L3 in order to extrapolate the findings to correlate with potential characteristics of ubiquitin binding events for UCH-L5N240 (Table 4, Supplementary Results). Under standard conditions, 1 atm and 25 °C, UCH-L1 exhibits the strongest binding affinity towards ubiquitin, with a dissociation constant (KD) of 0.26 μM, while that of UCH-L3 is significantly weaker (Table 4). While UCH-L1 and -L3 bind to ubiquitin with comparable free energy change, the binding event is entropically favourable at 25 °C in the case of UCH-L1 and entropically unfavourable for UCH-L3, resulting in marked entropy-enthalpy compensation (Supplementary Results). As a result, the changes in heat capacity (ΔCp) associated with ubiquitin binding is significantly larger for UCH-L3 (−703 cal mol−1 K−1) compared to that of UCH-L1 (−525 cal mol−1 K−1). ΔCp is associated with the changes in solvent accessible surface area (ΔSASA) upon binding. Indeed, UCH-L3 has a larger ΔSASA value (2691 Å2) compared to that of UCH-L1 (2372 Å2), according to their respective crystal structures in complex with ubiquitin (PDB ID: 1XD3 and 3KW5, respectively). This is accompanied by a large conformational change and cross-over loop ordering in UCH-L3 manifested in the large entropic penalty associated with ubiquitin binding. It is plausible to speculate that for UCH-L5N240, the entropic penalty associated with mono-ubiquitin binding is far greater that that for UCH-L3, hence the much weaker binding that is not accessible to ITC analysis.

Table 4
Thermodynamic parameters of ubiquitin binding to UCH-L1 and –L3.

Discussion

In this work, we have characterized in great detail the folding equilibrium and kinetics of the UCH domain of UCH-L5, i.e., UCH-L5N240, in comparison with previously reported results on UCH-L1 and -L3. Despite the apparent two-state folding equilibrium behaviour (Fig. 2), four well-resolved kinetic phases induced by urea could be observed by stopped-flow fluorescence measurements, indicating the presence of folding intermediates (Fig. 3). While the four kinetic folding phases have also been observed for UCH-L1 and -L3, and two parallel folding pathways were proposed to connect the native and denatured states, with two distinct folding intermediates being transiently populated along the two parallel pathways27,28, our interrupted refolding data suggest that the two fastest kinetic phases (k3 and k4) are connected sequentially rather than separately with two slower kinetic phases (k1 and k4) (Supplementary Results). Subsequent assignments of the connection(s) between the two fastest and two slowest phases were hindered because our stopped-flow instrumentation could not follow the large timescale separating the kinetic phases. The current data nevertheless suggest that the kinetic folding pathways of UCH-L5N240 are different from those of UCH-L1 and -L3.

Human UCHs are important DUBs for maintaining cellular proteostasis2. Their unique domain architecture – a single domain that binds to and hydrolyses the isopeptide bond of ubiqutinated substrates – with very slow off-rates for releasing hydrolysed ubiquitin products39 potentially increase the likelihood of encountering the proteasome machinery that actively unfolds ubiquitinated substrates for degradation. Indeed, it has been suggested by the recent single molecule mechanical unfolding study of UCH-L1 that the mechanical stability of the UCH domain may be beneficial for resisting UPS degradation40. Another important finding from the single molecule study is the consistent timescale of refolding triggered from mechanically and chemically unfolded states, implying that the intrinsic unfolding rates in the absence of chemical denaturant can be reliably derived from extrapolation of the linear dependency of chemical denaturant concentration (Fig. 3).

Our comparative kinetic analysis of UCH-L1, -L3 and -L5N240 highlighted the large dynamic range of the intrinsic unfolding rates (An external file that holds a picture, illustration, etc.
Object name is srep45174-m3.jpg) of the most dominant folding events, i.e., unfolding of the native state to one of the intermediate states, amongst the three UCHs. UCH-L5N240 has the longest cross-over loop and exhibits the slowest global unfolding rate, four orders of magnitude slower than those of UCH-L1 and -L3 (Fig. 3 and Table 3). The length of the cross-over loop is essential in substrate recognition for UCHs. The longer cross-over loop of UCH-L5 is necessary for binding and hydrolysing the isopeptide bond of K48-linked di-ubiquitin15. According to the crystal structure of UCH-L5 in complex with ubiquitin and Rpn136, UCH-L5 binds to mono-ubiquitin in the same pose as other UCHs, but with partially resolved electron density for the cross-over loop (Fig. 5a). However, the cross-over loop and the structural elements around the ubiquitin binding, site such as α-helix 7, need to undergo significant configuration rearrangement in order to accommodate a compact K48-linked di-ubiquitin without steric clashes between the cross-over loop and the substrate (Fig. 5b and c). We deduced an intrinsic structural plasticity from the NMR HDX analysis (Fig. 4), presenting the solution to this issue. Our findings also echo the conformational plasticity implicated in the regulation of the UCH-L5 DUB activity underscored by the recent crystallographic studies of UCH-L5 in complex with different auxiliary factors such as Rpn13 and INO80G for up- and down-regulating the DUB activity of UCH-L5, respectively6,7.

Figure 5
Conformational plasticity of UCH-L5 is required for hydrolysing isopeptide bond of K48-linked di-ubiquitin.

Generally speaking, longer loop lengths result in higher contact orders and thus slower folding rates41. However, the differences in the contact orders between UCH-L1, -L3, and -L5N240 are too small to account for the four orders of magnitudes difference in refolding and unfolding rates; the relative contact orders of UCH-L1, -L3, and -L5N240 are 0.137, 0.144, and 0.139, respectively (Fig. 3c). The differences in sequence compositions and the folding stabilities of individual secondary structure elements within different UCHs may have more significant contributions to the folding dynamics and kinetics than the contact orders per se. Indeed, despite the abundant internal dynamics within UCH-L5N240 in the absence of ubiquitin, it is highly stable thermodynamically and kinetically in terms of chemical stability compared to UCH-L1 and -L3 (Tables 1 and and3).3). The stabilization effect may be attributed to the higher configurational entropy of UCH-L5N240 with respect to that of UCH-L1 and -L3. Indeed, our NMR HDX analysis revealed that the individual secondary structural elements within UCH-L5N240 are much less protected from solvent exchange, indicative of abundant local fluctuations (Fig. 4). Rapid HDX in secondary structure elements (mostly peripheral helices) of UCH-L1 has been attributed to highly populated PUFs under native conditions28. It is therefore plausible that UCH-L5N240 also exhibits abundant PUFs under native conditions in solution, which help stabilize the folding of UCH-L5N240.

The extent to which local fluctuations occur within the PUFs of UCH-L5 may be far greater than what was observed in the crystal structure of UCH-L5, which shows structural disorder limited to the cross-over loop in all reported crystal structures and α-helix 2 in the apo form (Fig. 1)6,7,11. Similar loop dynamics was observed in UCH-L5 from Trichinella spiralis42. Indeed, cross-over loop flexibility was observed in apo UCH-L3, but it was fully resolved in the ubiquitin-form13. In contrast, the cross-over loop of UCH-L1 was resolved both in the apo- and ubiquitin-bound forms9,10. The entropic cost of structural ordering in UCH-L3 is confirmed in our ITC analysis (Table 4 and Supplementary Results). There is a strong correlation between the free energies of unfolding of UCH-L1, -L3 and -L5N240 with the lengths of the cross-over loops (Table 3). In addition, the degree of local fluctuations deduced from NMR HDX analysis is proportional to the length of the cross-over loop (Fig. 4).

Using model systems, the effects of loop insertions have been investigated, generally through inserting repeating glycine residues43,44,45. Based on the findings, an empirical function was proposed to probe residual structures in the denatured states with the assumptions that the inserted loop residues marginally destabilize the proteins, that the insertions do not interact with the remaining structure in the denatured, transition, or native states, and that the energetic changes are predominantly entropic46. Indeed, such assumptions were met for several other model systems47,48. Loop elongation however, is not always destabilizing for proteins. A recent study showed that polyglycine loop insertion could in fact significantly increase the entropy of the native state of an acylphosphatase, thereby increasing the native state stability without perturbing the transition and denatured states49. In contrast to the artificial manipulation of the loop lengths of model systems, our current results present a biologically relevant example, which has likely evolved to maximize the structural plasticity of UCH-L5 by a longer substrate recognition loop and peripheral helical elements. The structural plasticity therefore contributes to an entropic gain in thermodynamics while the higher flexibility is beneficial for recognizing a broad range of substrates with potentially different ubiquitin linkages without steric clashes with the cross-over loop (Fig. 5). Our NMR titration data (Supplementary Results) strongly suggest that UCH-L5N240 remains dynamic in complex with ubiquitin, which may be the underlying reason for its poor KM value (Supplementary Results). The entropically stabilized native state of UCH-L5N240 significantly reduces the unfolding rate that could be essential for withstanding the potential mechanical unfolding driven by the proteasome when it is associated with the regulatory subunit Rpn13 for DUB activity activation. It is therefore plausible that the functional contribution of the cross-over loops of UCHs that are part of the knotted elements is to balance the need for structural plasticity in substrate regulation and kinetic folding stability that is needed for withstanding proteasome-mediated degradation.

Methods

Recombinant protein purification

The pGEX-6P1 vector-based expression construct of UCH-L5N240 was a kind gift of Dr. Chittaranjan Das (Purdue University, USA). The resulting recombinant protein construct contains a glutathione S-transferase (GST) fused at the N-terminus of UCH-L5N240. The expression vector was transformed into a BL21 (DE3) E. coli strain and cultured in LB medium containing 0.1 mg/mL ampicillin for antibiotic selection. Recombinant protein over-expression was induced by adding 0.5 mM isopropyl thio-β-d-galactoside (IPTG) when the optical density of bacterial culture at 600 nm (OD600) was close to 0.8, followed by overnight incubation with shaking at 16 °C. The cells were harvested and purified by a standard protocol for GST-based affinity column purification, followed by fusion tag removal using PreScission protease (GE Biosciences, USA), as described previously11. The resulting recombinant protein was subjected to size-exclusion chromatography using a gel-filtration column (Superdex 75 26/60, GE Biosciences, USA) in 20 mM Tris-HCl (pH 7.6), 100 mM NaCl and 5 mM DTT. The purity of the recombinant protein was higher than 95% judging from SDS-PAGE by visual inspection.

Intrinsic fluorescence spectroscopy

The intrinsic fluorescence of UCH-L5N240 was monitored by exciting the samples at 260 nm and recording the emission spectra between 300 and 450 nm at 25 °C using a temperature-controlled fluorimeter (FP8500, JASCO, Japan). Urea-induced equilibrium unfolding was monitored by generating 41 aliquots of protein solution containing a gradient of urea (0–6 M) concentrations with a linear increment step of 2.5% made by a two-channel liquid dispenser (Hamilton, USA) to minimize manual pipetting errors.

Far-UV CD spectroscopy

Urea-induced equilibrium unfolding was also monitored by far-UV CD between 200 and 260 nm at 25 °C using a temperature-controlled CD spectrometer (J-815, JASCO, Japan). The spectra were collected with a bandwidth of 1 nm, a data interval of 0.5 nm, and an averaging time of one second. Thermal denaturation of UCH-L5N240 was monitored by analysing the CD signal changes at 215 nm between 25 and 80 °C while the changes at 222 nm were monitored for UCH-L1 and -L3. In all cases, 5 μM of protein was used and 1.2 molar equivalent of ubiquitin was added to assess the effect of ubiquitin binding. The resulting isotherms were normalized and fitted to a two-state equilibrium unfolding model to extract the thermodynamics parameters, i.e., melting temperature (Tm), enthalpy (ΔH), entropy (ΔS) and changes in heat capacity (ΔCp) as described previously50.

Thermodynamic parameters of equilibrium unfolding of UCH-L5N240 induced by urea

All equilibrium unfolding data were subjected to singular value decomposition (SVD) analysis using an in-house written script for MATLAB (MATLAB and Statistics Toolbox Release 2012b, The MathWorks, USA) as described previously29,30,31. The titration series were used to generate an m × n matrix, M, as inputs for SVD analysis, where m corresponds to the number of recorded wavelengths and n corresponds to the number of titration points. SVD reconstruction generates an m × m unitary rotation matrix Ur, an n × n unitary rotation matrix Vr, and an m × n diagonal singular value matrix S to reconstruct the input matrix M as follows

An external file that holds a picture, illustration, etc.
Object name is srep45174-m4.jpg

where Ur corresponds to the basis function of the observed spectral patterns, and Vr corresponds to the basis coefficient as a function of denaturant concentration. To determine the number of significant components necessary for reconstructing the input matrix, the normalized correlation coefficients of individual singular values were calculated and a threshold of 0.8 was used, leading to the conclusion that UCH-L5N240 unfolds in a two-state model.

The fluorescence and CD data were subjected to a two-state equilibrium-unfolding model as described previously28,29,30,31. The intrinsic fluorescence and far-UV CD data were fit separately and globally. Additionally, the changes of maximum fluorescence emission wavelength as a function of denaturant concentration were fitted to the same model for comparison.

Stopped-flow fluorescence spectroscopy

The folding kinetics of UCH-L5N240 was monitored using a stopped-flow spectrometer in fluorescence detection mode (SX18 stopped-flow spectrometer, Applied Photophysics, UK) as described previously28,29,30,31. Changes in total fluorescence were monitored using an excitation wavelength of 280 nm with a 320 nm cut-off filter. All experiments were carried out at 25 °C. 20 μM UCH-L5N240 was buffered in 20 mM Tris-HCl (pH 7.6), 100 mM NaCl with 8 M urea for unfolding or without urea for refolding measurements. For single-jump kinetic measurements, the folding reactions were triggered by rapidly mixing the folding or unfolding buffer with protein solution at an asymmetric mixing ratio of 10:1. After the 11-fold dilution, the final protein concentration was 1.8 μM. For double-jump kinetic measurements, three different modes were employed: interrupted refolding, interrupted unfolding, or sequential refolding kinetic measurement.

For the interrupted refolding measurements, 7.7 M urea-denatured protein solution was first mixed with refolding buffer to initiate refolding over different aging periods, followed by the second mixing step with denaturing buffer containing 8 M urea immediately before we took kinetic fluorescence measurements under denaturing conditions. The first mixing step was achieved by mixing protein solutions with unfolding or refolding buffer with a 1:5 mixing ratio. After a specific aging period, a second mixing step was carried out by mixing the aged protein solution with an equal amount of unfolding or refolding buffer (1:1 mixing ratio) immediately before kinetic fluorescence measurements. The resulting kinetic traces were subjected to global fitting in which the folding rates were shared amongst all the kinetic traces of different aging times, and the amplitudes corresponding to different kinetic phases were shown as a function of aging time, and the rates at which the corresponding amplitudes built up were extracted by single exponential fitting to the observed amplitude as a function of aging time (tage).

For all kinetic analyses, the observed reaction rates were extracted by fitting the kinetic traces to a single, double, or triple exponential function with an offset using the software package GraphPad Prism (GraphPad Software, USA). The choice of model, i.e., number of kinetic phases, was decided using the F-test statistics by Prism.

Chevron plot analysis

The observed reaction rates in the stopped-flow fluorescence measurements were fit to a simple two-state folding model to extract the associated kinetic and thermodynamic parameters. The observed reaction rates (kobs), with a linear refolding arm and a linear unfolding arm, were fit to the following equation:

An external file that holds a picture, illustration, etc.
Object name is srep45174-m5.jpg

where An external file that holds a picture, illustration, etc.
Object name is srep45174-m6.jpg and An external file that holds a picture, illustration, etc.
Object name is srep45174-m7.jpg are the folding and unfolding rate in the water, i.e., in the absence of denaturant, mf and mu are the m-values associated with folding and unfolding, R is the gas constant, and T is the sample temperature, which is set to 298 K. These parameters were subsequently used to calculate the free energy of unfolding associated with this reaction as follows:

An external file that holds a picture, illustration, etc.
Object name is srep45174-m8.jpg

NMR spectroscopy

Uniformly 13C- and 15N-labelled UCH-L5N240 was obtained by growing the transformed E. coli cells in M9 minimal medium containing 15N-labelled ammonium chloride (1 g/mL) and 13C-labelled D-glucose (2 g/mL) followed by the same purification procedure as described above. For backbone resonance assignments, a suite of triple resonance experiments were collected, including 15N-1H heteronuclear single quantum correlation (HSQC) spectroscopy, HNCO, HNCA, HNcoCA, HNCACB, and CBCAcoNH at 25 °C. All NMR data were collected using a Bruker AVIII NMR spectrometer, equipped with a cryogen-cooled TCI probe head, operating at 14.1 or 20.0 Tesla, corresponding to a proton Larmor frequency of 600 or 850 MHz, respectively. The data were collected using Topspin 3.0 (Bruker, Germany), processed by NMRPipe51, and visualized and analyzed by Sparky (https://www.cgl.ucsf.edu/home/sparky/). Near complete backbone (H, N, C, Cα and Cβ) assignments were achieved through the use of an iterative computer-aided procedure as described previously52,53. The backbone assignments have been deposited to the BMRB database under the accession number of 12009.

The ubiquitin titration experiment was carried out using 0.1 mM of uniformly 13C- and 15N-labelled UCH-L5N240 with unlabelled ubiquitin (5 mM stock solution concentration) in a 3 mm MATCH tube at 25 °C. and 20.0 Tesla. Two-dimensional 15N-1H SOFAST-HMQC and 13C-1H HSQC spectra were collected for individual titration points (0, 0.125, 0.25, 0.375, 0.5, 0.75, 1.0, 1.25, and 1.5 ubiquitin to UCH-L5N240 molar ratios) as described previously54,55.

The 15N spin relaxation dynamics of UCH-L5N240 was characterized by measuring the longitudinal (R1) and transverse (R2) relaxation rates and 15N-1H heteronuclear nuclear Overhauser effects (hetNOEs) of individual backbone amide groups at 37 °C and 14.1 Tesla as described previously56. The resulting relaxation rates and hetNOEs were used as inputs to determine the residue-specific order parameters (S2) using model-free analysis by TENSOR257. According to TENSOR2, the use of anisotropic diffusion model did not yield statistically significant improvement in the data analysis. We there used a simple isotropic model to analyse NMR relaxation data.

Hydrogen-deuterium exchange (HDX)

NMR HDX analyses of UCH-L1, -L3, and -L5N240 were carried out in 20 mM Tris-HCl (pH 7.6), 100 mM NaCl at 25 °C, as described previously28,33,58. The previously reported backbone assignments of UCH-L3 (BMRB accession number 15121)59 were used to confirm the assignments of the 15N-1H correlations. The amide 15N-1H cross-peak intensity as a function of HDX time was fit to a single exponential function using the built-in relaxation analysis module of the software package Sparky, and the rate constants were subsequently used to derive the protection factors (PFs) together with the protein sequences using the Excel spreadsheet with pre-defined functions from the Englander laboratory at University of Pennsylvania, USA (http://hx2.med.upenn.edu/download.html).

Isothermal titration calorimetry (ITC)

A Microcal 200 instrument (GE, USA) was used for the ITC measurements. The protein concentrations of UCHs and human ubiquitin were set to 27.5 and 275 μM, respectively. The latter was used as the titrant to be injected into the UCH solutions. The ITC experiments were carried out at 20, 25, 30, and 37 °C, and the resulting ITC data were individually analysed using NITPIC60 and subsequently exported to Sedphat61 for global fitting to extract the associated thermodynamics parameters.

DUB activity analysis

The enzyme kinetics of UCH variants were assessed by UbAMC as a model substrate as described previously16,17,18,19. The enzyme concentrations were set to 1 nM, and UbAMC substrate concentrations ranged from 1 to 1000 nM. The fluorescence as a function of time was recorded individually using a temperature-controlled fluorimeter (FP8500, JASCO, Japan). The results were analysed using the built-in Michaelis-Menten equation of the software package GraphPad Prism (GraphPad Software, USA).

Additional Information

How to cite this article: Lee, Y.-T. C. et al. Entropic stabilization of a deubiquitinase provides conformational plasticity and slow unfolding kinetics beneficial for functioning on the proteasome. Sci. Rep. 7, 45174; doi: 10.1038/srep45174 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supporting Information:

Acknowledgments

We could like to thank Jose Rizo-Rey (UT Southwestern Medical Center), Kuenphon Wu (St. Jude Children’s Hospital), Nunilo Cremades (University of Zaragoza), Yuan-Chao Lou (Academia Sinica), and Shu-Chuan (Chris) Jao (Academia Sinica) for their critical reading and comments for our manuscript. This work was supported by funding from the Ministry of Science and Technology (MOST104-2113-M-001-016) and Academia Sinica, Taiwan. S.T.D.H. was supported by a Career Development Award (CDA-00025/2010-C) from the International Human Frontier Science Program. The NMR data were collected at the High-field NMR Center of Academia Sinica, Taiwan.

Footnotes

The authors declare no competing financial interests.

Author Contributions Y.T.C.L. and S.T.D.H. designed the experiments. Y.T.C.L., C.Y.C., S.Y.C., and Y.R.P. performed the experiments. Y.T.C.L., M.R.H. and S.T.D.H. analysed the data. S.T.D.H. wrote the manuscript.

References

  • Hsu S.-T. D. Folding dynamics and structural basis of the enzyme mechanism of ubiquitin C-terminal hydroylases. In Understanding enzymes, Function, Design, Engineering and Analysis (ed. Svendsen A. A.) 167–202 (Pan Stanford Publishing, Singapore, 2016).
  • Reyes-Turcu F. E., Ventii K. H. & Wilkinson K. D. Regulation and cellular roles of ubiquitin-specific deubiquitinating enzymes. Annu. Rev. Biochem. 78, 363–97 (2009). [PMC free article] [PubMed]
  • Nijman S. M. et al. . A genomic and functional inventory of deubiquitinating enzymes. Cell 123, 773–86 (2005). [PubMed]
  • Hamazaki J. et al. . A novel proteasome interacting protein recruits the deubiquitinating enzyme UCH37 to 26S proteasomes. EMBO J. 25, 4524–36 (2006). [PubMed]
  • Qiu X. B. et al. . hRpn13/ADRM1/GP110 is a novel proteasome subunit that binds the deubiquitinating enzyme, UCH37. EMBO J. 25, 5742–53 (2006). [PubMed]
  • Sahtoe D. D. et al. . Mechanism of UCH-L5 activation and inhibition by DEUBAD domains in RPN13 and INO80G. Mol. Cell 57, 887–900 (2015). [PMC free article] [PubMed]
  • VanderLinden R. T. et al. . Structural basis for the activation and inhibition of the UCH37 deubiquitylase. Mol. Cell 57, 901–11 (2015). [PMC free article] [PubMed]
  • Johnston S. C., Larsen C. N., Cook W. J., Wilkinson K. D. & Hill C. P. Crystal structure of a deubiquitinating enzyme (human UCH-L3) at 1.8 Å resolution. EMBO J. 16, 3787–96 (1997). [PubMed]
  • Boudreaux D. A., Maiti T. K., Davies C. W. & Das C. Ubiquitin vinyl methyl ester binding orients the misaligned active site of the ubiquitin hydrolase UCHL1 into productive conformation. Proc. Natl. Acad. Sci. USA 107, 9117–22 (2010). [PubMed]
  • Das C. et al. . Structural basis for conformational plasticity of the Parkinson’s disease-associated ubiquitin hydrolase UCH-L1. Proc. Natl. Acad. Sci. USA 103, 4675–80 (2006). [PubMed]
  • Maiti T. K. et al. . Crystal structure of the catalytic domain of UCHL5, a proteasome-associated human deubiquitinating enzyme, reveals an unproductive form of the enzyme. FEBS J. 278, 4917–26 (2011). [PMC free article] [PubMed]
  • Nishio K. et al. . Crystal structure of the de-ubiquitinating enzyme UCH37 (human UCH-L5) catalytic domain. Biochem. Biophys. Res. Commun. 390, 855–60 (2009). [PubMed]
  • Misaghi S. et al. . Structure of the ubiquitin hydrolase UCH-L3 complexed with a suicide substrate. J. Biol. Chem. 280, 1512–20 (2005). [PubMed]
  • Lee Y. T. C. & Hsu S.-T. D. Familial mutations and post-translational modifications of UCH-L1 in Parkinson’s disease and neurodegenerative disorders. Curr. Prot. Pept. Sci., doi: , in press (2016).10.2174/1389203717666160217143721 [PubMed] [Cross Ref]
  • Zhou Z. R., Zhang Y. H., Liu S., Song A. X. & Hu H. Y. Length of the active-site crossover loop defines the substrate specificity of ubiquitin C-terminal hydrolases for ubiquitin chains. Biochem. J. 441, 143–9 (2012). [PubMed]
  • Case A. & Stein R. L. Mechanistic studies of ubiquitin C-terminal hydrolase L1. Biochemistry 45, 2443–52 (2006). [PubMed]
  • Luchansky S. J., Lansbury P. T. Jr. & Stein R. L. Substrate recognition and catalysis by UCH-L1. Biochemistry 45, 14717–25 (2006). [PubMed]
  • Boudreaux D. A., Chaney J., Maiti T. K. & Das C. Contribution of active site glutamine to rate enhancement in ubiquitin C-terminal hydrolases. FEBS J. 279, 1106–18 (2012). [PMC free article] [PubMed]
  • Dang L. C., Melandri F. D. & Stein R. L. Kinetic and mechanistic studies on the hydrolysis of ubiquitin C-terminal 7-amido-4-methylcoumarin by deubiquitinating enzymes. Biochemistry 37, 1868–79 (1998). [PubMed]
  • Lim N. C. & Jackson S. E. Molecular knots in biology and chemistry. J. Phys. Condens. Matter 27, 354101 (2015). [PubMed]
  • Virnau P., Mallam A. & Jackson S. Structures and folding pathways of topologically knotted proteins. J. Phys. Condens. Matter 23, 033101 (2011). [PubMed]
  • Jackson S. E., Suma A. & Micheletti C. How to fold intricately: using theory and experiments to unravel the properties of knotted proteins. Curr. Opin. Struct. Biol. 42, 6–14 (2017). [PubMed]
  • Mallam A. L. & Jackson S. E. Knot formation in newly translated proteins is spontaneous and accelerated by chaperonins. Nat. Chem. Biol. 8, 147–53 (2012). [PubMed]
  • Mallam A. L., Rogers J. M. & Jackson S. E. Experimental detection of knotted conformations in denatured proteins. Proc. Natl. Acad. Sci. USA 107, 8189–94 (2010). [PubMed]
  • Shih P.-M. et al. . Random-coil behavior of chemically denatured topologically knotted proteins revealed by small-angle X-ray scattering. J. Phys. Chem. B 119, 5437–5443 (2015). [PubMed]
  • Christian T. et al. . Methyl transfer by substrate signaling from a knotted protein fold. Nat. Struct. Mol. Biol. (2016). [PubMed]
  • Andersson F. I., Pina D. G., Mallam A. L., Blaser G. & Jackson S. E. Untangling the folding mechanism of the 52-knotted protein UCH-L3. FEBS J. 276, 2625–35 (2009). [PubMed]
  • Lou S. C. et al. . The knotted protein UCH-L1 exhibits partially unfolded forms under native conditions that share common structural features with its kinetic folding intermediates. J. Mol. Biol. 428, 2507–20 (2016). [PubMed]
  • Wang I., Chen S. Y. & Hsu S.-T. D. Unraveling the folding mechanism of the smallest knotted protein, MJ0366. J. Phys. Chem. B 119, 4359–70 (2015). [PubMed]
  • Wang L. W., Liu Y. N., Lyu P. C., Jackson S. E. & Hsu S.-T. D. Comparative analysis of the folding dynamics and kinetics of an engineered knotted protein and its variants derived from HP0242 of Helicobacter pylori. J. Phys. Condens. Matter 27, 354106 (2015). [PubMed]
  • Wang I., Chen S. Y. & Hsu S.-T. D. Folding analysis of the most complex Stevedore’s protein knot. Sci. Rep. 6, 31514 (2016). [PMC free article] [PubMed]
  • Hsu S.-T. D. Protein knotting through concatenation significantly reduces folding stability. Sci. Rep. 6, 39357 (2016). [PMC free article] [PubMed]
  • Andersson F. I. et al. . The effect of Parkinson’s-disease-associated mutations on the deubiquitinating enzyme UCH-L1. J. Mol. Biol. 407, 261–72 (2011). [PubMed]
  • Krishna M. M., Hoang L., Lin Y. & Englander S. W. Hydrogen exchange methods to study protein folding. Methods 34, 51–64 (2004). [PubMed]
  • Nishikawa K. et al. . Alterations of structure and hydrolase activity of parkinsonism-associated human ubiquitin carboxyl-terminal hydrolase L1 variants. Biochem. Biophys. Res. Commun. 304, 176–183 (2003). [PubMed]
  • Ohayon S., Spasser L., Aharoni A. & Brik A. Targeting deubiquitinases enabled by chemical synthesis of proteins. J. Am. Chem. Soc. 134, 3281–9 (2012). [PubMed]
  • Kranz J. K. & Schalk-Hihi C. Protein thermal shifts to identify low molecular weight fragments. Methods Enzymol. 493, 277–98 (2011). [PubMed]
  • Kumar S. M., Lyu P. C. & Hsu S.-T. D. Structural perturbation of the Parkinson’s disease-associated I93M mutation in human UCH-L1 revealed by solution state NMR spectroscopy. Chinese Journal of Magnetic Resonance 32, 329–340 (2015).
  • Larsen C. N., Krantz B. A. & Wilkinson K. D. Substrate specificity of deubiquitinating enzymes: ubiquitin C-terminal hydrolases. Biochemistry 37, 3358–68 (1998). [PubMed]
  • Ziegler F. et al. . Knotting and unknotting of a protein in single molecule experiments. Proc. Natl. Acad. Sci. USA (2016). [PubMed]
  • Plaxco K. W., Simons K. T. & Baker D. Contact order, transition state placement and the refolding rates of single domain proteins. J. Mol. Biol. 277, 985–94 (1998). [PubMed]
  • Morrow M. E. et al. . Stabilization of an unusual salt bridge in ubiquitin by the extra C-terminal domain of the proteasome-associated deubiquitinase UCH37 as a mechanism of its exo specificity. Biochemistry 52, 3564–78 (2013). [PMC free article] [PubMed]
  • Viguera A. R. & Serrano L. Loop length, intramolecular diffusion and protein folding. Nat. Struct. Biol. 4, 939–46 (1997). [PubMed]
  • Ladurner A. G. & Fersht A. R. Glutamine, alanine or glycine repeats inserted into the loop of a protein have minimal effects on stability and folding rates. J. Mol. Biol. 273, 330–7 (1997). [PubMed]
  • Nagi A. D. & Regan L. An inverse correlation between loop length and stability in a four-helix-bundle protein. Fold Des. 2, 67–75 (1997). [PubMed]
  • Sanchez I. E. Protein folding transition states probed by loop extension. Protein Sci. 17, 183–6 (2008). [PubMed]
  • Sondek J. & Shortle D. Accommodation of single amino acid insertions by the native state of staphylococcal nuclease. Proteins 7, 299–305 (1990). [PubMed]
  • Martinez J. C. et al. . Thermodynamic analysis of alpha-spectrin SH3 and two of its circular permutants with different loop lengths: discerning the reasons for rapid folding in proteins. Biochemistry 38, 549–59 (1999). [PubMed]
  • Dagan S. et al. . Stabilization of a protein conferred by an increase in folded state entropy. Proc. Natl. Acad. Sci. USA 110, 10628–33 (2013). [PubMed]
  • Wang I. et al. . Nucleotide contributions to the structural integrity and DNA replication initiation activity of noncoding Y RNA. Biochemistry 53, 5848–5863 (2014). [PubMed]
  • Delaglio F. et al. . NMRPipe: a multidimensional spectral processing system based on UNIX pipes. J. Biomol. NMR 6, 277–93 (1995). [PubMed]
  • Hsu S.-T. D., Behrens C., Cabrita L. D. & Dobson C. M. 1H, 15N and 13C assignments of yellow fluorescent protein (YFP) Venus. Biomol. NMR Assign. 3, 67–72 (2009). [PubMed]
  • Hsu S.-T. D., Cabrita L. D., Christodoulou J. & Dobson C. M. 1H, 15N and 13C assignments of domain 5 of Dictyostelium discoideum gelation factor (ABP-120) in its native and 8M urea-denatured states. Biomol. NMR Assign. 3, 29–31 (2009). [PubMed]
  • Hsu S.-T. D., Cabrita L. D., Fucini P., Christodoulou J. & Dobson C. M. Probing side-chain dynamics of a ribosome-bound nascent chain using methyl NMR spectroscopy. J. Am. Chem. Soc. 131, 8366–69 (2009). [PubMed]
  • Hsu S.-T. D. et al. . Structure and dynamics of a ribosome-bound nascent chain by NMR spectroscopy. Proc. Natl. Acad. Sci. USA 104, 16516–16521 (2007). [PubMed]
  • Hsu S.-T. D., Cabrita L. D., Fucini P., Dobson C. M. & Christodoulou J. Structure, dynamics and folding of an immunoglobulin domain of the gelation factor (ABP-120) from Dictyostelium discoideum. J. Mol. Biol. 388, 865–879 (2009). [PubMed]
  • Dosset P., Hus J. C., Blackledge M. & Marion D. Efficient analysis of macromolecular rotational diffusion from heteronuclear relaxation data. J. Biomol. NMR 16, 23–8 (2000). [PubMed]
  • Hsu S.-T. D. et al. . Folding study of Venus reveals a strong ion dependence of its yellow fluorescence under mildly acidic conditions. J. Biol. Chem. 285, 4859–69 (2010). [PMC free article] [PubMed]
  • Harris R. et al. . Backbone 1H, 15N and 13C resonance assignments for the 26-kD human de-ubiquitinating enzyme UCH-L3. Biomol. NMR Assign. 1, 51–3 (2007). [PubMed]
  • Keller S. et al. . High-precision isothermal titration calorimetry with automated peak-shape analysis. Anal. Chem. 84, 5066–73 (2012). [PMC free article] [PubMed]
  • Zhao H., Piszczek G. & Schuck P. SEDPHAT–a platform for global ITC analysis and global multi-method analysis of molecular interactions. Methods 76, 137–48 (2015). [PMC free article] [PubMed]
  • Xie D. & Freire E. Structure based prediction of protein folding intermediates. J. Mol Biol. 242, 62–80 (1994). [PubMed]
  • Robertson A. D. & Murphy K. P. Protein structure and the energetics of protein stability. Chem. Rev. 97, 1251–1268 (1997). [PubMed]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group