Short linear motifs - ex nihilo evolution of protein regulation

Recent advances in engineering synthetic biomolecular condensates

Biomolecular condensates are intriguing entities found within living cells. These structures possess the ability to selectively concentrate specific components through phase separation, thereby playing a crucial role in the spatiotemporal regulation of a wide range of cellular processes and metabolic activities. To date, extensive studies have been dedicated to unraveling the intricate connections between molecular features, physical properties, and cellular functions of condensates. This collective effort has paved the way for deliberate engineering of tailor-made condensates with specific applications. In this review, we comprehensively examine the underpinnings governing condensate formation. Next, we summarize the material states of condensates and delve into the design of synthetic intrinsically disordered proteins with tunable phase behaviors and physical properties. Subsequently, we review the diverse biological functions demonstrated by synthetic biomolecular condensates, encompassing gene regulation, cellular behaviors, modulation of biochemical reactions, and manipulation of endogenous protein activities. Lastly, we discuss future challenges and opportunities in constructing synthetic condensates with tunable physical properties and customized cellular functions, which may shed light on the development of new types of sophisticated condensate systems with distinct functions applicable to various scenarios.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known. Here, we use Saccharomyces cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, the resulting tyrosine phosphorylation is biologically spurious. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3500 proteins. The number of spurious pY sites generated correlates strongly with decreased growth, and we predict over 1000 pY events to be deleterious. However, we also find that many of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with tyrosine kinases. Our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

The Chlamydia effector IncE employs two short linear motifs to reprogram host vesicle trafficking

Chlamydia trachomatis, a leading cause of bacterial sexually transmitted infections, creates a specialized intracellular replicative niche by translocation and insertion of a diverse array of effectors (Incs [inclusion membrane proteins]) into the inclusion membrane. Here, we characterize IncE, a multifunctional Inc that encodes two non-overlapping short linear motifs (SLiMs) within its short cytosolic C terminus. The proximal SLiM, by mimicking just a small portion of an R-N-ethylmaleimide-sensitive factor adaptor protein receptor (SNARE) motif, binds and recruits syntaxin (STX)7- and STX12-containing vesicles to the inclusion. The distal SLiM mimics the sorting nexin (SNX)5 and SNX6 cargo binding site to recruit SNX6-containing vesicles to the inclusion. By simultaneously binding two distinct vesicle classes, IncE brings these vesicles in close apposition with each other at the inclusion to facilitate C. trachomatis intracellular development. Our work suggests that Incs may have evolved SLiMs to enable rapid evolution in a limited protein space to disrupt host cell processes.

Mitosis: An expanded view of mitotic mechanisms that arose in evolution

Mitosis exhibits astonishing evolutionary plasticity, with dividing eukaryotic cells differing in the organization of the mitotic spindle and the extent of nuclear envelope breakdown. A new study suggests that a multinucleated lifestyle may favor the evolution of closed nuclear division.

Unveiling epithelial plasticity regulation in lung cancer: Exploring the cross-talk among Tks4 scaffold protein partners

The epithelial-to-mesenchymal transition (EMT) represents a hallmark event in the evolution of lung cancer. This work aims to study a recently described EMT-regulating protein, Tks4, and to explore its potential as a prognostic biomarker in non-small cell lung cancer. In this study, we used CRISPR/Cas9 method to knockout (KO) Tks4 to study its functional roles in invadopodia formation, migration, and regulation of EMT marker expressions and we identified Tks4-interacting proteins. Tks4-KO A549 cells exhibited an EMT-like phenotype characterized by elongated morphology and increased expression of EMT markers. Furthermore, analyses of a large-scale lung cancer database and a patient-derived tissue array data revealed that the Tks4 mRNA level was decreased in more aggressive lung cancer stages. To understand the regulatory role of Tks4 in lung cancer, we performed a Tks4-interactome analysis via Tks4 immunoprecipitation-mass spectrometry on five different cell lines and identified CAPZA1 as a novel Tks4 partner protein. Thus, we propose that the absence of Tks4 leads to disruption of a connectome of multiple proteins and that the resulting undocking and likely mislocalization of signaling molecules impairs actin cytoskeleton rearrangement and activates EMT-like cell fate switches, both of which likely influence disease severity.

Elaboration of the Homer1 recognition landscape reveals incomplete divergence of paralogous EVH1 domains

Short sequences that mediate interactions with modular binding domains are ubiquitous throughout eukaryotic proteomes. Networks of short linear motifs (SLiMs) and their corresponding binding domains orchestrate many cellular processes, and the low mutational barrier to evolving novel interactions provides a way for biological systems to rapidly sample selectable phenotypes. Mapping SLiM binding specificity and the rules that govern SLiM evolution is fundamental to uncovering the pathways regulated by these networks and developing the tools to manipulate them. We used high-throughput screening of the human proteome to identify sequences that bind to the Enabled/VASP homology 1 (EVH1) domain of the postsynaptic density scaffolding protein Homer1. This expanded our understanding of the determinants of Homer EVH1 binding preferences and defined a new motif that can facilitate the discovery of additional Homer-mediated interactions. Interestingly, the Homer1 EVH1 domain preferentially binds to sequences containing an N-terminally overlapping motif that is bound by the paralogous family of Ena/VASP actin polymerases, and many of these sequences can bind to EVH1 domains from both protein families. We provide evidence from orthologous EVH1 domains in pre-metazoan organisms that the overlap in human Ena/VASP and Homer binding preferences corresponds to an incomplete divergence from a common Ena/VASP ancestor. Given this overlap in binding profiles, promiscuous sequences that can be recognized by both families either achieve specificity through extrinsic regulatory strategies or may provide functional benefits via multi-specificity. This may explain why these paralogs incompletely diverged despite the accessibility of further diverged isoforms.

Diverse genetic conflicts mediated by molecular mimicry and computational approaches to detect them

In genetic conflicts between intergenomic and selfish elements, driver and killer elements achieve biased survival, replication, or transmission over sensitive and targeted elements through a wide range of molecular mechanisms, including mimicry. Driving mechanisms manifest at all organismal levels, from the biased propagation of individual genes, as demonstrated by transposable elements, to the biased transmission of genomes, as illustrated by viruses, to the biased transmission of cell lineages, as in cancer. Targeted genomes are vulnerable to molecular mimicry through the conserved motifs they use for their own signaling and regulation. Mimicking these motifs enables an intergenomic or selfish element to control core target processes, and can occur at the sequence, structure, or functional level. Molecular mimicry was first appreciated as an important phenomenon more than twenty years ago. Modern genomics technologies, databases, and machine learning approaches offer tremendous potential to study the distribution of molecular mimicry across genetic conflicts in nature. Here, we explore the theoretical expectations for molecular mimicry between conflicting genomes, the trends in molecular mimicry mechanisms across known genetic conflicts, and outline how new examples can be gleaned from population genomic datasets. We discuss how mimics involving short sequence-based motifs or gene duplications can evolve convergently from new mutations. Whereas, processes that involve divergent domains or fully-folded structures occur among genomes by horizontal gene transfer. These trends are largely based on a small number of organisms and should be reevaluated in a general, phylogenetically independent framework. Currently, publicly available databases can be mined for genotypes driving non-Mendelian inheritance patterns, epistatic interactions, and convergent protein structures. A subset of these conflicting elements may be molecular mimics. We propose approaches for detecting genetic conflict and molecular mimicry from these datasets.

