Distinct stress conditions result in aggregation of proteins with similar properties

1Protein aggregation is the abnormal association of misfolded proteins into larger, often insoluble structures. Aggregation can be classified into two general categories: amyloid and amorphous. The amyloid state is a highly structured, insoluble, fibrillar deposit, usually consisting of many repeats of the same protein. This type of aggregation is central in the pathology of many neurodegenerative diseases including Alzheimer’s, Parkinson’s and Huntington’s disease. Amorphous protein aggregation can be best described as the apparently unordered aggregation of proteins, with each individual protein not generally associated with disease when aggregated. Protein misfolding that leads to aggregation can arise as a consequence of age, environmental stress, chemical modifications, destabilising mutations or lack of oligomeric assembly partners. Newly synthesized proteins appear to be particularly vulnerable to misfolding events and widespread protein aggregation is thought to be toxic, especially when the proteostasis network is compromised.

2Proteins are highly dynamic molecules, where various modifications or changes in the cellular environment can affect their native conformational fold. Whilst conformational flexibility is required for many proteins to function biologically, aberrant conformations, or misfolding, can lead to protein aggregation. Various stress conditions, such as high temperature, heavy metals and oxidative stress may cause protein misfolding and aggregation by shifting the conformational equilibrium towards more aggregation-prone states where exposed hydrophobic regions of misfolded proteins can interact with other exposed hydrophobic regions leading to aberrant protein-protein interactions. Therefore, to ensure protein homeostasis, the cell contains an arsenal of molecular chaperones that are able to detect non-native misfolded proteins and act upon them to prevent aggregation or amyloid formation. Generally, it is believed that chaperones are able to target unfolded proteins via recognition of hydrophobic stretches that would otherwise be buried within the native fold and therefore protected from the external environment. In addition, there are ATP-dependent chaperone classes that are involved in co-translational folding of nascent polypeptides and refolding of proteins (i.e. Hsp70, chaperonins, Hsp90). In the case of Hsp70 s, a co-factor (Hsp40 class chaperones) first assists in recruiting Hsp70 to substrates and then stimulates the ATPase activity of Hsp70 to drive substrate refolding. ATP-independent chaperones, such as the small Hsp class of chaperones, exhibit a ‘holdase’ function by binding to misfolded proteins, preventing their aggregation. Finally, certain chaperones contain a disaggregase function, such as the fungal specific Hsp104. Rather than preventing aggregation, these chaperones act to disassemble protein aggregates.

3In a previous study, we isolated aggregation-prone proteins under physiological conditions and arsenite stress and used bioinformatic analyses to identify characteristics that are linked to protein aggregation in living yeast cells. We found that these proteins have high translation rates and are substrates of ribosome-associated Hsp70 chaperones, indicating that they are susceptible to aggregation primarily during translation/folding. The toxic metalloid arsenite promotes protein aggregation by interfering with the folding of nascent polypeptides and by chaperone inhibition. Moreover, bioinformatic analysis of arsenite-induced aggregates suggests that arsenite stress lowers the general threshold for protein aggregation. Together, these studies suggest that protein aggregation is a normal physiological event, but conditions which perturb cellular homeostasis can increase the burden of protein aggregation.

4In this current study, we have extended our analysis of protein aggregation to include two additional stress conditions – azetidine-2-carboxylic acid (AZC) and hydrogen peroxide (H2O2) – that promote misfolding and aggregation through different mechanisms. AZC is a proline analogue which is competitively incorporated into proteins in place of proline. AZC incorporation alters the conformation of the polypeptide backbone, resulting in widespread protein misfolding and aggregation. H2O2 is a ubiquitous stress agent that is formed as a byproduct of aerobic respiration and following exposure to diverse biological and environmental factors. H2O2 gives rise to oxidative stress in cells which may in turn damage proteins and promote protein aggregation. Using computational approaches, we characterized the proteins that aggregate following AZC and H2O2 stress and compared them with proteins which aggregate during arsenite stress or during physiological conditions. Our data indicate that the three stress conditions, which work by distinct mechanisms, promote the aggregation of similar types of proteins, probably by lowering the threshold of protein aggregation. This suggests that the proteins in aggregates are intrinsically aggregation-prone, rather than being proteins which are affected in a stress-specific manner. Most proteins are susceptible for aggregation during synthesis/folding. In addition, certain proteins may aggregate post-translationally due to an imbalance between abundance and solubility.

