A Two-step Protein Quality Control Pathway for a Misfolded DJ-1 Variant in Fission Yeast

Synonymous mutations in the phosphoglycerate kinase 1 gene induce an altered response to protein misfolding in Schizosaccharomyces pombe

Background

Proteostasis refers to the processes that regulate the biogenesis, folding, trafficking, and degradation of proteins. Any alteration in these processes can lead to cell malfunction. Protein synthesis, a key proteostatic process, is highly-regulated at multiple levels to ensure adequate adaptation to environmental and physiological challenges such as different stressors, proteotoxic conditions and aging, among other factors. Because alterations in protein translation can lead to protein misfolding, examining how protein translation is regulated may also help to elucidate in part how proteostasis is controlled. Codon usage bias has been implicated in the fine-tuning of translation rate, as more-frequent codons might be read faster than their less-frequent counterparts. Thus, alterations in codon usage due to synonymous mutations may alter translation kinetics and thereby affect the folding of the nascent polypeptide, without altering its primary structure. To date, it has been difficult to predict the effect of synonymous mutations on protein folding and cellular fitness due to a scarcity of relevant data. Thus, the purpose of this work was to assess the effect of synonymous mutations in discrete regions of the gene that encodes the highly-expressed enzyme 3-phosphoglycerate kinase 1 (pgk1) in the fission yeast Schizosaccharomyces pombe.

Results

By means of systematic replacement of synonymous codons along pgk1, we found slightly-altered protein folding and activity in a region-specific manner. However, alterations in protein aggregation, heat stress as well as changes in proteasome activity occurred independently of the mutated region. Concomitantly, reduced mRNA levels of the chaperones Hsp9 and Hsp16 were observed.

Conclusion

Taken together, these data suggest that codon usage bias of the gene encoding this highly-expressed protein is an important regulator of protein function and proteostasis.

Disease-linked mutations cause exposure of a protein quality control degron

More than half of disease-causing missense variants are thought to lead to protein degradation, but the molecular mechanism of how these variants are recognized by the cell remains enigmatic. Degrons are stretches of amino acids that help mediate recognition by E3 ligases and thus confer protein degradation via the ubiquitin-proteasome system. While degrons that mediate controlled degradation of, for example, signaling components and cell-cycle regulators are well described, so-called protein-quality-control degrons that mediate the degradation of destabilized proteins are poorly understood. Here, we show that disease-linked dihydrofolate reductase (DHFR) missense variants are structurally destabilized and chaperone-dependent proteasome targets. We find two regions in DHFR that act as degrons, and the proteasomal turnover of one of these was dependent on the molecular chaperone Hsp70. Structural analyses by nuclear magnetic resonance (NMR) and hydrogen/deuterium exchange revealed that this degron is buried in wild-type DHFR but becomes transiently exposed in the disease-linked missense variants.

MutateX: an automated pipeline for in silico saturation mutagenesis of protein structures and structural ensembles

Mutations, which result in amino acid substitutions, influence the stability of proteins and their binding to biomolecules. A molecular understanding of the effects of protein mutations is both of biotechnological and medical relevance. Empirical free energy functions that quickly estimate the free energy change upon mutation (ΔΔG) can be exploited for systematic screenings of proteins and protein complexes. In silico saturation mutagenesis can guide the design of new experiments or rationalize the consequences of known mutations. Often software such as FoldX, while fast and reliable, lack the necessary automation features to apply them in a high-throughput manner. We introduce MutateX, a software to automate the prediction of ΔΔGs associated with the systematic mutation of each residue within a protein, or protein complex to all other possible residue types, using the FoldX energy function. MutateX also supports ΔΔG calculations over protein ensembles, upon post-translational modifications and in multimeric assemblies. At the heart of MutateX lies an automated pipeline engine that handles input preparation, parallelization and outputs publication-ready figures. We illustrate the MutateX protocol applied to different case studies. The results of the high-throughput scan provided by our tools can help in different applications, such as the analysis of disease-associated mutations, to complement experimental deep mutational scans, or assist the design of variants for industrial applications. MutateX is a collection of Python tools that relies on open-source libraries. It is available free of charge under the GNU General Public License from https://github.com/ELELAB/mutatex.

Chaperone-Facilitated Aggregation of Thermo-Sensitive Proteins Shields Them from Degradation during Heat Stress

