PRKN-linked familial Parkinson's disease: cellular and molecular mechanisms of disease-linked variants.

Parkinson's disease (PD) is a common and incurable neurodegenerative disorder that arises from the loss of dopaminergic neurons in the substantia nigra and is mainly characterized by progressive loss of motor function. Monogenic familial PD is associated with highly penetrant variants in specific genes, notably the PRKN gene, where homozygous or compound heterozygous loss-of-function variants predominate. PRKN encodes Parkin, an E3 ubiquitin-protein ligase important for protein ubiquitination and mitophagy of damaged mitochondria. Accordingly, Parkin plays a central role in mitochondrial quality control but is itself also subject to a strict protein quality control system that rapidly eliminates certain disease-linked Parkin variants. Here, we summarize the cellular and molecular functions of Parkin, highlighting the various mechanisms by which PRKN gene variants result in loss-of-function. We emphasize the importance of high-throughput assays and computational tools for the clinical classification of PRKN gene variants and how detailed insights into the pathogenic mechanisms of PRKN gene variants may impact the development of personalized therapeutics.

Destabilization and Degradation of a Disease-Linked PGM1 Protein Variant.

PGM1-linked congenital disorder of glycosylation (PGM1-CDG) is an autosomal recessive disease characterized by several phenotypes, some of which are life-threatening. Research focusing on the disease-related variants of the α-D-phosphoglucomutase 1 (PGM1) protein has shown that several are insoluble in vitro and expressed at low levels in patient fibroblasts. Due to these observations, we hypothesized that some disease-linked PGM1 protein variants are structurally destabilized and subject to protein quality control (PQC) and rapid intracellular degradation. Employing yeast-based assays, we show that a disease-associated human variant, PGM1 L516P, is insoluble, inactive, and highly susceptible to ubiquitylation and rapid degradation by the proteasome. In addition, we show that PGM1 L516P forms aggregates in S. cerevisiae and that both the aggregation pattern and the abundance of PGM1 L516P are chaperone-dependent. Finally, using computational methods, we perform saturation mutagenesis to assess the impact of all possible single residue substitutions in the PGM1 protein. These analyses identify numerous missense variants with predicted detrimental effects on protein function and stability. We suggest that many disease-linked PGM1 variants are subject to PQC-linked degradation and that our in silico site-saturated data set may assist in the mechanistic interpretation of PGM1 variants.

Lysine deserts prevent adventitious ubiquitylation of ubiquitin-proteasome components.

In terms of its relative frequency, lysine is a common amino acid in the human proteome. However, by bioinformatics we find hundreds of proteins that contain long and evolutionarily conserved stretches completely devoid of lysine residues. These so-called lysine deserts show a high prevalence in intrinsically disordered proteins with known or predicted functions within the ubiquitin-proteasome system (UPS), including many E3 ubiquitin-protein ligases and UBL domain proteasome substrate shuttles, such as BAG6, RAD23A, UBQLN1 and UBQLN2. We show that introduction of lysine residues into the deserts leads to a striking increase in ubiquitylation of some of these proteins. In case of BAG6, we show that ubiquitylation is catalyzed by the E3 RNF126, while RAD23A is ubiquitylated by E6AP. Despite the elevated ubiquitylation, mutant RAD23A appears stable, but displays a partial loss of function phenotype in fission yeast. In case of UBQLN1 and BAG6, introducing lysine leads to a reduced abundance due to proteasomal degradation of the proteins. For UBQLN1 we show that arginine residues within the lysine depleted region are critical for its ability to form cytosolic speckles/inclusions. We propose that selective pressure to avoid lysine residues may be a common evolutionary mechanism to prevent unwarranted ubiquitylation and/or perhaps other lysine post-translational modifications. This may be particularly relevant for UPS components as they closely and frequently encounter the ubiquitylation machinery and are thus more susceptible to nonspecific ubiquitylation.

HSP70-binding motifs function as protein quality control degrons.