5To extend our previous analyses of protein aggregation, protein aggregates were isolated and identified following H2O2 and AZC stress. Insoluble protein aggregates were prepared as previously described and identified using mass spectrometry following three independent experiments (see Materials and Methods). The H2O2-set and AZC-set were compared with our previously identified set of proteins which aggregate following arsenite stress (As-set). We noted that a considerable fraction of the proteins (~35%) aggregated in response to more than one stress condition (Fig. 1). Therefore, to facilitate a true comparative analysis of the stress-induced protein aggregates, we partitioned the identified proteins into non-overlapping datasets; 45 proteins uniquely identified in the As-set, 140 proteins within the AZC-set, 53 proteins within the H2O2-set, and a stress-set (Common-set) which contains 128 proteins that aggregate in at least two of the three stress conditions (Fig. 1).

6We first performed gene ontology analysis to examine what functional categories of proteins are enriched in the aggregate fractions following the three distinct stress conditions. For this and all subsequent analyses, our stress datasets were compared with an unstressed dataset, which is comprised of proteins that aggregate only under normal physiological conditions (Unstressed-set). Significantly enriched (5% FDR) functional categories were determined within the datasets using the MIPS Functional Catalogue. As we previously described, factors involved in protein synthesis including ribosomal and translation related proteins are strongly enriched within proteins that aggregate in the absence of stress (Fig. 2; Unstressed dataset). Additionally, proteins involved in energy and transport functions are enriched within these aggregates. More functional categories were enriched in the Common-set compared with the Unstressed-set; these include many protein synthesis related functions, as well as proteins involved in metabolism and energy related processes. There was also enrichment for proteins involved in protein folding, stabilisation and processing, as well as components of the unfolded protein response. These latter classes of proteins would be expected to constitute part of the cellular response to protein misfolding and aggregation. Stress-specific differences were found in the functional classes that are enriched under different stress conditions. The As-specific set was significantly enriched for proteins related to protein synthesis and translation (Fig. 2), in line with the notion that arsenite interferes with folding of nascent polypeptides. The AZC-specific set was enriched in a large number of categories, including metabolism, energy and protein synthesis-related functions, as well as cell rescue and defence proteins including a number of chaperones (Fig. 2). The large number of functional groups enriched in the AZC-specific data-set may be a reflection of its mode of action, as AZC will affect all proline-containing proteins. In contrast to the other sets, no functional groups were significantly enriched within the H2O2-specific set. Taken together, these data indicate that protein aggregates isolated from different conditions show enrichment for a number of similar functional categories.

7Our previous observations suggested that arsenite stress lowers the overall threshold for protein aggregation. To determine whether this is also true for other stresses that promote protein aggregation, we assessed a number of physicochemical properties of the proteins within our datasets. For comparison, a list of yeast proteins detectable by mass spectrometry in logarithmically growing cells was used to represent the properties of unaggregated proteins. Aggregated proteins in the Unstressed-set are more abundant (i.e. present in more molecules/cell), more highly expressed (indicated by a high codon adaptation index (CAI)), smaller in size (i.e. lower molecular weight (MW)), and have a higher isoelectric point (pI) than proteins in the Unaggregated set (Fig. 3a–d). Similar to the Unstressed-set, highly expressed and abundant proteins are significantly enriched in the aggregate fractions following all three stress conditions (Fig. 3a,b). However, the proteins which aggregate under stress conditions have considerably lower abundance and expression levels compared with the Unstressed-set (Fig. 3a,b). Proteins that aggregate following the three stress conditions also have a lower pI than the proteins in the Unstressed set; their pIs are similar to the Unaggregated-set (H2O2-set and As-set) or lower (Common-set and AZC-set) (Fig. 3d). Stress-specific differences are also found when the sizes of the aggregated proteins are compared (Fig. 3c). Whilst the proteins in the Unstressed-set are generally smaller in size than the unaggregated set, proteins which aggregate in the Common-set, H2O2-set and AZC-set are significantly larger than those in the Unstressed-set. Thus, the proteins that aggregate following stress conditions have lower expression levels, are less abundant, are more acidic and are larger than the proteins which aggregate during physiological non-stress conditions.