Cells have developed protein quality-control strategies to manage the accumulation of misfolded substrates during heat stress. Using a soluble reporter of misfolding in fission yeast, Rho1.C17R-GFP, we demonstrate that upon mild heat shock, the reporter collapses in protein aggregate centers (PACs). They contain and/or require several chaperones, such as Hsp104, Hsp16, and the Hsp40/70 couple Mas5/Ssa2. Stress granules do not assemble at mild temperatures and, therefore, are not required for PAC formation; on the contrary, PACs may serve as nucleation centers for the assembly of stress granules. In contrast to the general belief, the dominant fate of these PACs is not degradation, and the aggregated reporter can be disassembled by chaperones and recovers native structure and activity. Using mass spectrometry, we show that thermo-unstable endogenous proteins form PACs as well. In conclusion, formation of PACs during heat shock is a chaperone-mediated adaptation strategy.

Folliculin variants linked to Birt-Hogg-Dubé syndrome are targeted for proteasomal degradation

Germline mutations in the folliculin (FLCN) tumor suppressor gene are linked to Birt-Hogg-Dubé (BHD) syndrome, a dominantly inherited genetic disease characterized by predisposition to fibrofolliculomas, lung cysts, and renal cancer. Most BHD-linked FLCN variants include large deletions and splice site aberrations predicted to cause loss of function. The mechanisms by which missense variants and short in-frame deletions in FLCN trigger disease are unknown. Here, we present an integrated computational and experimental study that reveals that the majority of such disease-causing FLCN variants cause loss of function due to proteasomal degradation of the encoded FLCN protein, rather than directly ablating FLCN function. Accordingly, several different single-site FLCN variants are present at strongly reduced levels in cells. In line with our finding that FLCN variants are protein quality control targets, several are also highly insoluble and fail to associate with the FLCN-binding partners FNIP1 and FNIP2. The lack of FLCN binding leads to rapid proteasomal degradation of FNIP1 and FNIP2. Half of the tested FLCN variants are mislocalized in cells, and one variant (ΔE510) forms perinuclear protein aggregates. A yeast-based stability screen revealed that the deubiquitylating enzyme Ubp15/USP7 and molecular chaperones regulate the turnover of the FLCN variants. Lowering the temperature led to a stabilization of two FLCN missense proteins, and for one (R362C), function was re-established at low temperature. In conclusion, we propose that most BHD-linked FLCN missense variants and small in-frame deletions operate by causing misfolding and degradation of the FLCN protein, and that stabilization and resulting restoration of function may hold therapeutic potential of certain disease-linked variants. Our computational saturation scan encompassing both missense variants and single site deletions in FLCN may allow classification of rare FLCN variants of uncertain clinical significance.

Computational and cellular studies reveal structural destabilization and degradation of MLH1 variants in Lynch syndrome

Defective mismatch repair leads to increased mutation rates, and germline loss-of-function variants in the repair component MLH1 cause the hereditary cancer predisposition disorder known as Lynch syndrome. Early diagnosis is important, but complicated by many variants being of unknown significance. Here we show that a majority of the disease-linked MLH1 variants we studied are present at reduced cellular levels. We show that destabilized MLH1 variants are targeted for chaperone-assisted proteasomal degradation, resulting also in degradation of co-factors PMS1 and PMS2. In silico saturation mutagenesis and computational predictions of thermodynamic stability of MLH1 missense variants revealed a correlation between structural destabilization, reduced steady-state levels and loss-of-function. Thus, we suggest that loss of stability and cellular degradation is an important mechanism underlying many MLH1 variants in Lynch syndrome. Combined with analyses of conservation, the thermodynamic stability predictions separate disease-linked from benign MLH1 variants, and therefore hold potential for Lynch syndrome diagnostics.

Toward mechanistic models for genotype-phenotype correlations in phenylketonuria using protein stability calculations

Phenylketonuria (PKU) is a genetic disorder caused by variants in the gene encoding phenylalanine hydroxylase (PAH), resulting in accumulation of phenylalanine to neurotoxic levels. Here, we analyzed the cellular stability, localization, and interaction with wild-type PAH of 20 selected PKU-linked PAH protein missense variants. Several were present at reduced levels in human cells, and the levels increased in the presence of a proteasome inhibitor, indicating that proteins are proteasome targets. We found that all the tested PAH variants retained their ability to associate with wild-type PAH, and none formed aggregates, suggesting that they are only mildly destabilized in structure. In all cases, PAH variants were stabilized by the cofactor tetrahydrobiopterin (BH₄ ), a molecule known to alleviate symptoms in certain PKU patients. Biophysical calculations on all possible single-site missense variants using the full-length structure of PAH revealed a strong correlation between the predicted protein stability and the observed stability in cells. This observation rationalizes previously observed correlations between predicted loss of protein destabilization and disease severity, a correlation that we also observed using new calculations. We thus propose that many disease-linked PAH variants are structurally destabilized, which in turn leads to proteasomal degradation and insufficient amounts of cellular PAH protein.