PairK: Pairwise k-mer alignment for quantifying protein motif conservation in disordered regions

Protein-protein interactions are often mediated by a modular peptide recognition domain binding to a short linear motif (SLiM) in the disordered region of another protein. The ability to predict domain-SLiM interactions would allow researchers to map protein interaction networks, predict the effects of perturbations to those networks, and develop biologically meaningful hypotheses. Unfortunately, sequence database searches for SLiMs generally yield mostly biologically irrelevant motif matches or false positives. To improve the prediction of novel SLiM interactions, researchers employ filters to discriminate between biologically relevant and improbable motif matches. One promising criterion for identifying biologically relevant SLiMs is the sequence conservation of the motif, exploiting the fact that functional motifs are more likely to be conserved than spurious motif matches. However, the difficulty of aligning disordered regions has significantly hampered the utility of this approach. We present PairK (pairwise k-mer alignment), an MSA-free method to quantify motif conservation in disordered regions. PairK outperforms both standard MSA-based conservation scores and a modern LLM-based conservation score predictor on the task of identifying biologically important motif instances. PairK can quantify conservation over wider phylogenetic distances than MSAs, indicating that SLiMs may be more conserved than is implied by MSA-based metrics. PairK is available as open-source code at https://github.com/jacksonh1/pairk.

Reconstruction of Eriocheir sinensis Protein-Protein Interaction Network Based on DGO-SVM Method

Eriocheir sinensis is an economically important aquatic animal. Its regulatory mechanisms underlying many biological processes are still vague due to the lack of systematic analysis tools. The protein-protein interaction network (PIN) is an important tool for the systematic analysis of regulatory mechanisms. In this work, a novel machine learning method, DGO-SVM, was applied to predict the protein-protein interaction (PPI) in E. sinensis, and its PIN was reconstructed. With the domain, biological process, molecular functions and subcellular locations of proteins as the features, DGO-SVM showed excellent performance in Bombyx mori, humans and five aquatic crustaceans, with 92-96% accuracy. With DGO-SVM, the PIN of E. sinensis was reconstructed, containing 14,703 proteins and 7,243,597 interactions, in which 35,604 interactions were associated with 566 novel proteins mainly involved in the response to exogenous stimuli, cellular macromolecular metabolism and regulation. The DGO-SVM demonstrated that the biological process, molecular functions and subcellular locations of proteins are significant factors for the precise prediction of PPIs. We reconstructed the largest PIN for E. sinensis, which provides a systematic tool for the regulatory mechanism analysis. Furthermore, the novel-protein-related PPIs in the PIN may provide important clues for the mechanism analysis of the underlying specific physiological processes in E. sinensis.

Identification and characterization of specific motifs in effector proteins of plant parasites using MOnSTER

Plant pathogens cause billions of dollars of crop loss every year and are a major threat to global food security. Identifying and characterizing pathogens effectors is crucial towards their improved control. Because of their poor sequence conservation, effector identification is challenging, and current methods generate too many candidates without indication for prioritizing experimental studies. In most phyla, effectors contain specific sequence motifs which influence their localization and targets in the plant. Therefore, there is an urgent need to develop bioinformatics tools tailored for pathogen effectors. To circumvent these limitations, we have developed MOnSTER a specific tool that identifies clusters of motifs of protein sequences (CLUMPs). MOnSTER can be fed with motifs identified by de novo tools or from databases such as Pfam and InterProScan. The advantage of MOnSTER is the reduction of motif redundancy by clustering them and associating a score. This score encompasses the physicochemical properties of AAs and the motif occurrences. We built up our method to identify discriminant CLUMPs in oomycetes effectors. Consequently, we applied MOnSTER on plant parasitic nematodes and identified six CLUMPs in about 60% of the known nematode candidate parasitism proteins. Furthermore, we found co-occurrences of CLUMPs with protein domains important for invasion and pathogenicity. The potentiality of this tool goes beyond the effector characterization and can be used to easily cluster motifs and calculate the CLUMP-score on any set of protein sequences.

Modulation of biophysical properties of nucleocapsid protein in the mutant spectrum of SARS-CoV-2

Genetic diversity is a hallmark of RNA viruses and the basis for their evolutionary success. Taking advantage of the uniquely large genomic database of SARS-CoV-2, we examine the impact of mutations across the spectrum of viable amino acid sequences on the biophysical phenotypes of the highly expressed and multifunctional nucleocapsid protein. We find variation in the physicochemical parameters of its extended intrinsically disordered regions (IDRs) sufficient to allow local plasticity, but also observe functional constraints that similarly occur in related coronaviruses. In biophysical experiments with several N-protein species carrying mutations associated with major variants, we find that point mutations in the IDRs can have nonlocal impact and modulate thermodynamic stability, secondary structure, protein oligomeric state, particle formation, and liquid-liquid phase separation. In the Omicron variant, distant mutations in different IDRs have compensatory effects in shifting a delicate balance of interactions controlling protein assembly properties, and include the creation of a new protein-protein interaction interface in the N-terminal IDR through the defining P13L mutation. A picture emerges where genetic diversity is accompanied by significant variation in biophysical characteristics of functional N-protein species, in particular in the IDRs.

The Chlamydia effector IncE employs two short linear motifs to reprogram host vesicle trafficking

Chlamydia trachomatis, a leading cause of bacteria sexually transmitted infections, creates a specialized intracellular replicative niche by translocation and insertion of a diverse array of effectors (Incs) into the inclusion membrane. Here, we characterize IncE, a multi-functional Inc that encodes two non-overlapping short linear motifs (SLiMs) within its short cytosolic C-terminus. The proximal SLiM mimics an R-SNARE motif to recruit syntaxin (STX) 7 and 12-containing vesicles to the inclusion. The distal SLiM mimics the Sorting Nexin (SNX) 5 and 6 cargo binding site to recruit SNX6-containing vesicles to the inclusion. By simultaneously binding to two distinct vesicle classes, IncE reprograms host cell trafficking to promote the formation of a class of hybrid vesicles at the inclusion that are required for C. trachomatis intracellular development. Our work suggests that Incs may have evolved SLiMs to facilitate rapid evolution in a limited protein space to disrupt host cell processes.

Elaboration of the Homer1 Recognition Landscape Reveals Incomplete Divergence of Paralogous EVH1 Domains

