Acute unfolding of a single protein immediately stimulates recruitment of ubiquitin protein ligase E3C (UBE3C) to 26S proteasomes

The intracellular accumulation of aggregated misfolded proteins is a cytopathological hallmark of neurodegenerative diseases. However, the functional relationship between protein misfolding or aggregation and the cellular proteostasis network that monitors and maintains proteome health is poorly understood. Previous studies have associated translational suppression and transcriptional remodeling with the appearance of protein aggregates, but whether these responses are induced by aggregates or their misfolded monomeric or oligomeric precursors remains unclear. Because aggregation in cells is rapid, nonlinear, and asynchronous, it has not been possible to deconvolve these kinetically linked processes to determine the earliest cellular responses to misfolded proteins. Upon removal of the synthetic, biologically inert ligand shield-1 (S1), AgDD, an engineered variant FK506-binding protein (FKBP1A), rapidly (t½ ∼5 min) unfolds and self-associates, forming detergent-insoluble, microscopic cytoplasmic aggregates. Using global diglycine-capture (K-GG) proteomics, we found here that this solubility transition is associated with immediate increases in ubiquitylation of AgDD itself, along with that of endogenous proteins that are components of the ribosome and the 26S proteasome. We also found that the earliest cellular responses to acute S1 removal include recruitment of ubiquitin protein ligase E3C (UBE3C) to the 26S proteasome and ubiquitylation of two key proteasomal ubiquitin receptors, 26S proteasome regulatory subunit RPN10 (RPN10) and Rpn13 homolog (RPN13 or ADRM1). We conclude that these proteasomal responses are due to AgDD protein misfolding and not to the presence of detergent-insoluble aggregates.

To perform their biological functions, proteins must fold into and maintain a correct three-dimensional structure in a crowded intracellular environment (1). Moreover, the concentration of each protein must be precisely maintained by a carefully orchestrated balance of protein synthesis and degradation. This balance is maintained by the proteostasis network (PN), 5 composed of over 2000 proteins that comprise molecular chaperones, proteolytic systems, protein synthesis machinery, and regulators of these processes (2). Failure to maintain this homoeostatic balance leads to the deposition of aggregated proteins into intracellular inclusion bodies (IBs) in the nucleus and cytosol of diseased neurons (3). Although the nearly ubiquitous association of IBs in neurodegenerative diseases has long suggested a functional link between PN disruption and pathogenesis, the precise mechanistic relationship between protein misfolding and aggregation and neuropathology has been difficult to sort out. Indeed, there is no consensus whether IBs are toxic, neutral, or even cytoprotective (4).
Because the PN is a highly interconnected, multilayered network that controls all aspects of protein metabolism, its disruption has pleiotropic consequences that have been extremely challenging to deconvolve. Although conformationally compromised proteins accumulate in cells slowly and asynchronously and conformational diseases typically take decades to manifest, the conversion of a folding-compromised protein from a soluble to an insoluble state occurs rapidly and acutely both in cell culture (5) and in animal models (6) of disease. However, although genetically encoded conformational derangement (misfolding) of a single defined protein, such as huntingtin or tau, is sufficient to initiate a cascade of events leading ultimately to massive PN collapse (5,7,8), the asynchronous nature of protein misfolding and aggregation in animal and cell models of disease has hitherto precluded the ability to accurately identify the early (minute time scale) responses to conformationally induced PN challenge that precede the wellstudied transcriptional responses to protein aggregation.
In this study, we exploit a model of ligand-inducible protein misfolding and aggregation that permits highly synchronous, population-wide analysis of global PN perturbation at the proteome level with exquisite temporal precision (9). This approach uses cell lines that express variants of the FK506binding protein (FKBP1A) termed degradation domains, whose stability is absolutely dependent on the binding of a high-affinity synthetic ligand, Shield-1 (S1) (10,11). In the presence of S1, degradation domains are well-behaved, folded globular proteins. In the absence of ligand, degradation domains rapidly misfold and, when expressed in mammalian cells, are rapidly degraded by the ubiquitin proteasome system (UPS), together with any protein to which they are attached (10). This work exploits "AgDD," a variant degradation domain that differs from the widely used "DD" by the presence of an N-terminal decapeptide that is enriched in hydrophobic and aliphatic residues ( Fig. 1A) (9). Both DD and AgDD are conjugated to superfolder GFP (sfGFP), are monomeric, and are uniformly distributed in the cell when bound to their ligand, S1 (9,10). Upon acute removal of S1 from the culture medium ("washout" (w/o)), both proteins unfold, sampling different conformations and potentially transitioning dynamically among them (9,11). These partially unfolded proteins are targeted for ubiquitylation and, in the case of DD, rapid proteasomal degradation (t1 ⁄ 2 ϭ 90 min) (12,13). By contrast, S1 w/o causes AgDD to rapidly aggregate into puncta that are refractory to degradation and coalesce into cytoplasmic IBs (9). Our data indicate that one of the earliest cellular responses (Ͻ10 min) to acute protein misfolding is ubiquitylation of RPN13, one of the key Ub recep-tors on the apical cap of the 26S proteasome. We find that partially unfolded rather than aggregated AgDD induces recruitment of the Ub ligase, UBE3C to the proteasome as well as UBE3C-dependent multiubiquitylation of RPN13 and polyubiquitylation of proteasome-associated AgDD.