Protein stability and degradation in health and disease

The cellular proteome performs highly varied functions to sustain life. Since most of these functions require proteins to fold properly, they can be impaired by mutations that affect protein structure, leading to diseases such as Alzheimer's disease, cystic fibrosis, and Lynch syndrome. The cell has evolved an intricate protein quality control (PQC) system that includes degradation pathways and a multitude of molecular chaperones and co-chaperones, all working together to catalyze the refolding or removal of aberrant proteins. Thus, the PQC system limits the harmful consequences of dysfunctional proteins, including those arising from disease-causing mutations. This complex system is still not fully understood. In particular the structural and sequence motifs that, when exposed, trigger degradation of misfolded proteins are currently under investigation. Moreover, several attempts are being made to activate or inhibit parts of the PQC system as a treatment for diseases. Here, we briefly review the present knowledge on the PQC system and list current strategies that are employed to exploit the system in disease treatment.

Computational Structural Biology of S-nitrosylation of Cancer Targets

Nitric oxide (NO) plays an essential role in redox signaling in normal and pathological cellular conditions. In particular, it is well known to react in vivo with cysteines by the so-called S-nitrosylation reaction. S-nitrosylation is a selective and reversible post-translational modification that exerts a myriad of different effects, such as the modulation of protein conformation, activity, stability, and biological interaction networks. We have appreciated, over the last years, the role of S-nitrosylation in normal and disease conditions. In this context, structural and computational studies can help to dissect the complex and multifaceted role of this redox post-translational modification. In this review article, we summarized the current state-of-the-art on the mechanism of S-nitrosylation, along with the structural and computational studies that have helped to unveil its effects and biological roles. We also discussed the need to move new steps forward especially in the direction of employing computational structural biology to address the molecular and atomistic details of S-nitrosylation. Indeed, this redox modification has been so far an underappreciated redox post-translational modification by the computational biochemistry community. In our review, we primarily focus on S-nitrosylated proteins that are attractive cancer targets due to the emerging relevance of this redox modification in a cancer setting.

Amyotrophic lateral sclerosis as a protein level, non-genomic disease: Therapy with S2RM exosome released molecules

Amyotrophic lateral sclerosis (ALS) is a rapidly progressing neurodegenerative disease that leads to death. No effective treatments are currently available. Based on data from epidemiological, etiological, laboratory, and clinical studies, I offer a new way of thinking about ALS and its treatment. This paper describes a host of extrinsic factors, including the exposome, that disrupt the extracellular matrix and protein function such that a spreading, prion-like disease leads to neurodegeneration in the motor tracts. A treatment regimen is described using the stem cell released molecules from a number of types of adult stem cells to provide tissue dependent molecules that restore homeostasis, including proteostasis, in the ALS patient. Because stem cells themselves as a therapeutic are cumbersome and expensive, and when implanted in a host cause aging of the host tissue and often fail to engraft or remain viable, only the S2RM molecules are used. Rebuilding of the extracellular matrix and repair of the dysfunctional proteins in the ALS patient ensues.

A non-transcriptional role for the glucocorticoid receptor in mediating the cell stress response

The glucocorticoid receptor (GR) is essential for the stress response in mammals. We investigated potential non-transcriptional roles of GR in cellular stress response using fission yeast as a model.We surprisingly discovered marked heat stress resistance in yeast ectopically expressing human GR, which required expression of both the N-terminal transactivation domain, and the C-terminal ligand binding domain, but not the DNA-binding domain of the GR. This effect was not affected by GR ligand exposure, and occurred without significant GR nuclear accumulation. Mechanistically, the GR survival effect required Hsp104, and, indeed, GR expression increased Hsp104 expression. Proteomic analysis revealed GR binding to translasome components, including eIF3, a known partner for Sty1, a pattern of protein interaction which we confirmed using yeast two-hybrid studies.Taken together, we find evidence for a novel pathway conferring stress resistance in yeast that can be activated by the human GR, acting by protein-protein mechanisms in the cytoplasm. This suggests that in organisms where GR is natively expressed, GR likely contributes to stress responses through non-transcriptional mechanisms in addition to its well-established transcriptional responses.

Blocking protein quality control to counter hereditary cancers