Short sequences that mediate interactions with modular binding domains are ubiquitous throughout eukaryotic proteomes. Networks of Short Linear Motifs (SLiMs) and their corresponding binding domains orchestrate many cellular processes, and the low mutational barrier to evolving novel interactions provides a way for biological systems to rapidly sample selectable phenotypes. Mapping SLiM binding specificity and the rules that govern SLiM evolution is fundamental to uncovering the pathways regulated by these networks and developing the tools to manipulate them. We used high-throughput screening of the human proteome to identify sequences that bind to the Enabled/VASP homology 1 (EVH1) domain of the postsynaptic density scaffolding protein Homer1. In doing so, we expanded current understanding of the determinants of Homer EVH1 binding preferences and defined a new motif that can facilitate the discovery of additional Homer-mediated interactions. Interestingly, the Homer1 EVH1 domain preferentially binds to sequences containing an N-terminally overlapping motif that is bound by the paralogous family of Ena/VASP actin polymerases, and many of these sequences can bind to EVH1 domains from both protein families. We provide evidence from orthologous EVH1 domains in pre-metazoan organisms that the overlap in human Ena/VASP and Homer binding preferences corresponds to an incomplete divergence from a common Ena/VASP ancestor. Given this overlap in binding profiles, promiscuous sequences that can be recognized by both families either achieve specificity through extrinsic regulatory strategies or may provide functional benefits via multi-specificity. This may explain why these paralogs incompletely diverged despite the accessibility of further diverged isoforms.

Structured Tandem Repeats in Protein Interactions

Tandem repeats (TRs) in protein sequences are consecutive, highly similar sequence motifs. Some types of TRs fold into structural units that pack together in ensembles, forming either an (open) elongated domain or a (closed) propeller, where the last unit of the ensemble packs against the first one. Here, we examine TR proteins (TRPs) to see how their sequence, structure, and evolutionary properties favor them for a function as mediators of protein interactions. Our observations suggest that TRPs bind other proteins using large, structured surfaces like globular domains; in particular, open-structured TR ensembles are favored by flexible termini and the possibility to tightly coil against their targets. While, intuitively, open ensembles of TRs seem prone to evolve due to their potential to accommodate insertions and deletions of units, these evolutionary events are unexpectedly rare, suggesting that they are advantageous for the emergence of the ancestral sequence but are early fixed. We hypothesize that their flexibility makes it easier for further proteins to adapt to interact with them, which would explain their large number of protein interactions. We provide insight into the properties of open TR ensembles, which make them scaffolds for alternative protein complexes to organize genes, RNA and proteins.

The molecular basis for cellular function of intrinsically disordered protein regions

Intrinsically disordered protein regions exist in a collection of dynamic interconverting conformations that lack a stable 3D structure. These regions are structurally heterogeneous, ubiquitous and found across all kingdoms of life. Despite the absence of a defined 3D structure, disordered regions are essential for cellular processes ranging from transcriptional control and cell signalling to subcellular organization. Through their conformational malleability and adaptability, disordered regions extend the repertoire of macromolecular interactions and are readily tunable by their structural and chemical context, making them ideal responders to regulatory cues. Recent work has led to major advances in understanding the link between protein sequence and conformational behaviour in disordered regions, yet the link between sequence and molecular function is less well defined. Here we consider the biochemical and biophysical foundations that underlie how and why disordered regions can engage in productive cellular functions, provide examples of emerging concepts and discuss how protein disorder contributes to intracellular information processing and regulation of cellular function.

The logic of protein post-translational modifications (PTMs): Chemistry, mechanisms and evolution of protein regulation through covalent attachments

Protein post-translational modifications (PTMs) play a crucial role in all cellular functions by regulating protein activity, interactions and half-life. Despite the enormous diversity of modifications, various PTM systems show parallels in their chemical and catalytic underpinnings. Here, focussing on modifications that involve the addition of new elements to amino-acid sidechains, I describe historical milestones and fundamental concepts that support the current understanding of PTMs. The historical survey covers selected key research programmes, including the study of protein phosphorylation as a regulatory switch, protein ubiquitylation as a degradation signal and histone modifications as a functional code. The contribution of crucial techniques for studying PTMs is also discussed. The central part of the essay explores shared chemical principles and catalytic strategies observed across diverse PTM systems, together with mechanisms of substrate selection, the reversibility of PTMs by erasers and the recognition of PTMs by reader domains. Similarities in the basic chemical mechanism are highlighted and their implications are discussed. The final part is dedicated to the evolutionary trajectories of PTM systems, beginning with their possible emergence in the context of rivalry in the prokaryotic world. Together, the essay provides a unified perspective on the diverse world of major protein modifications.

Linear motifs regulating protein secretion, sorting and autophagy in Leishmania parasites are diverged with respect to their host equivalents

The pathogenic, tropical Leishmania flagellates belong to an early-branching eukaryotic lineage (Kinetoplastida) with several unique features. Unfortunately, they are poorly understood from a molecular biology perspective, making development of mechanistically novel and selective drugs difficult. Here, we explore three functionally critical targeting short linear motif systems as well as their receptors in depth, using a combination of structural modeling, evolutionary sequence divergence and deep learning. Secretory signal peptides, endoplasmic reticulum (ER) retention motifs (KDEL motifs), and autophagy signals (motifs interacting with ATG8 family members) are ancient and essential components of cellular life. Although expected to be conserved amongst the kinetoplastids, we observe that all three systems show a varying degree of divergence from their better studied equivalents in animals, plants, or fungi. We not only describe their behaviour, but also build models that allow the prediction of localization and potential functions for several uncharacterized Leishmania proteins. The unusually Ala/Val-rich secretory signal peptides, endoplasmic reticulum resident proteins ending in Asp-Leu-COOH and atypical ATG8-like proteins are all unique molecular features of kinetoplastid parasites. Several of their critical protein-protein interactions could serve as targets of selective antimicrobial agents against Leishmaniasis due to their systematic divergence from the host.

Design of functional intrinsically disordered proteins

Many proteins do not fold into a fixed three-dimensional structure, but rather function in a highly disordered state. These intrinsically disordered proteins pose a unique challenge to protein engineering and design: How can proteins be designed de novo if not by tailoring their structure? Here, we will review the nascent field of design of intrinsically disordered proteins with focus on applications in biotechnology and medicine. The design goals should not necessarily be the same as for de novo design of folded proteins as disordered proteins have unique functional strengths and limitations. We focus on functions where intrinsically disordered proteins are uniquely suited including disordered linkers, desiccation chaperones, sensors of the chemical environment, delivery of pharmaceuticals, and constituents of biomolecular condensates. Design of functional intrinsically disordered proteins relies on a combination of computational tools and heuristics gleaned from sequence-function studies. There are few cases where intrinsically disordered proteins have made it into industrial applications. However, we argue that disordered proteins can perform many roles currently performed by organic polymers, and that these proteins might be more designable due to their modularity.

B56δ long-disordered arms form a dynamic PP2A regulation interface coupled with global allostery and Jordan's syndrome mutations

