Mutant thermal proteome profiling for characterization of missense protein variants and their associated phenotypes within the proteome

Temperature-sensitive (TS) missense mutants have been foundational for characterization of essential gene function. However, an unbiased approach for analysis of biochemical and biophysical changes in TS missense mutants within the context of their functional proteomes is lacking. We applied MS-based thermal proteome profiling (TPP) to investigate the proteome-wide effects of missense mutations in an application that we refer to as mutant thermal proteome profiling (mTPP). This study characterized global impacts of temperature sensitivity – inducing missense mutations in two different subunits of the 26S proteasome. The majority of alterations identified by RNA-Seq and global proteomics were similar between the mutants, which could suggest that a similar functional disruption is occurring in both missense variants. Results from mTPP, how-ever, provide unique insights into the mechanisms that contribute to the TS phenotype in each mutant, revealing distinct changes that were not obtained using only steady-state transcriptome and proteome analyses. Computationally, multisite l -dynamics simulations add clear support for mTPP experimental findings. This shows that mTPP is a precise approach to measure changes in missense mutant – containing proteomes without the requirement for large amounts of starting material, specific antibodies against proteins of interest, and/or genetic manipulation of the biological system. Although experiments were performed under permissive conditions, mTPP provided insights into the underlying protein stability changes that cause dramatic cellular phenotypes observed at nonpermissive temperatures. Overall, mTPP provides unique mechanistic insights into missense mutation dysfunction and connection of genotype to phenotype in a rapid, nonbiased fashion.