Protein quality control (PQC) degrons are short protein segments that target misfolded proteins for proteasomal degradation, and thus protect cells against the accumulation of potentially toxic non-native proteins. Studies have shown that PQC degrons are hydrophobic and rarely contain negatively charged residues, features which are shared with chaperone-binding regions. Here we explore the notion that chaperone-binding regions may function as PQC degrons. When directly tested, we found that a canonical Hsp70-binding motif (the APPY peptide) functioned as a dose-dependent PQC degron both in yeast and in human cells. In yeast, Hsp70, Hsp110, Fes1, and the E3 Ubr1 target the APPY degron. Screening revealed that the sequence space within the chaperone-binding region of APPY that is compatible with degron function is vast. We find that the number of exposed Hsp70-binding sites in the yeast proteome correlates with a reduced protein abundance and half-life. Our results suggest that when protein folding fails, chaperone-binding sites may operate as PQC degrons, and that the sequence properties leading to PQC-linked degradation therefore overlap with those of chaperone binding.

Mapping the degradation pathway of a disease-linked aspartoacylase variant.

Canavan disease is a severe progressive neurodegenerative disorder that is characterized by swelling and spongy degeneration of brain white matter. The disease is genetically linked to polymorphisms in the aspartoacylase (ASPA) gene, including the substitution C152W. ASPA C152W is associated with greatly reduced protein levels in cells, yet biophysical experiments suggest a wild-type like thermal stability. Here, we use ASPA C152W as a model to investigate the degradation pathway of a disease-causing protein variant. When we expressed ASPA C152W in Saccharomyces cerevisiae, we found a decreased steady state compared to wild-type ASPA as a result of increased proteasomal degradation. However, molecular dynamics simulations of ASPA C152W did not substantially deviate from wild-type ASPA, indicating that the native state is structurally preserved. Instead, we suggest that the C152W substitution interferes with the de novo folding pathway resulting in increased proteasomal degradation before reaching its stable conformation. Systematic mapping of the protein quality control components acting on misfolded and aggregation-prone species of C152W, revealed that the degradation is highly dependent on the molecular chaperone Hsp70, its co-chaperone Hsp110 as well as several quality control E3 ubiquitin-protein ligases, including Ubr1. In addition, the disaggregase Hsp104 facilitated refolding of aggregated ASPA C152W, while Cdc48 mediated degradation of insoluble ASPA protein. In human cells, ASPA C152W displayed increased proteasomal turnover that was similarly dependent on Hsp70 and Hsp110. Our findings underscore the use of yeast to determine the protein quality control components involved in the degradation of human pathogenic variants in order to identify potential therapeutic targets.

Biophysical and Mechanistic Models for Disease-Causing Protein Variants.

The rapid decrease in DNA sequencing cost is revolutionizing medicine and science. In medicine, genome sequencing has revealed millions of missense variants that change protein sequences, yet we only understand the molecular and phenotypic consequences of a small fraction. Within protein science, high-throughput deep mutational scanning experiments enable us to probe thousands of variants in a single, multiplexed experiment. We review efforts that bring together these topics via experimental and computational approaches to determine the consequences of missense variants in proteins. We focus on the role of changes in protein stability as a driver for disease, and how experiments, biophysical models, and computation are providing a framework for understanding and predicting how changes in protein sequence affect cellular protein stability.

Rare Catechol-O-methyltransferase Missense Variants Are Structurally Unstable Proteasome Targets.

Catechol-O-methyltransferase (COMT) is a key enzyme in the metabolism of catecholamines. Substrates of the enzyme include neurotransmitters such as dopamine and epinephrine, and therefore, COMT plays a central role in neurobiology. Since COMT also metabolizes catecholamine drugs such as L-DOPA, variation in COMT activity could affect pharmacokinetics and drug availability. Certain COMT missense variants have been shown to display decreased enzymatic activity. Additionally, studies have shown that such missense variants may lead to loss of function induced by impaired structural stability, which results in activation of the protein quality control system and degradation by the ubiquitin-proteasome system. Here, we demonstrate that two rare missense variants of COMT are ubiquitylated and targeted for proteasomal degradation as a result of structural destabilization and misfolding. This results in strongly reduced intracellular steady-state levels of the enzyme, which for the L135P variant is rescued upon binding to the COMT inhibitors entacapone and tolcapone. Our results reveal that the degradation is independent of the COMT isoform as both soluble (S-COMT) and ER membrane-bound (MB-COMT) variants are degraded. In silico structural stability predictions identify regions within the protein that are critical for stability overlapping with evolutionarily conserved residues, pointing toward other variants that are likely destabilized and degraded.

