It is generally held that random-coil polypeptide chains undergo a barrier-less continuous collapse when the solvent conditions are changed to favor the fully folded native conformation. We test this hypothesis by probing intramolecular distance distributions during folding in one of the paradigms of folding reactions, that of cytochrome c. The Trp59-to-heme distance was probed by time-resolved Förster resonance energy transfer in the microsecond time range of refolding. Contrary to expectation, a state with a Trp59-heme distance close to that of the guanidinium hydrochloride (GdnHCl) denatured state is present after ~27 μs of folding. A concomitant decrease in the population of this state and an increase in the population of a compact high-FRET (Förster resonance energy transfer) state (efficiency>90%) show that the collapse is barrier limited. Small-angle X-ray scattering (SAXS) measurements over a similar time range show that the radius of gyration under native favoring conditions is comparable to that of the GdnHCl denatured unfolded state. An independent comprehensive global thermodynamic analysis reveals that marginally stable partially folded structures are also present in the nominally unfolded GdnHCl denatured state. These observations suggest that specifically collapsed intermediate structures with low stability in rapid equilibrium with the unfolded state may contribute to the apparent chain contraction observed in previous fluorescence studies using steady-state detection. In the absence of significant dynamic averaging of marginally stable partially folded states and with the use of probes sensitive to distance distributions, barrier-limited chain contraction is observed upon transfer of the GdnHCl denatured state ensemble to native-like conditions.

Accumulation of misfolded proteins in the lumen of the endoplasmic reticulum (ER) activates the unfolded protein response (UPR). Ire1, an ER-resident transmembrane kinase/RNase, senses the protein folding status inside the ER. When activated, Ire1 oligomerizes and trans-autophosphorylates, activating its RNase and initiating a nonconventional mRNA splicing reaction. Splicing results in production of the transcription factor Hac1 that induces UPR target genes; expression of these genes restores ER homeostasis by increasing its protein folding capacity and allows abatement of UPR signaling. Here, we uncouple Ire1's RNase from its kinase activity and find that cells expressing kinase-inactive Ire1 can regulate Ire1's RNase, splice HAC1 mRNA, produce Hac1 protein, and induce UPR target genes. Unlike wild-type IRE1, kinase-inactive Ire1 cells display defects in Ire1 deactivation. Failure to properly inactivate Ire1 causes chronic ER stress and reduces cell survival under UPR-inducing conditions. Thus, Ire1-catalyzed phosphoryl-transfer aids disassembly of Ire1 signaling complexes and is a critical component of the UPR homeostatic feedback loop.

The functional diversity of kinases enables specificity in cellular signal transduction. Yet how more than 500 members of the human kinome specifically receive regulatory inputs and convey information to appropriate substrates-all while using the common signaling output of phosphorylation-remains enigmatic. Here, we perform statistical co-evolution analysis, mutational scanning, and quantitative live-cell assays to reveal a hierarchical organization of the kinase domain that facilitates the orthogonal evolution of regulatory inputs and substrate outputs while maintaining catalytic function. We find that three quasi-independent "sectors"-groups of evolutionarily coupled residues-represent functional units in the kinase domain that encode for catalytic activity, substrate specificity, and regulation. Sector positions impact both disease and pharmacology: the catalytic sector is significantly enriched for somatic cancer mutations, and residues in the regulatory sector interact with allosteric kinase inhibitors. We propose that this functional architecture endows the kinase domain with inherent regulatory plasticity.

Genetic, epidemiological and experimental evidence implicate lysosomal dysfunction in Parkinson's disease (PD) and related synucleinopathies. Investigate several mouse models of lysosomal storage diseases (LSDs) and evaluate pathologies reminiscent of synucleinopathies. We obtained brain tissue from symptomatic mouse models of Gaucher, Fabry, Sandhoff, Niemann-Pick A (NPA), Hurler, Pompe and Niemann-Pick C (NPC) diseases and assessed for the presence of Lewy body-like pathology (proteinase K-resistant α-synuclein and tau aggregates) and neuroinflammation (microglial Iba1 and astrocytic GFAP) by immunofluorescence. All seven LSD models exhibited evidence of proteinopathy and/or inflammation in the central nervous system (CNS). However, these phenotypes were divergent. Gaucher and Fabry mouse models displayed proteinase K-resistant α-synuclein and tau aggregates but no neuroinflammation; whereas Sandhoff, NPA and NPC showed marked neuroinflammation and no overt proteinopathy. Pompe disease animals uniquely displayed widespread distribution of tau aggregates accompanied by moderate microglial activation. Hurler mice also demonstrated proteinopathy and microglial activation. The present study demonstrated additional links between LSDs and pathogenic phenotypes that are hallmarks of synucleinopathies. The data suggest that lysosomal dysregulation can contribute to brain region-specific protein aggregation and induce widespread neuroinflammation in the brain. However, only a few LSD models examined exhibited phenotypes consistent with synucleinopathies. While no model can recapitulate the complexity of PD, they can enable the study of specific pathways and mechanisms contributing to disease pathophysiology. The present study provides evidence that there are existing, previously unutilized mouse models that can be employed to study pathogenic mechanisms and gain insights into potential PD subtypes, helping to determine if they are amenable to pathway-specific therapeutic interventions.