Temperature-sensitive (TS) missense mutants have been foundational for characterization of essential gene function. However, an unbiased approach for analysis of biochemical and biophysical changes in TS missense mutants within the context of their functional proteomes is lacking. We applied MS-based thermal proteome profiling (TPP) to investigate the proteomewide effects of missense mutations in an application that we refer to as mutant thermal proteome profiling (mTPP). This study characterized global impacts of temperature sensitivityinducing missense mutations in two different subunits of the 26S proteasome. The majority of alterations identified by RNA-Seq and global proteomics were similar between the mutants, which could suggest that a similar functional disruption is occurring in both missense variants. Results from mTPP, however, provide unique insights into the mechanisms that contribute to the TS phenotype in each mutant, revealing distinct changes that were not obtained using only steady-state transcriptome and proteome analyses. Computationally, multisite l-dynamics simulations add clear support for mTPP experimental findings. This work shows that mTPP is a precise approach to measure changes in missense mutant-containing proteomes without the requirement for large amounts of starting material, specific antibodies against proteins of interest, and/or genetic manipulation of the biological system. Although experiments were performed under permissive conditions, mTPP provided insights into the underlying protein stability changes that cause dramatic cellular phenotypes observed at nonpermissive temperatures. Overall, mTPP provides unique mechanistic insights into missense mutation dysfunction and connection of genotype to phenotype in a rapid, nonbiased fashion.
The ability to systematically assess the function of individual genes and gene products is critical for gaining a full picture of how the cell works (1). Much of our advancement in knowledge on biological processes and the underlying mechanisms is owed to the use of genetic methodologies that perturb normal gene function (discussed in Refs. 2 and 3). One such approach to study the function of essential genes is to use conditional sequence variants, such as temperature-sensitive (TS) mutants. TS mutants have been key tools for decades in the study of many essential biological processes including but not limited to characterization of the cell cycle (4)(5)(6)(7)(8)(9), RNA polymerase II transcription (10)(11)(12)(13)(14)(15)(16)(17), and autophagy (18)(19)(20)(21)(22), among others. TS mutants have distinct advantages in that they provide the ability to negatively regulate protein function without having to change the transcriptional context of the gene, add chemicals, or change media (23). Furthermore, by simply altering the temperature at which cells are grown, TS mutant proteins can be studied in a context in which they retain function (permissive temperature), lose function (nonpermissive temperature), or have partial function (semipermissive temperature) (24,25).
Despite the long-term use of TS mutants, much is still left to discover in terms of the properties and mechanisms leading to temperature sensitivity. Many of the forward genetics screening processes that produce TS mutants result in multiple mutations, which can make determination of which mutations are causative and what mechanism is underlying the altered phenotype difficult. It has been predicted that the majority of phenotype-causing missense mutations in humans alter protein structure and lead to changes in protein stability or proteinprotein interactions (PPIs) (26). This suggests that biophysical details behind a phenotype can likely be revealed by analyzing the cell-wide impacts of missense mutations on the proteome. For instance, an unbiased proteome-wide method for characterizing the global impact of temperature-sensitive mutations on PPIs would increase our ability to link genotype to phenotype in existing, uncharacterized mutants. It would also aid in the development of TS mutants for mechanistic studies of essential genes from any model system. More importantly, any approach developed to examine TS mutants could have broad applicability to studies on phenotype and/or disease-causing protein sequence variants in complex biological systems.
Current widely used approaches, such as mRNA sequencing and global proteomics, provide valuable insight on expression levels; however, they provide no insights into potential changes in protein structure and/or PPIs. Large-scale PPI analysis continues to develop using MS-based workflows (27)(28)(29)(30)(31)(32). These methods have greatly increased our knowledge of PPIs but often require large amounts of starting material and involve other manipulations to the cell, such as introducing an epitope tag, or are reliant on antibodies that can differ widely in affinity and specificity (33,34). Additionally, large-scale interactome studies have to this point been focused on the normal cellular sequence variant for each target protein (known as the bait), leaving protein sequence variants out of consideration. Studies that have incorporated the use of human protein sequence variants into their analysis have been limited by the use of expression systems that are unlikely to reflect the native proteome (26). A nontargeted, proteome-wide MS methodology to evaluate how sequence variants in a protein of interest could affect PPI networks would improve our ability to predict the phenotypic outcomes of genomic mutations and would be more scalable than affinity purification-based approaches. Additionally, broad surveys of protein stability could provide unique biophysical insights, which can be directly compared with complementary computational methods, such as molecular dynamics and multisite l-dynamics simulations.
Nontargeted MS-based analyses have included approaches to characterize protein complex characteristics without affinity purification, including chromatography-based co-fractionation and thermal proximity co-aggregation (TPCA), which is a cellular thermal shift assay/TPP-related method (35)(36)(37)(38)(39)(40). These methods have shown a great deal of promise but have not been fully developed for the investigation of disease-causing protein variants. From these studies, we propose that TPP-based analysis of protein thermal stability could have sufficient specificity and sensitivity to detect both large and small changes in PPI dynamics. Prior TPP experiments have revealed both direct and indirect consequences of on-and off-target drug and ligand binding through shifts in protein melt, suggesting that it could have broad applicability to discover biophysical changes in proteins in varying conditions (41)(42)(43)(44)(45). The principle that changes in protein context can be implied from alterations in thermal stability has also been used to study protein-binding metabolites or nucleic acids (43,(46)(47)(48)(49), post-translational modifications (50)(51)(52), or changes in redox status (53). Additionally, TPP has been used to analyze protein complex behavior in steady-state conditions (40) and throughout different stages of the cell cycle (46,54). Whereas effects of genetic perturbations such as knockouts have been described (55,56), this technology has not yet been applied to study the impacts of missense mutations that lead to changes in cellular fitness. To explore the utility of thermal profiling-based methods for analysis of protein sequence variants that result in pronounced cellular fitness defects within the context of the proteome, we have performed mutant TPP analysis (mTPP) in TS mutant yeast using quantitative MS (57)(58)(59). Missense mutations could alter protein structure and lead to changes in PPIs within associated complexes and, in theory, affect the thermal stability of either the individual protein or the entire PPI network as a consequence of structural or PPI interface changes. Therefore, we hypothesized that mTPP could be applied to study the effects of missense mutations on the proteome to link phenotype-causing mutations to PPI networks. This work clearly shows that mutant TPP is a useful application of thermal profiling for studying genetic perturbations in diverse cellular systems, laying a foundation for novel characterization of the biophysical mechanisms within the context of native proteomes.

Results
Characterizing the global protein and RNA abundance effects of missense mutants in the 26S proteasome For our model system in these experiments, we chose TS yeast strains with missense mutations within individual subunits of the 26S proteasome complex, an essential protein complex responsible for the majority of selective proteolysis occurring in the cell (60,61). The yeast proteasome is a large multisubunit protein complex made up of 33 subunits that form two subcomplexes: a 19S regulatory particle and a 20S core complex. Subunits within the 19S regulatory particle recognize and selectively bind ubiquitylated substrates and deubiquitylate and unfold them (62). Substrates are then mechanically translocated into the proteolytic chamber of the 20S core for degradation (63,64). Because the proteasome plays a major role in protein turnover and homeostasis, the large-scale changes that could occur in terms of protein abundance levels and protein post-translational modifications make it an intriguing model for these studies. One TS mutant was chosen from each of the subcomplexes of the proteasome: one representing the 19S regulatory particle (rpn5-ts) and the other representing the 26S proteasome core particle (pup2-ts).
Yeast strains were obtained from a repository of temperature-sensitive mutants the Hieter laboratory developed as a resource for studying essential genes (23,24). Random missense mutations were introduced within an essential gene of interest with the resulting strains screened to identify mutations that caused a temperature-sensitive phenotype (23). The mutated essential genes are expressed from their native promoter and chromosomal context, allowing for analysis of missense mutation impact on the proteome without the confounding variable of altered protein expression levels. Additionally, the use of haploid yeast strains allows for the characterization of this approach in a genetic system that is not complicated by the presence of multiple alleles for a given gene of interest. These characteristics make these strains an ideal representative model for temperature sensitivity-inducing missense mutations. The initial global proteomics studies were performed at both the permissive (30°C) and restrictive temperatures (37°C) and revealed that there were extensive differences in the global proteome even at the permissive temperature. With these data in mind, all subsequent experiments were performed at the permissive temperature. Because the goal of these studies was to test whether mTPP was sufficient to identify changes in protein thermal stability of protein variants and their associated PPIs, it was reasoned that the less detrimental temperature would allow us to assess thermal stability changes without confounding changes in cellular biology. For instance, extensive changes in cellular signaling would be expected to occur as a consequence of proteasome dysfunction at restrictive temperatures in addition to the general changes that would occur as a consequence of heat stress.
Global proteomics and poly(A) RNA-Seq experiments were performed to gain a broad picture of the abundance changes occurring in these strains. Both mutant and WT cells were grown at permissive temperature (30°C) prior to protein or RNA extraction. For the global proteomics, cells were subjected to the workflow presented in Fig. S1, and protein was extracted using 8 M urea and mechanical disruption. RNA was extracted using an established hot acid phenol method (65) prior to library preparation and sequencing on an Illumina system. A total of 3,889 proteins and 5,750 transcripts were quantified (Tables S1 and S2) across all three genotypes. In the pup2-ts mutant, 1,731 proteins changed (Fig. 1A), and 702 mRNA transcripts changed in abundance (Fig. 1B). In the rpn5-ts mutant, 1,997 proteins (Fig. 1C) and 539 mRNA transcripts (Fig. 1D) changed in abundance. Comparisons between protein and mRNA transcript abundance show a low, but positive, correlation (Fig. S2). It is becoming increasingly evident that the interplay across the different levels of gene expression is complex and that changes in protein abundance cannot be explained by changes in mRNA levels alone (66). The major role of the proteasome is to maintain proteostasis, and mutations in subunits that disrupt its degradation activity are likely to result in protein abundance changes that are not explained by mRNA transcript abundance alone. Additionally, the proteasome has been shown to play critical roles, both proteolytically and nonproteolytically, in the regulation of transcription and chromatin organization (67)(68)(69). The numerous roles of the proteasome across multiple stages of gene expression can have various effects on gene expression, further complicating the relationship between protein and mRNA abundance levels in proteasome mutants.
As both rpn5-ts and pup2-ts are within the same protein complex, it was unsurprising that both global proteomics and RNA-Seq presented a large number of similar changes between both mutants. There was a 70% overlap in protein expression changes (Fig. S3) and a 52% overlap in mRNA expression  Thermal profiling of missense mutants changes ( Fig. S4) in the two mutants. Both pup2-ts and rpn5-ts have increased abundance of many of the proteasome subunits ( Fig. S5) as well as the transcription factor Rpn4 at both the mRNA and protein levels (Tables S1 and S2). The up-regulation at both the mRNA and protein level of the 26S proteasome is a known negative feedback mechanism for proteasome function that occurs through the transcription factor Rpn4. Rpn4 is a known substrate of the 26S proteasome and increases in abundance as a result of defective proteasome activity (70). The increase in Rpn4 abundance in both of these strains suggests that it is not being efficiently degraded, strongly indicating that the proteasome is not functioning properly. In an attempt to compensate for decreased proteasome activity in these strains, Rpn4 likely up-regulates the expression of the subunits of the proteasome to increase the number of proteasomes in the cell. Another strong indicator of proteasome malfunction is an increase of the proteasome chaperones Ump1 in both strains (Tables S1 and S2). Ump1 functions in assembly of the core proteasome and is degraded by the 20S proteasome prior to its binding the 19S regulatory particle and formation of the full 26S proteasome (71,72). The large increase in Ump1 protein accumulation, 12.7-fold increase in protein relative to a 2.8fold increase in mRNA, in pup2-ts may suggest that the TS mutations in Pup2 cause a proteasome assembly defect that exacerbates Ump1 accumulation (Tables S1 and S2). In rpn5-ts, the increase in Ump1 protein is comparable with the change in mRNA of UMP1 (2.5-fold increase in protein and 2.2-fold increase in mRNA), suggesting this may occur through a transcriptional response to the lack of proteasome activity. The increase in the chaperone Ump1 supports the likelihood of proteasome functional defects in these strains, consistent with previous work (73,74). However, the large increase in Ump1 protein in pup2-ts suggests a potential defect in proper proteasome assembly that differs from the changes observed in rpn5-ts cells.
Thermal proteome profiling can be applied to measure changes resulting from genetic perturbation MS methods are widely utilized to measure changes in protein abundance, protein turnover rates, post-translational modifications (PTMs), and PPIs (75). Thermal proteome profiling (TPP) is an emerging method that is revolutionizing how MS is able to gain insights into the proteome. TPP has been used to study such events as target engagement, post-translational modifications, and cell cycle variation (40)(41)(42)(43)(44)(45)(46)(47)(48)(49)(50)(51)(52)(53)(54)(55)(56)(57)(58)76); however, it has not yet been applied to biological questions regarding changes that occur as a consequence of single-gene missense mutations known to cause defects in cellular fitness through unknown mechanisms. We have applied TPP ( Fig. 2A) (57,77) to address the global effects of TS-inducing missense mutations on the proteome, which we refer to as mTPP.
Strains were grown at 30°C as performed for the global proteomics and mRNA expression analyses shown in Fig. 1 and lysed in a nondenaturing lysis buffer containing 0.1% Tween. Cell lysate from each genotype was subjected to the workflow described in Fig. 2A for thermal profiling. Three biological replicates were prepared and were subjected to high-pH basic frac-tionation, resulting in 72 LC-MS/MS runs for each genotype. We detected and quantified a total of 4,073 proteins across three biological replicates for each of the three genotypes ( Fig.  2B and Table S3). Data obtained from quantitative MS/MS experiments were used to make melt curves in two different ways (Fig. 2, C and D). Raw abundances from each temperature treatment were normalized to the ion abundance detected in the 45°C sample and were plotted to visualize melt curves of protein complexes in mutant and WT (Fig. 2C). Melt curves of three different protein complexes, the RNA exosome, the 40S ribosome, and the chaperonin-containing T-complex, are shown in Fig. 2C. From these, we can see that proteins within a complex melt in a similar fashion, consistent with what has been previously shown in TPCA studies (40). The TPP R package (78) was used to generate normalized protein melt curves and calculate melt temperatures (T m ) for ;3,400-3,600 proteins, depending on the replicate (numerical data for these in Table S4).
Average T m value distributions were plotted for the proteomes of WT as compared with individual mutants (Fig. S6). Observation of average T m in the range of 50-56°C indicated that the vast majority of proteins denature in this region, which is comparable with what has been previously observed by TPP (42,45,46) and with other methods measuring T m of the proteome (79). Similar to the global melt, average T m for the individual protein complexes in Fig. 2C was ;52-53°C in WT and in both mutants. The largest variance relative to WT in the average T m for the highlighted complexes was 0.32°C (the calculated difference relative to WT for the exosome in rpn5-ts and the ribosome in pup2-ts) and was not statistically significant. Together, these data suggest that the mutations in rpn5-ts and pup2-ts do not have a large overall effect on global melt or the stability of the exosome, ribosome, or chaperonin protein complexes. Additionally, for the majority of individual proteins, such as those shown in Fig. 2D, we saw very similar T m values between mutant and WT and tight fitting to a sigmoidal curve (as indicated by the r 2 value). These data show a successful adaptation of the TPP workflow for mTPP with the cost-effective change of using only six temperature points, which is of particular advantage for screening studies. All biological replicate plots for the proteins shown in Fig. 2D are provided in Fig. S7.  (Table S5). Further parsing was accomplished by limiting our data to melt curves with r 2 values 0.9 and then by proteins that were detected in at least two of the three replicates, providing us with a final analysis of DT m in ;2,050 proteins (Table S6). The 19S regulatory particle mutant, rpn5-ts, contains 13 missense mutations within Rpn5 (E347G, F402L, M171T, S200T, S416P, S54P, E294G, K147E, E214K, S295G, A328V, I353V, and Q342R), which confer temperature sensitivity at 37°C (Fig. S8) (24). Structurally, these mutations are found around the exterior of Rpn5 ( Fig. S9A) (80). The charge-change associated mutations are solventexposed, lessening their direct effects on the stability of rpn5-ts or the 19S regulatory particle, and many other mutations feature large-to-small changes in side-chain sizes, which collectively are expected to minimally perturb the secondary and tertiary folding of rpn5-ts, although they could impact proteinprotein interactions between Rpn5 and non-19S-binding partners. Whole-proteome analysis of thermal stability found significant (defined by a 95% confidence interval) decreases in the thermal stability of 70 proteins and increases in the stability of 40 proteins ( Fig. 3A and Table S7). Analysis of the melt curves of the subunits of the proteasome in rpn5-ts showed no changes in the 20S core relative to WT (Fig. 3B). At ;55-60°C, a flattening of the thermal melt curve was observed in both WT cells and rpn5-ts in raw melt curves. This observation could suggest a phased melting of the 20S core particle, perhaps representing dissociation of the subunits of the 20S complex followed by single subunit melting at slightly higher temperatures (Fig. 3B). Although this observation is interesting, following normalization, the curve flattening had minor impacts on sigmoidal curve fitting using the TPP package, which revealed that most of the 20 S subunits had average T m values close to 58.5-61.5°C with an S.D. of ;0.58 across all three replicates. Despite the numerous missense mutations within RPN5, neither the thermal stability of Rpn5 itself nor of any subunits of the 19S regulatory particle were significantly affected (Fig. 3C). All replicate data for the raw melt curves for each proteasome subunits are provided in Fig. S10. Following normalization, thermal melt profiles for the proteasome subunits showed tight reproducibility across all replicates and genotypes, indicating that no changes occur in proteasome protein subunit stability in rpn5-ts (Fig. 3, D and E; normalized melt curves for all proteasome subunits are provided in Fig. S11). A common mechanism that is used to explain a temperature sensitivity phenotype is that elevated temperature may disrupt the structural integrity of a mutant protein. In fact, this is a general concept that is thought to explain temperature sensitivity in many conditional mutant organisms (81). However, the lack of an observed shift in stability of Rpn5 suggests that thermal sensitivity of this strain is not necessarily due to structural changes in Rpn5 at restrictive temperatures, but likely because of other PPI effects, changes in Rpn5 function, or compensation pathways that adjust for mutant protein function at the permissive temperatures.
The 20S core particle mutant, pup2-ts, contains four missense mutations within Pup2 (C76R, L204Q, T113M, and S260P)   Thermal profiling of missense mutants and is unable to grow properly at 37°C (Fig. S8) (24). In contrast to rpn5-ts, the pup2-ts mutations, with the exception of S260P, occur within the hydrophobic core of the subunit. These mutations include small-to-large side-chain size changes and a chargechange mutation, which suggests that these amino acid substitutions could be more disruptive to the stability of pup2-ts and, perhaps, to the larger 20S core particle (Fig. S9B) (82). mTPP analysis showed that a total of 22 proteins were thermally destabilized in the pup2-ts mutant ( Fig. 4 and Table S7). Remarkably, of those 22 proteins, all 14 of the 20S core subunits were significantly destabilized, which was clearly observed in melt curve analysis (Fig. 5A). None of the regulatory particle subunits had altered thermal stability in pup2-ts (Fig. 5B), suggesting that the mutations in this protein likely affect proteasome activity by disrupting the core proteasome and not the regulatory particle. Normalization of the data via the TPP package resulted in similar findings showing significant destabilization of the 14 20S core particle subunits (Fig. 5, C and D). The specificity of these findings is striking, with nearly 64% of the destabilized proteins occurring within the same protein complex as the mutant protein. The normalized melt curves for each proteasome subunit in pup2-ts are provided in Fig. S11. Gene ontology analysis (83) of the destabilized protein group identified the gene ontology term proteasomal ubiquitin-independent protein catabolic process as the most enriched, with a -fold change enrichment value .100 and a false discovery rate of 4.25e 225 . These findings clearly show that mTPP has a very high degree of selectivity to identify proteins with altered thermal properties and, in this case, their close PPI partners. Interestingly, 103 proteins were significantly stabilized in pup2-ts (Fig. 4). The pup2-ts-stabilized proteins could be either directly impacted by the pup2-ts mutation or due to secondary impacts of the mutations, such as disruption of proteasome function, leading to an accumulation of polyubiquitylation on these proteins or changes in their PPI network.