Lynch syndrome, molecular mechanisms and variant classification.

Patients with the heritable cancer disease, Lynch syndrome, carry germline variants in the MLH1, MSH2, MSH6 and PMS2 genes, encoding the central components of the DNA mismatch repair system. Loss-of-function variants disrupt the DNA mismatch repair system and give rise to a detrimental increase in the cellular mutational burden and cancer development. The treatment prospects for Lynch syndrome rely heavily on early diagnosis; however, accurate diagnosis is inextricably linked to correct clinical interpretation of individual variants. Protein variant classification traditionally relies on cumulative information from occurrence in patients, as well as experimental testing of the individual variants. The complexity of variant classification is due to (1) that variants of unknown significance are rare in the population and phenotypic information on the specific variants is missing, and (2) that individual variant testing is challenging, costly and slow. Here, we summarise recent developments in high-throughput technologies and computational prediction tools for the assessment of variants of unknown significance in Lynch syndrome. These approaches may vastly increase the number of interpretable variants and could also provide important mechanistic insights into the disease. These insights may in turn pave the road towards developing personalised treatment approaches for Lynch syndrome.

Mammalian 26S proteasomes remain intact during protein degradation.

It has been suggested that degradation of polyubiquitylated proteins is coupled to dissociation of 26S proteasomes. In contrast, using several independent types of experiments, we find that mammalian proteasomes can degrade polyubiquitylated proteins without disassembling. Thus, immobilized, (35)S-labeled 26S proteasomes degraded polyubiquitylated Sic1 and c-IAP1 without releasing any subunits. In addition, it is predicted that if 26S proteasomes dissociate into 20S proteasomes and regulatory complexes during a degradation cycle, the reassembly rate would be limiting at low proteasome concentrations. However, the rate with which each proteasome degraded polyubiquitylated Sic1 was independent of the proteasome concentration. Likewise, substrate-dependent dissociation of 26S proteasomes could not be detected by nondenaturing electrophoresis. Lastly, epoxomicin-inhibited 20S proteasomes can trap released regulatory complexes, forming inactive 26S proteasomes, but addition of epoxomicin-inhibited 20S proteasomes had no effect on the degradation of either polyubiquitylated Sic1 or UbcH10 by 26S proteasomes or of endogenous substrates in cell extracts.

A context-dependent and disordered ubiquitin-binding motif.

Ubiquitin is a small, globular protein that is conjugated to other proteins as a posttranslational event. A palette of small, folded domains recognizes and binds ubiquitin to translate and effectuate this posttranslational signal. Recent computational studies have suggested that protein regions can recognize ubiquitin via a process of folding upon binding. Using peptide binding arrays, bioinformatics, and NMR spectroscopy, we have uncovered a disordered ubiquitin-binding motif that likely remains disordered when bound and thus expands the palette of ubiquitin-binding proteins. We term this motif Disordered Ubiquitin-Binding Motif (DisUBM) and find it to be present in many proteins with known or predicted functions in degradation and transcription. We decompose the determinants of the motif showing it to rely on features of aromatic and negatively charged residues, and less so on distinct sequence positions in line with its disordered nature. We show that the affinity of the motif is low and moldable by the surrounding disordered chain, allowing for an enhanced interaction surface with ubiquitin, whereby the affinity increases ~ tenfold. Further affinity optimization using peptide arrays pushed the affinity into the low micromolar range, but compromised context dependence. Finally, we find that DisUBMs can emerge from unbiased screening of randomized peptide libraries, featuring in de novo cyclic peptides selected to bind ubiquitin chains. We suggest that naturally occurring DisUBMs can recognize ubiquitin as a posttranslational signal to act as affinity enhancers in IDPs that bind to folded and ubiquitylated binding partners.

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.

Understanding the Origins of Loss of Protein Function by Analyzing the Effects of Thousands of Variants on Activity and Abundance.