The heat shock response (HSR) controls expression of molecular chaperones to maintain protein homeostasis. Previously, we proposed a feedback loop model of the HSR in which heat-denatured proteins sequester the chaperone Hsp70 to activate the HSR, and subsequent induction of Hsp70 deactivates the HSR (Krakowiak et al., 2018; Zheng et al., 2016). However, recent work has implicated newly synthesized proteins (NSPs) - rather than unfolded mature proteins - and the Hsp70 co-chaperone Sis1 in HSR regulation, yet their contributions to HSR dynamics have not been determined. Here, we generate a new mathematical model that incorporates NSPs and Sis1 into the HSR activation mechanism, and we perform genetic decoupling and pulse-labeling experiments to demonstrate that Sis1 induction is dispensable for HSR deactivation. Rather than providing negative feedback to the HSR, transcriptional regulation of Sis1 by Hsf1 promotes fitness by coordinating stress granules and carbon metabolism. These results support an overall model in which NSPs signal the HSR by sequestering Sis1 and Hsp70, while induction of Hsp70 - but not Sis1 - attenuates the response.

Alpha-synuclein (αSyn) protein levels correlate with the risk and severity of Parkinson's disease and related neurodegenerative diseases. Lowering αSyn is being actively investigated as a therapeutic modality. Here, we systematically map the regulatory network that controls endogenous αSyn using sequential CRISPR-knockout and -interference screens in an αSyn gene (SNCA)-tagged cell line and induced pluripotent stem cell-derived neurons (iNeurons). We uncover αSyn modifiers at multiple regulatory layers, with amino-terminal acetyltransferase B (NatB) enzymes being the most potent endogenous αSyn modifiers in both cell lines. Amino-terminal acetylation protects the cytosolic αSyn from rapid degradation by the proteasome in a Ube2w-dependent manner. Moreover, we show that pharmacological inhibition of methionyl-aminopeptidase 2, a regulator of NatB complex formation, attenuates endogenous αSyn in iNeurons carrying SNCA triplication. Together, our study reveals several gene networks that control endogenous αSyn, identifies mechanisms mediating the degradation of nonacetylated αSyn, and illustrates potential therapeutic pathways for decreasing αSyn levels in synucleinopathies.

Despite the vast number of modification sites mapped within mRNAs, known examples of consequential mRNA modifications remain rare. Here, we provide multiple lines of evidence to show that Ime4p, an N6-methyladenosine (m6A) methyltransferase required for meiosis in yeast, acts by methylating a site in the 3' UTR of the mRNA encoding Rme1p, a transcriptional repressor of meiosis. Consistent with this mechanism, genetic analyses reveal that IME4 functions upstream of RME1. Transcriptome-wide, RME1 is the primary message that displays both increased methylation and reduced expression in an Ime4p-dependent manner. In yeast strains for which IME4 is dispensable for meiosis, a natural polymorphism in the RME1 promoter reduces RME1 transcription, obviating the requirement for methylation. Mutation of a single m6A site in the RME1 3' UTR increases Rme1p repressor production and reduces meiotic efficiency. These results reveal the molecular and physiological consequences of a modification in the 3' UTR of an mRNA.

Heat shock factor 1 is the master transcriptional regulator of molecular chaperones and binds to the same cis-acting heat shock element (HSE) across the eukaryotic lineage. In budding yeast, Hsf1 drives the transcription of ∼20 genes essential to maintain proteostasis under basal conditions, yet its specific targets and extent of inducible binding during heat shock remain unclear. Here we combine Hsf1 chromatin immunoprecipitation sequencing (seq), nascent RNA-seq, and Hsf1 nuclear depletion to quantify Hsf1 binding and transcription across the yeast genome. We find that Hsf1 binds 74 loci during acute heat shock, and these are linked to 46 genes with strong Hsf1-dependent expression. Notably, Hsf1's induced DNA binding leads to a disproportionate (∼7.5-fold) increase in nascent transcription. Promoters with high basal Hsf1 occupancy have nucleosome-depleted regions due to the presence of "pioneer factors." These accessible sites are likely critical for Hsf1 occupancy as the activator is incapable of binding HSEs within a stably positioned, reconstituted nucleosome. In response to heat shock, however, Hsf1 accesses nucleosomal sites and promotes chromatin disassembly in concert with the Remodels Structure of Chromatin (RSC) complex. Our data suggest that the interplay between nucleosome positioning, HSE strength, and active Hsf1 levels allows cells to precisely tune expression of the proteostasis network.