8Given the increased abundance of aggregated proteins during both unstressed and stressed conditions compared to the Unaggregated set, we examined whether this correlates with protein stability. For this analysis, we used a data-set of protein half-lives as determined by measuring protein abundance over time after inhibition of protein biosynthesis. The proteins which aggregate under non-stress (Unstressed-set) or arsenite stress (As-set) have on the average a longer half-life than the proteins in the Unaggregated set (Fig. 3e). However, no significant differences in protein stability are observed for the proteins in the Common-, H2O2- or AZC-set compared with the Unaggregated set (Fig. 3e). Thus, increased protein stability does not appear to account for the increased abundance of the proteins which tend to aggregate.

9Finally, we investigated the hydrophobicity (GRAVY score) of the aggregated proteins. Proteins in the Unstressed-set, the Common-set and the AZC-set are generally more hydrophobic than proteins in the Unaggregated-set (Fig. 3f). However, the proteins that aggregate under stress are less hydrophobic than those that aggregate in the absence of stress. We conclude that proteins that aggregate during physiological conditions and stress share several features, and that stress conditions shift the criteria for protein aggregation propensity.

10We next examined the relative amino acid content of the proteins enriched in our aggregate fractions. Proteins in the Unstressed-set are enriched in aliphatic amino acids including Ala, Gly and Val compared with the Unaggregated-set (Fig. 4), in agreement with their hydrophobic character (Fig. 3f). Furthermore, basic amino acids (Lys and Arg) are strongly enriched, whereas acidic amino acids (Asp, Glu) are underrepresented in the Unstressed-set (Fig. 4), which accounts for their higher pI compared with the Unaggregated-set (Fig. 3d). A number of amino acids are underrepresented in the Unstressed-set including Pro, Ser, His and sulphur-containing amino acids (Met, Cys) (Fig. 4). We also found that Gln and Asn are underrepresented in the Unstressed-set suggesting that the proteins in these aggregates are distinct from the well-known amyloid forming proteins. These amino acids are normally thought to underlie amyloid formation and have been linked to prion formation in yeast and mammalian neurological disorders including Huntington’s disease.

11The stress-aggregated proteins also show differences in their amino acid content compared to the Unaggregated protein set (Fig. 4). However, we did not observe any strong correlations with amino acids that are known to be targeted by the different stress conditions. For example, the AZC-set is not enriched in proline residues suggesting that they are not simply proteins where excess AZC is incorporated. Arsenite and H2O2 are known to target cysteine-containing proteins as part of their mode of toxicity, but the relative cysteine content of proteins within the As- and H2O2-set is not enriched (Fig. 4). This suggests that the proteins identified during stress conditions represent intrinsically aggregation-prone proteins. In agreement with this hypothesis, and similar to the Unstressed-set, the Common-set is enriched for hydrophobic proteins (Fig. 3f) and aliphatic amino acids including Ala, Gly, Ile, Val (Fig. 4). However, stress-specific differences are seen in the content of hydrophobic proteins (and aliphatic amino acids) as the As- and H2O2-specific sets are not significantly enriched for hydrophobic proteins compared with the Unaggregated-set (Fig. 3f). We note though that some aliphatic amino acids (Ala, Val) are enriched in the As-set (Fig. 4), indicating a weak correlation between hydrophobicity and As-specific aggregation. A number of amino acids that are underrepresented in the Unstressed-set are also underrepresented in the Common-set including Asn, Gln, Ser, His and Met. Thus, while a high content of aliphatic amino acids is positively correlated with protein aggregation, a high content of Asn, Gln, Ser, His and Met might be negatively correlated with physiological and stress-induced aggregation.