Understanding and predicting how amino acid substitutions affect proteins are keys to our basic understanding of protein function and evolution. Amino acid changes may affect protein function in a number of ways including direct perturbations of activity or indirect effects on protein folding and stability. We have analyzed 6,749 experimentally determined variant effects from multiplexed assays on abundance and activity in two proteins (NUDT15 and PTEN) to quantify these effects and find that a third of the variants cause loss of function, and about half of loss-of-function variants also have low cellular abundance. We analyze the structural and mechanistic origins of loss of function and use the experimental data to find residues important for enzymatic activity. We performed computational analyses of protein stability and evolutionary conservation and show how we may predict positions where variants cause loss of activity or abundance. In this way, our results link thermodynamic stability and evolutionary conservation to experimental studies of different properties of protein fitness landscapes.

Dss1 is a 26S proteasome ubiquitin receptor.

The ubiquitin-proteasome system is the major pathway for protein degradation in eukaryotic cells. Proteins to be degraded are conjugated to ubiquitin chains that act as recognition signals for the 26S proteasome. The proteasome subunits Rpn10 and Rpn13 are known to bind ubiquitin, but genetic and biochemical data suggest the existence of at least one other substrate receptor. Here, we show that the phylogenetically conserved proteasome subunit Dss1 (Sem1) binds ubiquitin chains linked by K63 and K48. Atomic resolution data show that Dss1 is disordered and binds ubiquitin by binding sites characterized by acidic and hydrophobic residues. The complementary binding region in ubiquitin is composed of a hydrophobic patch formed by I13, I44, and L69 flanked by two basic regions. Mutations in the ubiquitin-binding site of Dss1 cause growth defects and accumulation of ubiquitylated proteins.

DSS1/Sem1, a Multifunctional and Intrinsically Disordered Protein.

DSS1/Sem1 is a versatile intrinsically disordered protein. Besides being a bona fide subunit of the 26S proteasome, DSS1 associates with other protein complexes, including BRCA2-RPA, involved in homologous recombination; the Csn12-Thp3 complex, involved in RNA splicing; the integrator, involved in transcription; and the TREX-2 complex, involved in nuclear export of mRNA and transcription elongation. As a subunit of the proteasome, DSS1 functions both in complex assembly and possibly as a ubiquitin receptor. Here, we summarise structural and functional aspects of DSS1/Sem1 with particular emphasis on its multifunctional and disordered properties. We suggest that DSS1/Sem1 can act as a polyanionic adhesive to prevent nonproductive interactions during construction of protein assemblies, uniquely employing different structures when associating with the diverse multisubunit complexes.

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.

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 (BH4 ), 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.

Random Mutagenesis Analysis of the Influenza A M2 Proton Channel Reveals Novel Resistance Mutants.

The influenza M2 proton channel is a major drug target, but unfortunately, the acquisition of resistance mutations greatly reduces the functional life span of a drug in influenza treatment. New M2 inhibitors that inhibit mutant M2 channels otherwise resistant to the early adamantine-based drugs have been reported, but it remains unclear whether and how easy resistance could arise to such inhibitors. We have combined a newly developed proton conduction assay with an established method for selection and screening, both Escherichia coli-based, to enable the study of M2 function and inhibition. Combining this platform with two groups of structurally different M2 inhibitors allowed us to isolate drug resistant M2 channels from a mutant library. Two groups of M2 variants emerged from this analysis. A first group appeared almost unaffected by the inhibitor, M_089 (N13I, I35L, and F47L) and M_272 (G16C and D44H), and the single-substitution variants derived from these (I35L, L43P, D44H, and L46P). Functionally, these resemble the known drug resistant M2 channels V27A, S31N, and swine flu. In addition, a second group of tested M2 variants were all still inhibited by drugs but to a lesser extent than wild type M2. Molecular dynamics simulations aided in distinguishing the two groups where drug binding to the wild type and the less resistant M2 group showed a stable positioning of the ligand in the canonical binding pose, as opposed to the drug resistant group in which the ligand rapidly dissociated from the complex during the simulations.

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.

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.