Clonal populations of cells exhibit cell-to-cell variation in the transcription of individual genes. In addition to this noise in gene expression, heterogeneity in the proteome and the proteostasis network expands the phenotypic diversity of a population. Heat shock factor 1 (Hsf1) regulates chaperone gene expression, thereby coupling transcriptional noise to proteostasis. Here we show that cell-to-cell variation in Hsf1 activity is an important determinant of phenotypic plasticity. Budding yeast cells with high Hsf1 activity were enriched for the ability to acquire resistance to an antifungal drug, and this enrichment depended on Hsp90, a known phenotypic capacitor and canonical Hsf1 target. We show that Hsf1 phosphorylation promotes cell-to-cell variation, and this variation, rather than absolute Hsf1 activity, promotes antifungal resistance. We propose that Hsf1 phosphorylation enables differential tuning of the proteostasis network in individual cells, allowing populations to access a range of phenotypic states.

Models for regulation of the eukaryotic heat shock response typically invoke a negative feedback loop consisting of the transcriptional activator Hsf1 and a molecular chaperone. Previously we identified Hsp70 as the chaperone responsible for Hsf1 repression and constructed a mathematical model that recapitulated the yeast heat shock response (Zheng et al., 2016). The model was based on two assumptions: dissociation of Hsp70 activates Hsf1, and transcriptional induction of Hsp70 deactivates Hsf1. Here we validate these assumptions. First, we severed the feedback loop by uncoupling Hsp70 expression from Hsf1 regulation. As predicted by the model, Hsf1 was unable to efficiently deactivate in the absence of Hsp70 transcriptional induction. Next, we mapped a discrete Hsp70 binding site on Hsf1 to a C-terminal segment known as conserved element 2 (CE2). In vitro, CE2 binds to Hsp70 with low affinity (9 µM), in agreement with model requirements. In cells, removal of CE2 resulted in increased basal Hsf1 activity and delayed deactivation during heat shock, while tandem repeats of CE2 sped up Hsf1 deactivation. Finally, we uncovered a role for the N-terminal domain of Hsf1 in negatively regulating DNA binding. These results reveal the quantitative control mechanisms underlying the heat shock response.

The unfolded protein response (UPR) is an intracellular signaling pathway that counteracts variable stresses that impair protein folding in the endoplasmic reticulum (ER). As such, the UPR is thought to be a homeostat that finely tunes ER protein folding capacity and ER abundance according to need. The mechanism by which the ER stress sensor Ire1 is activated by unfolded proteins and the role that the ER chaperone protein BiP plays in Ire1 regulation have remained unclear. Here we show that the UPR matches its output to the magnitude of the stress by regulating the duration of Ire1 signaling. BiP binding to Ire1 serves to desensitize Ire1 to low levels of stress and promotes its deactivation when favorable folding conditions are restored to the ER. We propose that, mechanistically, BiP achieves these functions by sequestering inactive Ire1 molecules, thereby providing a barrier to oligomerization and activation, and a stabilizing interaction that facilitates de-oligomerization and deactivation. Thus BiP binding to or release from Ire1 is not instrumental for switching the UPR on and off as previously posed. By contrast, BiP provides a buffer for inactive Ire1 molecules that ensures an appropriate response to restore protein folding homeostasis to the ER by modulating the sensitivity and dynamics of Ire1 activity.

Alpha-synuclein (αS) is a conformationally plastic protein that reversibly binds to cellular membranes. It aggregates and is genetically linked to Parkinson's disease (PD). Here, we show that αS directly modulates processing bodies (P-bodies), membraneless organelles that function in mRNA turnover and storage. The N terminus of αS, but not other synucleins, dictates mutually exclusive binding either to cellular membranes or to P-bodies in the cytosol. αS associates with multiple decapping proteins in close proximity on the Edc4 scaffold. As αS pathologically accumulates, aberrant interaction with Edc4 occurs at the expense of physiologic decapping-module interactions. mRNA decay kinetics within PD-relevant pathways are correspondingly disrupted in PD patient neurons and brain. Genetic modulation of P-body components alters αS toxicity, and human genetic analysis lends support to the disease-relevance of these interactions. Beyond revealing an unexpected aspect of αS function and pathology, our data highlight the versatility of conformationally plastic proteins with high intrinsic disorder.

