A functionally divergent intrinsically disordered region underlying the conservation of stochastic signaling.

Stochastic signaling dynamics expand living cells' information processing capabilities. An increasing number of studies report that regulators encode information in their pulsatile dynamics. The evolutionary mechanisms that lead to complex signaling dynamics remain uncharacterized, perhaps because key interactions of signaling proteins are encoded in intrinsically disordered regions (IDRs), whose evolution is difficult to analyze. Here we focused on the IDR that controls the stochastic pulsing dynamics of Crz1, a transcription factor in fungi downstream of the widely conserved calcium signaling pathway. We find that Crz1 IDRs from anciently diverged fungi can all respond transiently to calcium stress; however, only Crz1 IDRs from the Saccharomyces clade support pulsatility, encode extra information, and rescue fitness in competition assays, while the Crz1 IDRs from distantly related fungi do none of the three. On the other hand, we find that Crz1 pulsing is conserved in the distantly related fungi, consistent with the evolutionary model of stabilizing selection on the signaling phenotype. Further, we show that a calcineurin docking site in a specific part of the IDRs appears to be sufficient for pulsing and show evidence for a beneficial increase in the relative calcineurin affinity of this docking site. We propose that evolutionary flexibility of functionally divergent IDRs underlies the conservation of stochastic signaling by stabilizing selection.

Identifying molecular features that are associated with biological function of intrinsically disordered protein regions.

In previous work, we showed that intrinsically disordered regions (IDRs) of proteins contain sequence-distributed molecular features that are conserved over evolution, despite little sequence similarity that can be detected in alignments (Zarin et al., 2019). Here, we aim to use these molecular features to predict specific biological functions for individual IDRs and identify the molecular features within them that are associated with these functions. We find that the predictable functions are diverse. Examining the associated molecular features, we note some that are consistent with previous reports and identify others that were previously unknown. We experimentally confirm that elevated isoelectric point and hydrophobicity, features that are positively associated with mitochondrial localization, are necessary for mitochondrial targeting function. Remarkably, increasing isoelectric point in a synthetic IDR restores weak mitochondrial targeting. We believe feature analysis represents a new systematic approach to understand how biological functions of IDRs are specified by their protein sequences.

Proteome-wide signatures of function in highly diverged intrinsically disordered regions.

Intrinsically disordered regions make up a large part of the proteome, but the sequence-to-function relationship in these regions is poorly understood, in part because the primary amino acid sequences of these regions are poorly conserved in alignments. Here we use an evolutionary approach to detect molecular features that are preserved in the amino acid sequences of orthologous intrinsically disordered regions. We find that most disordered regions contain multiple molecular features that are preserved, and we define these as 'evolutionary signatures' of disordered regions. We demonstrate that intrinsically disordered regions with similar evolutionary signatures can rescue function in vivo, and that groups of intrinsically disordered regions with similar evolutionary signatures are strongly enriched for functional annotations and phenotypes. We propose that evolutionary signatures can be used to predict function for many disordered regions from their amino acid sequences.

A Noisy Analog-to-Digital Converter Connects Cytosolic Calcium Bursts to Transcription Factor Nuclear Localization Pulses in Yeast.

Several examples of transcription factors that show stochastic, unsynchronized pulses of nuclear localization have been described. Here we show that under constant calcium stress, nuclear localization pulses of the transcription factor Crz1 follow stochastic variations in cytosolic calcium concentration. We find that the size of the stochastic calcium bursts is positively correlated with the number of subsequent Crz1 pulses. Based on our observations, we propose a simple stochastic model of how the signaling pathway converts a constant external calcium concentration into a digital number of Crz1 pulses in the nucleus, due to the time delay from nuclear transport and the stochastic decoherence of individual Crz1 molecule dynamics. We find support for several additional predictions of the model and suggest that stochastic input to nuclear transport may produce noisy digital responses to analog signals in other signaling systems.

Short linear motifs in intrinsically disordered regions modulate HOG signaling capacity.

BACKGROUND: The effort to characterize intrinsically disordered regions of signaling proteins is rapidly expanding. An important class of disordered interaction modules are ubiquitous and functionally diverse elements known as short linear motifs (SLiMs). RESULTS: To further examine the role of SLiMs in signal transduction, we used a previously devised bioinformatics method to predict evolutionarily conserved SLiMs within a well-characterized pathway in S. cerevisiae. Using a single cell, reporter-based flow cytometry assay in conjunction with a fluorescent reporter driven by a pathway-specific promoter, we quantitatively assessed pathway output via systematic deletions of individual motifs. We found that, when deleted, 34% (10/29) of predicted SLiMs displayed a significant decrease in pathway output, providing evidence that these motifs play a role in signal transduction. Assuming that mutations in SLiMs have quantitative effects on mechanisms of signaling, we show that perturbations of parameters in a previously published stochastic model of HOG signaling could reproduce the quantitative effects of 4 out of 7 mutations in previously unknown SLiMs. CONCLUSIONS: Our study suggests that, even in well-characterized pathways, large numbers of functional elements remain undiscovered, and that challenges remain for application of systems biology models to interpret the effects of mutations in signaling pathways.