12Given that the aggregated proteins identified under both non-stress and stress conditions are highly expressed and abundant, we compared their translation rates with unaggregated proteins. For this analysis we compared the aggregated proteins with a genome-wide estimate of translation rates. Proteins that aggregate under unstressed conditions are significantly enriched for proteins which show high rates of mRNA translation compared with unaggregated proteins (Fig. 5a). Similarly, in agreement with their higher protein abundances and expression levels (Fig. 3a,b), translation rates are also significantly increased in the Common-set as well as the As- and AZC-specific sets (Fig. 5A). Nevertheless, these proteins have considerably lower translation rates compared with the proteins in the Unstressed-set (Fig. 5a). In contrast, translation rates are not increased for the H2O2-set despite their enrichment for higher abundance and expression proteins (Fig. 5a).

13It is now well established that many proteins are subject to co-translational folding, predominantly mediated by Hsp70 family chaperones. For example, Ssb1 and Ssb2 are ribosome-associated chaperones that are important in the folding of nascent polypeptide chains. We therefore examined whether aggregated proteins are enriched for co-translational substrates of Ssb chaperones using available Ssb2 data. In accordance with their high translation rates, proteins in the As-, AZC-, Common- and Unstressed-set are significantly enriched in proteins that are co-translational Ssb2 substrates compared to Unaggregated proteins (Fig. 5B). In contrast, the H2O2-specific set is not enriched for co-translational Ssb2 substrates. Taken together, these findings suggest that proteins are primarily susceptible for aggregation during translation/folding, both during physiological conditions and during stress. The exception appears to be H2O2 (see Discussion).

14We previously reported that proteins which aggregate following arsenite stress are significantly enriched in proteins with more chaperone interactions per protein compared with proteins in the proteome. This analysis was performed using a chaperone-protein interaction atlas for 63 chaperones present in yeast. When this analysis was repeated here using the non-overlapping datasets, we found that proteins that aggregate following stress conditions (Common-set) are significantly enriched for multiple chaperone interactions (Fig. 5c,d). The As-, H2O2-, and AZC-specific sets and the Unstressed-set were not enriched, nor depleted, for chaperone-interacting proteins compared to the Unaggregated-set.

15Due to the importance of Hsp70 s in co-translational folding of nascent polypeptides and refolding of proteins, we next asked whether the aggregated protein sets are enriched for interactions with specific members of the Hsp70 family. For this, we chose seven cytoplasmic (Ssa1–4, Sse1–2, Ssz1), two ribosome-associated (Ssb1, Ssb2), two ER (Kar2, Lhs1) and three mitochondrial (Ssc1, Ssq1, Ecm10) chaperones (Fig. 5e) from the chaperone-protein interaction atlas. The Common-set is enriched for Ssa1, Ssa2, Ssb1, Ssb2, Sse1, Ssq1, and Ssz1 interacting proteins compared to the Unaggregated-set (Fig. 5E). The Unstressed- and AZC-specific sets were both enriched for Ssb2 interactions whilst the H2O2-set was enriched for Ssz1 interactions. The Unstressed-set was also significantly underrepresented for Ssa1, Ssa2, and Ssb1 interactions (Fig. 5e), possibly because the Unstressed-set contains many ribosomal proteins that require specialized chaperones for proper folding. Thus, Hsp70 chaperone-interacting proteins therefore appear to be significantly enriched in stress-induced aggregates compared with unstressed protein aggregates. Few stress-specific interactions were observed, except for the ribosome associated Ssz1-interacting proteins which is enriched in the H2O2-set. None of the aggregated protein sets were significantly enriched for interactions with Ecm10, Kar2, Lhs1, Ssa3, Ssa4, Ssc1 and Sse2 compared with the Unaggregated set. Taken together, stress-aggregated proteins are enriched for multiple chaperone interactions.