Transcriptional induction of heat shock protein (HSP) genes is accompanied by dynamic changes in their 3D structure and spatial organization, yet the molecular basis for these phenomena remains unknown. Using chromosome conformation capture and single-cell imaging, we show that genes transcriptionally activated by Hsf1 specifically interact across chromosomes and coalesce into diffraction-limited intranuclear foci. Genes activated by the alternative stress regulators Msn2/Msn4, in contrast, do not interact among themselves nor with Hsf1 targets. Likewise, constitutively expressed genes, even those interposed between HSP genes, show no detectable interaction. Hsf1 forms discrete subnuclear puncta when stress activated, and these puncta dissolve in concert with transcriptional attenuation, paralleling the kinetics of HSP gene coalescence and dissolution. Nuclear Hsf1 and RNA Pol II are both necessary for intergenic HSP gene interactions, while DNA-bound Hsf1 is necessary and sufficient to drive heterologous gene coalescence. Our findings demonstrate that Hsf1 can dynamically restructure the yeast genome.

In Parkinson's disease (PD), α-synuclein (αS) pathologically impacts the brain, a highly lipid-rich organ. We investigated how alterations in αS or lipid/fatty acid homeostasis affect each other. Lipidomic profiling of human αS-expressing yeast revealed increases in oleic acid (OA, 18:1), diglycerides, and triglycerides. These findings were recapitulated in rodent and human neuronal models of αS dyshomeostasis (overexpression; patient-derived triplication or E46K mutation; E46K mice). Preventing lipid droplet formation or augmenting OA increased αS yeast toxicity; suppressing the OA-generating enzyme stearoyl-CoA-desaturase (SCD) was protective. Genetic or pharmacological SCD inhibition ameliorated toxicity in αS-overexpressing rat neurons. In a C. elegans model, SCD knockout prevented αS-induced dopaminergic degeneration. Conversely, we observed detrimental effects of OA on αS homeostasis: in human neural cells, excess OA caused αS inclusion formation, which was reversed by SCD inhibition. Thus, monounsaturated fatty acid metabolism is pivotal for αS-induced neurotoxicity, and inhibiting SCD represents a novel PD therapeutic approach.

Although tRNA modifications have been well catalogued, the precise functions of many modifications and their roles in mediating gene expression are still being elucidated. Whereas tRNA modifications were long assumed to be constitutive, it is now apparent that the modification status of tRNAs changes in response to different environmental conditions. The URM1 pathway is required for thiolation of the cytoplasmic tRNAs tGlu(UUC), tGln(UUG), and tLys(UUU) in Saccharomyces cerevisiae. We demonstrate that URM1 pathway mutants have impaired translation, which results in increased basal activation of the Hsf1-mediated heat shock response; we also find that tRNA thiolation levels in wild-type cells decrease when cells are grown at elevated temperature. We show that defects in tRNA thiolation can be conditionally advantageous, conferring resistance to endoplasmic reticulum stress. URM1 pathway proteins are unstable and hence are more sensitive to changes in the translational capacity of cells, which is decreased in cells experiencing stresses. We propose a model in which a stress-induced decrease in translation results in decreased levels of URM1 pathway components, which results in decreased tRNA thiolation levels, which further serves to decrease translation. This mechanism ensures that tRNA thiolation and translation are tightly coupled and coregulated according to need.

Cell signaling systems transmit information by post-translationally modifying signaling proteins, often via phosphorylation. While thousands of sites of phosphorylation have been identified in proteomic studies, the vast majority of sites have no known function. Assigning functional roles to the catalog of uncharacterized phosphorylation sites is a key research challenge. Here we present a general approach to address this challenge and apply it to a prototypical signaling pathway, the pheromone response pathway in Saccharomyces cerevisiae. The pheromone pathway includes a mitogen activated protein kinase (MAPK) cascade activated by a G-protein coupled receptor (GPCR). We used published mass spectrometry-based proteomics data to identify putative sites of phosphorylation on pheromone pathway components, and we used evolutionary conservation to assign priority to a list of candidate MAPK regulatory sites. We made targeted alterations in those sites, and measured the effects of the mutations on pheromone pathway output in single cells. Our work identified six new sites that quantitatively tuned system output. We developed simple computational models to find system architectures that recapitulated the quantitative phenotypes of the mutants. Our results identify a number of putative phosphorylation events that contribute to adjust the input-output relationship of this model eukaryotic signaling system. We believe this combined approach constitutes a general means not only to reveal modification sites required to turn a pathway on and off, but also those required for more subtle quantitative effects that tune pathway output. Our results suggest that relatively small quantitative influences from individual phosphorylation events endow signaling systems with plasticity that evolution may exploit to quantitatively tailor signaling outcomes.