Inhibitors of molecular chaperones and the ubiquitin-proteasome system have already been clinically implemented to counter certain cancers, including multiple myeloma and mantle cell lymphoma. The efficacy of this treatment relies on genomic alterations in cancer cells causing a proteostatic imbalance, which makes them more dependent on protein quality control (PQC) mechanisms than normal cells. Accordingly, blocking PQC, e.g. by proteasome inhibitors, may cause a lethal proteotoxic crisis in cancer cells, while leaving normal cells unaffected. Evidence, however, suggests that the PQC system operates by following a better-safe-than-sorry principle and is thus prone to target proteins that are only slightly structurally perturbed, but still functional. Accordingly, implementing PQC inhibitors may also, through an entirely different mechanism, hold potential for other cancers. Several inherited cancer susceptibility syndromes, such as Lynch syndrome and von Hippel-Lindau disease, are caused by missense mutations in tumor suppressor genes, and in some cases, the resulting amino acid substitutions in the encoded proteins cause the cellular PQC system to target them for degradation, although they may still retain function. As a consequence of this over-meticulous PQC mechanism, the cell may end up with an insufficient amount of the abnormal, but functional, protein, which in turn leads to a loss-of-function phenotype and manifestation of the disease. Increasing the amounts of such proteins by stabilizing with chemical chaperones, or by targeting molecular chaperones or the ubiquitin-proteasome system, may thus avert or delay the disease onset. Here, we review the potential of targeting the PQC system in hereditary cancer susceptibility syndromes.

The exocyst subunit Sec3 is regulated by a protein quality control pathway

Exocytosis involves fusion of secretory vesicles with the plasma membrane, thereby delivering membrane proteins to the cell surface and releasing material into the extracellular space. The tethering of the secretory vesicles before membrane fusion is mediated by the exocyst, an essential phylogenetically conserved octameric protein complex. Exocyst biogenesis is regulated by several processes, but the mechanisms by which the exocyst is degraded are unknown. Here, to unravel the components of the exocyst degradation pathway, we screened for extragenic suppressors of a temperature-sensitive fission yeast strain mutated in the exocyst subunit Sec3 (sec3-913). One of the suppressing DNAs encoded a truncated dominant-negative variant of the 26S proteasome subunit, Rpt2, indicating that exocyst degradation is controlled by the ubiquitin-proteasome system. The temperature-dependent growth defect of the sec3-913 strain was gene dosage-dependent and suppressed by blocking the proteasome, Hsp70-type molecular chaperones, the Pib1 E3 ubiquitin-protein ligase, and the deubiquitylating enzyme Ubp3. Moreover, defects in cell septation, exocytosis, and endocytosis in sec3 mutant strains were similarly alleviated by mutation of components in this pathway. We also found that, particularly under stress conditions, wild-type Sec3 degradation is regulated by Pib1 and the 26S proteasome. In conclusion, our results suggest that a cytosolic protein quality control pathway monitors folding and proteasome-dependent turnover of an exocyst subunit and, thereby, controls exocytosis in fission yeast.

Predicting the impact of Lynch syndrome-causing missense mutations from structural calculations

Accurate methods to assess the pathogenicity of mutations are needed to fully leverage the possibilities of genome sequencing in diagnosis. Current data-driven and bioinformatics approaches are, however, limited by the large number of new variations found in each newly sequenced genome, and often do not provide direct mechanistic insight. Here we demonstrate, for the first time, that saturation mutagenesis, biophysical modeling and co-variation analysis, performed in silico, can predict the abundance, metabolic stability, and function of proteins inside living cells. As a model system, we selected the human mismatch repair protein, MSH2, where missense variants are known to cause the hereditary cancer predisposition disease, known as Lynch syndrome. We show that the majority of disease-causing MSH2 mutations give rise to folding defects and proteasome-dependent degradation rather than inherent loss of function, and accordingly our in silico modeling data accurately identifies disease-causing mutations and outperforms the traditionally used genetic disease predictors. Thus, in conclusion, in silico biophysical modeling should be considered for making genotype-phenotype predictions and for diagnosis of Lynch syndrome, and perhaps other hereditary diseases.

The protein quality control system targets misfolded proteins for degradation. So far it has not been possible from sequence or structural data to predict the biological stability of a misfolded protein, or the effect of mutations on intracellular protein levels. Here we demonstrate that in silico saturation mutagenesis and biophysical calculations of the structural stability of the human mismatch repair protein MSH2 correlate with cellular protein levels, turnover and function. Of 24 different MSH2 variants, some of which are linked to Lynch syndrome, a destabilization of as little as 3 kcal/mol is sufficient to cause rapid degradation via the ubiquitin-proteasome pathway. Thus, biophysical modeling can, to a large extent, predict the metabolic stability of proteins. We also show that the same biophysical calculations can be used to distinguish with high accuracy neutral sequence variation from pathogenic variants, and that the calculations outperform several traditionally used disease predictors. We therefore suggest the method to be of potential value for patient stratification in Lynch syndrome, and perhaps other hereditary diseases.