Intrinsically disordered regions (IDR) and short linear motifs (SLiMs) play pivotal roles in the intricate signaling networks governed by phosphatases and kinases. B56δ (encoded by PPP2R5D) is a regulatory subunit of protein phosphatase 2A (PP2A) with long IDRs that harbor a substrate-mimicking SLiM and multiple phosphorylation sites. De novo missense mutations in PPP2R5D cause intellectual disabilities (ID), macrocephaly, Parkinsonism, and a broad range of neurological symptoms. Our single-particle cryo-EM structures of the PP2A-B56δ holoenzyme reveal that the long, disordered arms at the B56δ termini fold against each other and the holoenzyme core. This architecture suppresses both the phosphatase active site and the substrate-binding protein groove, thereby stabilizing the enzyme in a closed latent form with dual autoinhibition. The resulting interface spans over 190 Å and harbors unfavorable contacts, activation phosphorylation sites, and nearly all residues with ID-associated mutations. Our studies suggest that this dynamic interface is coupled to an allosteric network responsive to phosphorylation and altered globally by mutations. Furthermore, we found that ID mutations increase the holoenzyme activity and perturb the phosphorylation rates, and the severe variants significantly increase the mitotic duration and error rates compared to the normal variant.

Diversity of short linear interaction motifs in SARS-CoV-2 nucleocapsid protein

IMPORTANCE

Short linear motifs (SLiMs) are 3-10 amino acid long binding motifs in intrinsically disordered protein regions (IDRs) that serve as ubiquitous protein-protein interaction modules in eukaryotic cells. Through molecular mimicry, viruses hijack these sequence motifs to control host cellular processes. It is thought that the small size of SLiMs and the high mutation frequencies of viral IDRs allow rapid host adaptation. However, a salient characteristic of RNA viruses, due to high replication errors, is their obligate existence as mutant swarms. Taking advantage of the uniquely large genomic database of SARS-CoV-2, here, we analyze the role of sequence diversity in the presentation of SLiMs, focusing on the highly abundant, multi-functional nucleocapsid protein. We find that motif mimicry is a highly dynamic process that produces an abundance of motifs transiently present in subsets of mutant species. This diversity allows the virus to efficiently explore eukaryotic motifs and evolve the host-virus interface.

Exploiting the endogenous yeast nuclear proteome to identify short linear motifs in vivo

Peptide-domain interactions mediated by short linear motifs (SLiMs) play crucial roles in cellular biology. The simplicity of SLiMs poses challenges in their computational identification. Existing high-throughput methods for discovering SLiMs lack cellular context as they are typically performed in vitro. We developed a functional selection method using yeast to identify peptides that interact with the endogenous yeast nuclear proteome. Remarkably, peptides selected for in yeast also mediated nuclear import in human cells. Notably, the identified peptides did not resemble classical nuclear localization sequences. This platform has the potential to identify and investigate motifs that interact with the nuclear proteome of yeast and human and to aid in the identification and understanding of alternative protein nuclear import mechanisms.

The fitness cost of spurious phosphorylation

The fidelity of signal transduction requires the binding of regulatory molecules to their cognate targets. However, the crowded cell interior risks off-target interactions between proteins that are functionally unrelated. How such off-target interactions impact fitness is not generally known, but quantifying this is required to understand the constraints faced by cell systems as they evolve. Here, we use the model organism S. cerevisiae to inducibly express tyrosine kinases. Because yeast lacks bona fide tyrosine kinases, most of the resulting tyrosine phosphorylation is spurious. This provides a suitable system to measure the impact of artificial protein interactions on fitness. We engineered 44 yeast strains each expressing a tyrosine kinase, and quantitatively analysed their phosphoproteomes. This analysis resulted in ~30,000 phosphosites mapping to ~3,500 proteins. Examination of the fitness costs in each strain revealed a strong correlation between the number of spurious pY sites and decreased growth. Moreover, the analysis of pY effects on protein structure and on protein function revealed over 1000 pY events that we predict to be deleterious. However, we also find that a large number of the spurious pY sites have a negligible effect on fitness, possibly because of their low stoichiometry. This result is consistent with our evolutionary analyses demonstrating a lack of phosphotyrosine counter-selection in species with bona fide tyrosine kinases. Taken together, our results suggest that, alongside the risk for toxicity, the cell can tolerate a large degree of non-functional crosstalk as interaction networks evolve.

Computational analysis of regulatory regions in human protein kinases

Eukaryotic proteins often feature modular domain structures comprising globular domains that are connected by linker regions and intrinsically disordered regions that may contain important functional motifs. The intramolecular interactions of globular domains and nonglobular regions can play critical roles in different aspects of protein function. However, studying these interactions and their regulatory roles can be challenging due to the flexibility of nonglobular regions, the long insertions separating interacting modules, and the transient nature of some interactions. Obtaining the experimental structures of multiple domains and functional regions is more difficult than determining the structures of individual globular domains. High-quality structural models generated by AlphaFold offer a unique opportunity to study intramolecular interactions in eukaryotic proteins. In this study, we systematically explored intramolecular interactions between human protein kinase domains (KDs) and potential regulatory regions, including globular domains, N- and C-terminal tails, long insertions, and distal nonglobular regions. Our analysis identified intramolecular interactions between human KDs and 35 different types of globular domains, exhibiting a variety of interaction modes that could contribute to orthosteric or allosteric regulation of kinase activity. We also identified prevalent interactions between human KDs and their flanking regions (N- and C-terminal tails). These interactions exhibit group-specific characteristics and can vary within each specific kinase group. Although long-range interactions between KDs and nonglobular regions are relatively rare, structural details of these interactions offer new insights into the regulation mechanisms of several kinases, such as HASPIN, MAPK7, MAPK15, and SIK1B.

Ehrlichia Wnt SLiM ligand mimic deactivates the Hippo pathway to engage the anti-apoptotic Yap-GLUT1-BCL-xL axis

Ehrlichia chaffeensis TRP120 effector has evolved short linear motif (SLiM) ligand mimicry to repurpose multiple evolutionarily conserved cellular signaling pathways, including Wnt, Notch, and Hedgehog. In this investigation, we demonstrate that E. chaffeensis and recombinant TRP120 deactivate Hippo signaling, resulting in the activation of Hippo transcription coactivator Yes-associated protein (Yap). Moreover, a homologous 6 amino acid (QDVASH) SLiM shared by TRP120 and Wnt3a/5a ligands phenocopied Yap and β-catenin activation induced by E. chaffeensis, rTRP120, and Wnt5a. Similar Hippo gene expression profiles were also stimulated by E. chaffeensis, rTRP120, SLiM, and Wnt5a. Single siRNA knockdown of Hippo transcription co-activator/factors, Yap, and transcriptional enhanced associate domain (TEAD) significantly decreased E. chaffeensis infection. Yap activation was abolished in THP-1 Wnt Frizzled-5 (Fzd5) receptor knockout cells (KO), demonstrating Fzd5 receptor dependence. In addition, the TRP120-Wnt-SLiM antibody blocked Hippo deactivation (Yap activation). Expression of anti-apoptotic Hippo target gene SLC2A1 (encodes glucose transporter 1; GLUT1) was upregulated by E. chaffeensis and corresponded to increased levels of GLUT1. Conversely, siRNA knockdown of SLC2A1 significantly inhibited infection. Higher GLUT1 levels correlated with increased B cell lymphoma-extra large (BCL-xL) and decreased BCL2-associated X, apoptosis regulator (Bax) levels. Moreover, blocking Yap activation with the inhibitor Verteporfin induced apoptosis that corresponded to significant reductions in GLUT1 and BCL-xL levels and activation of Bax and Caspase-3 and -9. This study identifies a novel shared Wnt/Hippo SLiM ligand mimic and demonstrates that E. chaffeensis deactivates the Hippo pathway to engage the anti-apoptotic Yap-GLUT1-BCL-xL axis.