The assembly of signaling components and transcription factors in ordered subcellular structures is increasingly implicated as an important regulatory strategy for modulating the activity of cellular pathways. Here, we document the inducible formation of subnuclear foci formed by two mitogen-activated protein kinases (MAPKs) in Saccharomyces cerevisiae upon hyperosmotic stress. Specifically, we demonstrate that activation of the hyperosmotic stress response pathway induces the mating pathway MAPK Fus3 and the filamentation pathway MAPK Kss1 to form foci in the nucleus that are organized by their shared downstream transcription factor Ste12. Foci formation of colocalized Ste12, Fus3, and Kss1 requires the kinase activity of the hyperosmotic response MAPK Hog1 and correlates with attenuated signaling in the mating pathway. Conversely, activation of the mating pathway prevents foci formation upon subsequent hyperosmotic stress. These results suggest that Hog1-mediated spatial localization of Fus3 and Ste12 into subnuclear foci could contribute to uncoupling the pheromone and osmolarity pathways, which share signaling components, under high-osmolarity conditions.

Amino acid replacements at dozens of positions in the dimeric protein human, Cu,Zn superoxide dismutase (SOD1) can cause amyotrophic lateral sclerosis (ALS). Although it has long been hypothesized that these mutations might enhance the populations of marginally-stable aggregation-prone species responsible for cellular toxicity, there has been little quantitative evidence to support this notion. Perturbations of the folding free energy landscapes of metal-free versions of five ALS-inducing variants, A4V, L38V, G93A, L106V and S134N SOD1, were determined with a global analysis of kinetic and thermodynamic folding data for dimeric and stable monomeric versions of these variants. Utilizing this global analysis approach, the perturbations on the global stability in response to mutation can be partitioned between the monomer folding and association steps, and the effects of mutation on the populations of the folded and unfolded monomeric states can be determined. The 2- to 10-fold increase in the population of the folded monomeric state for A4V, L38V and L106V and the 80- to 480-fold increase in the population of the unfolded monomeric states for all but S134N would dramatically increase their propensity for aggregation through high-order nucleation reactions. The wild-type-like populations of these states for the metal-binding region S134N variant suggest that even wild-type SOD1 may also be prone to aggregation in the absence of metals.

Protein aggregation is a hallmark of many diseases but also underlies a wide range of positive cellular functions. This phenomenon has been difficult to study because of a lack of quantitative and high-throughput cellular tools. Here, we develop a synthetic genetic tool to sense and control protein aggregation. We apply the technology to yeast prions, developing sensors to track their aggregation states and employing prion fusions to encode synthetic memories in yeast cells. Utilizing high-throughput screens, we identify prion-curing mutants and engineer "anti-prion drives" that reverse the non-Mendelian inheritance pattern of prions and eliminate them from yeast populations. We extend our technology to yeast RNA-binding proteins (RBPs) by tracking their propensity to aggregate, searching for co-occurring aggregates, and uncovering a group of coalescing RBPs through screens enabled by our platform. Our work establishes a quantitative, high-throughput, and generalizable technology to study and control diverse protein aggregation processes in cells.

Mutations in GBA, the gene encoding the lysosomal enzyme glucocerebrosidase (GCase), represent the greatest genetic risk factor for developing synucleinopathies including Parkinson's disease (PD). Additionally, PD patients harboring a mutant GBA allele present with an earlier disease onset and an accelerated disease progression of both motor and non-motor symptoms. Preclinical studies in mouse models of synucleinopathy suggest that modulation of the sphingolipid metabolism pathway via inhibition of glucosylceramide synthase (GCS) using a CNS-penetrant small molecule may be a potential treatment for synucleinopathies. Here, we aim to alleviate the lipid storage burden by inhibiting the de novo synthesis of the primary glycosphingolipid substrate of GCase, glucosylceramide (GlcCer). We have previously shown that systemic GCS inhibition reduced GlcCer and glucosylsphingosine (GlcSph) accumulation, slowed α-synuclein buildup in the hippocampus, and improved cognitive deficits. Here, we studied the efficacy of a brain-penetrant clinical candidate GCS inhibitor, venglustat, in mouse models of GBA-related synucleinopathy, including a heterozygous Gba mouse model which more closely replicates the typical GBA-PD patient genotype. Collectively, these data support the rationale for modulation of GCase-related sphingolipid metabolism as a therapeutic strategy for treating GBA-related synucleinopathies.