16Given that proteins with multiple chaperone interactions are enriched within stress-aggregated sets, we examined whether chaperones were isolated in our aggregate fractions. It should be emphasised that our analysis does not allow us to differentiate between chaperones which are functional components of the aggregates, versus chaperones which are themselves aggregation-prone. Of the 63 known chaperones in S. cerevisiae, we identified 30 chaperones distributed between all the datasets: 11 Hsp70 s, five Hsp40 s, seven chaperonin subunits, three AAA+ family members, two Hsp90 s, one Hsp60 and one small Hsp chaperone (Fig. 6). A total of seven chaperones are present in the Unstressed-set: Ssb1, Ssa1, Ssa4, Ssc1, Hsp82, Hsc82 and Sec63. We also identified chaperones within the non-overlapping stress sets. The Common-set included 19 chaperones spanning all six chaperone classes (Fig. 6). All of these chaperones except Ssc1, Hsp78, Hsp60 and Sec63, are present in the cytoplasm. Their inclusion within the Common-set may indicate that these chaperones are part of a general cytoplasmic response to stress and that they are associated with their client proteins. Seven chaperones were identified within the AZC-specific set; four of these (Ssa2, Ssa3, Sse2, Kar2) belong to the Hsp70 family, whilst one (Mdj1) is a known co-factor of proteins of the Hsp70 family (Fig. 6). This may be indicative of the mechanism of action of AZC, as Hsp70 family proteins are important for co-translational folding. Two chaperones were present in the As-specific set (Zuo1 and Ssb2), and both are ribosome-associated chaperones (Fig. 6). This is in agreement with the notion that arsenite primarily targets nascent proteins for aggregation. The H2O2-specific set contained two chaperones (Cct6 and Mcx1) but it is currently unclear how these chaperones relate to H2O2’s mode of action. The H2O2-set was enriched for proteins that interact with Ssz1 – however, this chaperone is absent from the aggregates themselves. This may suggest that H2O2 might inhibit Ssz1 function.

17Protein aggregation is a hallmark of many ageing related diseases. A recent study identified 480 proteins that aggregate during postmitotic ageing in yeast. When we compared these proteins with our datasets we found no significant overlap with our Unstressed-set. This makes sense given that our unstressed proteins were identified in exponential phase yeast cells rather than in an aged population. Interestingly, a significant overlap (1.4-fold, p = <0.001) was found between the ageing dataset and the Common-set, although there was no significant overlap for any single stress-specific condition (Fig. 7). Therefore, it appears that proteins which generally aggregate in response to stress are also likely to aggregate in aged yeast cells.

18Finally, we wanted to compare our datasets with proteins that aggregate in another organism. A previous study identified 461 proteins in aged C. elegans protein aggregates. Of these 461 proteins, 120 have a recognisable yeast orthologue which we used for comparative analysis with our datasets (referred to as CE-set). Within our protein aggregate datasets, we found that 126 (Unstressed-set), 30 (As-set), 25 (H2O2-set), 83 (AZC-set) and 84 (Common-set) proteins have recognisable C. elegans orthologues. Our analysis revealed a significant overlap between the yeast stress datasets (Common-set) and the CE-set. Overall, 69 (57.5%) of the proteins in our stress dataset are present in the CE-set (Fig. 8). 26 of the 126 orthologous proteins in the Unstressed-set overlap with the CE-set (22% of CE-set; 2.3 fold above the expected overlap; p = <0.001) (Fig. 8). Of the 26 overlapping proteins, 18 are ribosomal proteins and two are Hsp70 family chaperones. Significant overlaps are also found with the AZC- and Common-sets. Notable proteins, in which their presence in aggregates is conserved between our stress datasets and the CE-set, include components of the essential chaperonin Ctt ring complex. Subunits, Tcp1, Cct3, Cct4 and Cct8 are present in all our stress sets and have like-for-like orthologues in C.elegans, which are present in the CE-set. Another essential protein, Cdc48, is also conserved in all stress induced and CE-set aggregates. C. elegans has two orthologous of Cdc48 (cdc-48.1 and cdc-48.2), of which both forms are present in CE-set.