Phosphorylation of Schizosaccharomyces pombe Dss1 mediates direct binding to the ubiquitin-ligase Dma1 in vitro

Intrinsically disordered proteins (IDPs) are often multifunctional and frequently posttranslationally modified. Deleted in split hand/split foot 1 (Dss1-Sem1 in budding yeast) is a highly multifunctional IDP associated with a range of protein complexes. However, it remains unknown if the different functions relate to different modified states. In this work, we show that Schizosaccharomyces pombe Dss1 is a substrate for casein kinase 2 in vitro, and we identify three phosphorylated threonines in its linker region separating two known disordered ubiquitin-binding motifs. Phosphorylations of the threonines had no effect on ubiquitin-binding but caused a slight destabilization of the C-terminal α-helix and mediated a direct interaction with the forkhead-associated (FHA) domain of the RING-FHA E3-ubiquitin ligase defective in mitosis 1 (Dma1). The phosphorylation sites are not conserved and are absent in human Dss1. Sequence analyses revealed that the Txx(E/D) motif, which is important for phosphorylation and Dma1 binding, is not linked to certain branches of the evolutionary tree. Instead, we find that the motif appears randomly, supporting the mechanism of ex nihilo evolution of novel motifs. In support of this, other threonine-based motifs, although frequent, are nonconserved in the linker, pointing to additional functions connected to this region. We suggest that Dss1 acts as an adaptor protein that docks to Dma1 via the phosphorylated FHA-binding motifs, while the C-terminal α-helix is free to bind mitotic septins, thereby stabilizing the complex. The presence of Txx(D/E) motifs in the disordered regions of certain septin subunits may be of further relevance to the formation and stabilization of these complexes.

Updated protein domain annotation of the PARP protein family sheds new light on biological function

AlphaFold2 and related computational tools have greatly aided studies of structural biology through their ability to accurately predict protein structures. In the present work, we explored AF2 structural models of the 17 canonical members of the human PARP protein family and supplemented this analysis with new experiments and an overview of recent published data. PARP proteins are typically involved in the modification of proteins and nucleic acids through mono or poly(ADP-ribosyl)ation, but this function can be modulated by the presence of various auxiliary protein domains. Our analysis provides a comprehensive view of the structured domains and long intrinsically disordered regions within human PARPs, offering a revised basis for understanding the function of these proteins. Among other functional insights, the study provides a model of PARP1 domain dynamics in the DNA-free and DNA-bound states and enhances the connection between ADP-ribosylation and RNA biology and between ADP-ribosylation and ubiquitin-like modifications by predicting putative RNA-binding domains and E2-related RWD domains in certain PARPs. In line with the bioinformatic analysis, we demonstrate for the first time PARP14's RNA-binding capability and RNA ADP-ribosylation activity in vitro. While our insights align with existing experimental data and are probably accurate, they need further validation through experiments.

Substrate and phosphorylation site selection by phosphoprotein phosphatases

Dynamic protein phosphorylation and dephosphorylation are essential regulatory mechanisms that ensure proper cellular signaling and biological functions. Deregulation of either reaction has been implicated in several human diseases. Here, we focus on the mechanisms that govern the specificity of the dephosphorylation reaction. Most cellular serine/threonine dephosphorylation is catalyzed by 13 highly conserved phosphoprotein phosphatase (PPP) catalytic subunits, which form hundreds of holoenzymes by binding to regulatory and scaffolding subunits. PPP holoenzymes recognize phosphorylation site consensus motifs and interact with short linear motifs (SLiMs) or structural elements distal to the phosphorylation site. We review recent advances in understanding the mechanisms of PPP site-specific dephosphorylation preference and substrate recruitment and highlight examples of their interplay in the regulation of cell division.

Distinct Clades of Protein Phosphatase 2A Regulatory B'/B56 Subunits Engage in Different Physiological Processes