Ubiquitylation of key PN nodes in response to acute misfolding of a single protein
To study the earliest global changes to the PN that occur in response to acute protein folding stress, we exploited the ability of AgDD to unfold and aggregate following acute removal of the stabilizing ligand, S1 (9). To evaluate the suitability of AgDD for such global analysis, the spatial distribution of sfGFP fluorescence was monitored by time-lapse imaging of HEK293 cells stably expressing AgDD-sfGFP (heretofore, AgDD) following acute w/o of S1 (Fig. 1B). In the presence of ligand, AgDD fluorescence was diffuse, indistinguishable from that typically observed for unmodified GFP. Upon S1 w/o, AgDD rapidly shifted to a highly granular, punctate distribution observable at the earliest imaged time point, 2.5 min (Fig. 1C). Importantly, this pattern was observed in all cells at this time point, indicating that acute removal of S1 leads to simultaneous changes in AgDD distribution across an entire population of cells. Over the next 60 min, the puncta continued to increase in size and decrease in number, as reported previously (9), suggesting that individual AgDD IBs grow by coalescence of puncta. Coincident with this, diffuse AgDD fluorescence declined and was  Figure 1. Synchronous aggregation of AgDD induced by S1 w/o. A, AgDD-sfGFP and DD-sfGFP protein domains ("AgDD" and "DD" throughout for brevity). B, time course of AgDD partitioning into the Triton X-100 -insoluble fraction of cell lysates following S1 w/o. C, quantification of the data in C (n ϭ 4). D, time course of AgDD transition from diffuse to punctate subcellular distribution following S1 w/o. Error bars, S.D.
EDITORS' PICK: Unfolding-induced recruitment of UBE3C to 26S largely undetectable at later time points, suggesting that newly synthesized AgDD molecules, in the absence of S1, efficiently add onto pre-existing IBs. This transition of AgDD fluorescence from a diffuse to punctate appearance was accompanied by a rapid decline in the Triton X-100 solubility of AgDD, as assessed by sedimentation analysis (Fig. 1, C and B). The halftime for this solubility transition was ϳ5.2 min, with 63% insoluble after 10 min of S1 w/o. These data indicate that AgDD rapidly, reproducibly, and synchronously transitions into detergent-insoluble aggregates upon acute S1 w/o. This approach is, therefore, well-suited to conduct a systems-level interrogation of temporally resolved PN responses to acute unfolding and aggregation of a single, defined cytosolic protein.
To assess the immediate (i.e. Յ10 min) effects of acute AgDD unfolding on the PN that occur prior to transcriptional responses, we used an unbiased quantitative K-GG proteomic strategy to comprehensively assess the impact of S1 w/o on the cellular ubiquitylome ( Fig. 2A). Because the UPS tags conformationally altered proteins with Ub chains that direct them for degradation by the 26S proteasome (14), changes to the stability or ubiquitylation status of cellular proteins constitute a sensitive measure of proteome stress (15). K-GG proteomics exploits an antibody against the diglycine signature remnant on peptides generated from ubiquitylated proteins after complete digestion with trypsin (due to cleavage of the C-terminal Arg-Gly-Gly sequence of ubiquitin) (16,17). Thus, it allows global specific profiling of ubiquitin post-translational modification on proteins present in a cell population at a given time or under a given stress. By combining K-GG immunocapture with tandem mass tagging (TMT) technology (a combination of highly multiplexed isobaric tags), it is possible to obtain precise timeresolved assessments of changes to the ubiquitylome in a single experiment. Lysates from AgDD cells subjected to S1 w/o for 0, 2.5, 5, 7.5, and 10 min ( Fig. 2A) were digested with trypsin, and K-GG immunoprecipitates from each time point were isolated, labeled with 11-plex isobaric tags, pooled, fractionated, and analyzed by nLC-MS/MS as described previously (18). In total, we identified and quantified 4,549 common unique sites of ubiquitylation across all 11 experimental conditions in two independent experiments (Table S1). Parallel analyses of control samples from WT HEK293 cells not expressing AgDD cultured in the presence or absence of S1 or subjected to S1 w/o were performed to ensure that in the absence of AgDD, S1 and S1 w/o have no significant biological impact on protein ubiquitylation ( Fig. S1 (B and C) and Table S2).
Unsurprisingly, in AgDD cells, the earliest and most substantial changes in ubiquitylation (3.8 -5.2-fold increase after 10-min S1 w/o) were detected on AgDD itself, with three sites (Lys-35, Lys-48, and Lys-53) significantly modified ( Fig. 2B; note that ubiquitylated AgDD lysines are numbered according to their position in native FKBP1A). This result is consistent with previous publications showing that the DD domain undergoes partial unfolding and becomes rapidly ubiquitylated following S1 w/o (11).
In addition to ubiquitylation of AgDD itself, significant changes to the endogenous ubiquitylome were also observed following the acute 10-min S1 w/o, with some sites exhibiting increased or decreased levels of ubiquitylation and many sites trending monotonically over the time course ( Fig. 2C and Fig.  S1A). At 10 min, we identified 57 ubiquitylation sites corresponding to 48 proteins exhibiting statistically significant changes, including 42 sites with increased, and 15 sites with decreased modification (Fig. 2, B and C). The sites exhibiting significantly decreased ubiquitylation correspond to a group of cytosolic and organellar proteins with no obvious functional relationship to one another or to known stress-response pathways. By contrast, the 42 sites exhibiting statistically significant increased ubiquitylation in response to S1 w/o are strikingly enriched in functionally related PN proteins, including chaperones, proteasome subunits, ribosomal proteins, and ribosome biogenesis factors, suggesting that acute unfolding and aggregation of a single abundant protein leads to rapid (Յ10-min) Ub modification of both protein degradation and protein synthesis machinery. In addition to proteins with documented roles in the PN, which will be discussed below, our analysis identified a significant increase in PARP4 ubiquitylation that increased monotonically across all time points following S1 w/o. PARP4 is a cytosolic poly-ADP ribosylase and one of three protein components of mammalian vault complexes. Purified vaults demonstrate ADP-ribosylase activity, confirming that PARP4 is enzymatically active in vaults (19). Whereas many PARP proteins are active in stress-response pathways (20), the relationship of PARP4 and potentially vault-associated PARP4 to the PN remains unknown.
The most striking early changes to the ubiquitylome in response to S1 w/o were increased Ub modification of components of the UPS itself. We observed a significant monotonic increase in the modification of Ub at Lys-29 and Lys-33, suggesting that increased production or decreased turnover of these noncanonical Ub chains is an early and progressive response to acute unfolding and aggregation of AgDD. K-GG proteomic profiling also revealed that S1 w/o induces immediate 2-3-fold increases in Ub modification of the proteasome itself, notably on the Ub receptor, RPN13 (ADRM1), at two sites (Lys-21 and Lys-34) that were second in magnitude only to modification of AgDD itself. A single site in RPN10 (PSMD4, Lys-40) also exhibited a modest (1.5-fold) increase in ubiquitylation, reaching statistical significance only at 10 min following S1 w/o. These data reveal that acute unfolding and/or aggregation of a single highly expressed cytosolic protein leads to immediate and progressive increases in ubiquitylation of key nodes in the PN network, particularly the ribosome and the 26S proteasome.

