Plasticity of the proteasome-targeting signal Fat10 enhances substrate degradation
Abstract
The proteasome controls levels of most cellular proteins, and its activity is regulated under stress, quiescence, and inflammation. However, factors determining the proteasomal degradation rate remain poorly understood. Proteasome substrates are conjugated with small proteins (tags) like ubiquitin and Fat10 to target them to the proteasome. It is unclear if the structural plasticity of proteasome-targeting tags can influence substrate degradation. Fat10 is upregulated during inflammation, and its substrates undergo rapid proteasomal degradation. We report that the degradation rate of Fat10 substrates critically depends on the structural plasticity of Fat10. While the ubiquitin tag is recycled at the proteasome, Fat10 is degraded with the substrate. Our results suggest significantly lower thermodynamic stability and faster mechanical unfolding in Fat10 compared to ubiquitin. Long-range salt bridges are absent in the Fat10 structure, creating a plastic protein with partially unstructured regions suitable for proteasome engagement. Fat10 plasticity destabilizes substrates significantly and creates partially unstructured regions in the substrate to enhance degradation. NMR-relaxation-derived order parameters and temperature dependence of chemical shifts identify the Fat10-induced partially unstructured regions in the substrate, which correlated excellently to Fat10-substrate contacts, suggesting that the tag-substrate collision destabilizes the substrate. These results highlight a strong dependence of proteasomal degradation on the structural plasticity and thermodynamic properties of the proteasome-targeting tags.
Editor's evaluation
This manuscript probes the ways in which a protein tag might influence the structure, dynamics and stability of a covalently-attached substrate protein. These important findings are of significance to several fields, particularly in understanding how these influences control the abundance of proteins within a cell. The evidence provided is solid and the manuscript will be of interest to scientists working on protein folding and cellular degradation.
https://doi.org/10.7554/eLife.91122.sa0Introduction
Proteasome substrates must be conjugated with other small proteins, known as proteasome-targeting tags, to target substrates to the proteasome (Kwon and Ciechanover, 2017). Ubiquitin is a proteasome-targeting tag conjugated to substrates by posttranslational modification. It interacts with the 19 S proteasome subunit, which harbors ubiquitin receptors. Consequently, 19 S aligns its base subunits, consisting of AAA +ATPases, with the 20 S Core Particle (CP) to create a channel for substrate entry. Deubiquitinase enzymes cleave the ubiquitin, while the substrate enters the ATPases to unfold and translocate to CP, where it is cleaved into short peptides by proteases (de la Peña et al., 2018; Wehmer et al., 2017; Worden et al., 2017). Proteasomal degradation regulates most nuclear and cytosolic protein levels (Zhao et al., 2015). Global cellular protein degradation rate is modulated in the quiescent state, under nutrient stress, or during inflammation (Zhao et al., 2015; VerPlank et al., 2019; Liu et al., 2021; Zhang et al., 2017). However, mechanisms that regulate proteasome degradation in such conditions remain poorly understood. Partially unstructured regions at the N/C-termini or within the substrate are critical to proteasome engagement and degradation. Global disorder and topology of substrates also influence the degradation (Correa Marrero and Barrio-Hernandez, 2021). It is unclear if the plasticity of the proteasome-targeting tag can regulate the substrate degradation rate.
During stress or inflammation, a swift change in proteasomal degradation activity by changing the global disorder in proteins or its biased sequences may be challenging. An alternate mechanism is to upregulate a proteasome-targeting tag that rapidly degrades substrate proteins. The human leukocyte antigen-F adjacent transcript 10 (Fat10) is a proteasome-targeting tag that directly targets substrates for proteasomal degradation (Schmidtke et al., 2014). Fat10 expression is restricted to immune cells (Buerger et al., 2015), and other cell types express Fat10 upon induction by proinflammatory cytokines (Buchsbaum et al., 2012; Hipp et al., 2005; Raasi et al., 1999). Fat10 includes two ubiquitin-like domains attached by a flexible linker and is conjugated to substrate proteins by posttranslational modification. Fat10’s surface properties and interactions are distinct from ubiquitin (Aichem et al., 2018). Fat10 has a lower half-life of 1 hr, while ubiquitin has a longer half-life of ~24 hr (Haas and Bright, 1987; Aichem and Groettrup, 2020). Fat10’s shorter half-life prevents a prolonged inflammatory response that can potentially lead to apoptosis (Aichem and Groettrup, 2020). Mechanisms underlying the rapid degradation of Fat10 are poorly understood.
The substrates of Fat10 overlap significantly with ubiquitin substrates (Aichem and Groettrup, 2020), suggesting that Fat10 is an auxiliary proteasomal-targeting signal activated during inflammation. However, several mechanistic differences exist between the ubiquitin- and Fat10-proteasome pathways. Ubiquitin is cleaved by deubiquitinating enzymes at the proteasome before the substrate enters ATPases, but Fat10 remains uncleaved and is degraded along with the substrate (Hipp et al., 2005). Well-folded substrates must be globally or partially unfolded for engagement with the proteasome. Accessory unfoldases like Cdc48 assist the process by unfolding substrates (Olszewski et al., 2019) before they interact with the proteasome. Substrate unfolding by accessory unfoldases is necessary for the ubiquitin substrates but not for the Fat10 substrates (Aichem et al., 2018), suggesting a profound impact of Fat10 on the substrate’s structure and energetics. However, Fat10’s effect on the substrate’s structure and thermodynamics is unknown.
While ubiquitin has a rigid structure, Fat10 has a ductile fold, providing an opportunity to probe if the plasticity of proteasome-targeting tags affects substrate degradation rate. In this work, we correlate the thermodynamic properties and conformational plasticity of Fat10 and ubiquitin with the proteasomal degradation rate of their substrates. Our results suggest that the Fat10 free energy barrier is substantially lower than ubiquitin. Consequently, the Fat10 unfolding kinetics are sevenfold faster than ubiquitin, and Fat10’s mechanical resistance to unfolding is weaker than ubiquitin’s. The absence of long-range salt bridges in Fat10 creates partially unstructured regions, leading to efficient proteasome engagement and degradation. Furthermore, Fat10’s structural plasticity reduces the thermodynamic stability of substrate proteins in cellular and in vitro conditions and creates local partially unstructured regions in the substrate. The substrate reciprocally reduces Fat10 stability by thermodynamic coupling, creating more disorder in the conjugate and enhancing degradation. NMR experiments measured enthalpy and conformational entropy to reveal Fat10-induced sites of local disorder in the substrate. Our experimental and computational data suggest that nonspecific collisions with the proteasome-targeting tag destabilize the substrate. These results highlight that the proteasome-targeting tags’s plasticity regulates proteasomal degradation.
Results
Fat10 undergoes rapid proteasomal degradation
Although UBLs are structurally similar, they are diverse in sequence landscape. The N-terminal and C-terminal Fat10 domains (Fat10D1 and Fat10D2) share 29% and 36% of their sequence identity with ubiquitin, respectively (Figure 1A–C). Moreover, the two domains of Fat10 share only 18% of their sequence identity. Important regions like the L8-I44-V70 hydrophobic patch, an interaction hotspot in ubiquitin, are absent in Fat10. A phylogenetic tree analysis based on available ULM structures showed that among the two Fat10 domains, the C-terminal Fat10D2 is structurally closer to ubiquitin (Figure 1—figure supplement 1A–D). Given the poor sequence identity between ubiquitin and Fat10, their cellular levels could be distinctly regulated.
The degradation rate of ubiquitin and Fat10 was compared in cells using cycloheximide chase assay. Non-conjugable variants UbK0-GV, UbK0-AV and Fat10-AV were expressed in HEK293T3 cells, and the protein levels were measured against time (Figure 1D, Figure 1—figure supplement 1E). Fat10 levels dropped sharply compared to ubiquitin, suggesting that Fat10 undergoes rapid degradation (Figure 1E). Treatment with a proteasomal inhibitor restores Fat10 (Figure 1F), indicating that Fat10 undergoes proteasomal degradation. The protein levels of individual Fat10 domains were also studied. Both domains degraded faster than the full-length Fat10 (Figure 1E, G and H), suggesting that the Fat10 is more stable than the individual domains. The isolated Fat10 domains were also degraded via the proteasomal pathway (Figure 1—figure supplement 1F and G). Protein cofactors specific to the full-length Fat10 may bind and bury its partially disordered regions to reduce proteasomal degradation. Transient interactions between the isolated domains may also bury the disordered region and provide a more compact structure to the full-length Fat10, reducing its degradation.
Thermodynamic characterization of Fat10
To study the correlation between Fat10 degradation and its structural and thermodynamic properties, Fat10 unfolding was studied in a purified system. Fat10 was recombinantly expressed and purified. The far-UV Circular Dichroism (CD) spectra and two-dimensional 15N Heteronuclear Single Quantum Coherence (15N-HSQC) NMR spectra of Fat10 suggested a well-folded tertiary structure (Figure 1—figure supplement 2A, B). The thermodynamics of unfolding Fat10 and other ULMs were measured by chemical denaturation using guanidinium hydrochloride (GdnHCl). Gibbs free energy of unfolding (ΔGunfolding) for ubiquitin was 8 kcal/mol as shown in Figure 1I–J. ΔGunfolding for another ULM, SUMO1 was lower than ubiquitin (3.7 kcal/mol). However, ΔGunfolding was lowest for Fat10 (2.3 kcal/mol). Ubiquitin also has high thermodynamic stability, and its melting temperature (Tm) is greater than 95 °C (Wintrode et al., 1994), whereas Fat10 unfolded with a Tm of 55 °C (Figure 1—figure supplement 2C), suggesting Fat10 has a malleable structure compared to ubiquitin. Differential scanning fluorimetry measurements reported similar differences in Tm between Fat10 and ubiquitin (Aichem et al., 2018).
Transient interactions between the Fat10 domains may bury partially disordered regions and increase the thermodynamic stability of the full-length protein. In that case, the unfolding energies of individual domains shall be lower than the full-length protein (Batey et al., 2008). The two Fat10 domains were isolated, and their energetics were measured. Far-UV CD and 15N-HSQC spectra confirmed that the purified isolated Fat10 domains are folded (Figure 1—figure supplement 3A–C). GdnHCl-induced and temperature-induced denaturation of the Fat10 domains were carried out (Figure 1—figure supplement 3D–F). ΔGunfolding of individual domains is lower than full-length Fat10 (Figure 1J), suggesting transient interactions between the Fat10 domains stabilize the full-length protein. ΔGunfolding and Tm values indicate that Fat10 has substantially lower thermodynamic stability than ubiquitin.
Fat10 has low resistance to mechanical unfolding
Protein mechanical unfolding by molecular motors, such as the proteasomal ATPases, can be simulated by steered molecular dynamics (Lu et al., 1998). The mechanical unfolding of Fat10 was studied in the explicit solvent by Adaptive Steered Molecular Dynamics (ASMD) (Figure 2A, Figure 2—video 1, Figure 2—figure supplement 1, Mücksch and Urbassek, 2016). A linear di-ubiquitin (Ub2), whose size is similar to Fat10, was used for comparison. While the work required to unfold and linearize Ub2 was 735 kcal/mol, the same for Fat10 was 581 kcal/mol (Figure 2B). The lower work required to unfold Fat10 corroborated the lower value of ΔGunfolding measured during Fat10 denaturation (Figure 1J). The simulation was repeated for each isolated Fat10 domain and compared to monoubiquitin. The work required to unfold the isolated Fat10 domains was lower than ubiquitin by 30–60 kcal/mol (Figure 2B). Resistance during mechanical unfolding arises from cooperative packing interactions between buried sidechain atoms (Brockwell et al., 2005). The lower resistance in Fat10 suggests weaker cooperative packing interactions at the protein’s buried core.
Weak cooperative intramolecular interactions should result in faster unfolding kinetics. Fat10 unfolding kinetics was studied by MD simulations at high temperatures. Three μs (300ns * 10 replica) simulations were performed at 450 K for Ub2 and Fat10 (Figure 3—video 1). Native intramolecular contacts, including Van der Waals (VdW) interactions and hydrogen bonds (hbonds), decrease over time as the protein unfolds at high temperatures. While Ub2 unfolded gradually in 300 ns, Fat10 unfolded within 40 ns (Figure 2C). The native contacts and hbonds in the individual isolated Fat10 domains also disrupted faster than ubiquitin (Figure 2D). The ULM proteins adopt a β-grasp fold consisting of a five-strand β-sheet. The β-sheet hbonds were disrupted promptly in Fat10 as opposed to ubiquitin (Figure 2E and F), confirming weak cooperative interactions in Fat10.
The fraction of native β-sheet hbonds in Fat10 and ubiquitin were plotted in Figure 2G–I. Distinct sequences of inter β-strand hbonds disruption between the Fat10 domains and ubiquitin suggest distinct unfolding pathways. A comprehensive analysis of the secondary structure hbonds during unfolding showed that the 310 helix α2 unfolds first in ubiquitin, followed by β3-β4, α1, β3-β5, and β1-β5 contacts (Figure 2J). The β1-β2 contacts break at the last step. The unfolding pathway is similar to that observed in equilibrium simulations of ubiquitin (Piana et al., 2013). In contrast, the first unfolding event in Fat10D1 is the disruption of β3-β4 contacts, followed by α2, α1, and β1-β5 contacts (Figure 2K). The β3-β5 and β1-β2 contacts are disrupted simultaneously at the last step. Similarly, the unfolding pathway of the second domain in Fat10 is distinct from ubiquitin (Figure 2L). Altogether, the resistance to mechanical unfolding is lower in Fat10 than in ubiquitin, the kinetics of Fat10 unfolding is faster, and the Fat10 unfolding pathway is distinct from the ubiquitin unfolding pathway.
The absence of key interactions creates partially unstructured regions in Fat10
The flexible Fat10 structure may sample higher energy states with partially unstructured regions, which was investigated by MD simulations at room temperature. Root Mean Square Fluctuations (RMSF) of Fat10 domains are high at the β1-β2 loop, the C-terminal end of α1, and the α2 loop, indicating regions with high flexibility (Figure 3—figure supplement 1A). 2D free energy of Fat10 conformations was calculated as a function of the negative logarithm of root mean square deviation (RMSD) and the radius of gyration (Rgyr) populations. Ubiquitin was used as a control in these simulations. The 2D free energy plot of ubiquitin reflected a single highly populated state (Figure 3A), whose representative structure superimposes with the native ubiquitin x-ray structure (Figure 3B). In contrast, the native states of Fat10 domains are in equilibrium with higher energy partially unfolded states (Figure 3A). Local interactions between β1-β2 and helix α1 are disrupted, and loops at the N- and C-termini of α1 are disordered in these partially unfolded states.
The radius of gyration (Rog) reflects the nature and extent of packing interactions in the protein. Although ubiquitin and the Fat10 domains are similar in size, the Rog of the Fat10 ground state is higher than ubiquitin, suggesting more flexibility and lower packing (Figure 3C). VdW and electrostatic interaction energies were calculated between the interacting secondary structures in Fat10 and ubiquitin to study the basis of reduced packing interactions. The VdW interaction energy was comparable between Fat10 domains and ubiquitin (Figure 3D). The electrostatic interaction energy was also similar in all regions except between β2 and α1, where Fat10 has a significantly lower negative value than ubiquitin, suggesting fewer electrostatic contacts (Figure 3D). Two hydrogen bonds and a salt bridge between β2 and α1 are present in ubiquitin but absent in Fat10, which reduces the interaction between β2 and α1 (Figure 3—figure supplement 1B). Another salt bridge between the α2 loop and α1 is also exclusive to ubiquitin (Figure 3—figure supplement 1C). These salt bridges may be crucial to the stability and compactness of ubiquitin, and their absence makes Fat10 flexible.
When these salt bridges were disrupted in ubiquitin by substitution, the free energy landscape changed, giving rise to partially unfolded forms similar to Fat10 (Figure 3—figure supplement 1D). The RMSF values of mutant ubiquitin increase at the C-terminal end of α1 and the β4-α2 loop. Principle component analysis (PCA) and the RMSF values suggested that the ubiquitin mutants have increased dynamics in the regions of broken salt bridges (Figure 3—figure supplement 1E-G). When the unfolding kinetics between mutants and wt protein were compared by simulations at 450 K, the mutants unfolded much faster than ubiquitin but were similar to Fat10 domains (Figure 3E). The thermodynamic stability of ubiquitin mutants was measured experimentally. Disrupting the salt bridges reduces the melting point of ubiquitin by 50 °C, suggesting the importance of the salt bridges for ubiquitin stability (Figure 3F-G). The salt bridges were engineered in the N-terminal Fat10 domain by appropriate substitutions. The protein with engineered salt bridges populated the ground state but not the partially unfolded states (Figure 3—figure supplement 2). RMSF and PCA analysis suggested that engineered salt bridges reduced the fluctuations and increased rigidity in the domain. Overall, critical electrostatic interactions present between α1 and β1β2 and between α1 and α2 loop stabilize the ubiquitin fold. These interactions are absent in Fat10, creating a flexible structure in dynamic equilibrium with partially unfolded forms.
Weak hydrogen bonds and high conformational entropy detected in Fat10
Protein thermodynamic stability depends on the energy of enthalpic interactions and conformational entropy. The temperature dependence of backbone amide proton chemical shifts is an excellent reporter of backbone hbond strength, which contributes to enthalpic interactions (Doyle et al., 2016). The temperature dependence of the amide proton chemical shifts is the temperature coefficient Tc, where Tc = ΔδNH/ΔT (Cierpicki and Otlewski, 2001). Higher negative Tc values suggest weaker hbonds and disorder propensity, while lower negative Tc values suggest stronger hbonds and structural rigidity. Tc’s were measured in Fat10 and plotted against each residue in Figure 4A, and the fitting errors are provided in Figure 4—figure supplement 1A. In the N-terminal Fat10 domain, residues in the region β1 to β2 and the loop between including α2 had high negative Tc values, suggesting weaker hbonds and disorder propensity in these regions (Figure 4A). Similarly, several residues in the β1, β2, α1, and α2 loop in the C-terminal domain had high negative Tc values.
Tc’s were measured in ubiquitin and plotted against each residue in Figure 4B (fitting errors in Figure 4—figure supplement 1B). The loop between β1 and β2 had high negative Tc values in ubiquitin, which correlates with the lack of interactions and disorder propensity in this region (Peters and de Groot, 2012). Due to the low sequence similarity between Fat10 domains and ubiquitin, comparing the Tc values between individual residues of the two proteins is challenging. Instead, Tc values were averaged over the different protein segments and compared (Figure 4—figure supplement 1C). The difference in averaged Tc values between the N-terminal domain and ubiquitin indicates that Fat10 hbonds are weaker in the β1β2 region and the α2 loop (Figure 4C and D). Similarly, the β1β2 regions and α1 have weaker hbonds in the C-terminal domain. The Tc measurements of the isolated Fat10 domains yielded similar results (Figure 4—figure supplement 1D-J).
To estimate the conformational entropy in Fat10, the backbone order parameters (S2) were measured by standard NMR relaxation experiments. The NMR spin-lattice relaxation parameter R1, spin-spin relaxation parameter R2, and heteronuclear NOEs (hetNOE) were measured in Fat10 and ubiquitin (Figure 4—figure supplement 2). S2 values were calculated using the Lipari-Szabo Model-Free Analysis method and averaged over the protein segments. The difference in averaged S2 between Fat10 and ubiquitin was converted to conformational entropy –TΔS (Sharp et al., 2015; Zeng et al., 2017) and plotted in Figure 4F–H. The higher value of -TΔS suggests higher dynamics and flexibility. The β1β2 and α2 loops had higher conformational entropy in the first Fat10 domain than ubiquitin (Figure 4F and G). In the second domain, β1β2, β3β4, and helix α1 were more entropic compared to ubiquitin, the most significant being α1 (Figure 4G and H). Overall, salt bridges increase interactions between α1 and β1β2 and between α1 and the α2 loop in ubiquitin. In their absence, these regions in Fat10 have weaker hbonds and higher conformational entropy.
Fat10 increases substrate unfolding and degradation in cells
We then studied the degradation of Fat10 conjugated substrates in cells by a cycloheximide assay using GFP as the model. Fat10-conjugated GFP degraded by 60% within 4 hr after cycloheximide treatment, much faster than diubiquitin-conjugated GFP (Figure 5A, Figure 5—figure supplement 1A). To study the effect of individual Fat10 domains, we compared the degradation of GFP in Fat10D1-GFP and Ub-GFP (Figure 5B, Figure 5—figure supplement 1B). The Fat10D1-GFP degraded by 40% within 4 hr, whereas the levels of Ub-GFP remained unchanged. It is unclear if the higher degradation rate is exclusive because the proteasome-targeting tag Fat10 degrades rapidly, thereby accelerating substrate degradation, or if it also affects the substrate’s thermodynamic stability and induces partially unstructured regions. To study changes in the substrate thermodynamic stability due to Fat10 in cellular conditions, a heterologous system has to be used where Fat10 substrates are not degraded. The bacterial Pup-proteasome system functions by a distinct mechanism (Striebel et al., 2014) that does not recognize or degrade the eukaryotic Ub/Fat10 conjugated proteins and can be used as a heterologous host. The CRABP1 protein was chosen as the model substrate. CRABP1 can be engineered to bind a fluorescent dye that is quenched in the native state but not in the unfolded state (Ignatova and Gierasch, 2004; Figure 5C). Cells expressing CRABP1 were treated with various urea concentrations. The difference in fluorescent signals between 0 M and 3 M Urea was lowest in apo CRABP1, higher in Ub-CRABP1, and highest in Fat10-CRABP1, suggesting the CRABP1 stability is lowest when conjugated to Fat10 (Figure 5—figure supplement 1C-E). ΔGunfolding was 5.9 kcal/mole for CRABP1, 4.8 kcal/mole for Ub-CRABP1, and 2.4 kcal/mole for Fat10-CRABP1 (Figure 5D–E, Figure 5—figure supplement 1F). The difference of 2.4 kcal/mol between Fat10-CRABP1 and Ub-CRABP1 suggested a significant increase in substrate unfolding by Fat10 conjugation.
Fat10 induces partially unfolded forms in the substrate
In highly crowded cellular environments, proteins experience various physical interactions, such as the exclusion volume effect and nonspecific transient interactions. To study whether the substrate destabilization is solely due to Fat10 and not the above effects, the impact of Fat10 was studied in purified proteins. Producing purified Fat10 isopeptide conjugated substrates is technically challenging owing to the lack of a suitable Fat10 E3-substrate pair active under in vitro conditions. Hence, Fat10 was covalently conjugated at the N-terminus of a model ultra-stable substrate protein like the Cyan Fluorescent Protein (CFP) to mimic the N-terminal isopeptide conjugation (Ciechanover and Ben-Saadon, 2004; Figure 6A). The free energy of unfolding CFP was calculated by measuring CFP fluorescence during chemical denaturation. CFP is a well-folded substrate with high thermodynamic stability and ΔGunfolding = 11 kcal/mol. The change in unfolding energy of CFP is modest ΔΔGunfolding = 0.2 kcal/mol when conjugated with ubiquitin (Figure 6B, Figure 6—figure supplement 1A). The modest decrease in CFP stability correlates with the finding that ubiquitin conjugation is insufficient for the direct degradation of well-folded substrates and requires Cdc48 (Olszewski et al., 2019). Interestingly, the ΔΔGunfolding was 3.0 kcal/mol and 3.8 kcal/mol when CFP is conjugated to Fat10D1 and Fat10D2, respectively (Figure 6B, Figure 6—figure supplement 1B), suggesting Fat10 domains destabilize the substrate 15-fold more than ubiquitin. Change in free energy could be due to the nonspecific or specific interactions between Fat10 and CFP (Bigman and Levy, 2020). Specific interactions should depend on the complementary surface properties of Fat10 and CFP and should differentiate between the full-length Fat10 and the smaller individual domains. CFP’s thermodynamic stability after conjugation to the full-length Fat10 reduced by 3kcal/mol, similar to the isolated domains, suggesting that nonspecific interactions due to conjugation reduce CFP stability.
Fat10 may increase the partially unfolded regions in the substrate, which is critical to interact with the proteasome. The susceptibility of partially unfolded conformations to a typical protease thermolysin was measured in CFP, ubiquitin-CFP, and Fat10-CFP (Figure 6B–D, Figure 6—figure supplement 1C-F). Although CFP was stable in the presence of protease, ubiquitin-CFP and Fat10-CFP were proteolyzed. The Fat10-CFP proteolysis rate (kobs) was significantly more than ubiquitin-CFP at a given thermolysin concentration (Figure 6D). The change in observed proteolysis rates with different thermolysin concentrations (mproteolysis) was sevenfold higher in Fat10-CFP than in ubiquitin-CFP (Figure 6D). The values of mproteolysis can be correlated to the free energy of proteolysis (Park and Marqusee, 2004). The free energy of proteolysis decreased from 8.7 kcal/mol for ubiquitin-CFP to 7.5 kcal/mol for Fat10-CFP (Figure 6B), indicating that Fat10 substantially increases the partially unstructured regions in CFP. Our results indicate that Fat10 has a more drastic effect than ubiquitin on substrate stability, suggesting that Fat10 substrates may be degraded independently of the unfoldases (Aichem et al., 2018). Degradation of Fat10-GFP was monitored in the HEK293T cells in the presence and absence of the Cdc48 inhibitor DBeQ (Figure 6E). Fat10-GFP degradation was unperturbed by DBeQ, suggesting that Fat10-modified substrates are degraded independently of Cdc48. Altogether, Fat10 significantly destabilizes substrates and creates partially unstructured regions for direct proteasomal degradation.
Fat10 reduces local stability in the substrate
To probe the Fat10-induced disorder further, we used a combination of computational and experimental frameworks to detect changes in the substrate’s local stability. Ubiquitin was chosen as a well-folded small substrate, as covalent interactions between Fat10 and ubiquitin have been reported previously (Buchsbaum et al., 2012). MD simulations studied the impact of Fat10 on the kinetics of substrate ubiquitin unfolding at 450 K. As a control, we used diubiquitin, where the first ubiquitin serves as the tag and the second ubiquitin is the substrate. The substrate ubiquitin in diubiquitin unfolded similarly to its free form (Figure 7—figure supplement 1A, C). However, the substrate ubiquitin conjugated to Fat10 domains unfolded faster (Figure 7—figure supplement 1B, D), suggesting that Fat10 domains increase the rate of substrate unfolding. Fat10 also enhanced fluctuations in the helix α1 and the loop near α2 in the substrate (Figure 7A). It increased the radius of gyration in substrate ubiquitin from 11 Å to 12 Å and created a broad RMSD profile starting from 2 Å to 3 Å, indicating a less compact substrate (Figure 7B). The domains had a similar effect and created partially disordered states in the substrate (Figure 7B).
Fat10’s effect on substrate stability was experimentally measured using Ub45F/W mutant as a substrate, where the intrinsic tryptophan fluorescence reports the thermodynamics of substrate unfolding. Ub45F/W is functionally active, and its thermodynamic stability is similar to ubiquitin (Khorasanizadeh et al., 1993). A single tryptophan in the N-terminal Fat10 domain is mutated to Phenylalanine (F) to not interfere with the fluorescence spectra (Figure 7C). The free energy of unfolding the substrate reduced from 11.7 kcal/mol to 7 kcal/mol when covalently bound to Fat10 (Figure 7D). Individual domains of Fat10 had a similar effect on the substrate (Figure 7E). C-terminal conjugation has a higher impact than at the N-termini (Figure 7F and G, Ub-Fat10D2 versus Fat10D2-Ub), suggesting that the conjugation site may regulate substrate stability.
To further understand substrate destabilization, we analyzed the tag-substrate and intra-substrate contacts in MD simulations of the Fat10~Ub conjugate. Many tag-substrate contacts were formed between the D2 domain and the β1β2 and β4-α2-β5 regions in the substrate ubiquitin (Figure 8—figure supplement 1A, B). We compared the long-range intra-substrate contacts in the free substrate and the Fat10-substrate conjugate (Figure 8—figure supplement 1C). The regions where intra-substrate contacts are disrupted overlapped with those where Fat10 contacts the substrate, suggesting that Fat10 interaction disrupts the intra-substrate contacts. To study the higher destabilizing effect of Fat10 compared to ubiquitin, we noted that Fat10 has a higher exposed hydrophobic surface area (23%) than ubiquitin (13%) and can form more hydrophobic interactions with the substrate. Moreover, Fat10 domains are labile and sample partially unfolded forms, exposing their buried hydrophobic residues to enhance hydrophobic interactions with the substrate. We compared the intra-substrate long-range contacts between a linear diubiquitin molecule and Fat10-ubiquitin conjugate. The first ubiquitin serves as the tag in diubiquitin, and the second ubiquitin is the substrate. The number of intra-substrate contacts disrupted in the Fat10-substrate conjugate was higher than diubiquitin, confirming the higher destabilization effect of Fat10 (Figure 8—figure supplement 2A and B).
Most intermolecular contacts were observed between the Fat10 D2 domain and the substrate, and few between the D1 domain and substrate. Hence, we conjugated the D2 domain to the N-terminal or C-terminal end of the substrate ubiquitin to study the significance of the conjugation site (Figure 8A and B). The C-terminal tail in substrate ubiquitin is dynamic and explores a larger conformational space, providing greater flexibility for the tag to interact with the substrate (Figure 8B). Consequently, the intermolecular contacts between the tag and substrate were more evenly distributed (Figure 8—figure supplement 2C, D). In addition, a more significant number of the intra-substrate contacts were disrupted (Figure 8—figure supplement 2A, E, F), which corroborates the reduced substrate stability upon C-terminal conjugation (Figure 7F, Ub-Fat10D2 versus Fat10D2-Ub). These results highlight that conformational dynamics at the conjugation site increase tag-substrate contacts and decrease substrate stability.
We used NMR spectroscopy to monitor the Fat10-induced changes in the substrate’s local stability. Since the destabilizing effect of Fat10 and the D2 domain are similar, and to avoid overlap in the NMR spectra, we studied the conjugates where the D2 domain is conjugated to substrate ubiquitin at its N- and C-termini (Fat10D2-Ub and Ub-Fat10D2). 15N-edited HSQC spectra of the conjugated proteins showed well-dispersed backbone amide resonances, confirming that they are folded (Figure 8—figure supplement 3A and B). Chemical shift perturbations (CSPs) of amide NMR resonances report the change in the chemical environment in a protein. High NMR CSPs were observed in additional regions in Ub when Fat10D2 is conjugated to the C-terminus (Figure 8C and D), as expected from MD simulations. We then measured Tc values to estimate Fat10D2’s effect on the strength of hbonds in the substrate (Figure 8—figure supplement 3C). Differences in averaged Tc values between Fat10D2-ubiquitin and free ubiquitin show that the hbonds in the strands β1β2, helix α1, and the α2 loop in the substrate ubiquitin are weakened (Figure 8E). Additional hbonds are destabilized in β3β4 and β5 when Fat10 is conjugated at the C-terminus, which is commensurate with its lower ΔGunfolding. Experimental NMR Tc measurements correlate excellently with the changes in local intra-substrate contacts observed in the MD simulations (Figure 8E and Figure 8—figure supplement 2E, F). The substrate’s conformational entropy was quantified from NMR relaxation order parameters (Figure 8—figure supplement 4). Conjugation at the N-terminal increases entropy at the β1β2 loop, β2, and α1 in the substrate (Figure 8F). Entropy increases further in these regions when D2 is conjugated at the C-terminus. In addition, the α2 loop and β3β4 regions become more entropic. Overall, the Fat10 domain contacts the substrate at multiple regions to weaken its local stability and increase disorder. The intrinsic flexibility of the conjugation site regulates Fat10’s impact on substrate stability.
The substrate may also influence the stability of the proteasome-targeting tag. The Fat10D2 tag is more compact (low Rog) in free form than in the tag-substrate conjugate (Figure 8—figure supplement 5A). Moreover, the Fat10D2 tag unfolds faster when conjugated with the substrate (Figure 8—figure supplement 5B), suggesting a thermodynamic coupling between the tag and the substrate. The tag-substrate contacts are greater when the tag is conjugated to ubiquitin C-termini (Figure 8—figure supplement 1A). Consequently, the tag unfolded faster when conjugated to the substrate C-termini (Figure 8—figure supplement 5C), suggesting that the tag stability also depends on the conjugation site.
Discussion
Posttranslational modification with proteasome-targeting tags is necessary but insufficient for protein degradation. The substrate/proteasome engagement and the substrate unfolding are rate-limiting steps for degradation (Bard et al., 2019; Matyskiela et al., 2013). Global disorder, topology, local regions with a high disorder, and biased sequences in the substrate are critical factors that regulate proteasomal degradation (Correa Marrero and Barrio-Hernandez, 2021). Whether the structural flexibility of the proteasome-targeting tags impacts substrate degradation is unclear. We report that the proteasome-targeting tag Fat10 has a malleable structure that samples multiple partially unfolded forms at physiological temperature. The lack of strong long-range electrostatic interactions between the central helix α1, the β1β2 strands, and the α2 loop, creates the flexible fold that provides weak resistance to mechanical unfolding. These properties of Fat10 expedite its unfolding, degradation, and degradation of Fat10 substrates, highlighting the role of structural plasticity of proteasome-targeting tags in regulating protein degradation rate. Our data correlate well with the previous observations that substituting Fat10 Ubl domains with ubiquitin impedes proteasomal degradation (Aichem et al., 2018).
Fat10’s impact on the substrate structure and dynamics is noteworthy, as it destabilized various substrates in both in-vitro and cellular conditions. We chose two ultra-stable proteins, CFP and Ub, whose melting points are above 90 °C. Fat10 reduced the unfolding energies of both substrates by 3–5 kcal/mol, which was considerably higher than the effect of ubiquitin on these substrates. Our simulations suggested that Fat10 makes several nonspecific intermolecular contacts with the substrate to perturb the intra-substrate contacts. Competition between inter- and intra-molecular interactions has been recently highlighted, where critical salt-bridge interactions at protein-protein interfaces are “stolen” by new posttranslational modifications within the interacting proteins (Skinner et al., 2017; Mohanty et al., 2019). A similar mechanism may disrupt critical intra-substrate contacts by substrate-tag collision. Intriguingly, the substrate reciprocally increases partially unfolded regions in Fat10. Thermodynamic coupling between the substrate and tag functions in tandem to increase partially unfolded regions in the substrate-tag conjugate and accelerate its degradation.
Increased enthalpy and conformational entropy suggest local order-to-disorder transitions in proteins, giving rise to partially unfolded forms. We have measured the Tc values to estimate hbond strength, where higher ΔTc values indicate reduced hbond strength and increased enthalpy. Entropy was calculated from NMR relaxation order parameters. The increase in local enthalpy and entropy values identified Fat10-induced local disorder in the substrate, which agreed well with the computational data. An intriguing question is whether the chemical environment at the Fat10 conjugation site can regulate substrate stability. Our results show that conjugation at disordered regions allows greater conformational space for Fat10, increasing its collisions with the substrate and reducing substrate stability. We explored further by studying the effect of Fat10 conjugation on several substrates with varied electrostatic surface potentials (Ravichandran and Das, 2024) and found that Fat10 conjugation on positively charged surfaces has a lower/insignificant effect on substrate energies because the negatively charged Fat10 tag forms stable interactions with the substrate, which reduces its conformational flexibility. These studies suggest that the local environment around the conjugation site can regulate Fat10’s effect on the substrate energetics.
Posttranslational modification with ubiquitin destabilizes substrates (Carroll et al., 2020; Hagai and Levy, 2010; Carroll et al., 2021; Morimoto et al., 2016). However, the proteasomal degradation pathways of ubiquitin and Fat10 substrates are distinct. Post cleavage by proteasome deubiquitinases, ubiquitin’s impact on the substrate is unsustained. Hence, few substrates may not interact with the proteasome ATPases and escape degradation (Correa Marrero and Barrio-Hernandez, 2021). Fat10 remains conjugated to the substrates, and its effect on the substrate persists until the conjugate translocates to the ATPases. Since Fat10 is degraded along with the substrate, its malleability directly influences the substrate degradation rate. While the N-termini of ubiquitin is rigid, a short N-terminal disordered region in Fat10 contributes to its degradation (Aichem et al., 2018). Moreover, Fat10 and ubiquitin’s impact on the substrate structure and energetics are distinct. Fat10 domains have a greater exposed hydrophobic surface than ubiquitin, which can induce more nonspecific substrate-tag collisions. Fat10’s ground state is in equilibrium with partially unfolded forms with exposed hydrophobic patches, further enhancing interactions with the substrate. Together, these structural and thermodynamic properties of Fat10 create a higher disorder in the substrate than in ubiquitin. These mechanistic differences explain the unfoldase(Cdc48)-independent rapid degradation of Fat10 substrates. In hypoxic conditions, however, the 20 S proteasome degrades ubiquitin conjugated to disordered substrates (Sahu et al., 2021). Since ubiquitin is uncleaved from the substrate, the mechanistic differences between the Fat10 and ubiquitin for 20 S proteasomal degradation is an intriguing question for further investigation.
The study has a few caveats. Fat10 was conjugated to the N- or C-terminus of the substrates here and not to specific lysines. Without an identified FAT10 E3-substrate pair that is active under in vitro conditions, it is technically challenging to conjugate substrate lysines site-specifically. The discovery of novel FAT10 E3s will facilitate such studies in the future. Previous computational research suggested that conjugation between two proteins can destabilize their folded state, stabilize their unfolded state, or both (Bigman and Levy, 2020; Sokolovski et al., 2015). While our NMR experiments have exclusively studied Fat10’s effect on the substrate’s folded state, further studies are required to reveal its impact on the substrate’s unfolded state.
The proteasome degrades 80–90% of cellular proteins (Jang, 2018), and understanding its activity has broad implications for understanding cellular response mechanisms to stress and infection. The immune cells undergo major reprogramming of signaling and antigen presentation during inflammation, and the proteasome plays an essential role in this process (Çetin et al., 2021). An important question is how proteasome degradation is modulated in these cells and why Fat10 is upregulated during infection. Our work provides a plausible understanding that the upregulation of a malleable proteasome-targeting signal Fat10 helps immune cells modulate the proteasomal degradation rate in response to infection. Fat10 is designed as an efficient proteasomal degradation signal activated as an inflammatory response. Intriguingly, Fat10 also presents scope for designing therapeutics by controlled protein degradation. Understanding the physical and thermodynamic effects of Fat10 modification is essential to engineer it for therapeutics.
Materials and methods
Cloning
Request a detailed protocolA Fat10 gene (Life Technologies) was constructed where C7, C9, C160, and C162 were mutated to alanine. Such cysteine mutants of Fat10 have been used previously for in vitro studies (Aichem et al., 2018; Truongvan et al., 2022). The thermal stability of cysteine mutant Fat10 is similar to wild-type Fat10 (Aichem et al., 2018). The gene was then sub-cloned in a pET3a vector tagged with 6xHis using NdeI and BamHI restriction sites and was used in the in vitro experiments. The N-terminal domain (Fat10D1: 1-83aa with the GG motif at the C-terminal end) and the C-terminal domain (Fat10D2: 84-165aa) were PCR amplified and further cloned into pET3a and pET14b vector, respectively. Fat10D1 was tagged with 6xHis, whereas Fat10D2 was tagged with 6xHis-scSMT. A precision protease site was introduced between the scSMT tag and Fat10D2.
The CFP plasmid was a generous gift from Dr. Aakash Gulyani, which was further sub-cloned into a pET3a vector using NdeI and BamHI restriction sites with a 6xHis tag at its N-terminus. CFP chimeric constructs Ub-CFP, Fat10-CFP, Fat10D1-CFP, and Fat10D2-CFP were cloned using overlapping PCR (Bryksin and Matsumura, 2010), where the amplified product of CFP was inserted at the C-terminus of the genes. For the ubiquitin chimeric constructs, Fat10D2 and Fat10D1 were swapped by ubiquitin to make Fat10D1-Ub and Ub-Fat10D2 constructs, respectively. The Fat10D2-Ub was first created by sub-cloning Fat10D2 in the pGEX6P1 vector with the GST tag and then fusing the PCR product of ubiquitin at the C-termini of Fat10D2. Similarly, ubiquitin was inserted at the C-terminus of Fat10 to create the Fat10-Ub construct. For the tryptophan-based melt studies, phenylalanine (F) at the 45th position of ubiquitin was mutated to tryptophan (W) in Ub, Fat10-Ub, Fat10D2-Ub, and Ub-Fat10D2. The tryptophan residue at the 17th position in Fat10 was mutated to phenylalanine in Fat10-UbF45W and Fat10D1-UbF45W by site-directed mutagenesis. All the clones were verified using the Sanger sequencing method.
3xFLAG-wtFat10 (wild type 1–165) was constructed from Life Technologies for the cellular experiments. The gene was further sub-cloned into pcDNA3.1/hyg(+) vector (a kind gift from Dr. Apurva Sarin, NCBS) using HindIII and BamH1 sites. 3xFLAG-wtFat10D1 (aa:1–81) and 3xFLAG-wtFat10D2 (aa:82–165) were further made from 3xFLAG-wtFat10 in pcDNA3.1. After confirming the 3xFLAG-wtFat10D1 clone, a C-terminal tail sequence CYCIGG (same as the C-terminal tail of wtFat10) was inserted. The GG residues in 3xFLAG-wtFat10-GG, 3xFLAG-wtFat10D1-CYCIGG, and 3xFLAG-wtFat10D2-GG were also mutated to non-conjugal AV residues to create 3xFLAG-wtFat10-AV, 3xFLAG-wtFat10D1-CYCIAV, and 3xFLAG-wtFat10D2-AV in the pcDNA3.1 vector. Mutations in Fat10 were performed by site-directed mutagenesis. The Ubk0GV-GFP construct was obtained from Addgene (#11932) (Dantuma et al., 2006) and further sub-cloned into pcDNA3.1 to get 3xFLAG-Ub-k0GV. For the cycloheximide assay, GFP chimeric constructs were made by inserting the PCR product of GFP at the C-terminus end of 1xFLAG-UbGV and 3xFLAG-wtFat10AV by PCR. For the in-cell protein stability assay, a pseudo wild type CRABP1 gene (with R131Q stabilizing mutation and tetra-Cysteine motif by replacing G102C, D103C substitution followed by further insertion of 2 Cysteine residues between 105P and 106K; Ignatova and Gierasch, 2004) was constructed from Life Technologies and sub-cloned in pET3a vector with 6xHis tag. The PCR-amplified product of CRABP1 was inserted at the C-terminus end of ubiquitin, and Fat10 constructs were used to make Ub-CRABP1 and Fat10-CRABP1. All primer sequences will be provided upon request to authors by email.
Protein purification
Request a detailed protocolAll the clones were expressed in E. coli BL21 (DE3) cells. Post-induction with 1 mM IPTG, the bacterial culture transformed with Fat10, Fat10D1, Fat10D2, Ub-CFP, Fat10-CFP, Fat10D1-CFP, and Fat10D2-CFP were grown at 18 °C for 16 hr, lysed, and resuspended in lysis buffer (Buffer A: 50 mM Tris buffer pH 7.5, 500 mM NaCl, 20 mM Imidazole and 1 mM DTT). The supernatant was loaded onto IMAC (GE) pre-packed column and eluted with imidazole gradient using Buffer A and Buffer B (Buffer A with 500 mM Imidazole). The purified fractions were applied on a size exclusion chromatography on a Superdex 16/600 75 pg column (GE Healthcare) pre-equilibrated with SEC buffer (Buffer C: 50 mM Tris buffer pH 7.5, 250 mM NaCl, with or without 1 mM DTT). All the ubiquitin proteins were also purified, as given above. GST-ppx-Fat10D2-Ub was purified using a GST pre-packed column and eluted with 10–20 mM reduced Glutathione. The purified proteins His-scSMT-Fat10D2, His-scSMT-ppx-Fat10D2-CFP, and GST-ppx-Fat10D2-Ub, were further treated with GST-tagged Precision protease enzyme followed by reverse IMAC. The digested Fat10D2, Fat10D2-CFP, and Fat10D2-Ub were purified by size exclusion chromatography. Ubiquitin and SUMO1 proteins were purified using protocols published elsewhere (Chatterjee et al., 2019; Negi et al., 2020). For the NMR experiments, Fat10, Fat10D1, Fat10D2, Fat10-Ub, Fat10D1-Ub, Ub-Fat10D2, and Fat10D2-Ub were grown in M9 media containing isotopic Ammonium chloride (15NH4Cl) and/or isotopic Glucose (13C-Glucose). The purification protocol was the same as given above.
Unfolding studies
Request a detailed protocolChemical denaturation studies by CD spectroscopy were carried out using 20 µM of protein (Ub, SUMO1, Fat10, Fat10D1, and Fat10D2) incubated overnight with various concentrations of Guanidium Chloride (GdnCl, 0–6 M) made in the native buffer at 25 °C. The change in the far-UV signal at 222 nm was monitored using the Jasco J-1500 spectropolarimeter. CFP and its variants were studied similarly by incubating 15 nM of protein with GdnCl. CFP fluorescence was measured using Horiba Fluromax-4 with an excitation wavelength of 434 nm. The emission signal was collected at 474 nm. UbF45W and its variants were similarly incubated with GdnCl. The signal at 340 nm was monitored for the change in intrinsic tryptophan fluorescence. Chemical denaturation data were normalized to a two-state unfolding equation using SIGMA plot software. The final curve fitting used a monomeric two-parameter melt equation (Privalov, 1979). The calculated free energy of unfolding, ΔGunfolding, and slope of the transition curve, m, were used to calculate free energy (ΔG0) at 0 M [GdnCl] using the Equation 1.
Proteins were dialyzed to phosphate buffer (50 mM, pH 6.5, 250 mM NaCl) with or without DTT (1 mM) for thermal equilibrium unfolding studies. The unfolding transition was monitored by observing the change in the CD signal at 222 nm on the JASCO J-1500 spectropolarimeter connected to a Peltier. Data were collected from 20°C to 95°C after every 1 °C per minute rise in temperature with 32 s of data integration. Raw data of the CD signal were normalized to a two-state unfolding equation and plotted against temperature. The curve was fitted with a sigmoidal equation to yield Tm. Thermal equilibrium unfolding data were analyzed and processed using SIGMA plot software.
Native-state proteolysis
Request a detailed protocolThe protocol for native-state proteolysis has been described previously (Carroll et al., 2020). A total of 100 μl aliquots of Ub-CFP and Fat10-CFP were treated with multiple thermolysin protease concentrations (stock concentration 10 mg ml–1). At different times, 10 μl of the reaction mixture was taken, quenched with 6x-SDS loading dye, and run on SDS-PAGE gels. An I-Bright imaging instrument captured the in-gel fluorescence image (Life Technologies). Images were inverted to greyscale and quantified by normalizing them to no protease control sample. The quantified data were fitted to the first-order exponential equation and observed proteolysis kinetics (kobs) were calculated at different thermolysin concentrations. The mean of kobs (n=3) at different concentrations was fitted to a linear equation against thermolysin concentration. ΔGproteolysis was calculated using Equation 3 by using the slope of the linear fit (Equation 2). The kcat/KM of Thermolysin is 99,000 M–1 s–1 (Park and Marqusee, 2004).
NMR studies
Request a detailed protocol13C,15N isotopically labeled Fat10D1, and Fat10D2 were prepared in 50 mM Tris, 250 mM NaCl buffer with pH 7.5. Standard triple resonance multidimensional NMR experiments 15N-HSQC, HNCA, HN(CO)CA, CBCACONH, and HNCACB were recorded at 298 K on an 800MHz Bruker Avance III HD spectrometer equipped with a cryoprobe head. These experiments yielded the backbone assignment of Fat10D1 and Fat10D2. For proteins Fat10D1-Ub and Ub-Fat10D2, only the ubiquitin residue peaks were assigned using standard triple resonance experiments HNCA, HNCOCA, CBCACONH, and HNCACB. NMR data were processed in NMRpipe (Delaglio et al., 1995) and analyzed by NMRFAM-SPARKY (Lee et al., 2015) software. After peak-picking of the backbone experimental data in SPARKY, the peaks were assigned by the PINE server (Lee et al., 2009) and confirmed manually.
Cell culture experiments
Request a detailed protocolApproximately 0.1 million cells/well of a 12-well plate were seeded with HEK293T cell line and transfected after 70% confluency using a Lipofectamine reagent (Promega). For the cycloheximide assay, cells were transfected with 3xFLAG-wtFat10-AV, 3xFLAG-wtFat10D1-CYCIAV, 3xFLAG-wtFat10D2-AV, 3xFLAG-UbK0GV and 3xFLAG-UbK0AV. The transfected cells were incubated at 37 °C for 18–20 hr, treated with Cycloheximide (CHX, final concentration 50 µg/mL), and lysed at different time points. For 6 hr, these cells were treated with MG132 (final concentration 10 µM). 20 mM Tris (pH 7.5), 150OD600 = mM NaCl, 1 X protease inhibitor cocktail, and 1% NP40 were used to lyse the cell pellet. The BCA kit (Thermo) estimated the total protein amount. Similar transfections were carried out for FLAG-UbGV-GFP, FLAG-Fat10AV-GFP, FLAG-diUbGV-GFP, and FLAG-Fat10D1-GFP. After 18 hr of incubation, the cells were treated with CHX (100 μg/mL) in the presence and absence of p97 inhibitor (DBeQ; final concentration 10 µM) and further incubated for different time points. All the Immunoblots were probed with Mouse anti-FLAG antibody (SIGMA, 1 in 10,000 dilutions) and Mouse anti-β-actin antibody (Santa Cruz, 1 in 5000 dilutions) and further probed with HRP conjugated anti-mouse secondary antibody at 1:10,000 dilutions. The blots were developed using clarity ECL (Bio-Rad) staining reagent, and the chemiluminescence signal was observed in Image Quant LAS 4000 (GE). The Immunoblots for GFP variants were probed similarly with the mouse anti-FLAG antibody. Tubulin was used as the loading control for normalization and probed with a Mouse anti-tubulin antibody (SIGMA, 1 in 3000 dilutions). All the immunoblots of GFP variants were developed, and the chemiluminescence signal was observed in the Amersham Imager 400 (GE). Quantification was done using ImageJ software with three biological replicates for each experiment.
In-cell protein stability assay
Request a detailed protocolThe stability of CRABP1 with/without ubiquitin and Fat10 in the cellular environment was monitored by pre-labeling the E. coli BL21 (DE3) cells with FIAsH-EDT2 fluorescent dye (final concentration 150 µM) at OD600 = 0.5, followed by incubation till OD600 reaches 1. IPTG was added to cultures to induce protein synthesis. After two hours, the culture solution was aliquoted (150 µl) and treated with various concentrations of Urea ranging from 0M to 3M, keeping the final volume at 400 µl. The culture solutions were incubated for 30 min before pelleting and washing with native buffer (10 mM Tris buffer, pH 7.5), followed by the fluorescent measurements (excitation 500 nm and emission 531 nm). The rest of the protocol is the same as described previously (Ignatova and Gierasch, 2004).
High-temperature unfolding simulations
Request a detailed protocolSystems were prepared in a cubic box of TIP3p water, with a minimum distance of at least 50 Å between solute atoms and the box edge. Counter ions were added to neutralize the system. The system setup and equilibration procedures were similar to the unrestrained simulations, with 450 K as the equilibration temperature. Protein unfolding simulations were performed at 450 K for 300 ns. The simulations were performed with ten replicas to obtain an average folding pathway.
ASMD simulations
Request a detailed protocolConstant force explicit solvent ASMD simulations of Fat10D1, Fat10D2, Ub, full-length Fat10, and Ub2 were carried out. The reaction coordinate is the end-to-end distance between the Cα atoms of the first amino acid and the last amino acid of the respective protein. For simulations consisting of Fat10D1, Fat10D2, and Ub, the proteins were pulled at 1 Å/ns velocity, and for full-length Fat10 and Ub2, the pulling velocity was 5 Å/ns. The reaction coordinate was partitioned into ten equal segments, each with ten trajectories. The system was energy-minimized equilibrated, and the resulting coordinates & velocities were used as starting points for ASMD simulations. The simulations were performed with periodic boundary conditions in the NPT ensemble; electrostatic interactions were computed using the PME (Particle Mesh Ewald) method. Nonbonded interactions were treated with a cutoff of 8 Å.
Structure-based Phylogenetic analysis
Request a detailed protocolThe structures of 10 available UBL PDBs were retrieved from RCSB PDB as tabulated in Fig. S1A. Crystal structures were cleaned by removing co-crystal structures, multiple domains, ions, and water. For NMR structures, the best representative model was used. For Fat10 and ISG15, individual UBL domains were isolated. RMSD matrix and structure-based sequence alignment were obtained for all structures based on multiple superimpositions using the MUSTANG tool (Konagurthu et al., 2006). Further, the phylogenetic tree was plotted using aligned sequences in MEGA X (Kumar et al., 2018).
Molecular dynamics simulations
Request a detailed protocolStructures for individual domains of Fat10D1 & Fat10D2 were retrieved from PDB IDs 6GF1 and 6GF2, respectively. The full-length Fat10 structure was modeled using Swiss-Modeller with 6GF1 and 6GF2 as templates (Waterhouse et al., 2018). Chain B of 6GF1 was further processed for simulation and called Fat10D1. All water and sulfate molecules were removed from the Fat10D1 crystal structure. 6GF2 was an NMR ensemble structure, and the best representative structure, Fat10D2, was further considered for simulation. For chimeric models, the individual domains of Fat10D1, Fat10D2, and ubiquitin were fused between the N-terminal of one domain and the C-terminal of the other using UCSF Chimera (Pettersen et al., 2004). These models were further simulated for 500 ns before using the final structure for unfolding simulations. The salt bridge mutants were designed using the Rotamer function in UCSF Chimera. Systems for Fat10D1, Fat10D2, Fat10, and ubiquitin were prepared with the LEaP program of Ambertools18. Systems were prepared in a cubic box of TIP3p water, with a minimum distance of at least 12 Å between solute atoms and the box edge. Counter ions were added to neutralize the system.
Parameters describing system topology were based on the Amber ff99SBildn force field (Lindorff‐Larsen et al., 2010). The systems were first relaxed by energy minimization in two stages using the Sander module of Amber18. In the first stage, water molecules were minimized with restraint on protein, and then the entire system was minimized. The respective systems were then heated incrementally in NVT from 0K to 300K for 5 ns with positional restraints (20 kcal/mol/Å2) on protein atoms. Further, system density was equilibrated for 5 ns in the NPT ensemble with positional restraints (20 kcal/mol/Å2) on protein atoms. Further, four subsequent equilibration stages reduced the restraints on the backbone atom from 20 to 0 through a series of molecular dynamics simulations in an NPT ensemble for 400 ps each. The final production run was performed for 2.5 µs in NPT with three replicas. The distance cutoff for short-range nonbonded interactions was set to 1 nm. The particle mesh Ewald (PME) method was used to treat long-range electrostatic interactions. The SHAKE algorithm (Ryckaert et al., 1977) was applied to constrain all bonds involving hydrogen atoms. Temperature was set to 300 K using a Langevin thermostat, and pressure was maintained at 1 bar using the Berendsen barostat. Using the hydrogen mass repartitioning (HMR) scheme (Hopkins et al., 2015), the integration time step was set to 4 fs. Dynamics were propagated using the leapfrog integrator. Snapshots were saved every 40 ps, giving 65200 conformations from a single run. A total of 7.5 µs (3*2.5 µs) data was pooled for further analysis.
The trajectory analyses were performed using the AMBER suite’s CPPTRAJ module (Hipp et al., 2005). The averages and standard errors were calculated using in-house scripts and were plotted using the R program. Native contacts were calculated using the native-contacts method in CPPTRAJ with a 7 Å distance cutoff. For backbone native contacts across secondary structure pairs, the native contacts were calculated, defining a 3.5 Å distance cutoff on backbone atoms. The data was averaged across ten individual simulations at 450 K to calculate the mean and standard error. The Free Energy plots obtained from bin populations of the 2-dimensional histograms obtained from binning Root Mean Squared Deviation (RMSD) and Radius of Gyration values using the formula:
where ΔGi is the free energy value of bin i, kβ is Boltzmann’s constant in kcal/mol*K, T is the temperature in K, Ni is the population of bin i, and N0 is the population of the most populated bin.
RMSD and Rgyr were calculated using rms and rog methods available in CPPTRAJ with their respective crystal structure as a reference. The data were obtained from three individual 2.5 µs repeats at room temperature. Rgyr represents the Radius of Gyration for the protein without loops. The loops were omitted during the Free Energy calculation to avoid false positives. Rog represents the Radius of Gyration for the entire protein. Principle Component Analysis was performed on Ub, Fat10D1, Fat10D2, and their respective salt bridge mutants. PCA was also performed for Fat10D1-Ub, Fat10D2-Ub, and Ub-Fat10D2 using the Bio3D package Grant et al., 2006 in R Development Core Team, 2021. PDB format trajectory was further produced that interpolates between the most dissimilar structures in the distribution along PC1. To obtain contact maps, per-residue native and non-native contacts were calculated from trajectories, normalized, and plotted using GNUplot.
Far-UV Circular Dichroism (CD) measurements
Request a detailed protocolThe Far-UV CD data was collected using Jasco J-1500 spectropolarimeter in a 1 mm cylindrical quartz cuvette, using a 1 nm/s scan speed from 200 to 250 nm with a digital integration time of 4 s. Five scans at 25 °C were averaged and plotted using the SIGMA plot after converting the mdeg to MRE (Mean Residual Ellipticity). Appropriate buffer scans were subtracted from the protein’s scan.
NMR relaxation measurements
Request a detailed protocolFor the studies of the Fat10 dynamics, uniformly labeled 15N-Fat10, 15N-Fat10D1, 15N-Fat10D1-Ub, 15N-Fat10D2-Ub, 15N-Ub-Fat10D2, and 15N-ubiquitin were dialyzed in 50 mM Sodium phosphate buffer, 250 mM NaCl with pH 7.5. The NMR relaxation data were collected at 298 K on an 800 MHz and 600 MHz NMR spectrometer (Bruker). Longitudinal (T1), Transverse (T2) time constraints, and hetero-nuclear Overhauser Enhancement (hetNOE) experiments were carried out using standard pulse sequences in Bruker. For T1 measurements, data were recorded at the following relaxation delays: 0.004, 0.03, 0.06, 0.1, 0.15, 0.2, 0.4, and 0.8 s. For T2 measurements, the relaxation delays were set to 0.004, 0.017, 0.035, 0.051, 0.068, 0.086, 0.102, 0.119, and 0.136 s. A 5 s pre-saturation was used in the 15N-Heteronuclear NOE experiment. The reference experiment was carried out with 10 s delay without pre-saturation. The T1-relaxation delays for 15N-Fat10 were: 0.01, 0.025, 0.05, 0.075, 0.1, 0.2, 0.3 s, whereas the T2 delay time was: 0.01, 0.02, 0.4, 0.06, 0.08, 0.1 s and the hetNOE was recorded using standard pulse sequence in Bruker with 5 s as mixing time. The hetNOE was calculated as the ratio of intensities Isat/Iref. The hetNOE error was calculated as (NOE err/NOE) = [(Iref err/Iref)2+(Isat err/Isat)2]1/2, where the error of each intensity measurement is the rmsd noise of each plane. Order parameter (S2) was calculated using T1, T2, and hetNOE experiments in RELAX software (d’Auvergne and Gooley, 2008a; d’Auvergne and Gooley, 2008b) using the Lipari Szabo model-free approach where the S2, correlation time (τe), and Rex were obtained. The S2 values were converted to backbone conformational entropy using the entropy meter approach (Sharp et al., 2015; Zeng et al., 2017).
The temperature dependence of chemical shifts was measured for uniformly labeled 15N-Ub, 15N-Fat10D1, and 15N-Fat10D2 in 50 mM PO4 buffer with 250 mM NaCl at pH 7.5 on the 600MHz spectrometer. The temperature varied from 283K to 313K. The same was measured for 15N-Fat10, 15N-Fat10-Ub, and 15N-ubiquitin at the 800MHz spectrometer. Similar experiments were carried out in 50 mM PO4 buffer with 250 mM NaCl at pH 6.5 for the proteins 15N-Fat10D1-Ub, 15N-Fat10D2-Ub, and 15N-Ub-Fat10D2. The temperature varied from 283K to 323K. The chemical shift of water was used for reference and calibrated using a temperature-independent 4,4-dimethyl-4-silapentane-1-sulfonic acid (DSS) signal. All the NMR HSQC experiments were processed by NMRpipe and analyzed in SPARKY. The chemical shifts in the 1H-dimension were analyzed using linear regression by MATLAB. The linear fit between the ΔδNH and temperature provided the temperature coefficient of a residue, Tc = (ΔδNH/ΔT). The Residual Sum Square (RSS) was calculated as , where yi is the ith measured value and f(xi) is the fitted value of yi.
Appendix 1
Data availability
The manuscript and supporting file include all data generated or analyzed during this study.
References
-
The ubiquitin-like modifier FAT10 - much more than a proteasome-targeting signalJournal of Cell Science 133:jcs246041.https://doi.org/10.1242/jcs.246041
-
Studying the folding of multidomain proteinsHFSP Journal 2:365–377.https://doi.org/10.2976/1.2991513
-
Entropic contributions to protein stabilityIsrael Journal of Chemistry 60:705–712.https://doi.org/10.1002/ijch.202000032
-
Mechanically unfolding the small, topologically simple protein LBiophysical Journal 89:506–519.https://doi.org/10.1529/biophysj.105.061465
-
FAT10 is a proteasomal degradation signal that is itself regulated by ubiquitinationMolecular Biology of the Cell 23:225–232.https://doi.org/10.1091/mbc.E11-07-0609
-
Site-specific ubiquitination affects protein energetics and proteasomal degradationNature Chemical Biology 16:866–875.https://doi.org/10.1038/s41589-020-0556-3
-
A conserved and buried edge-to-face aromatic interaction in small ubiquitin-like modifier (SUMO) has A role in SUMO stability and functionThe Journal of Biological Chemistry 294:6772–6784.https://doi.org/10.1074/jbc.RA118.006642
-
N-terminal ubiquitination: more protein substrates join inTrends in Cell Biology 14:103–106.https://doi.org/10.1016/j.tcb.2004.01.004
-
Amide proton temperature coefficients as hydrogen bond indicators in proteinsJournal of Biomolecular NMR 21:249–261.https://doi.org/10.1023/a:1012911329730
-
A dynamic ubiquitin equilibrium couples proteasomal activity to chromatin remodelingThe Journal of Cell Biology 173:19–26.https://doi.org/10.1083/jcb.200510071
-
NMRPipe: A multidimensional spectral processing system based on UNIX pipesJournal of Biomolecular NMR 6:277–293.https://doi.org/10.1007/BF00197809
-
The dynamics of ubiquitin pools within cultured human lung fibroblastsThe Journal of Biological Chemistry 262:345–351.
-
FAT10, a ubiquitin-independent signal for proteasomal degradationMolecular and Cellular Biology 25:3483–3491.https://doi.org/10.1128/MCB.25.9.3483-3491.2005
-
Long-time-step molecular dynamics through hydrogen mass repartitioningJournal of Chemical Theory and Computation 11:1864–1874.https://doi.org/10.1021/ct5010406
-
Regulation of protein degradation by proteasomes in cancerJournal of Cancer Prevention 23:153–161.https://doi.org/10.15430/JCP.2018.23.4.153
-
MEGA X: molecular evolutionary genetics analysis across computing platformsMolecular Biology and Evolution 35:1547–1549.https://doi.org/10.1093/molbev/msy096
-
The ubiquitin code in the ubiquitin-proteasome system and autophagyTrends in Biochemical Sciences 42:873–886.https://doi.org/10.1016/j.tibs.2017.09.002
-
Conformational switching of the 26S proteasome enables substrate degradationNature Structural & Molecular Biology 20:781–788.https://doi.org/10.1038/nsmb.2616
-
Ubiquitylation directly induces fold destabilization of proteinsScientific Reports 6:39453.https://doi.org/10.1038/srep39453
-
Accelerating steered molecular dynamics: toward smaller velocities in forced unfolding simulationsJournal of Chemical Theory and Computation 12:1380–1384.https://doi.org/10.1021/acs.jctc.5b01024
-
A novel polyubiquitin chain linkage formed by viral Ubiquitin is resistant to host deubiquitinating enzymesThe Biochemical Journal 477:2193–2219.https://doi.org/10.1042/BCJ20200289
-
Probing the high energy states in proteins by proteolysisJournal of Molecular Biology 343:1467–1476.https://doi.org/10.1016/j.jmb.2004.08.085
-
UCSF Chimera--A visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Stability of proteins: small globular proteinsAdvances in Protein Chemistry 33:167–241.https://doi.org/10.1016/S0065-3233(08)60460-X
-
SoftwareR: a language and environment for statistical computingR Foundation for Statistical Computing, Vienna, Austria.
-
Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanesJournal of Computational Physics 23:327–341.https://doi.org/10.1016/0021-9991(77)90098-5
-
The 20S as a stand-alone proteasome in cells can degrade the ubiquitin tagNature Communications 12:26427-0.https://doi.org/10.1038/s41467-021-26427-0
-
FAT10ylation as a signal for proteasomal degradationBiochimica et Biophysica Acta (BBA) - Molecular Cell Research 1843:97–102.https://doi.org/10.1016/j.bbamcr.2013.01.009
-
Thermodynamic protein destabilization by GFP tagging: a case of interdomain allosteryBiophysical Journal 109:1157–1162.https://doi.org/10.1016/j.bpj.2015.04.032
-
Pupylation as a signal for proteasomal degradation in bacteriaBiochimica et Biophysica Acta (BBA) - Molecular Cell Research 1843:103–113.https://doi.org/10.1016/j.bbamcr.2013.03.022
-
Structures of UBA6 explain its dual specificity for ubiquitin and FAT10Nature Communications 13:4789.https://doi.org/10.1038/s41467-022-32040-6
-
SWISS-MODEL: homology modelling of protein structures and complexesNucleic Acids Research 46:W296–W303.https://doi.org/10.1093/nar/gky427
-
Thermodynamic contribution of backbone conformational entropy in the binding between SH3 domain and proline-rich motifBiochemical and Biophysical Research Communications 484:21–26.https://doi.org/10.1016/j.bbrc.2017.01.089
Article and author information
Author details
Funding
Tata Institute of Fundamental Research (RTI4006)
- Ranabir Das
Science and Engineering Research Board (CRG/2021/006032)
- Ranabir Das
Department of Biotechnology, Ministry of Science and Technology, India (dbt/pr12422/med/31/287/2014)
- Ranabir Das
- Aravind Ravichandran
Council of Scientific and Industrial Research, India
- Hitendra Negi
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
This work was supported by the Tata Institute of Fundamental Research, Department of Atomic Energy, Government of India, under project identification no RTI 4006, and Science and Engineering Research Board, under project identification no CRG/2021/006032. The NMR data were acquired at the NCBS-TIFR NMR Facility, supported by the Department of Atomic Energy, Government of India, under project no RTI 4006. The NMR facility is also partially supported by the Department of Biotechnology, India, under project number dbt/pr12422/med/31/287/2014. The authors would like to thank Andy Byrd (NCI), Jess Li (NCI), and Kylie Walters (NCI) for reagents and helpful discussions. HN acknowledges a fellowship from UGC-CSIR, India. AR and PD acknowledge scholarships from the Department of Biotechnology, India.
Copyright
© 2024, Negi et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 531
- views
-
- 88
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Neuroscience
- Structural Biology and Molecular Biophysics
More than 20 recurrent missense gain-of-function (GOF) mutations have been identified in the sodium-activated potassium (KNa) channel gene KCNT1 in patients with severe developmental and epileptic encephalopathies (DEEs), most of which are resistant to current therapies. Defining the neuron types most vulnerable to KCNT1 GOF will advance our understanding of disease mechanisms and provide refined targets for precision therapy efforts. Here, we assessed the effects of heterozygous expression of a Kcnt1 GOF variant (Kcnt1Y777H) on KNa currents and neuronal physiology among cortical glutamatergic and GABAergic neurons in mice, including those expressing vasoactive intestinal polypeptide (VIP), somatostatin (SST), and parvalbumin (PV), to identify and model the pathogenic mechanisms of autosomal dominant KCNT1 GOF variants in DEEs. Although the Kcnt1Y777H variant had no effects on glutamatergic or VIP neuron function, it increased subthreshold KNa currents in both SST and PV neurons but with opposite effects on neuronal output; SST neurons became hypoexcitable with a higher rheobase current and lower action potential (AP) firing frequency, whereas PV neurons became hyperexcitable with a lower rheobase current and higher AP firing frequency. Further neurophysiological and computational modeling experiments showed that the differential effects of the Kcnt1Y777H variant on SST and PV neurons are not likely due to inherent differences in these neuron types, but to an increased persistent sodium current in PV, but not SST, neurons. The Kcnt1Y777H variant also increased excitatory input onto, and chemical and electrical synaptic connectivity between, SST neurons. Together, these data suggest differential pathogenic mechanisms, both direct and compensatory, contribute to disease phenotypes, and provide a salient example of how a pathogenic ion channel variant can cause opposite functional effects in closely related neuron subtypes due to interactions with other ionic conductances.
-
- Structural Biology and Molecular Biophysics
Polycystin-1 (PC1) is the protein product of the PKD1 gene whose mutation causes autosomal dominant Polycystic Kidney Disease (ADPKD). PC1 is an atypical G protein-coupled receptor (GPCR) with an autocatalytic GAIN domain that cleaves PC1 into extracellular N-terminal and membrane-embedded C-terminal (CTF) fragments. Recently, activation of PC1 CTF signaling was shown to be regulated by a stalk tethered agonist (TA), resembling the mechanism observed for adhesion GPCRs. Here, synthetic peptides of the first 9- (p9), 17- (p17), and 21-residues (p21) of the PC1 stalk TA were shown to re-activate signaling by a stalkless CTF mutant in human cell culture assays. Novel Peptide Gaussian accelerated molecular dynamics (Pep-GaMD) simulations elucidated binding conformations of p9, p17, and p21 and revealed multiple specific binding regions to the stalkless CTF. Peptide agonists binding to the TOP domain of PC1 induced close TOP-putative pore loop interactions, a characteristic feature of stalk TA-mediated PC1 CTF activation. Additional sequence coevolution analyses showed the peptide binding regions were consistent with covarying residue pairs identified between the TOP domain and the stalk TA. These insights into the structural dynamic mechanism of PC1 activation by TA peptide agonists provide an in-depth understanding that will facilitate the development of therapeutics targeting PC1 for ADPKD treatment.