Multiomics intersection analysis of mTPP with global proteome abundance and transcriptome data provides unique mechanistic insights into mutant dysfunction
Intersection analysis was used to visualize the total overlap between the mTPP, RNA-Seq, and global protein abundance data sets through the generation of Upset plots using the UpSetR package in R (84) and manual data set intersection analysis (Fig. 6). Integration of mTPP with proteomic and transcriptomic data can provide a clearer picture of the mechanisms linking genotype to phenotype and aid in mechanistic hypothesis generation. For example, a possible explanation for an increase in thermal stability in proteins in these mutants is through a change in the actions, abundance, and/or function of cellular chaperone proteins. There was an increase in the   Thermal profiling of missense mutants to increase in protein abundance 1.7-and 1.4-fold in pup2-ts and rpn5-ts, respectively (86). The up-regulation of cellular chaperones could be a compensation mechanism for defects in 26S proteasome function and could help explain why these cells are able to survive at permissive temperature despite altered proteasomal activity. Multiomics intersection analysis of the three -omics data sets revealed subcategories of interest (Fig. 6). There was little overlap between gene products that increased or decreased in parallel, and in some cases, there were some that changed in opposite ways (e.g. protein up/RNA down or protein up/destabilized), providing potential interesting candidates for further mechanistic studies. In the pup2-ts mutant, 14 of the 15 proteins seen to be both thermally destabilized and have increased mRNA. Intersects indicated in Fig. 6A by arrows are the subunits of the proteasome core, showing the strength of using mTPP in multiomics intersection analysis for identifying candidate protein subsets involved in the mechanism behind an observed phenotype. Aside from the core proteasome subunits, the additional member within this subset is the DNA damageinducible protein 1 (Ddi1), a ubiquitin-binding protein that functions as a proteasome shuttle (reviewed in Ref. 87) that has recently been characterized as a ubiquitin-dependent protease (88). Ddi1 is destabilized in rpn5-ts as well (Table S7) and has been shown to interact with the proteasome via affinity purification MS (89). It is possible that the disruption in proteasome activity is impacting critical interactions between Ddi1 and the proteasome or its substrates, leading to Ddi1 destabilization. In rpn5-ts, one particular intersection that stood out is a gene that was thermally destabilized and decreased in abundance but had an increase in mRNA abundance, Tma17 (also known as Adc17), highlighted by an arrow in Fig. 6B. Tma17 is an important proteasome assembly chaperone that has been shown to help enhance lid assembly in response to stress in a regulatory mechanism independent of Rpn4 (90). An increase in TMA17 mRNA suggests that the cell may be trying to compensate for the destabilization and decreased abundance of this protein by up-regulating its mRNA transcript.
Analyzing the overlap of changes occurring in pup2-ts and rpn5-ts has the potential for discovery of general proteasome interactors as well as highlighting subunit-specific interactors. For example, the intersect in proteins that are increasing in abundance but whose mRNA levels are not increased in both mutants could identify candidate proteasome substrates. An Upset plot presenting the overlaps of like changes occurring in pup2-ts and rpn5-ts is shown in Fig. 7A. Almost 72% (1,117 of 1,554) of proteins that were found to increase in abundance using global proteomics were found to increase in both pup2-ts and rpn5-ts (Fig. 7B, abundance increased). A total of 363 mRNA transcripts were up-regulated in both pup2-ts and rpn5-ts; however, this was only ;50% of the total increases observed in mRNA levels in both mutants (Fig. 7). Strikingly, Overlap between datasets  Figure 7. Comparing changes between rpn5-ts and pup2-ts. Shown is an upset plot of the overlap in changes between rpn5-ts (orange) and pup2-ts (teal). Gene products that are changing in the same way in both genotypes are indicated in gray. B, a stacked bar plot shows the percentage of each data set that is changing in the same direction in both mutants (gray) or changes unique to rpn5-ts (orange) or pup2-ts (teal).
there was a low percentage of proteins that had changes in thermal stability in the same direction in both mutants, six stabilized and three destabilized (Fig. 7), underscoring the remarkable selectivity and reproducibility of the mTPP approach. Compared with the commonly applied approaches of transcriptome and proteome abundance analysis, measuring the thermal stability of proteins using mTPP provided the most unique insights into the underlying impact of each proteasome subunit mutant on the biophysical state of the proteasome and the mechanism of missense mutation dysfunction.
Molecular modeling quantifies the contributions of individual missense mutations to destabilization of pup2-ts and the 20S proteasome Following the observation that pup2-ts cells show destabilization of all of the proteins within the 20S core particle, molecular dynamics (MD) and multisite l-dynamic (MSlD) simulations were performed to further investigate the structural and thermodynamic mechanisms of pup2-ts destabilization. MD simulations enable the dynamic motion of atoms in molecules to be modeled as a function of time, yielding atomistic insights into biophysical observables (91). From all-atom MD simulations of Pup2 and pup2-ts, isolated in solution, differences in root mean square fluctuation (DRMSF) of individual residues over the course of the simulation revealed that the C76R, T113M, and L204Q mutations induce several dynamic changes in the structure of Pup2 (Table S8 and Fig. S12). The S260P mutation was not included in the molecular model; this residue is positioned on the flexible C-terminal tail of pup2-ts, and no electron density for this region was captured in the crystal structure (82). Notably, large increases of dynamic fluctuations of ;0.2-1.0 Å are observed in the two a-helices composed of residues 84-139, which contain the T113M mutation and sit directly adjacent to the C76R mutation. It is clear that these two small-to-large side-chain perturbations destabilize this region of Pup2 dramatically. In contrast, the L204Q mutation induces very few fluctuation changes. Extrapolating these findings onto the 1RYP crystal structure of the yeast 20S proteasome illustrates that Pup2 residues with the largest DRMSF line the internal cavity of the 20S core and form important PPIs with neighboring 20S a subunits (Fig. S12, C and D) (82). Inducing large dynamic fluctuations in this region of the 20S core particle is thus expected to weaken important interactions to neighboring subunits and disrupt the rigid structure of the core proteasome.
To further investigate the individual contributions of each mutation to Pup2 destabilization, MSlD calculations were performed. MSlD is a rigorous free-energy simulation method able to calculate changes in protein-binding affinity in response to point mutations (92,93). Compared with other free-energy simulation methods, MSlD is highly efficient, requiring an order of magnitude less computation than standard approaches without loss of precision, and it facilitates sampling of different perturbations at multiple sites in a single simulation (92,(94)(95)(96)(97). This offers the advantage of exploring nonadditive effects between multiple point mutations within a protein simultaneously (95). This study also represents the first application of MSlD to investigate protein-protein interactions. The mutations in pup2-ts, specifically C76R, T113M, and L204Q, were explored in a combinatorial fashion in both insolated and complexed states of Pup2. Relative free energies of binding (DDG bind ) between WT Pup2 and eight different combinations of the mutations, including the combination of all three found in pup2-ts, were computed using a standard thermodynamic cycle at 30 and 50°C (Table 1 and Fig. S13) (98). The freeenergy data in Table 1 fully support the experimental mTPP findings that the complete 20S core particle is destabilized by the incorporation of pup2-ts, with a computed DDG bind of ;3.8 kcal/mol. The most destabilizing perturbation is C76R, which introduces a large, charged residue into the hydrophobic core of Pup2 and causes several structural deformations in Pup2.
Analysis of the MSlD trajectories shows that the Arg-76 can either protrude between residues 84-101 and 109-120 (referred to as "inward") or between residues 66-74 and 84-101 (referred to as "outward"), as shown in Fig. 8. In both conformations, arginine disrupts the hydrophobic network of nearby valine, leucine, and isoleucine residues, opens up the core of the protein, and facilitates the penetration of nearby water molecules into the core of Pup2 (Fig. 8B). Interestingly, the T113M and L204Q mutations, alone, do not have a significant effect on the free energy of binding of Pup2 to the 20S core particle, with DDG bind results of 20.20 kcal/mol or less. However, these residues play an important complementary role toward stabilizing the C76R mutation. Clear stabilizing, nonadditive, DDG bind effects of 20.7 to 20.8 kcal/mol per residue are observed to reduce the singular C76R mutation DDG bind from 5.15 to 3.79 kcal/mol in the full pup2-ts mutant. Structural observations help explain these effects. When Arg-76 protrudes inward, Met-113 can help fill the additional space introduced by Arg-76 as well as hydrogen-bond with water molecules that enter the core (Fig. 8A); when Arg-76 protrudes outward, Gln-205 is able to hydrogen-bond to the backbone of Gln-218 and thus stabilize large fluctuations in the bottom half of Pup2 (Fig. 8D). As temperature is increased to 50°C, however, this nonadditive complementarity is reduced, and Pup2 is more fully destabilized by the C76R mutation. These effects likely contribute to the mTPP observed decrease in T m for pup2-ts compared with WT Pup2. Collectively, these findings support that the C76R, T113M, and L204Q temperature-sensitive mutations in Pup2 destabilize the structure of Pup2 and weaken its ability to bind to the 20S core particle. They also reveal that the majority of Thermal profiling of missense mutants these effects stem from the C76R mutation, with nonadditive, but complementary and stabilizing, effects from T113M and L204Q.

Discussion
Temperature sensitivity is generally thought to lead to loss of function of a mutant protein at restrictive temperatures due to decreased protein stability at higher temperatures. Our results suggest that this is not always the case and that the biophysical properties of TS proteins are likely mutation-specific, with unique mechanisms resulting in temperature sensitivity in pup2-ts and rpn5-ts. In the pup2-ts mutant, we calculated DT m values for over 2,000 proteins and remarkably only saw destabilization of 22 proteins, 14 of which are the 20S core proteasome itself (Fig. 4). These findings suggest that mutations in pup2-ts thermally destabilize Pup2 and that elevated temperatures can lead to a disruption of the stability of the entire proteasome core. This was well-supported by MD and MSlD simulations, which quantified the individual contributions of three pup2-ts mutations and demonstrated loss of binding affinity of pup2-ts to the 20S proteasome. These data also show the highly com-plementary nature of mTPP, MD, and MSlD, with mTPP providing experimental evidence in support of MD and MSlD findings. In summary, we conclude that the proteasome dysfunction in this mutant is most likely due to changes in folding and incorporation of Pup2 into the 20S core particle, leading to a defect in or a decreased efficiency of overall proteasome core assembly. As assembly of the regulatory particle occurs independently of the core (99), it is reasonable that mutations in the hydrophobic interior of Pup2 do not destabilize the 19S regulatory particle subunits alongside the 20S core. Interestingly, despite presenting with similar protein and RNA abundance alterations, the mutations in rpn5-ts do not appear to affect the thermal stability of Rpn5 or the 26S proteasome complex.
Mutations in rpn5-ts instead lead to temperature sensitivity via a different mechanism, likely disruption of other proteinprotein interactions essential to proteasome function. The fact that many rpn5-ts mutations are close to the surface of the protein supports this notion. Other studies of the rpn5-ts mutant found that these mutations in Rpn5 cause nuclear mislocalization of the regulatory particle as well as chromosomal instability (73,74). Of the 110 proteins mTPP identified to have stability changes in rpn5-ts, 50 are found within the nucleus, so it could be that the observed phenotype is due to the loss of essential nuclear interactions of the proteasome. Another possibility could be that the mutations in Rpn5 do not allow for proper maturation and function of the 26S proteasome (80, 100-108). The regulatory particle chaperone Tma17 was significantly destabilized and decreased in protein abundance, despite having an increase in mRNA abundance. Therefore, the mutations in rpn5-ts could prevent the proper interaction of Tma17 with the proteasome lid, leading to its destabilization. Due to the important role for Tma17 in increasing proteasome assembly when proteolysis capacity has been exceeded (109), the destabilization of Tma17 could be further crippling the proteasome defects in this mutant, and the increase in TMA17 transcript suggests the cell is trying to compensate for the destabilization and decreased abundance of this protein. Finally, it is possible that direct and/or indirect effects on protein PTMs, such as phosphorylation, could contribute to changes in stability, but phospho-TPP experiments suggest that these effects would likely be subtle (52).
Global protein and RNA abundance measurements from rpn5-ts and pup2-ts revealed a large number of possible candidates that could be causing phenotypic changes observed in the mutant strains. These experiments were performed at permissive temperature, indicating that even at temperatures that allow for normal cell growth, the mutations in rpn5-ts and pup2-ts are perturbing proteasome function. It remains a possibility that a number of the proteins that increase in abundance in both mutants as a consequence of proteasome perturbation are direct substrates of the proteasome. Additional studies measuring changes in PTMs could shine further light on these effects as well. For example, ubiquitylation in these mutant strains may help determine whether some of the changes in protein expression and/or stability in either of the proteasome mutants are due to accumulation of ubiquitylated proteins. Many of the alterations in global protein and RNA abundance were the same across the two mutants, which could lead to the conclusion that they have similar functional disruption. However, the additional data gained through mTPP on protein stability and potential PPI changes suggest that the rpn5-ts and pup2-ts mutants are altered through different mechanisms. These missense mutants contain relatively few proteins that are changed in thermal stability, most of which are unique to the mutant subunit studied, suggesting mTPP to be sensitive enough to measure mutation dependent shifts in thermal stability and protein-specific associations, even between two proteins within the same protein complex. As the mTPP protocol measures protein thermal stability at temperatures outside of physiological range, it allows for the exaggeration of biophysical effects that may be too small to cause catastrophic protein unfolding in live cells. Use of this application allows for the determination of mechanistic differences in the phenotypes of two mutant strains in the context of their native proteome that would be difficult, or perhaps impossible, to discover through other existing experimental methods in a high-throughput manner.
Our study produced quantitative melt curves for over 3,000 proteins across each of the three genotypes evaluated. Although these mutations affect overall proteostasis, only ;100-130 proteins had a change in thermal stability, and there was not a global effect on protein thermal stability, with the average global melt temperature remaining similar to WT. There were few changes in thermal stability relative to abundance changes, which suggests that proteasome targets (many of them likely polyubiquitylated) may not have major changes in overall stability that could assist in the unfolding necessary for degradation by the 26S proteasome. mTPP identified possible PPIs that are unique to Pup2 and Rpn5 individually, as seen through proteins that change in thermal stability in only one of the two mutants, providing an avenue through which one can study the unique biophysical mechanisms that are occurring due to mutations in different proteins that lead to similar phenotypes. Finally, these data show that mTPP can be used to reproducibly measure global protein abundance from the lowtemperature treatment conditions in yeast, removing the need for separate global proteomics experiments, consistent with what has previously been shown in TPP experiments in HeLa cells (46). Despite being isolated in two different buffers, correlation values were greater than 0.8 between experiments (Fig.  S14A), and there was an 82% overlap in protein identification (Fig. S14B).mTPP is capable of measuring changes in the thermal stability of proteomes as a result of temperature sensitivity-inducing missense mutations without the need for any additional genetic manipulation or antibody-based interrogation. Considering the inherent challenges with antibody-based methods because of specificity concerns (34), methods that avoid dependence on antibodies are likely to increase the rigor and reproducibility of mutant protein interrogation. Additionally, the lack of a requirement of protein subunit epitope tagging, a genetic manipulation, allows for rapid interrogation of proteome-transcriptome-phenotype characterization. Epitope tagging is a key molecular biology tool for affinity purification followed by MS (AP-MS). Although AP-MS is a powerful approach for protein-protein interaction identification, it requires a large amount of starting material, which is challenging to obtain from many biological systems, including patient samples (110)(111)(112), and can require a substantial amount of time for optimization of the purification. Comparison of mTPP thermal stability changes with MD and MSlD was extremely useful to explain the experimental changes in protein thermal stability using existing three-dimensional structures for the 26S proteasome. However, mTPP can also provide biophysical insight in the absence of existing structural data for proteins of interest.
Not only are mutant proteins important molecular biology tools, but more than 100,000 genomic mutations have been linked to human diseases, of which ;60% are predicted to affect protein stability and/or PPIs (26). Whereas genomic and transcriptomic sequencing data are abundant and readily accessible, changes in mRNA have been shown to explain only ;40% of protein level changes (113). A full understanding of disease relies not only on genomic data, but also on defining the status of cellular proteins in the context of their abundance, post-translational modification status, PPIs, stability, and turnover rate. Previous studies that have obtained high-throughput data on disease-causing mutations, for instance, have relied on approaches, such as yeast two-hybrid and chaperone Thermal profiling of missense mutants engagement studies, that are carried out in vitro (26). Although these approaches can provide unique molecular insights, they also rely on gene cloning and overexpression, which will artificially alter the balance of allele-specific (in diploid cells) and/or PPI partner expression, leading to confounding changes throughout the proteome. In addition to the capability of mTPP for the characterization of TS mutants for studies of normal biological mechanism, use of this method with disease-associated mutant proteins will help lay the foundation for a bioinformatic resource increasing understanding of disease mechanisms and addressing the pathogenicity of missense mutations.

Yeast strains and growth conditions for proteomics experiments
Temperature-sensitive yeast strains were obtained from the Hieter laboratory (23,24). Two strains with missense mutations within the 26S proteasome, rpn5-ts and pup2-ts, were chosen for these experiments. Both mutant and WT cells were inoculated at an A 600 = 0.3 and grown to an A 600 = 0.8 in yeast extract, peptone, dextrose (YPD) medium at permissive temperature (30°C). YPD was removed by filtration through a nitrocellulose membrane. Cells were flash-frozen with liquid nitrogen and stored at 280°C to be used in subsequent sample preparation steps.

mTPP sample preparation
Cells for each genotype and biologically distinct replicate (i.e. three genotypes and three biological replicates, for a total of nine), were removed from nitrocellulose membrane using water, pelleted, and reconstituted in lysis buffer containing 40 mM Hepes-KOH (pH 7.5), 10% glycerol, 250 mM NaCl, 0.1% Tween 20, and fresh yeast protease inhibitors. Resulting suspensions were lysed using a mini bead beater and a vortex genie. After removing the beads, lysates were subjected to centrifugation (14,000 3 g for 20 min) to pellet out cellular debris and to obtain a clear lysate. Protein concentration values for each resulting lysate were determined using the Bio-Rad protein assay (Bio-Rad), employing vendor-provided protocols. All lysates were then diluted to a protein concentration of 5 mg/ml for subsequent heat treatment. Six 50-ml aliquots for each genotype and replicate were then distributed across six PCR tubes and equilibrated at six temperature points-35.0, 45.3, 50.1, 55.2, 60.7, and 74.9°C for 3 min in a thermocycler (Mastercycler Pro, Eppendorf, Hamburg, Germany) system as described elsewhere (57). Following heat treatment, lysates were centrifuged for 20 min at 4°C to pellet out insoluble protein and then to decant the soluble fraction. A 20-ml aliquot from each heattreated sample was subject to TCA-based protein precipitation overnight. Following centrifugation and acetone-based wash steps, resulting protein pellets were reconstituted in 30 ml of 8 M urea in Tris-HCl (pH 8.0). Samples were subjected to reduction of Cys-Cys bonds with 5 mM tris(2-carboxyethyl)phosphine hydrochloride and alkylation with 10 mM chloroacetamide to protect the reduced Cys residues. Samples were diluted to 2 M urea and were digested in-solution overnight using tryp-sin (Promega Corp., Madison, WI, USA) as described previously to derive peptides (111,112,114). Samples were desalted using a Sep-Pak ® Vac 1 cc C18 Cartridge (Waters Corp., Milford, MA, USA) according to the vendor-provided vacuum manifold protocol that involves 1) priming the column, 2) loading the peptides onto the column, 3) washing the column to remove any buffers from the immobilized peptides, and 4) elution of immobilized peptides. This was performed to remove buffer additives and excess reagents as many lysis buffers are known to interfere with the TMT labeling reaction. Resulting peptide elutions were dried using a speed vacuum and subjected to TMT labeling using vendor-provided protocols. A TMTsixplex TM kit (Thermo Fisher Scientific, Waltham, MA, USA) was employed for the purpose, and, more specifically, channels TMT126, TMT127, TMT128, TMT129, TMT130, and TMT131 were used to label peptide solutions derived from the 35.0, 45.3, 50.1, 55.2, 60.7, and 74.9°C temperature treatments, respectively. Labeling reactions were quenched, and the six labeled samples were mixed and were dried using a speed vacuum system. Dried samples were fractionated using reversed-phase fractionation columns (eight fractions) employing vendor-provided protocols (Pierce). The resulting eight fractions were dried using a speed vacuum system and resuspended in 0.1% formic acid (30 ml) prior to nano-LC-MS/MS analysis as described below.