AgDD unfolding/aggregation do not impair global 26S proteasome activity
Acute unfolding and aggregation of AgDD could impair proteasome function by overwhelming substrate binding capacity or by "choking" the proteasome with hard-to-unfold aggregates (5,21,22). Moreover, ubiquitylation of 19S Ub receptors RPN13 and RPN10 could interfere with binding and degradation of UPS substrates, leading to their accumulation. Indeed, nearly stoichiometric modification of RPN13 through in vitro ubiquitylation of isolated proteasomes has been reported to impair substrate binding and degradation (23,24), although it is unclear whether those effects are entirely due to modification  of RPN13 as opposed to other proteasome subunits. To determine whether the observed increase in endogenous protein ubiquitylation induced by AgDD aggregation is the consequence of inhibited proteasome activity, we evaluated cleavage of a bioluminescent peptide reporter of chymotrypsin-like proteasome activity in cleared lysates from AgDD cells subjected to S1 w/o. Intriguingly, instead of a decrease, we observed a 40% increase in chymotrypsin-like activity relative to mock-treated controls after a 10-min S1 w/o (Fig. 3A), suggesting that unfolding and aggregation of AgDD stimulates, rather than inhibits, peptidase activity of the 20S proteasome. Despite this curious effect on 20S peptidase activity, S1 w/o impacted neither the steady-state level nor the degradation kinetics of core-glycosylated CD147 (CD147 (CG); Fig. 3B), an endogenous proteasome substrate whose degradation is strictly dependent on Ub conjugation and on 26S proteasome function (25). To assess whether acute AgDD unfolding causes global proteasome impairment, we compared the impact of acute S1 w/o on the K-GG ubiquitylome of AgDD cells with the published (17) effects of proteasome inhibition by bortezimib (Btz) in HCT116 cells. Overall, there was very poor (r 2 Ͻ 0.11) correlation between the ubiquitylation sites altered at 10 min following S1 w/o and published Btz-responsive sites (Fig. 3C). Among the 57 significantly altered ubiquitylation sites identified in our K-GG proteomic analysis, only five-nuclear complex protein 2 ho-molog (NOCL2, Lys-516), 40s ribosomal protein S25 (RPS25, Lys-66), 40s ribosomal protein S16 (RPS16, Lys-50), catenin ␤-1 (CTNNB1, Lys-233), and RPN13 (ADRM1, Lys-21)-were also increased following direct proteasome inhibition. The poor correlation between the S1 w/o and Btz-induced changes to the K-GG ubiquitylome suggests that acute S1 w/o does not lead to globally impaired proteasome function. Together, these data indicate that the altered ubiquitylation status of proteins following S1 w/o might reflect activation of a cellular stress response and not a secondary consequence of proteasome impairment.