Integrating images from multiple microscopy screens reveals diverse patterns of change in the subcellular localization of proteins.

The evaluation of protein localization changes on a systematic level is a powerful tool for understanding how cells respond to environmental, chemical, or genetic perturbations. To date, work in understanding these proteomic responses through high-throughput imaging has catalogued localization changes independently for each perturbation. To distinguish changes that are targeted responses to the specific perturbation or more generalized programs, we developed a scalable approach to visualize the localization behavior of proteins across multiple experiments as a quantitative pattern. By applying this approach to 24 experimental screens consisting of nearly 400,000 images, we differentiated specific responses from more generalized ones, discovered nuance in the localization behavior of stress-responsive proteins, and formed hypotheses by clustering proteins that have similar patterns. Previous approaches aim to capture all localization changes for a single screen as accurately as possible, whereas our work aims to integrate large amounts of imaging data to find unexpected new cell biology.

Parallel reorganization of protein function in the spindle checkpoint pathway through evolutionary paths in the fitness landscape that appear neutral in laboratory experiments.

Regulatory networks often increase in complexity during evolution through gene duplication and divergence of component proteins. Two models that explain this increase in complexity are: 1) adaptive changes after gene duplication, such as resolution of adaptive conflicts, and 2) non-adaptive processes such as duplication, degeneration and complementation. Both of these models predict complementary changes in the retained duplicates, but they can be distinguished by direct fitness measurements in organisms with short generation times. Previously, it has been observed that repeated duplication of an essential protein in the spindle checkpoint pathway has occurred multiple times over the eukaryotic tree of life, leading to convergent protein domain organization in its duplicates. Here, we replace the paralog pair in S. cerevisiae with a single-copy protein from a species that did not undergo gene duplication. Surprisingly, using quantitative fitness measurements in laboratory conditions stressful for the spindle-checkpoint pathway, we find no evidence that reorganization of protein function after gene duplication is beneficial. We then reconstruct several evolutionary intermediates from the inferred ancestral network to the extant one, and find that, at the resolution of our assay, there exist stepwise mutational paths from the single protein to the divergent pair of extant proteins with no apparent fitness defects. Parallel evolution has been taken as strong evidence for natural selection, but our results suggest that even in these cases, reorganization of protein function after gene duplication may be explained by neutral processes.

Local statistics allow quantification of cell-to-cell variability from high-throughput microscope images.

MOTIVATION: Quantifying variability in protein expression is a major goal of systems biology and cell-to-cell variability in subcellular localization pattern has not been systematically quantified. RESULTS: We define a local measure to quantify cell-to-cell variability in high-throughput microscope images and show that it allows comparable measures of variability for proteins with diverse subcellular localizations. We systematically estimate cell-to-cell variability in the yeast GFP collection and identify examples of proteins that show cell-to-cell variability in their subcellular localization. CONCLUSIONS: Automated image analysis methods can be used to quantify cell-to-cell variability in microscope images.

Detecting functional divergence after gene duplication through evolutionary changes in posttranslational regulatory sequences.

Gene duplication is an important evolutionary mechanism that can result in functional divergence in paralogs due to neo-functionalization or sub-functionalization. Consistent with functional divergence after gene duplication, recent studies have shown accelerated evolution in retained paralogs. However, little is known in general about the impact of this accelerated evolution on the molecular functions of retained paralogs. For example, do new functions typically involve changes in enzymatic activities, or changes in protein regulation? Here we study the evolution of posttranslational regulation by examining the evolution of important regulatory sequences (short linear motifs) in retained duplicates created by the whole-genome duplication in budding yeast. To do so, we identified short linear motifs whose evolutionary constraint has relaxed after gene duplication with a likelihood-ratio test that can account for heterogeneity in the evolutionary process by using a non-central chi-squared null distribution. We find that short linear motifs are more likely to show changes in evolutionary constraints in retained duplicates compared to single-copy genes. We examine changes in constraints on known regulatory sequences and show that for the Rck1/Rck2, Fkh1/Fkh2, Ace2/Swi5 paralogs, they are associated with previously characterized differences in posttranslational regulation. Finally, we experimentally confirm our prediction that for the Ace2/Swi5 paralogs, Cbk1 regulated localization was lost along the lineage leading to SWI5 after gene duplication. Our analysis suggests that changes in posttranslational regulation mediated by short regulatory motifs systematically contribute to functional divergence after gene duplication.

Endoproteolytic processing of the mammalian prion glycoprotein family.