mTPP data acquisition: Nano-LC-MS/MS
Nano-LC-MS/MS analyses were performed on a Q-Exactive Plus TM mass spectrometer (Thermo Fisher Scientific) coupled to an EASY-nLC TM HPLC system (Thermo Fisher Scientific). A 10-ml equivalent volume of the resuspended fractions from above was loaded using 300 bars as applied maximum pressure onto an in-house prepared reversed-phase column. Each reversed-phase column was prepared by pulling a 100 mM fused-silica column to carry a ;5-mm tip for the nanospray using a P-2000 laser puller and then packing the capillary with C18 reverse-phase resin (particle size, 3-mm diameter; Dr. Maisch HPLC GmbH, Ammerbuch, Germany). The peptides were eluted using a varying mobile phase (MP) gradient from 95% phase A (formic acid/H 2 O, 0.1:99.9 (v/v)) to 24% phase B (formic acid/acetonitrile, 0.4:99.6 (v/v)) for 150 min, from 24% phase B to 35% phase B for 25 min, and then keeping the same MP composition for an additional 5 min at 400 nl/min to ensure elution of all peptides. Nano-LC mobile phase was introduced into the mass spectrometer using a Nanospray Flex Source (Proxeon Biosystems A/S, Odense, Denmark). The heated capillary temperature was kept at 275°C, and ion spray voltage was kept at 2.5 kV. During peptide elution, the mass spectrometer method was operated in positive ion mode for 180 min, programmed to select the most intense ions from the full MS scan using a top 20 method. Additional parameters were as follows: microscans, 1; resolution, 70,000; AGC target, 3E6; maximum IT, 50 ms; number of scan ranges, 1; scan range, 400-1,600 m/z; and spectrum data type, "profile." Then datadependent MS/MS scans were performed with the following parameters: microscans, 1; resolution, 35,000; AGC target, 1E5; maximum IT, 64 ms; loop count, 20; MSX count, 1; isolation window, 0.7 m/z; fixed first mass, 100 m/z; NCE, 38.0; spectrum data type, "centroid." The respective data-dependent settings were set with the following parameters: minimum AGC target, 1.00e3; intensity threshold, 1.6e4; apex trigger, "2"; charge exclusion, "1, 7, 8, .8"; multiple charge states, "all"; peptide match, "preferred"; exclude isotopes, "on"; dynamic exclusion, 30.0 s; if idle "pick others." The data were recorded using Thermo Xcalibur (version 4.1.31.9) software (Thermo Fisher Scientific).
Sample preparation for global quantitative proteomics comparison of mutant proteomes relative to WT Three biological replicates of WT, pup2-ts, and rpn5-ts were prepared as explained above. Unlike the previous preparation, cells were lysed in 8 M urea using a mini bead beater and a vortex genie. A 50-mg equivalent of protein from each lysate was then subjected to reduction, alkylation, proteolytic digestion, and desalting as described above. Peptide samples were labeled using a 10plex kit (Thermo Fisher Scientific, Waltham, MA, USA). Specific labeling reagents corresponding to TMT127N, TMT127C, TMT128N, TMT128C, TMT129N, TMT129C, TMT130N, TMT130C, and TMT131 were respectively used to label the peptide solutions derived from three replicates from each genotype: WT, PUP2-TS, and RPN5-TS. Following reaction quenching, mixing, and the subsequent drying step, the resulting peptide mixture was fractionated as described, resulting in eight fractions. Each fraction was dried completely using a speed vacuum system and resuspended in 0.1% formic acid (30 ml) prior to the nano-LC-MS/MS analysis explained below.