Polyglutamine expansion of ataxin-3 alters its degree of ubiquitination and phosphorylation at specific sites

Ubiquitination and phosphorylation of proteins represent post translational modifications (PTMs) capable of regulating a variety of cellular processes. In the neurodegenerative disorder spinocerebellar ataxia type 3 (SCA3), the disease causing protein ataxin-3 carries an expanded polyglutamine (polyQ) stretch causing it to aggregate in nuclear inclusions. These inclusions are decorated with ubiquitin suggestive of a malfunction in the clearance of the mutant protein. Differences in ubiquitin chain topology between normal and polyQ expanded ataxin-3 could be involved in the differential clearance of the two proteins and play a role in SCA3 pathogenesis. Likewise, changes in phosphorylation patterns between the two variants could contribute to pathogenic processes involved in SCA3. We therefore determined the ubiquitination and phosphorylation patterns, together with the ubiquitin-linkage types associated with normal and polyQ expanded ataxin-3 by mass spectrometry (MS). This analysis revealed a similar ubiquitin linkage pattern on normal and expanded ataxin-3. However, the distribution of ubiquitinated lysine residues was altered in polyQ expanded ataxin-3, with increased ubiquitination at the new identified ubiquitination site lysine-8. MS analysis of phosphorylation also revealed novel phosphorylation sites in ataxin-3, and an increase in phosphorylation of expanded ataxin-3 at several positions. The study suggests that differences in clearance of normal and expanded ataxin-3 are not attributed to differences in ubiquitin-linkage pattern. However, the observed differences between the normal and polyQ expanded ataxin-3 with respect to the degree of ubiquitination and phosphorylation on specific sites may have an impact on ataxin-3 function and SCA3 pathogenesis.

UBL/BAG-domain co-chaperones cause cellular stress upon overexpression through constitutive activation of Hsf1

As a result of exposure to stress conditions, mutations, or defects during synthesis, cellular proteins are prone to misfold. To cope with such partially denatured proteins, cells mount a regulated transcriptional response involving the Hsf1 transcription factor, which drives the synthesis of molecular chaperones and other stress-relieving proteins. Here, we show that the fission yeast Schizosaccharomyces pombe orthologues of human BAG-1, Bag101, and Bag102, are Hsp70 co-chaperones that associate with 26S proteasomes. Only a subgroup of Hsp70-type chaperones, including Ssa1, Ssa2, and Sks2, binds Bag101 and Bag102 and key residues in the Hsp70 ATPase domains, required for interaction with Bag101 and Bag102, were identified. In humans, BAG-1 overexpression is typically observed in cancers. Overexpression of bag101 and bag102 in fission yeast leads to a strong growth defect caused by triggering Hsp70 to release and activate the Hsf1 transcription factor. Accordingly, the bag101-linked growth defect is alleviated in strains containing a reduced amount of Hsf1 but aggravated in hsp70 deletion strains. In conclusion, we propose that the fission yeast UBL/BAG proteins release Hsf1 from Hsp70, leading to constitutive Hsf1 activation and growth defects.

The Mutational Landscape of the Oncogenic MZF1 SCAN Domain in Cancer

SCAN domains in zinc-finger transcription factors are crucial mediators of protein-protein interactions. Up to 240 SCAN-domain encoding genes have been identified throughout the human genome. These include cancer-related genes, such as the myeloid zinc finger 1 (MZF1), an oncogenic transcription factor involved in the progression of many solid cancers. The mechanisms by which SCAN homo- and heterodimers assemble and how they alter the transcriptional activity of zinc-finger transcription factors in cancer and other diseases remain to be investigated. Here, we provide the first description of the conformational ensemble of the MZF1 SCAN domain cross-validated against NMR experimental data, which are probes of structure and dynamics on different timescales. We investigated the protein-protein interaction network of MZF1 and how it is perturbed in different cancer types by the analyses of high-throughput proteomics and RNASeq data. Collectively, we integrated many computational approaches, ranging from simple empirical energy functions to all-atom microsecond molecular dynamics simulations and network analyses to unravel the effects of cancer-related substitutions in relation to MZF1 structure and interactions.