19In this study, we used computational approaches to characterize proteins that aggregate in the absence of stress and under three distinct stress conditions (arsenite, H2O2 and AZC). We found that proteins that aggregate during physiological conditions and during stress are generally more abundant, highly expressed and translated at faster rates, when compared to the wider proteome. This extends our previous finding that high protein abundance positively correlates with an increased aggregation propensity. Previous studies have predicted that aggregation propensity is not correlated with protein expression and abundance, i.e. highly expressed and abundant proteins have evolved to be both highly soluble and resistant to aggregation. However, it was noted that, although proteins are expressed at a level to allow functionality in balance with their intrinsic aggregation propensity, there is almost no flexibility with this equilibrium. Hence, any factors that would decrease solubility or increase the concentration of a protein would result in unavoidable aggregation. Indeed, a recent study has identified a number of proteins which are maintained at a high concentration relative to their solubility. It was proposed that these ‘supersaturated’ proteins are highly dependent on the proteostasis network to maintain their solubility, and that any perturbations in this network will shift supersaturated proteins towards aggregation. For example, the proteostasis network is believed to decline during ageing and protein aggregates isolated from aged C. elegans were found to be enriched for supersaturated proteins. Protein abundance alone can therefore be a good indicator of protein aggregation in vivo. Our results are in striking agreement with this hypothesis; we found that proteins that aggregate during physiological unstressed conditions and during stress are generally more abundant, highly expressed and translated at faster rates, when compared to the wider proteome.

20Our data indicate that the three stress conditions, which each work by distinct mechanisms, promote the aggregation of similar types of proteins; these proteins are abundant, highly expressed, and translated at high rates. Stress-aggregated proteins tend to be somewhat hydrophobic, more acidic and larger in size compared to unaggregated proteins. As they share many features, these proteins are likely intrinsically aggregation-prone, rather than being proteins which are affected in a stress-specific manner. Our current study further revealed that cellular stress conditions act to lower the threshold of protein aggregation. Notably, the abundance, expression levels and translation rates of the stress-aggregated proteins were not as high as for the proteins that aggregate during physiological conditions. Likewise, proteins that aggregate during environmental stress were less hydrophobic than those aggregating during physiological conditions. Thus, proteins that are normally not susceptible to aggregation become aggregation-prone during stress.

21A large fraction of the aggregated proteins are co-translational substrates of ribosome-associated Ssb2. Moreover, the stress aggregated proteins (Common-set) are enriched for chaperone interactions including several Hsp70-interacting proteins (Ssb1, Ssb2, Ssa1, Ssa2, Sse1, Ssq1, Ssz1). Thus, these proteins appear susceptible for aggregation primarily during synthesis/folding, indicating that stress conditions promote co-translational aggregation. The exception appears to be H2O2 which may be because oxidative stress strongly inhibits translation and hence translation rate is no longer differentiated in the proteins which aggregate following H2O2 stress. Alternatively, H2O2-stress may additionally target native proteins for aggregation. The unstressed dataset is translated at higher rates than the unaggregated and stressed datasets and is also enriched for proteins that interact with ribosome-bound Ssb2 indicating co-translational aggregation. However, we found that the Unstressed set is generally depleted of Hsp70 chaperones (Ssb1-, Ssa1- and Ssa2-interacting), but not for chaperones in general. Thus, physiological aggregates which form under non-stressed conditions may also contain a fraction of proteins that are less prone to co-translational aggregation. In line with the supersaturated protein theory, these proteins may aggregate when their abundance exceeds their solubility.