Global proteomics data acquisition: Nano-LC-MS/MS
Nano-LC-MS/MS analyses were performed on an Orbitrap Fusion TM Lumos TM mass spectrometer (Thermo Fisher Scientific) coupled to an EASY-nLC TM HPLC system (Thermo Fisher Scientific). An 18-ml equivalent volume of the resuspended fractions was loaded onto an in-house prepared reversed phase column using 600 bars as the applied maximum pressure. Each reversed-phase column was prepared by pulling a 100-mm fused-silica column to carry a ;5-mm tip for the nanospray using a P-2000 laser puller and then packing the capillary with a C18 reversed-phase resin (particle size, 3-mm diameter; Dr. Maisch HPLC GmbH). The peptides were eluted using a 180-min gradient increasing from 95% buffer A (0.1% formic acid in water) and 5% buffer B (0.1% formic acid in acetonitrile) to 25% buffer B at a flow rate of 400 nl/min. The peptides were eluted using a 180-min gradient increasing from 95% buffer A (0.1% formic acid in water) and 5% buffer B (0.1% formic acid in acetonitrile) to 25% buffer B at a flow rate of 400 nl/ min. Nano-LC mobile phase was introduced into the mass spectrometer using a Nanospray Flex Source (Proxeon Biosystems A/S). During peptide elution, the heated capillary temperature was kept at 275°C, and ion spray voltage was kept at 2.6 kV. The mass spectrometer method was operated in positive ion mode for 180 min, having a cycle time of 3 s for MS/MS acquisition. MS data were acquired using a data-dependent acquisition method that was programmed to have two data-dependent scan events following the first survey MS scan. During MS1, using a wide quadrupole isolation, survey scans were obtained with an Orbitrap resolution of 120,000 with vendordefined parameters-m/z scan range, 375-1,500; maximum injection time, 50; AGC target, 4E5; microscans, 1; RF lens (%), 30; data type, "profile"; polarity, positive with no source fragmentation and to include charge states 2-7 for fragmentation. Dynamic exclusion for fragmentation was kept at 60 s. During MS2, the following vendor-defined parameters were assigned to isolate and fragment the selected precursor ions: isolation mode = quadrupole; isolation offset = off; isolation window = 0.7; multi-notch isolation = false; scan range mode = auto normal; FirstMass = 120; activation type = CID; collision energy (%) = 35; activation time = 10 ms; activation Q = 0.25; multistage activation = false; detector type = IonTrap; IonTrap scan rate = Turbo; maximum injection time = 50 ms; AGC target = 1E4; microscans = 1; data type = centroid. During MS3, daughter ions selected from neutral losses (e.g. H 2 O or NH 3 ) of precursor ion CID during MS2 were subjected to further fragmentation using higher-energy C-trap dissociation to obtain TMT reporter ions and peptide-specific fragment ions using the following vendor-defined parameters: isolation mode = quadrupole; isolation window = 2; multi-notch isolation = true; MS2 isolation window (m/z) = 2; number of notches = 3; collision energy (%) = 65; orbitrap resolution = 50,000; scan range (m/z) = 100-500; maximum injection time = 105 ms; AGC target = 1E5; data type = centroid. The data were recorded using Thermo Scientific Xcalibur (version 4.1.31.9) software (Thermo Fisher Scientific).