Protein phosphatase 2A (PP2A) is a strongly conserved and major protein phosphatase in all eukaryotes. The canonical PP2A complex consists of a catalytic (C), scaffolding (A), and regulatory (B) subunit. Plants have three groups of evolutionary distinct B subunits: B55, B' (B56), and B''. Here, the Arabidopsis B' group is reviewed and compared with other eukaryotes. Members of the B'α/B'β clade are especially important for chromatid cohesion, and dephosphorylation of transcription factors that mediate brassinosteroid (BR) signaling in the nucleus. Other B' subunits interact with proteins at the cell membrane to dampen BR signaling or harness immune responses. The transition from vegetative to reproductive phase is influenced differentially by distinct B' subunits; B'α and B'β being of little importance, whereas others (B'γ, B'ζ, B'η, B'θ, B'κ) promote transition to flowering. Interestingly, the latter B' subunits have three motifs in a conserved manner, i.e., two docking sites for protein phosphatase 1 (PP1), and a POLO consensus phosphorylation site between these motifs. This supports the view that a conserved PP1-PP2A dephosphorelay is important in a variety of signaling contexts throughout eukaryotes. A profound understanding of these regulators may help in designing future crops and understand environmental issues.

Predicting evolutionary outcomes through the probability of accessing sequence variants

Natural selection can only operate on available genetic variation. Thus, determining the probability of accessing different sequence variants from a starting sequence can help predict evolutionary trajectories and outcomes. We define the concept of "variant accessibility" as the probability that a set of genotypes encoding a particular protein function will arise through mutations before subject to natural selection. This probability is shaped by the mutational biases of nucleotides and the structure of the genetic code. Using the influenza A virus as a model, we discuss how a more accessible but less fit variant can emerge as an adaptation rather than a more fit variant. We describe a genotype-accessibility landscape, complementary to the genotype-fitness landscape, that informs the likelihood of a starting sequence reaching different parts of genotype space. The proposed framework lays the foundation for predicting the emergence of adaptive genotypes in evolving systems such as viruses and tumors.

The nexus between RNA-binding proteins and their effectors

RNA-binding proteins (RBPs) regulate essentially every event in the lifetime of an RNA molecule, from its production to its destruction. Whereas much has been learned about RNA sequence specificity and general functions of individual RBPs, the ways in which numerous RBPs instruct a much smaller number of effector molecules, that is, the core engines of RNA processing, as to where, when and how to act remain largely speculative. Here, we survey the known modes of communication between RBPs and their effectors with a particular focus on converging RBP-effector interactions and their roles in reducing the complexity of RNA networks. We discern the emerging unifying principles and discuss their utility in our understanding of RBP function, regulation of biological processes and contribution to human disease.

SLiM-binding pockets: an attractive target for broad-spectrum antivirals

Short linear motif (SLiM)-mediated interactions offer a unique strategy for viral intervention due to their compact interfaces, ease of convergent evolution, and key functional roles. Consequently, many viruses extensively mimic host SLiMs to hijack or deregulate cellular pathways and the same motif-binding pocket is often targeted by numerous unrelated viruses. A toolkit of therapeutics targeting commonly mimicked SLiMs could provide prophylactic and therapeutic broad-spectrum antivirals and vastly improve our ability to treat ongoing and future viral outbreaks. In this opinion article, we discuss the therapeutic relevance of SLiMs, advocating their suitability as targets for broad-spectrum antiviral inhibitors.

Extended regulation interface coupled to the allosteric network and disease mutations in the PP2A-B56δ holoenzyme

An increasing number of mutations associated with devastating human diseases are diagnosed by whole-genome/exon sequencing. Recurrent de novo missense mutations have been discovered in B56δ (encoded by PPP2R5D), a regulatory subunit of protein phosphatase 2A (PP2A), that cause intellectual disabilities (ID), macrocephaly, Parkinsonism, and a broad range of neurological symptoms. Single-particle cryo-EM structures show that the PP2A-B56δ holoenzyme possesses closed latent and open active forms. In the closed form, the long, disordered arms of B56δ termini fold against each other and the holoenzyme core, establishing dual autoinhibition of the phosphatase active site and the substrate-binding protein groove. The resulting interface spans over 190 Å and harbors unfavorable contacts, activation phosphorylation sites, and nearly all residues with ID-associated mutations. Our studies suggest that this dynamic interface is close to an allosteric network responsive to activation phosphorylation and altered globally by mutations. Furthermore, we found that ID mutations perturb the activation phosphorylation rates, and the severe variants significantly increase the mitotic duration and error rates compared to the wild variant.

Ehrlichia Wnt short linear motif ligand mimetic deactivates the Hippo pathway to engage the anti-apoptotic Yap-GLUT1-BCL-xL axis

Ehrlichia chaffeensis TRP120 effector has evolved short linear motif (SLiM) ligand mimicry to repurpose multiple evolutionarily conserved cellular signaling pathways including Wnt, Notch and Hedgehog. In this investigation, we demonstrate that E. chaffeensis and recombinant TRP120 deactivate Hippo signaling resulting in activation of Hippo transcription coactivator Yap and target gene expression. Moreover, a homologous 6 amino acid (QDVASH) SLiM shared by TRP120 and Wnt3a/5a ligands phenocopied Yap and β-catenin activation induced by E. chaffeensis, rTRP120 and Wnt5a. Similar Hippo gene expression profiles were also stimulated by E. chaffeensis, rTRP120, SLiM and Wnt5a. Single siRNA knockdown of Hippo transcription co-activator/factors (Yap and TEAD) significantly decreased E. chaffeensis infection. Yap activation was abolished in THP-1 Wnt Frizzled-5 (Fzd5) receptor knockout cells (KO), demonstrating Fzd5 receptor dependence. In addition, TRP120 Wnt-SLiM antibody blocked Hippo deactivation (Yap activation). Expression of anti-apoptotic Hippo target gene SLC2A1 (encodes glucose transporter 1; GLUT1) was upregulated by E. chaffeensis and corresponded to increased levels of GLUT1. Conversely, siRNA knockdown of SLC2A1 significantly inhibited infection. Higher GLUT1 levels correlated with increased BCL-xL and decreased Bax levels. Moreover, blocking Yap activation with the inhibitor Verteporfin induced apoptosis that corresponded to significant reductions in levels of GLUT1 and BCL-xL, and activation of Bax and Caspase-3 and -9. This study identifies a novel shared Wnt/Hippo SLiM ligand mimetic and demonstrates that E. chaffeensis deactivates the Hippo pathway to engage the anti-apoptotic Yap-GLUT1-BCL-xL axis.

Evolutionary Conserved Short Linear Motifs Provide Insights into the Cellular Response to Stress

Short linear motifs (SLiMs) are evolutionarily conserved functional modules of proteins composed of 3 to 10 residues and involved in multiple cellular functions. Here, we performed a search for SLiMs that exert sequence similarity to two segments of alpha-fetoprotein (AFP), a major mammalian embryonic and cancer-associated protein. Biological activities of the peptides, LDSYQCT (AFP_14-20) and EMTPVNPGV (GIP-9), have been previously confirmed under in vitro and in vivo conditions. In our study, we retrieved a vast array of proteins that contain SLiMs of interest from both prokaryotic and eukaryotic species, including viruses, bacteria, archaea, invertebrates, and vertebrates. Comprehensive Gene Ontology enrichment analysis showed that proteins from multiple functional classes, including enzymes, transcription factors, as well as those involved in signaling, cell cycle, and quality control, and ribosomal proteins were implicated in cellular adaptation to environmental stress conditions. These include response to oxidative and metabolic stress, hypoxia, DNA and RNA damage, protein degradation, as well as antimicrobial, antiviral, and immune response. Thus, our data enabled insights into the common functions of SLiMs evolutionary conserved across all taxonomic categories. These SLiMs can serve as important players in cellular adaptation to stress, which is crucial for cell functioning.

Evolution of SLiM-mediated hijack functions in intrinsically disordered viral proteins

Viruses and their hosts are involved in an ‘arms race’ where they continually evolve mechanisms to overcome each other. It has long been proposed that intrinsic disorder provides a substrate for the evolution of viral hijack functions and that short linear motifs (SLiMs) are important players in this process. Here, we review evidence in support of this tenet from two model systems: the papillomavirus E7 protein and the adenovirus E1A protein. Phylogenetic reconstructions reveal that SLiMs appear and disappear multiple times across evolution, providing evidence of convergent evolution within individual viral phylogenies. Multiple functionally related SLiMs show strong coevolution signals that persist across long distances in the primary sequence and occur in unrelated viral proteins. Moreover, changes in SLiMs are associated with changes in phenotypic traits such as host range and tropism. Tracking viral evolutionary events reveals that host switch events are associated with the loss of several SLiMs, suggesting that SLiMs are under functional selection and that changes in SLiMs support viral adaptation. Fine-tuning of viral SLiM sequences can improve affinity, allowing them to outcompete host counterparts. However, viral SLiMs are not always competitive by themselves, and tethering of two suboptimal SLiMs by a disordered linker may instead enable viral hijack. Coevolution between the SLiMs and the linker indicates that the evolution of disordered regions may be more constrained than previously thought. In summary, experimental and computational studies support a role for SLiMs and intrinsic disorder in viral hijack functions and in viral adaptive evolution.

An Evolutionary Perspective on the Origin, Conservation and Binding Partner Acquisition of Tankyrases

Tankyrases are poly-ADP-ribosyltransferases that regulate many crucial and diverse cellular processes in humans such as Wnt signaling, telomere homeostasis, mitotic spindle formation and glucose metabolism. While tankyrases are present in most animals, functional differences across species may exist. In this work, we confirm the widespread distribution of tankyrases throughout the branches of multicellular animal life and identify the single-celled choanoflagellates as earliest origin of tankyrases. We further show that the sequences and structural aspects of TNKSs are well-conserved even between distantly related species. We also experimentally characterized an anciently diverged tankyrase homolog from the sponge Amphimedon queenslandica and show that the basic functional aspects, such as poly-ADP-ribosylation activity and interaction with the canonical tankyrase binding peptide motif, are conserved. Conversely, the presence of tankyrase binding motifs in orthologs of confirmed interaction partners varies greatly between species, indicating that tankyrases may have different sets of interaction partners depending on the animal lineage. Overall, our analysis suggests a remarkable degree of conservation for tankyrases, and that their regulatory functions in cells have likely changed considerably throughout evolution.

Molecular determinants of TRAF6 binding specificity suggest that native interaction partners are not optimized for affinity

TRAF6 is an adaptor protein involved in signaling pathways that are essential for development and the immune system. It participates in many protein-protein interactions, some of which are mediated by the C-terminal MATH domain, which binds to short peptide segments containing the motif PxExx[FYWHDE], where x is any amino acid. Blocking MATH domain interactions is associated with favorable effects in various disease models. To better define TRAF6 MATH domain binding preferences, we screened a combinatorial library using bacterial cell-surface peptide display. We identified 236 of the best TRAF6-interacting peptides and a set of 1,200 peptides that match the sequence PxE but do not bind TRAF6 MATH. The peptides that were most enriched in the screen bound TRAF6 tighter than previously measured native peptides. To better understand the structural basis for TRAF6 interaction preferences, we built all-atom structural models of the MATH domain in complex with high-affinity binders and nonbinders identified in the screen. We identified favorable interactions for motif features in binders as well as negative design elements distributed across the motif that can disfavor or preclude binding. Searching the human proteome revealed that the most biologically relevant TRAF6 motif matches occupy a different sequence space from the best hits discovered in combinatorial library screening, suggesting that native interactions are not optimized for affinity. Our experimentally determined binding preferences and structural models support the design of peptide-based interaction inhibitors with higher affinities than endogenous TRAF6 ligands.

The evolution and polymorphism of mono-amino acid repeats in androgen receptor and their regulatory role in health and disease

Androgen receptor (AR) is a key member of nuclear hormone receptors with the longest intrinsically disordered N-terminal domain (NTD) in its protein family. There are four mono-amino acid repeats (polyQ1, polyQ2, polyG, and polyP) located within its NTD, of which two are polymorphic (polyQ1 and polyG). The length of both polymorphic repeats shows clinically important correlations with disease, especially with cancer and neurodegenerative diseases, as shorter and longer alleles exhibit significant differences in expression, activity and solubility. Importantly, AR has also been shown to undergo condensation in the nucleus by liquid-liquid phase separation, a process highly sensitive to protein solubility and concentration. Nonetheless, in prostate cancer cells, AR variants also partition into transcriptional condensates, which have been shown to alter the expression of target gene products. In this review, we summarize current knowledge on the link between AR repeat polymorphisms and cancer types, including mechanistic explanations and models comprising the relationship between condensate formation, polyQ1 length and transcriptional activity. Moreover, we outline the evolutionary paths of these recently evolved amino acid repeats across mammalian species, and discuss new research directions with potential breakthroughs and controversies in the literature.

Functional Tuning of Intrinsically Disordered Regions in Human Proteins by Composition Bias

Intrinsically disordered regions (IDRs) in protein sequences are flexible, have low structural constraints and as a result have faster rates of evolution. This lack of evolutionary conservation greatly limits the use of sequence homology for the classification and functional assessment of IDRs, as opposed to globular domains. The study of IDRs requires other properties for their classification and functional prediction. While composition bias is not a necessary property of IDRs, compositionally biased regions (CBRs) have been noted as frequent part of IDRs. We hypothesized that to characterize IDRs, it could be helpful to study their overlap with particular types of CBRs. Here, we evaluate this overlap in the human proteome. A total of 2/3 of residues in IDRs overlap CBRs. Considering CBRs enriched in one type of amino acid, we can distinguish CBRs that tend to be fully included within long IDRs (R, H, N, D, P, G), from those that partially overlap shorter IDRs (S, E, K, T), and others that tend to overlap IDR terminals (Q, A). CBRs overlap more often IDRs in nuclear proteins and in proteins involved in liquid-liquid phase separation (LLPS). Study of protein interaction networks reveals the enrichment of CBRs in IDRs by tandem repetition of short linear motifs (rich in S or P), and the existence of E-rich polar regions that could support specific protein interactions with non-specific interactions. Our results open ways to pin down the function of IDRs from their partial compositional biases.

PP1 phosphatases control PAR-2 localization and polarity establishment in C. elegans embryos

Cell polarity relies on the asymmetric distribution of the conserved PAR proteins, which is regulated by phosphorylation/dephosphorylation reactions. While the kinases involved have been well studied, the role of phosphatases remains poorly understood. In Caenorhabditis elegans zygotes, phosphorylation of the posterior PAR-2 protein by the atypical protein kinase PKC-3 inhibits PAR-2 cortical localization. Polarity establishment depends on loading of PAR-2 at the posterior cortex. We show that the PP1 phosphatases GSP-1 and GSP-2 are required for polarity establishment in embryos. We find that codepletion of GSP-1 and GSP-2 abrogates the cortical localization of PAR-2 and that GSP-1 and GSP-2 interact with PAR-2 via a PP1 docking motif in PAR-2. Mutating this motif in vivo, to prevent binding of PAR-2 to PP1, abolishes cortical localization of PAR-2, while optimizing this motif extends PAR-2 cortical localization. Our data suggest a model in which GSP-1/-2 counteracts PKC-3 phosphorylation of PAR-2, allowing its cortical localization at the posterior and polarization of the one-cell embryo.

A functional interaction between liprin-α1 and B56γ regulatory subunit of protein phosphatase 2A supports tumor cell motility

Scaffold liprin-α1 is required to assemble dynamic plasma membrane-associated platforms (PMAPs) at the front of migrating breast cancer cells, to promote protrusion and invasion. We show that the N-terminal region of liprin-α1 contains an LxxIxE motif interacting with B56 regulatory subunits of serine/threonine protein phosphatase 2A (PP2A). The specific interaction of B56γ with liprin-α1 requires an intact motif, since two point mutations strongly reduce the interaction. B56γ mediates the interaction of liprin-α1 with the heterotrimeric PP2A holoenzyme. Most B56γ protein is recovered in the cytosolic fraction of invasive MDA-MB-231 breast cancer cells, where B56γ is complexed with liprin-α1. While mutation of the short linear motif (SLiM) does not affect localization of liprin-α1 to PMAPs, localization of B56γ at these sites specifically requires liprin-α1. Silencing of B56γ or liprin-α1 inhibits to similar extent cell spreading on extracellular matrix, invasion, motility and lamellipodia dynamics in migrating MDA-MB-231 cells, suggesting that B56γ/PP2A is a novel component of the PMAPs machinery regulating tumor cell motility. In this direction, inhibition of cell spreading by silencing liprin-α1 is not rescued by expression of B56γ binding-defective liprin-α1 mutant. We propose that liprin-α1-mediated recruitment of PP2A via B56γ regulates cell motility by controlling protrusion in migrating MDA-MB-231 cells.

Structural basis for tunable affinity and specificity of LxCxE-dependent protein interactions with the retinoblastoma protein family

The retinoblastoma protein (Rb) and its homologs p107 and p130 are critical regulators of gene expression during the cell cycle and are commonly inactivated in cancer. Rb proteins use their "pocket domain" to bind an LxCxE sequence motif in other proteins, many of which function with Rb proteins to co-regulate transcription. Here, we present binding data and crystal structures of the p107 pocket domain in complex with LxCxE peptides from the transcriptional co-repressor proteins HDAC1, ARID4A, and EID1. Our results explain why Rb and p107 have weaker affinity for cellular LxCxE proteins compared with the E7 protein from human papillomavirus, which has been used as the primary model for understanding LxCxE motif interactions. Our structural and mutagenesis data also identify and explain differences in Rb and p107 affinities for some LxCxE-containing sequences. Our study provides new insights into how Rb proteins bind their cell partners with varying affinity and specificity.

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.

Evolution of short linear motifs and disordered proteins Topic: yeast as model system to study evolution

Evolutionary preservation of protein structure had a major influence on the field of molecular evolution: changes in individual amino acids that did not disrupt protein folding would either have no effect or subtly change the 'lock' so that it could fit a new 'key'. Homology of individual amino acids could be confidently assigned through sequence alignments, and models of evolution could be tested. This view of molecular evolution excluded large regions of proteins that could not be confidently aligned, such as intrinsically disordered regions (IDRs) that do not fold into stable structures. In the last decade, major progress has been made in understanding the evolution of IDRs, much of it facilitated by new experimental and computational approaches in yeast. Here, we review this progress as well as several still outstanding questions.

Coupling to short linear motifs creates versatile PME-1 activities in PP2A holoenzyme demethylation and inhibition

Protein phosphatase 2A (PP2A) holoenzymes target broad substrates by recognizing short motifs via regulatory subunits. PP2A methylesterase 1 (PME-1) is a cancer-promoting enzyme and undergoes methylesterase activation upon binding to the PP2A core enzyme. Here, we showed that PME-1 readily demethylates different families of PP2A holoenzymes and blocks substrate recognition in vitro. The high-resolution cryoelectron microscopy structure of a PP2A-B56 holoenzyme–PME-1 complex reveals that PME-1 disordered regions, including a substrate-mimicking motif, tether to the B56 regulatory subunit at remote sites. They occupy the holoenzyme substrate-binding groove and allow large structural shifts in both holoenzyme and PME-1 to enable multipartite contacts at structured cores to activate the methylesterase. B56 interface mutations selectively block PME-1 activity toward PP2A-B56 holoenzymes and affect the methylation of a fraction of total cellular PP2A. The B56 interface mutations allow us to uncover B56-specific PME-1 functions in p53 signaling. Our studies reveal multiple mechanisms of PME-1 in suppressing holoenzyme functions and versatile PME-1 activities derived from coupling substrate-mimicking motifs to dynamic structured cores.

Comparative Analysis of Structural Features in SLiMs from Eukaryotes, Bacteria, and Viruses with Importance for Host-Pathogen Interactions

Protein-protein interactions drive functions in eukaryotes that can be described by short linear motifs (SLiMs). Conservation of SLiMs help illuminate functional SLiMs in eukaryotic protein families. However, the simplicity of eukaryotic SLiMs makes them appear by chance due to mutational processes not only in eukaryotes but also in pathogenic bacteria and viruses. Further, functional eukaryotic SLiMs are often found in disordered regions. Although proteomes from pathogenic bacteria and viruses have less disorder than eukaryotic proteomes, their proteins can successfully mimic eukaryotic SLiMs and disrupt host cellular function. Identifying important SLiMs in pathogens is difficult but essential for understanding potential host-pathogen interactions. We performed a comparative analysis of structural features for experimentally verified SLiMs from the Eukaryotic Linear Motif (ELM) database across viruses, bacteria, and eukaryotes. Our results revealed that many viral SLiMs and specific motifs found across viruses and eukaryotes, such as some glycosylation motifs, have less disorder. Analyzing the disorder and coil properties of equivalent SLiMs from pathogens and eukaryotes revealed that some motifs are more structured in pathogens than their eukaryotic counterparts and vice versa. These results support a varying mechanism of interaction between pathogens and their eukaryotic hosts for some of the same motifs.

Use of viral motif mimicry improves the proteome-wide discovery of human linear motifs

Linear motifs have an integral role in dynamic cell functions, including cell signaling. However, due to their small size, low complexity, and frequent mutations, identifying novel functional motifs poses a challenge. Viruses rely extensively on the molecular mimicry of cellular linear motifs. In this study, we apply systematic motif prediction combined with functional filters to identify human linear motifs convergently evolved also in viral proteins. We observe an increase in the sensitivity of motif prediction and improved enrichment in known instances. We identify >7,300 non-redundant motif instances at various confidence levels, 99 of which are supported by all functional and structural filters. Overall, we provide a pipeline to improve the identification of functional linear motifs from interactomics datasets and a comprehensive catalog of putative human motifs that can contribute to our understanding of the human domain-linear motif code and the associated mechanisms of viral interference.

Ehrlichia SLiM ligand mimetic activates Hedgehog signaling to engage a BCL-2 anti-apoptotic cellular program

Ehrlichia chaffeensis (E. chaffeensis) has evolved eukaryotic ligand mimicry to repurpose multiple cellular signaling pathways for immune evasion. In this investigation, we demonstrate that TRP120 has a novel repetitive short linear motif (SLiM) that activates the evolutionarily conserved Hedgehog (Hh) signaling pathway to inhibit apoptosis. In silico analysis revealed that TRP120 has sequence and functional similarity with Hh ligands and a candidate Hh ligand SLiM was identified. siRNA knockdown of Hh signaling and transcriptional components significantly reduced infection. Co-immunoprecipitation and surface plasmon resonance demonstrated that rTRP120-TR interacted directly with Hh receptor Patched-2 (PTCH2). E. chaffeensis infection resulted in early upregulation of Hh transcription factor GLI-1 and regulation of Hh target genes. Moreover, soluble recombinant TRP120 (rTRP120) activated Hh and induced gene expression consistent with the eukaryotic Hh ligand. The TRP120-Hh-SLiM (NPEVLIKD) induced nuclear translocation of GLI-1 in THP-1 cells and primary human monocytes and induced a rapid and expansive activation of Hh pathway target genes. Furthermore, Hh activation was blocked by an α-TRP120-Hh-SLiM antibody. TRP120-Hh-SLiM significantly increased levels of Hh target, anti-apoptotic protein B-cell lymphoma 2 (BCL-2), and siRNA knockdown of BCL-2 dramatically inhibited infection. Blocking Hh signaling with the inhibitor Vismodegib, induced a pro-apoptotic cellular program defined by decreased mitochondria membrane potential, significant reductions in BCL-2, activation of caspase 3 and 9, and increased apoptotic cells. This study reveals a novel E. chaffeensis SLiM ligand mimetic that activates Hh signaling to maintain E. chaffeensis infection by engaging a BCL-2 anti-apoptotic cellular program.