AgDD and RPN13 are ubiquitylated in response to S1 w/o
To directly assess the effects of S1 w/o on RPN13 ubiquitylation, we used a mixture of Halo-tagged Ub-association domain (UBA) reagents that bind Ub conjugates without preference for chain linkage topology (26,27) to enrich ubiquitylated proteins from detergent lysates of AgDD cells (Fig. 4). GFP immunoblot analysis of Halo-UBA-captured proteins from mock-treated AgDD cells revealed a faint "ladder" of bands migrating more slowly than unmodified AgDD (Fig. 4A). S1 w/o caused a substantial increase in the intensity of this ladder and led to the appearance of a prominent high-molecular weight (HMW) GFP-immunoreactive smear. All of these species, both in mock-treated and S1 w/o cells, were eradicated by pretreat- Figure 2. S1 w/o in AgDD cells rapidly alters the ubiquitylation of PN components. A, experimental design and workflow: TMT labeling of S1 w/o time course sample and replicate. B, volcano plot (Ϫlog 10 p value versus log 2 ratio to untreated cells) of K-GG sites at 10 min following S1 w/o. Proteins with log 2 ratios ϽϪ0.585 or Ͼ0.585 (p Ͻ 0.01) are indicated as colored empty circles, and filled colored circles indicate statistically significant hits (Welch's t test (S0 ϭ 0.5), corrected for multiple comparison by permutation-based FDR (1%)). C, heat map of log 2 scaled changes in K-GG sites at all time points for sites determined to be statistically altered at 10 min. Gene names are listed with select protein names in parentheses. Ubiquitylated AgDD lysines are numbered according to their position in native FKBP1A. The numbering of these sites in AgDD is shifted relative to FKBP1A by the N-terminal decapeptide, becoming Lys-45, Lys-58, and Lys-63.   ). B, S1 w/o has no detectable impact on degradation of core-glycosylated CD147 (CD147 (CG)). Lysates of AgDD cells treated with emetine for the indicated times were assessed by immunoblot analysis with an antibody to endogenous CD147 (left) and quantified (right). C, site-specific ubiquitylation patterns induced by 2-, 4-, and 8-h Btz treatment in HCT116 cells (17) as compared with 10-min S1 w/o in AgDD HEK293 cells. Sites significantly (0.2% FDR (17) and 1% FDR (this work)) altered by both S1 w/o and Btz treatment are colored red and labeled.
EDITORS' PICK: Unfolding-induced recruitment of UBE3C to 26S ment of lysates with the deubiquitylase USP2-cc (28) (Fig. 4B), demonstrating that they correspond to multiubiquitylated and polyubiquitylated forms of AgDD. As the proteomic analysis revealed substantially increased abundance of the same three K-GG peptides following S1 w/o, rather than the appearance of additional K-GG peptides, it is likely that the HMW ubiquitylated species produced in response to S1 w/o are produced by extension of existing oligoubiquitin species on AgDD. These data demonstrate that AgDD, like its nonaggregating counterpart, DD (11,13), is rapidly and extensively polyubiquitylated in response to removal of the synthetic stabilizing ligand S1. Immunoblot analysis of Halo-UBA-captured proteins with RPN13 antibody also identified a ladder of slower migrating species in AgDD cells that increased in intensity following S1 w/o (Fig. 4C) and were lost following pretreatment with USP2-cc (Fig. 4D). Unlike AgDD, we did not observe evidence of HMW polyubiquitin chains on RPN13 in response to S1 w/o.
Together with the K-GG proteomic analysis, we conclude that basal ubiquitylation of RPN13 is augmented by acute unfolding and/or aggregation of AgDD. In contrast to AgDD and RPN13 ubiquitylation, we were unable to detect a slower-migrating form of RPN10 by immunoblotting among affinity-captured ubiquitylated conjugates or associated with affinity-captured proteasomes, suggesting that the stoichiometry of this modification is below the threshold for immunodetection (Fig. S2A).
Previous studies have suggested a role for the HECT-domain Ub ligase UBE3C in ubiquitylating proteasome-bound RPN13 in response to a variety of acute stresses (24). Like RPN13, UBE3C is a nonstoichiometric 26S proteasome subunit that partitions between cytosolic and proteasome-bound pools (31). The yeast UBE3C ortholog, Hul5, was previously shown to play a central role in ubiquitin-dependent degradation of low-solubility cytosolic proteins (32) and to increase the processivity of stalled proteasomes (33). To assess the role of UBE3C in ubiquitylating RPN13 in response to S1 w/o, we used siRNA to deplete UBE3C from AgDD cells expressing RPN11-HTBH and assessed both total (Halo-UBA) and proteasome-associated (RPN11-HTBH PD) RPN13 by immunoblot analysis of the affinity-captured eluates. UBE3C-targeted siRNAs reduced UBE3C expression by ϳ95% in Triton X-100 lysates and reduced the amount of UBE3C captured by RPN11-HTBH PDs by ϳ75% (Fig. 5C). UBE3C knockdown blocked S1 w/o-induced increases of total and proteasome-associated RPN13-Ub n conjugates but did not reduce basal RPN13 ubiquitylation (Fig. 5A). These data suggest that S1 w/o activates UBE3C-dependent RPN13 multiubiquitylation as reported previously for pleiotropic stresses, such as heat shock and pharmacological proteasome inhibition (23,24).
The effect of UBE3C knockdown on AgDD ubiquitylation differed markedly from that observed for RPN13 (Fig. 5B). As with RPN13, AgDD-Ub n conjugates were co-enriched with RPN11-HTBH-captured proteasomes. Whereas nonubiquitylated AgDD binds nonspecifically to streptavidin beads, Ub n -AgDD capture by RPN11-HTBH depends upon association with 26S proteasomes (Fig. S4A). However, unlike with RPN13, UBE3C knockdown resulted in increased levels of AgDD multiubiquitylation in the presence of S1, suggesting that UBE3C is dispensable for basal ubiquitylation of folded AgDD (Fig. 5B). Importantly, UBE3C deletion prevented the enhanced multiubiquitylation and HMW polyubiquitylation in response to S1 w/o (Fig. 5B). These data suggest that UBE3C can extend oligoubiquitin chains on AgDD, perhaps by virtue of the E4 activity reported for Hul5 (34).

UBE3C is recruited to the 26S proteasome in response to AgDD unfolding
Increased RPN13 ubiquitylation in response to S1 w/o could result from either globally increased UBE3C activity or from recruitment of UBE3C to the proteasome from a cytosolic pool. To distinguish between these possibilities, we used RPN11-HTBH capture to assess the impact of S1 w/o on UBE3C proteasome association (Fig. 5C). The amount of UBE3C that coprecipitated with RPN11-HTBH increased ϳ4-fold following S1 w/o, closely matching the magnitude of increase in total RPN13 ubiquitylation (Ub 1-4 ) under the same conditions ( EDITORS' PICK: Unfolding-induced recruitment of UBE3C to 26S 5C). At the same time, S1 w/o increased the fraction of proteasome-associated UBE3C from ϳ25 to 60% (Fig. S3, A and B). Thus, S1 w/o in AgDD cells stimulates UBE3C recruitment to the proteasome from the cytosol, leading to increased RPN13 ubiquitylation and AgDD polyubiquitylation.
In principle, UBE3C recruitment and consequent ubiquitylation events in response to S1 w/o could be due to either unfolding or aggregation of AgDD. Because AgDD aggregation occurs rapidly upon S1 w/o, it is not possible to kinetically determine whether these events are the consequence of the sudden appearance of a large bolus of a partially unfolded, aggregation-prone protein or its subsequent aggregation. To distinguish between these possibilities, we assessed the effects on the proteasome of S1 w/o in HEK293 cells stably expressing DD, a variant of FKBP1A that differs from AgDD only by the absence of an N-terminal decapeptide (Fig. 1A). In the absence of S1, DD adopts a partially folded "molten globule"-like state (11) that is efficiently degraded by the UPS (10) without generating fluorescent puncta, detergent-insoluble aggregates, or IBs (10,11,13). Remarkably, we found that S1 w/o from DD cells recapitulated all of the proteasome-associated effects that we observed in AgDD cells, including UBE3C recruitment (Fig. 6, A and B), RPN13 ubiquitylation (Fig. 6, A and C), enhanced 20S peptidase activity (Fig. 6D), and formation of HMW Ub conjugates (Fig. 6, A and E), although these events proceed with slower kinetics in DD compared with AgDD cells. Together, these data suggest that the sudden generation of partially unfolded AgDD and DD proteins is sufficient to drive UBE3C recruitment and activity increases at the proteasome.
Because these effects appear to occur more rapidly in AgDD compared with DD cells, it is possible that detergent-insoluble AgDD microaggregates could contribute to these responses. Alternatively, it is possible that the enhanced kinetics in AgDD cells simply reflect AgDD's aggregation propensity or the more efficient degradation of ubiquitylated DD. To distinguish among these possibilities, we used an S1 addback approach (9) to arrest AgDD unfolding by exchanging cells that had been subjected to 10-min S1 w/o into a fresh medium "chase" containing S1 and incubating for an additional 15 min (Fig. 7). Without S1 re-addition, this additional 15-min chase substantially increased the fraction of detergent-insoluble AgDD (Fig.  7A; compare lanes 5 and 6 with lanes 3 and 4). In contrast, the fraction of insoluble AgDD in cells that had been chased in the presence of S1 was comparable with the fraction in cells subjected to 10-min S1 w/o alone (Fig. 7A; compare lanes 7 and 8 Figure 5. UBE3C-catalyzed ubiquitylation of RPN13 and AgDD at the proteasome. A, RPN13-Ub n captured by UBA PD (left) or co-precipitated with affinity-captured proteasomes (right) increases ϳ4-fold following 10-min S1 w/o. UBE3C knockdown blocks S1 w/o-stimulated increase in RPN13-Ub n . B, AgDD-Ub n captured by UBA PD (left) or co-precipitated with affinity-captured proteasomes (right). UBE3C knockdown prevents the formation of only HMW-AgDD-Ub n that associates with the proteasome. C, proteasome-associated UBE3C increases ϳ4-fold following S1 w/o. Error bars, S.D.

EDITORS' PICK: Unfolding-induced recruitment of UBE3C to 26S
with lanes 3 and 4). Blocking ubiquitylation with an inhibitor of the Ub-activating enzyme E1 did not influence either the kinetics of aggregation or the effect of S1 re-addition (Fig. 7A, lanes  9 -12); nor did inclusion of a proteasome inhibitor (Fig. 7A,  lanes 13 and 14). These data confirm the previous report (9) that S1 re-addition prevents new aggregate formation but does not cause dissolution of already-formed aggregates. S1 re-addition also prevented the formation of additional proteasomeassociated RPN13-Ub n and HMW AgDD-Ub n (Fig. 7, B and C), suggesting that these events are responses to the presence of partially unfolded AgDD and not to the presence of AgDD aggregates (7, B and C). By contrast, inhibiting E1 during the chase led to a decrease in multiubiquitylated RPN13 and to an increase in the mobility, but not the abundance, of polyubiquitylated HMW AgDD. These observations suggest that proteasome-associated RPN13 and AgDD are subject to rapid de-ubiquitylation, a conclusion that is supported by the lack of effect of MG132. Taken together, these data suggest that ubiquitylation of RPN13 and UBE3C recruitment to the proteasome are the consequence of acute and rapid unfolding of a single abundant cytosolic protein and not a response to aggregation.

Discussion
Much of our understanding of cellular protein homeostasis networks comes from studies in which cells are subjected to insults like heat shock (35), Hsp90 inhibitors (36), or chemical denaturants or oxidants (37,38) that cause immediate disruption to global protein conformation. Although such approaches have been instrumental in revealing the existence of homeostatic networks that promote stereotypical, temporally orchestrated restoration of homeostasis, animal cells are rarely, if ever, subjected to such acute, massive, proteome-wide insults. Instead, most conformational diseases arise from the production of a single, conformationally impaired, aggregation-prone protein, such as the mutant proteins linked to genetic diseases (39). Although cell culture models expressing such mutant proteins have been val-uable tools in examining the impact on the PN of producing a single conformationally defective protein, the slow and asynchronous nature of protein expression and aggregation has limited the utility of these models to investigate the temporal orchestration of the PN responses. Consequently, there is a substantial gap in our understanding of the relationship between acute global insults like heat shock and the more physiologically relevant but slow PN responses to the presence of a single conformationally impaired protein. In this study, we have combined a systems-level analysis of PN perturbation with a chemical biological approach to define the earliest cellular responses to acute unfolding and aggregation of a single, defined protein. Our data reveal that cells respond to acute misfolding of a single protein with broad changes to protein synthesis and protein degradation  The most striking response of the PN to the induced unfolding and aggregation of AgDD was the ubiquitylation of two proteasomal Ub receptors, RPN10 and RPN13. Both proteins are located near the apex of the 19S regulatory particle (RP) and are proposed to facilitate cooperative binding of polyubiquitylated conjugates in a manner that optimally orients them for downstream processing and degradation (40 -43). Although the identification of RPN10 is intriguing, because of the borderline statistical significance of this identification in our K-GG analysis and our inability to biochemically confirm this modification, we focused our attention on ubiquitylation of RPN13, a protein that binds to the 19S regulatory particle via an interaction between its N-terminal domain and RPN2 (44). RPN13 also binds and activates the deubiquitylase UCH37 via its C-terminal domain (45)(46)(47). The two K-GG-modified lysines on RPN13 map to its Ubbinding surface (48). Whereas ubiquitylation of both RPN10 and RPN13 has been previously reported to occur in response to a number of stressors, including heat shock, acute oxidation, and proteasome inhibition (23,24), our data indicate that the 19S proteasome is a major target of early ubiquitylation events stimulated by protein unfolding.
Our study establishes a key role for UBE3C in modifying both RPN13 and AgDD. It is noteworthy that yeast RPN10 is monoubiquitylated by Hul5 (49). UBE3C and Hul5 possess E4 activities that can extend Ub chains on substrates bound to 26S proteasomes, perhaps antagonistically with proteasome-bound deubiquitylases (34); this antagonism is proposed to facilitate degradation of difficult-to-degrade proteins (13). UBE3C was previously identified in an RNAi screen for DD degradation (13). However, the modest stabilization of DD observed upon UBE3C depletion suggests that it does not act alone and, together with the findings presented here, suggests a model wherein partially unfolded clients like DD and AgDD are first multiubiquitylated in the cytosol by a yet unidentified E3 and subsequently polyubiquitylated by the E4 activity of 26S-associated UBE3C.
UBE3C can extend Lys-48 -linked Ub chains on substrates with Lys-29 linkages (50,51). In this regard, it is noteworthy that the K-GG analysis identified a significant increase in doubly modified Lys-29, Lys-33 Ub-Ub linkages in response to S1 w/o. UBE3C-dependent RPN13 ubiquitylation is one of the most robust and immediate changes to global ubiquitylation following S1 w/o in both AgDD and DD cells (Fig. 2, B and C), occurring simultaneously with proteasomal recruitment of UBE3C (Fig. 6B). These rapid kinetics pose intriguing questions about the mechanism by which UBE3C "senses" protein unfolding and the means by which it rapidly redistributes to the proteasome to aid in degradation.
Our finding that the changes to the proteasome in response to S1 w/o in AgDD cells are fully recapitulated in cells expressing aggregation-resistant DD strongly suggests that protein unfolding is sufficient to induce a proteasome stress response. This conclusion is supported by the observation that S1 re-addition following 10-min S1 w/o results in rapid deubiquitylation of both RPN13 and proteasome-associated AgDD, despite the persistence of detergent-insoluble AgDD aggregates. Moreover, these data suggest that the ubiquitylation states of both RPN13 and AgDD are dynamic, maintained by antagonism between 26S-associated Ub conjugation and deubiquitylation activities. Whereas these observations are consistent with previous data suggesting antagonism between Hul5 and the proteasome-associated deubiquitylase, Ubp6 (34), additional studies will be needed to identify the deubiquitylases that maintain this equilibrium for proteins like AgDD, which are refractory to degradation.

Experimental procedures
Cell culture and S1 w/o HEK293 cells stably expressing AgDD-sfGFP or DD-sfGFP were cultured in DMEM supplemented with 10% heat-inactivated FBS and 1 M S1. S1 w/o was performed by adding purified 500 M recombinant FKBP1A F36V protein (rDD, diluted in PBS ϩ 4% glycerol) (9) to the media of cultured DD or AgDD cells to create a final concentration of 5 M. Mock washouts were performed by adding an equal volume of PBS ϩ 4% glycerol to cell culture media.

UBA affinity capture
Purified Halo-TUBE Dsk2 and Halo-TUBE ubiquilin-1 (52) were separately conjugated to Magne HaloTag beads (Promega) by rotating at room temperature. Following mock or S1 w/o, cells were washed in ice-cold PBS, harvested by scraping, and stored at Ϫ80°C. For UBA pulldowns, cells were lysed in UBA lysis buffer (50 mM Tris, pH 7.5, 150 mM NaCl, 1% Nonidet P-40, 50 mM chloroacetamide, 2 mM NEM, 2.5 mM EDTA, 1 mM PMSF, cOmplete TM EDTA-free protease inhibitor mixture (Roche Applied Science)). 45 l of both 10% bead-TUBE slurries was added to 500 g of cleared AgDD cell lysates and rotated overnight at 4°C. Beads were then washed four times with 1 ml of wash buffer (50 mM Tris, pH 7.5, 500 mM NaCl, 1% Nonidet P-40, 50 mM chloroacetamide, 2 mM NEM, 2.5 mM EDTA, 1 mM PMSF, cOmplete TM EDTA-free protease inhibitor mixture (Roche Applied Science)). Protein was eluted from these beads with two sequential 10-min incubations at 95°C with 25 l of 2ϫ protein sample buffer, resolved on SDS-PAGE gels, and analyzed by Western blotting using the indicated antibodies.

Proteasome affinity capture
Cells were lysed in proteasome PD lysis buffer (50 mM sodium phosphate, pH 7.5, 100 mM NaCl, 1% Nonidet P-40, 10% glycerol, 5 mM ATP, 5 mM MgCl 2 , cOmplete TM EDTA-free protease inhibitor mixture (Roche Applied Science)). Lysates were cleared by centrifugation at 13,000 ϫ g for 15 min at 4°C. Cleared lysates were then combined with magnetic streptavidin resin (10 g of total protein/1 l of resin) (Pierce) and rotated for 30 min at 4°C. Resin was washed twice in lysis buffer (vortexing to suspend resin in each wash). Protein was eluted from resin with two sequential incubations with 25 l of 2ϫ protein sample buffer at 95°C for 10 min, resolved on SDS-polyacrylamide gels, and analyzed by Western blotting using the indicated antibodies.

20S proteasome activity assays
Cleared lysates from AgDD cells were normalized by total protein. Samples of cleared lysates were combined in a 1:1 ratio with reconstituted Proteasome Glo TM (Promega) solution, and bioluminescence was detected in a BioTek TM Synergy TM NEO HTS multimode microplate reader.
Solubility shift assays S1 w/o was performed in AgDD cells for the indicated time periods. Cells were harvested by aspirating media and immediately washing cells with ice-cold PBS. Cells were held on ice while PBS was aspirated and cells were collected by scraping. Cells were lysed in ice-cold buffer containing 50 mM Tris, pH 7.5, 150 mM NaCl, 1% Triton X-100, 2 g/ml aprotinin, 1 g/ml pepstatin, 10 mM NEM. The lysate was rotated at 4°C for 15 min and cleared by centrifugation at 20,000 ϫ g for 30 min at 4°C. The resulting pellet was washed two times with lysis buffer and repelleted after each wash at 20,000 ϫ g. After all supernatant was removed, the pellet was solubilized in buffer containing 50 mM Tris, 1% SDS, 150 mM NaCl, 2 g/ml aprotinin, 1 g/ml pepstatin, 10 mM NEM. S1 w/o pulse and S1 addback chase S1 w/o was performed on AgDD cells stably expressing Rpn11-HTBH. At 10 min following S1 w/o, medium was aspirated and replaced with DMEM alone or DMEM containing either 2 M S1, 10 M E1 inhibitor (TAK-243, MedChem Express), both together, or both together in combination with 10 M MG132. All replacement medium was warmed to 37°C prior to the exchange. Following medium replacement, cells were incubated at 37°C for an additional 15 min and then harvested and processed according to the proteasome affinity capture protocol (see above).

Live-cell imaging
Live-cell imaging was performed on AgDD cells, plated in DMEM with 10% heat-inactivated FBS, and incubated at 37°C and 5% CO 2 in a glass-bottom 8-well chambered coverglass (Nunc). For time-lapse imaging experiments, S1 w/o was performed immediately prior to the start of imaging. Cells were imaged by time-lapse fluorescence video microscopy using an inverted microscope (Zeiss Axiovert 200M) encased in a Perspex chamber heated to 37°C. The chambered coverglass was placed in a secondary chamber perfused with humidified 5% CO 2 . A motorized stage allowed imaging of multiple fields. Images were acquired with a cooled charge-coupled device (Cool-SNAP HQ) and UV light source (X-Cite; Lumen Dynamics), filters for visualization of GFP, and a ϫ20 air objective. MetaMorph software was used to control image acquisition, and ImageJ (version 1.6, National Institutes of Health) was used for image analysis.

Emetine chase
AgDD cells were left untreated (0-h time point) or treated with final concentrations of either 5 M rDD (S1 w/o) or vehicle and 1 M emetine for the indicated times before collection in cold PBS. Cells were pelleted with a tap spin and lysed in 1% Triton X-100, 50 mM Tris, pH 7.5, 150 mM NaCl, 1 M pepstatin A, and 1 mM PMSF, and lysates were cleared by centrifugation at 20,000 ϫ g. Protein concentration was measured using a BCA Protein Assay Kit (Thermo Fisher Scientific), and concentrations were normalized prior to SDS-PAGE and immunoblotting. Band intensities were quantified using Image Studio Lite software (LI-COR Biosciences). The protein remaining was calculated as a percentage of untreated, and one-phase exponential decay curves were fit using Prism version 7 (GraphPad Software).

K-GG proteomics
Protein isolation-S1 w/o was performed in AgDD cells for the indicated time period. At the indicated times, cells were washed twice with ice-cold PBS, collected, and snap-frozen. Cell pellets were lysed in lysis buffer (50 mM EPPS, pH 8.0, 8 M urea, 0.25% (v/v) Triton X-100, 1 mM 4-(2-aminoethyl)benzenesulfonyl fluoride hydrochloride, 10 M PR-619, 0.5 mM tris(2-carboxyethyl)phosphine, 50 mM chloroacetamide), to produce whole-cell extracts. Whole-cell extracts were sonicated and clarified by centrifugation (16,000 ϫ g for 10 min at 4°C) followed by filtration through a 0.45-m filter, and protein concentrations were determined by the Bradford assay.
General sample preparation-Protein extracts (1.5 mg) were subjected to disulfide bond reduction with 5 mM tris(2-carboxyethyl)phosphine (room temperature, 10 min) and alkylation with 25 mM chloroacetamide (room temperature, 20 min). Methanol-chloroform precipitation was performed prior to protease digestion. In brief, four parts undiluted methanol were added to each sample and vortexed, one part chloroform was then added to the sample and vortexed, and finally three parts water were added to the sample and vortexed. The sample was centrifuged at 6,000 rpm for 2 min at room temperature and subsequently washed twice with 100% methanol. Samples were resuspended in 100 mM EPPS, pH 8.5, containing 6 M urea and digested at 37°C for 2 h with LysC protease at a 200:1 protein/ protease ratio. Samples were then diluted with 100 mM EPPS, pH 8.5, to a final urea concentration of 0.75 M. Trypsin was then added at a 100:1 protein-to-protease ratio, and the reaction was incubated for a further 6 h at 37°C. Samples were acidified with 1% formic acid for 15 min and subjected to C18 solid-phase EDITORS' PICK: Unfolding-induced recruitment of UBE3C to 26S extraction (Sep-Pak, Waters). The Pierce quantitative colorimetric peptide assay (catalog no. 23275) was used to quantify the digest and to accurately aliquot 1 mg of peptides for diGLY enrichment.
Immunoprecipitation of diGLY-containing peptides-diGLY capture was performed largely as described (18). The diGly mAb (Cell Signaling Technology; D4A7 clone) (32 g of antibody/1 mg of peptide) was coupled to Protein A Plus Ultralink resin (1 l slurry to 1 g of antibody) (Thermo Fisher Scientific) overnight at 4°C prior to its chemical cross-linking reaction. Dried peptides (1 mg) were resuspended in 1.5 ml of ice-cold IAP buffer (50 mM MOPS (pH 7.2), 10 mM sodium phosphate, and 50 mM NaCl) and centrifuged at maximum speed for 5 min at 4°C to remove any insoluble material. Supernatants (pH ϳ7.2) were incubated with the antibody beads for 2 h at 4°C with gentle end-over-end rotation. After centrifugation at 215 ϫ g for 2 min, beads were washed three more times with ice-cold IAP buffer and twice with ice-cold PBS. The diGLY peptides were eluted twice with 0.15% TFA, desalted using homemade StageTips, and dried via vacuum centrifugation, prior to TMT labeling.
diGLY proteomics analysis using TMT-TMT-labeled diGLY peptides were fractionated according to the manufacturer's instructions using a High pH reversed-phase peptide fractionation kit (Thermo Fisher Scientific) for a final six fractions and subjected to C18 StageTip desalting prior to MS analysis. Mass spectrometry data were collected using an Orbitrap Fusion Lumos mass spectrometer (Thermo Fisher Scientific, San Jose, CA) coupled to a Proxeon EASY-nLC1200 liquid chromatography (LC) pump (Thermo Fisher Scientific). Peptides were separated on a 100-m inner diameter microcapillary column packed in-house with ϳ35 cm of Accucore150 resin (2.6 m, 150 Å, Thermo Fisher Scientific) with a gradient consisting of 5-26% (0 -155 min), 26 -32% (155-170 min) (acetonitrile, 0.1% formic acid) over a total 180-min run at ϳ600 nl/min. For analysis, we loaded one-half of each fraction onto the column. Each analysis used the Multi-Notch MS 3 -based TMT method (53). The scan sequence began with an MS 1 spectrum (Orbitrap analysis; resolution 120,000 at 200 Th; mass range 400 -1250 m/z; automatic gain control (AGC) target 1 ϫ 10 6 ; maximum injection time 100 ms). Precursors for MS 2 analysis were selected using a top 4 second method. MS 2 analysis consisted of collision-induced dissociation (quadrupole Orbitrap analysis; AGC 1 ϫ 10 5 ; isolation window 0.8 Th; normalized collision energy 35; maximum injection time 300 ms; resolution was 15,000 at 200 Th). Monoisotopic peak assignment was used, and previously interrogated precursors were excluded using a dynamic window (120 s Ϯ 10 ppm). As described previously, only precursors with a charge state between 3 and 6 were selected for downstream analysis (18). Following acquisition of each MS 2 spectrum, a synchronous precursor selection MS 3 scan was collected on the top 10 most intense ions in the MS 2 spectrum (53). MS 3 precursors were fragmented by high-energy collision-induced dissociation and analyzed using the Orbitrap (normalized collision energy 65; AGC 2 ϫ 10 5 ; maximum injection time 500 ms; resolution was 50,000 at 200 Th).
Data analysis-Mass spectra were processed using a Cometbased (2018.01 rev. 2) in-house software pipeline (54,55). Spectra were converted to mzXML using a modified version of ReAdW.exe. Database searching included all entries from the Human Reference Proteome (2017-05) UniProt database, as well as an in-house curated list of contaminants. This database was concatenated with one composed of all protein sequences in the reverse order. Searches were performed using a 20-ppm precursor ion tolerance for total protein level analysis. The product ion tolerance was set to 0.03 Da. TMT tags on lysine residues and peptide N termini (ϩ229.163 Da) and carbamidomethylation of cysteine residues (ϩ57.021 Da) were set as static modifications, whereas oxidation of methionine residues (ϩ15.995 Da) and GlyGly modification (ϩ114.0429 Da) were set as variable modification. Peptide-spectrum matches (PSMs) were adjusted to a 1% false discovery rate (FDR) (56). PSM filtering was performed using a linear discriminant analysis, as described previously (55), while considering the following parameters: Comet Log Expect, Diff Seq. Delta Log Expect, XCorr, missed cleavages, peptide length, charge state, and precursor mass accuracy. The presence of at least one GlyGly modification was required. For TMT-based reporter ion quantitation, we extracted the summed signal/noise ratio for each TMT channel and found the closest matching centroid to the expected mass of the TMT reporter ion (integration tolerance of 0.003 Da). Ubiquitylation site localization was determined using the AScore algorithm (57). AScore is a probability-based approach for high-throughput protein phosphorylation site localization. Specifically, a threshold of 13 corresponded to 95% confidence in site localization. Ubiquitylated peptides were quantified by summing reporter ion counts across all matching PSMs using in-house software, as described previously (55). PSMs with poor quality, MS3 spectra with isolation specificity Ͻ0.5 or with TMT reporter summed signal/noise ratios that were less than 110, or PSMs that had no MS3 spectra were excluded from quantification (53).
Peptide quantification values were exported for further analysis in Microsoft Excel and Perseus (58), and the statistical test and parameters used are indicated in the corresponding figure legends. Briefly, Welch's t test analysis was performed to compare two data sets, using the s0 parameter (in essence a minimal -fold change cut-off), and correction for multiple comparison was achieved by the permutation-based FDR method, both functions that are built-in in the Perseus software (58). Tables S1 and S2 list all quantified peptides used as well as the associated TMT reporter ratio to control channels used for quantitative analysis.