Protein identification and quantification
Resulting RAW files were analyzed using Proteome Discoverer TM 2.2 (Thermo Fisher Scientific). The SEQUEST HT search engine was used to search against a yeast protein database from the UniProt sequence database (December 2015) containing 6,279 yeast protein and common contaminant sequences (FASTA file used available on ProteomeXchange under accession number PXD017223). Specific search parameters used were: trypsin as the proteolytic enzyme, peptides with a maximum of two missed cleavages, precursor mass tolerance of 10 ppm, and a fragment mass tolerance of 0.6 Da. Static modifications used for the search were 1) carbamidomethylation on cysteine residues and 2) TMTsixplex label on lysine residues and the N termini of peptides. Dynamic modifications used for the search were oxidation of methionine and acetylation of N termini. The percolator false discovery rate was set to a strict setting of 0.01 and a relaxed setting of 0.05. Values from both unique and razor peptides were used for quantification. For the mTPP experiments, no normalization setting was used for protein quantification; for the global protein abundance measurements, peptides were normalized by "total peptide amount." Data are available via ProteomeXchange with identifiers PXD017222 and PXD017223.

TPP analysis and waterfall plots
The TPP package (version 3.12.0) (78) in R was used to generate normalized melt curves and to determine protein melt temperatures as described previously (77). Changes in T m for Thermal profiling of missense mutants each protein were calculated utilizing the function WT (T m )mutant (T m ), where WT (T m ) represents the median melt temperature of a protein determined in WT samples, and mutant (T m ) represents the median melt temperature of a protein determined in mutant samples. Resulting data processing was performed in R Studio (R Studio for Mac, version 1.1.456). Data were parsed down to proteins that were detected in at least two of the three replicates and had melt curves with r 2 values .0.9. Proteins were ranked according to median change in T m and ordered from the largest change (proteins that were destabilized in the mutant) to smallest change (proteins that were stabilized in the mutant). From the distribution of T m changes in a single comparison, T m changes that are out of the range, mean (T m change) 6 2s (where s is the S.D.), were considered statistically significant T m changes and identified as proteins destabilized or stabilized due to the mutation. Waterfall plots were created using ggplot2 (83).