Cellular prion protein (PrP(C)) misfolds to form infectivity-associated scrapie prion protein and generates C-terminal fragments C1 and C2 in healthy and prion-infected animals. C1 cleavage occurs N-terminally of PrP(C)'s hydrophobic domain, whereas the larger C2 fragment is generated by cleavage at the end of the octarepeat region. As the PrP-like proteins Doppel and Shadoo (Sho) have been reported to inhabit similar membrane environments as PrP(C), we investigated endoproteolysis by using a panel of mutant alleles. Doppel undergoes efficient in vivo cleavage at a C1 site mapped to the start of the globular domain, which is a structurally similar cleavage site to that in PrP(C). Sho is processed to C1 and C2 fragments, and proved refractory to mutagenesis to inactivate C1 cleavage. As a reciprocal product of C1 cleavage, Sho also engenders a metabolically stable N1 fragment with a C-terminus after its hydrophobic domain, an observation that may account for N1's association with membrane and/or cellular fractions in vitro and in vivo. Our data indicate that glycosylation status and yet to be identified proteases modulate internal C1 and C2 proteolysis events within the mammalian prion protein family.

The CNS glycoprotein Shadoo has PrP(C)-like protective properties and displays reduced levels in prion infections.

The cellular prion protein, PrP(C), is neuroprotective in a number of settings and in particular prevents cerebellar degeneration mediated by CNS-expressed Doppel or internally deleted PrP ('DeltaPrP'). This paradigm has facilitated mapping of activity determinants in PrP(C) and implicated a cryptic PrP(C)-like protein, 'pi'. Shadoo (Sho) is a hypothetical GPI-anchored protein encoded by the Sprn gene, exhibiting homology and domain organization similar to the N-terminus of PrP. Here we demonstrate Sprn expression and Sho protein in the adult CNS. Sho expression overlaps PrP(C), but is low in cerebellar granular neurons (CGNs) containing PrP(C) and high in PrP(C)-deficient dendritic processes. In Prnp(0/0) CGNs, Sho transgenes were PrP(C)-like in their ability to counteract neurotoxic effects of either Doppel or DeltaPrP. Additionally, prion-infected mice exhibit a dramatic reduction in endogenous Sho protein. Sho is a candidate for pi, and since it engenders a PrP(C)-like neuroprotective activity, compromised neuroprotective activity resulting from reduced levels may exacerbate damage in prion infections. Sho may prove useful in deciphering several unresolved facets of prion biology.

Dissociated phenotypes in presenilin transgenic mice define functionally distinct gamma-secretases.

Gamma-secretase depends on presence of presenilins (PS), Nct, Aph-1, and PEN-2 within a core complex. This endoproteolytic activity cleaves within transmembrane domains of amyloid-beta precursor protein (APP) and Notch, and familial Alzheimer's disease (FAD) mutations in PS1 or PS2 genes shift APP cleavage from production of amyloid-beta (Abeta) 40 peptide to greater production of Abeta42. Although studies in PS1/PS2-deficient embryonic cells define overlapping activities for these proteins, in vivo complementation of PS1-deficient animals described here reveals an unexpected spectrum of activities dictated by PS1 and PS2 alleles. Unlike PS1 transgenes, wild-type PS2 transgenes expressed in the mouse CNS support little Abeta40 or Abeta42 production, and FAD PS2 alleles support robust production of only Abeta42. Although wild-type PS2 transgenes failed to rescue Notch-associated skeletal defects in PS1 hypomorphs, a "gained" competence in this regard was apparent for FAD alleles of PS2. The range of discrete and divergent processing activities in mice reconstituted with different PS genes and alleles argues against gamma-secretase being a single enzyme with intrinsically relaxed substrate and cleavage site specificities. Instead, our studies define functionally distinct gamma-secretase variants. We speculate that extrinsic components, in combination with core complexes, may tailor functional variants of this enzyme to their preferred substrates.

Genetic mapping of activity determinants within cellular prion proteins: N-terminal modules in PrPC offset pro-apoptotic activity of the Doppel helix B/B' region.

The PrP-like Doppel (Dpl) protein causes apoptotic death of cerebellar neurons in transgenic mice, a process prevented by expression of the wild type (wt) cellular prion protein, PrP(C). Internally deleted forms of PrP(C) resembling Dpl such as PrPDelta32-121 produce a similar PrP(C)-sensitive pro-apoptotic phenotype in transgenic mice. Here we demonstrate that these phenotypic attributes of wt Dpl, wt PrP(C), and PrPDelta132-121 can be accurately recapitulated by transfected mouse cerebellar granule cell cultures. This system was then explored by mutagenesis of the co-expressed prion proteins to reveal functional determinants. By this means, neuroprotective activity of wt PrP(C) was shown to be nullified by a deletion of the N-terminal charged region implicated in endocytosis and retrograde axonal transport (PrPDelta23-28), by deletion of all five octarepeats (PrPDelta51-90), or by glycine replacement of four octarepeat histidine residues required for selective binding of copper ions (Prnp"H/G"). In the case of Dpl, overlapping deletions defined a requirement for the gene interval encoding helices B and B' (DplDelta101-125). These data suggest contributions of copper binding and neuronal trafficking to wt PrP(C) function in vivo and place constraints upon current hypotheses to explain Dpl/PrP(C) antagonism by competitive ligand binding. Further implementation of this assay should provide a fuller understanding of the attributes and subcellular localizations required for activity of these enigmatic proteins.