22It is well established that certain amino acids influence protein aggregation, and that proteins with solvent exposed stretches of high hydrophobicity and a low net charge are aggregation-prone. Indeed, aliphatic amino acids along with basic amino acids were over-represented in physiological aggregates (Unstressed set). Aliphatic amino acids were also enriched within stress-induced aggregates. However, in contrast to the Unstressed-set, the Common-set is enriched for acidic amino acids (Asp and Glu) rather than basic amino acids. Thus, proteins that aggregate under stress conditions are generally acidic whilst proteins that aggregate under physiological conditions are generally basic. It has been proposed that hydrophobic and electrostatic forces that promote the formation of functional protein complexes can also cause abnormal protein-protein interactions. Indeed, aggregation prone regions enriched in aliphatic amino acids appear to be important in protein-protein interactions. Interestingly, positively charged residues (including Arg and Lys) that flank aggregation prone regions strongly oppose close packing and aggregation. Hence, it is thought that these aggregation resistant flanks provide specificity to hydrophobic-mediated protein-protein interactions. In contrast, Glu and Asp appear to be ineffective in disrupting aggregation when compared to Lys and Arg. A large fraction of the aggregated proteins identified in this study are complex forming proteins (e.g. ribosomal proteins). Moreover, we previously showed that proteins aggregating during physiological conditions and arsenite stress are engaged in a significantly higher number of protein-protein interactions per protein than the proteins in the proteome. Thus, our data are in agreement with the notion that protein complex formation and aggregation may be governed by similar principles.

23We found a significant overlap between proteins that aggregate during yeast ageing and during stress (Common-set). Our analysis also revealed significant overlap between stress aggregating yeast proteins (Common-set) and proteins that aggregate in ageing C. elegans. There is also an overlap between ageing-dependent protein aggregation in C. elegans and disease-dependent aggregation in mammals. Several aggregation-prone yeast proteins have human homologues that are implicated in protein misfolding diseases. Thus, similar mechanisms may apply in disease- and non-disease settings and the factors and components that control protein aggregation may be evolutionary conserved.

24Following SD (untreated), AZC (5 mM) or H2O2 (1 mM) treatment for 2 hours, insoluble protein aggregates were isolated as previously described. Insoluble protein extracts were separated by reducing SDS/PAGE (12% gels) and visualized by silver staining with the Bio-Rad silver stain plus kit. Aggregated proteins were identified by mass spectrometry (performed by the Biomolecular Analysis Core Facility, Faculty of Life Sciences, The University of Manchester) in triplicate for each condition. For protein identification, protein samples were run a short distance into SDS-PAGE gels and stained using colloidal Coomassie blue (Sigma). Total proteins were excised, trypsin digested, and identified using liquid chromatography-mass spectrometry (LC-MS). Proteins were identified using the Mascot mass fingerprinting programme (www.matrixscience.com) to search the NCBInr and Swissprot databases. Final datasets for each condition were determined by selecting proteins that were identified in at least two of the three replicates.

25Datasets for each condition were assessed for functional enrichment (p-value > 0.01; 0.05 FDR) of functional categories (MIPS database) using FunCat (available at http://www.helmholtz-muenchen.de/en/ibis. Protein abundance data was retrieved from the PaxDB integrated dataset (available at http://pax-db.org). The proteins that aggregate in aged-yeast or C. elegans were retrieved from published sources and statistical significance of the overlap was evaluated with a hypergeometric test. Enrichment and significance of aggregation propensity under one, two or three stress conditions were estimated with Monte Carlo simulations, using as background a list of proteins which are expected to be identifiable with LC-MS. Venn diagrams and visualization of the distribution of protein hits between conditions were made using Venny (http://bioinfogp.cnb.csic.es/tools/venny/). Physicochemical data, translation rates and chaperone interactions were evaluated with pair wise Mann-Whitney-Wilcoxon U-tests, followed by Holm-Bonferroni adjustments to avoid inflation of Type I error rate. Amino acid content was assessed statistically with Mann-Whitney-Wilcoxon U-test, and p-values were filtered on 0.05 FDR.

26How to cite this article: Weids, A. J. et al. Distinct stress conditions result in aggregation of proteins with similar properties. Sci. Rep. 6, 24554; doi: 10.1038/srep24554 (2016).