RNA extraction
Both mutant and WT cells were inoculated at an A 600 = 0.5 in YPD medium and permitted to recover for 90 min at permissive temperature (30°C). Cells were pelleted, washed with water, and resuspended in 10 ml of acetate-EDTA buffer (50 mM sodium acetate, pH 5.2, 10 mM EDTA). RNA was isolated using a hot acid phenol method as described previously (65). In brief, 800 ml of 20% SDS and 10 ml of cold acid phenol were added to each sample in a Nalgene phenol-resistant tube. Samples were incubated for 10 min with rotation in a hybridization oven at 65°C and cooled on ice for 5 min. The samples were then centrifuged at 10,000 rpm for 15 min, and the top phase was transferred to a 50-ml 5 PRIME Phase Lock Gel Tube (reference no. 2302870). 13 ml of chloroform was added, and samples were centrifuged at 3,000 rpm for 10 min; the remaining top phase was poured into a new phenol-resistant tube. Onetenth volume sodium acetate at pH 5.2 and an equal volume of isopropyl alcohol were added, and precipitated RNA was collected by centrifugation at 12,000 rpm for 45 min. The RNA pellet was washed with 70% EtOH and dried in a fume hood. RNA was resuspended with molecular biology grade water. Following purification, RNA was treated using the TURBO DNA-free kit (Thermo Fisher Scientific) to remove DNA contamination.

Library preparation and sequencing
The concentration and quality of total RNA samples were first assessed using an Agilent 2100 Bioanalyzer. A RIN (RNA integrity number) of five or higher was required to pass the quality control. Then 100 ng of RNA per sample were used to prepare dual-indexed strand-specific cDNA libraries using the KAPA mRNA HyperPrep kit (Roche Applied Science). The resulting libraries were assessed for quantity and size distribution using Qubit and an Agilent 2100 Bioanalyzer. Pooled libraries were utilized for clustering amplification on cBot using a HiSeq 3000/4000 PE Cluster Kit and sequenced with a 2 3 75bp paired-end configuration on HiSeq4000 (Illumina) using a HiSeq 3000/4000 PE SBS Kit. A Phred quality score (Q score) was used to measure the quality of sequencing. More than 94% of the sequencing reads reached Q30 (99.9% base call accuracy).

Sequence alignment and gene counts
The sequencing data were first assessed using FastQC (version 0.11.5, Babraham Bioinformatics, Cambridge, UK) for quality control. All sequenced libraries were mapped to the yeast genome (UCSC sacCer3) using STAR RNA-seq aligner (version 2.5) (115) with the following parameter: "2outSAMmapqUnique 60." The reads distribution across the genome was assessed using bamutils (from ngsutils version 0.5.9) (116). Uniquely mapped sequencing reads were assigned to S288C R64-2-1 20150113 annotated genes using featureCounts (subread version 1.5.1) (117) with the following parameters: "-s 2 -p -Q 10." Genes with read count per million ,0.5 in 4 of the samples were removed. The data were normalized using the TMM (trimmed mean of M values) method. Differential expression analysis was performed using edgeR (version 3.12.1) (118,119). False discovery rate (FDR) was computed from p values using the Benjamini-Hochberg procedure. All files from the RNA-Seq experiments performed for this study have been uploaded to the Gene Expression Omnibus (GEO) under the accession number GSE143236.

Global abundance plots
Volcano plots and correlation plots were created in R Studio using ggplot2. Abundance values, -fold change ratios, and p values for the proteomics data were calculated by Proteome Discover TM 2.2 (Thermo Fisher Scientific). FPKM values, -fold change ratios and FDR values were calculated as described above. p value and FDR cutoff for significance was set at ,0.05. Venn diagram was created using Venny version 2.1 (120). Abundance dot plots were created using GraphPad Prism version 6.00 for Windows. Upset plots were created using the UpSetR package in R Studio (84).

Molecular modeling computational details
System preparation-Initial starting coordinates for Pup2 and the 20S proteasome complex were obtained from the Protein Data Bank file 1RYP (82). The protein complex was set up similarly to Vilseck et al. (96), with only protein residues within a 27-Å radius around Pup2 included in the complex simulations. All protein backbone atoms beyond 17 Å from Pup2 were harmonically restrained to their crystallographic positions with a force constant of 10 kcal/mol·Å 2 . All side chains and protein residues within a 17-Å radius from Pup2, and including Pup2, remained fully flexible. Starting coordinates for isolated, unbound Pup2 were generated by extracting Pup2 (chain E) from the 1RYP complex. Residue flips for histidine, glutamine, and asparagine were assessed using the MolProbity webserver (121,122), and protonation states of titratable residues were determined with the assistance of PROPKA (123). After examining the structures and considering the MolProbity and PROPKA results, the side chains of Glu-56 of chain L, His-16 of chain D, His-73 of chain E, and His-67 of chain M were protonated, and Lys-62 of chain F was deprotonated. All protein residues were modeled using the CHARMM36 all-atom force field for proteins (124,125). Standard CHARMM36 patches of GLYP, NTER, and CTER were applied to native N-and C-terminal residues, and peptide fragments resulting from the truncation of the full 20S complex were capped with ACE and CT3 patches, respectively. The truncated complex and unbound Pup2 were then solvated with the TIP3P water model (126) using the CHARMM-GUI webserver (127) with a minimum of 10 Å of solvent from each face of the structure. This produced cubic water boxes with edge lengths of 125.0 Å (truncated complex) and 85 Å (Pup2). Na 1 and Cl 2 ions were added to neutralize the net charge of the system and achieve an ionic strength of 100 mM NaCl. Prior to running MD or MSlD, the CHARMM molecular simulation package (version 44b1) was used to minimize each solvated system obtained from the CHARMM-GUI (128,129). First, all protein atoms were harmonically restrained, whereas all ions and water molecules were relaxed with 200 steps of steepest-descent minimization, after which all restraints were removed and the entire system was subject to an additional 250-500 steps of steepest-descent minimization.
Molecular dynamics-MD simulations were run with OpenMM on graphic-processing units (130). For Pup2 and pup2-ts systems, four independent MD duplicates were run for 100 ns each. MD simulations were run in the isobaricisothermal (NPT) ensemble at 30°C and 1 atmosphere using Langevin dynamics, with a friction coefficient of 10 ps 21 and a Monte Carlo barostat (131). Periodic boundary conditions were employed with nonbonded cutoffs of 12 Å to truncate all long-range interactions. An integration time step of 2 fs was used, and trajectory frames were saved every 5,000 steps. Trajectories were then superimposed onto the minimized starting structure using VMD (132). Root mean square fluctuations were calculated as root mean square average distances between an atom and its mean position in a set of structures, using Equation 1.
Over the course of N s frames in each of the MD trajectories, x i coordinates of the backbone a carbon for each k residue in Pup2 were analyzed to calculate an RMSF for each residue (133). RMSF values were then averaged over the four duplicate runs for Pup2 and pup2-ts, respectively, and plotted in Fig.  S12A. Root mean square deviations of all backbone a carbon atoms were then calculated with respect to the minimized starting structure to show global protein fluctuations as a function of time sampled.

RMSD
In this analysis, an RMSD is calculated for every frame by looping over all N a atoms with x 0 reference and x i target coordinates (133). RMSDs between duplicate MD runs were again averaged and plotted in Fig. S12B.
Multisite l-dynamics-MSlD calculations were performed in CHARMM using the domain decomposition (DOMDEC) module to enable the use of graphic-processing units (92,93,128,129,134). Simulations were run with a Nose-Hoover thermostat and Langevin pressure piston, with a friction coefficient of 20 ps 21 , to maintain sampling within the NPT ensemble at 30°C and 1 atmosphere (135)(136)(137). Trajectory frames were saved every 10,000 steps, and l states were saved every 10 steps for computing relative free-energy differences. The SHAKE algorithm was employed to constrain all hydrogen-heavy atom bond lengths, and periodic boundary conditions were employed (138). The particle mesh Ewald method was used to compute all long-range electrostatic interactions (139)(140)(141), and all nonbonded Lennard-Jones long-range interactions were truncated at 10 Å, with force-switching to zero between 9 and 10 Å (142). With MSlD, the C76R, T113M, and L204Q pup2-ts mutations were simulated as a multiple-topology model, where all atomic components of native and mutant side chains were explicitly included, using the BLOCK facility within CHARMM. Following standard reported procedures for using MSlD (92)(93)(94)(95)(96)(97), holonomic restraints were used to constrain l between 0 and 1, thus maintaining physically relevant end points during the simulation (143). For calculating free-energy differences, a threshold of l 0.99 was used as an approximation of the physical end state l = 1.0. A soft-core Lennard Jones potential was used to scale all nonbonded interactions for the alchemical substituents by l (94). Substituent dihedral angles were scaled by l, whereas bonds, angles, and improper dihedral angles were not, as this was found previously to yield better sampling without causing conformations to become trapped in local energy minima (95,97). The adaptive landscape flattening (ALF) algorithm was used to identify appropriate biasing potentials for MSlD (94,95). ALF was run in two stages: first with several 100-ps simulations, followed by additional 1-ns simulations for bias refinement. The unbound Pup2 at 303.15 K ran 176 100-ps simulations followed by 18 1-ns simulations. ALF biases for the complex at 303.15 K were initialized with estimates from preliminary single-mutant simulations, and then 221 100-ps simulations and 14 1-ns simulations were run. Unfortunately, despite these great lengths of sampling with ALF and MSlD, converged biases that yielded consistent sampling of all alchemical end states were elusive. This can be attributed to the difficulty of exploring small-to-large perturbations within the hydrophobic core of Pup2. To circumvent this difficulty, the best biases from unbound and bound Pup2 were then used in one 5-ns production run followed by three sequential 25ns production runs employing variable-bias replica exchange (VB-REX), an enhanced sampling algorithm (95). The unbound Pup2 at 323.15 K ran 500 100-ps simulations and 44 1-ns simulations, all using VB-REX for bias refinement. The bound Pup2 at 323.15 K ran 145 100-ps simulations using VB-REX for bias refinement. The best biases from unbound and bound Pup2 at 323.15 K were then used in three sequential 25-ns production runs employing VB-REX. Three to five duplicate runs with three replicas each were used in the 25-ns production iterations, and ALF was used in between iterations to compute biases for the next set of calculations. Ultimately, the last set of 25-ns production runs were used for computing differences in binding free energies between WT and mutant variants of Pup2. Thus, for the Thermal profiling of missense mutants final free energies reported in Table 1, a total of 225 ns of sampling was performed for both solvated and complexed states of Pup2 at 303.15 K, and a total of 375 ns of sampling was performed for both the solvated and complexed state of Pup2 at 323.15 K. MSlD trajectories were then analyzed with the assistance of VMD and PyMOL (132,144).

Data availability
All analyzed data sets used for this study are referenced in the article and included in the supporting information. The MS proteomics data have been deposited to the ProteomeXchange Consortium (145) via the PRIDE (146) partner repository with the data set identifiers PXD017223 (global) and PXD017222 (mTPP). The RNA-Seq data have been deposited to Gene Expression Omnibus (GEO) under the accession number GSE143236. Any code used for data analysis is available as referenced. Thermo Proteome Discoverer software is a commercial product; however, a free viewer is available for processed data files at https://thermo.flexnetoperations.com/.
Funding and additional information-This work was supported in part by National Institutes of Health Grant T32 HL007910 (to S. A. P. J.) and by the Showalter Research Trust (to A. L. M.). This work was also supported in part by the Indiana Clinical and Translational Sciences Institute, which is funded by National Institutes of Health Grant UL1TR002529, a National Center for Advancing Translational Sciences, Clinical, and Translational Sciences Award. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Conflict of interest-The authors declare that they have no conflicts of interest with the contents of this article.