Attributes of short linear motifs

Computational Prediction of Linear Interacting Peptides

Intrinsically disordered protein regions, IDRs, are observed in many eukaryotic proteins. They play critical roles in essentially all cellular processes because segments of these regions, known as linear interacting peptides (LIPs), are heavily involved in regulatory protein interactions across proteomes. This chapter presents an integrated summary of the results from the last two Critical Assessments of protein Intrinsic Disorder predictions, known as CAID events, on the computational prediction of LIP segments. Because the CAID community questioned the quality of the test dataset used by the first CAID event, we reannotated this dataset using more accurate annotations from the latest DisProt database release. Then, we compared the results of the first CAID with the updated data and the results of the second CAID event. Our comparison highlights the importance of data annotation on the evaluation outcome and provides recommendations for users of LIP predictors.

Liquid-liquid phase separation in microorganisms: Insights into existence, functions, and applications

Liquid-liquid phase separation (LLPS) is a universal mechanism essential for maintaining cellular integrity and function in microorganisms, facilitating the organization of biomolecules into dynamic compartments. Although extensively studied in mammalian cells, research on LLPS formation and regulation in microorganisms remains limited. This review integrates insights from diverse studies exploring LLPS across microorganisms. We discuss the role of intrinsic disorders in microbial proteins and their relationship with environmental adaptation. Additionally, we examine how microorganisms utilize LLPS to sense changes in environmental parameters, such as temperature, pH, and nutrient levels, enabling them to respond to stresses and regulate cellular processes, such as cell division, protein synthesis, and metabolic flux. We highlight that LLPS is a promising target for synthetic biology and therapeutic intervention against pathogenic microorganisms. We also explore the research landscape of LLPS in microorganisms and address challenges associated with the techniques used in LLPS research. Further research is needed to focus on the detailed molecular regulatory mechanisms of condensates, biotechnological and synthetic biology applications, facilitating improved manipulation of microorganisms, and the identification of novel therapeutic targets.

Sp140L Is a Novel Herpesvirus Restriction Factor

Herpesviruses, including the oncogenic Epstein-Barr Virus (EBV), must bypass host DNA sensing mechanisms to establish infection. The first viral latency protein expressed, EBNA-LP, is essential for transformation of naïve B cells, yet its role in evading host defenses remains unclear. Using single-cell RNA sequencing of EBNA-LP-Knockout (LPKO)-infected B cells, we reveal an antiviral response landscape implicating the 'speckled proteins' as key restriction factors countered by EBNA-LP. Specifically, loss of SP100 or the primate-specific SP140L reverses the restriction of LPKO, suppresses a subset of canonically interferon-stimulated genes, and restores viral gene transcription and cellular proliferation. Notably, we also identify Sp140L as a restriction target of the herpesvirus saimiri ORF3 protein, implying a role in immunity to other DNA viruses. This study reveals Sp140L as a restriction factor that we propose links sensing and transcriptional suppression of viral DNA to an IFN-independent innate immune response, likely relevant to all nuclear DNA viruses.

Phosphorylation of P-stalk proteins defines the ribosomal state for interaction with auxiliary protein factors

Ribosomal action is facilitated by the orchestrated work of trans-acting factors and ribosomal elements, which are subject to regulatory events, often involving phosphorylation. One such element is the ribosomal P-stalk, which plays a dual function: it activates translational GTPases, which support basic ribosomal functions, and interacts with the Gcn2 kinase, linking the ribosomes to the ISR pathway. We show that P-stalk proteins, which form a pentamer, exist in the cell exclusively in a phosphorylated state at five C-terminal domains (CTDs), ensuring optimal translation (speed and accuracy) and may play a role in the timely regulation of the Gcn2-dependent stress response. Phosphorylation of the CTD induces a structural transition from a collapsed to a coil-like structure, and the CTD gains conformational freedom, allowing specific but transient binding to various protein partners, optimizing the ribosome action. The report reveals a unique feature of the P-stalk proteins, indicating that, unlike most ribosomal proteins, which are regulated by phosphorylation in an on/off manner, the P-stalk proteins exist in a constantly phosphorylated state, which optimizes their interaction with auxiliary factors.

Enhancing prediction of short linear protein motifs with Wregex 3.0

Short linear motifs (SLiMs) play an important role in protein-protein interactions. However, SLiM patterns are intrinsically permissive and result into many matches that occur just by chance, specially when targeting large datasets. To prioritize these matches as candidates for functional testing, we developed Wregex (Weighted regular expression), which uses a position-specific scoring matrix (PSSM) to order a list of regular expression matches according to a PSSM-derived score. Here we present Wregex 3.0, an improved version with new functionalities such as the support for a second auxiliary motif to help refining prediction of a primary SLiM, and post-translational modifications (PTMs) enrichment taking into account that many regulatory SLiM-mediated interactions are modulated by one or more PTMs. This version also incorporates a number of new features such as a convenient use of subproteomes, showing UniProt annotations such as disordered regions, searching for all known motifs and generating decoy databases for enrichment analysis. We provide case studies to illustrate how these new Wregex functionalities enhance prediction of short linear protein motifs. The Wregex 3.0 server is freely accessible at https://ehubio.ehu.eus/wregex3/.

A quantitative intracellular peptide binding assay reveals recognition determinants and context dependence of short linear motifs

Transient protein-protein interactions play key roles in controlling dynamic cellular responses. Many examples involve globular protein domains that bind to peptide sequences known as Short Linear Motifs (SLiMs), which are enriched in intrinsically disordered regions of proteins. Here we describe a novel functional assay for measuring SLiM binding, called Systematic Intracellular Motif Binding Analysis (SIMBA). In this method, binding of a foreign globular domain to its cognate SLiM peptide allows yeast cells to proliferate by blocking a growth arrest signal. A high-throughput application of the SIMBA method involving competitive growth and deep sequencing provides rapid quantification of the relative binding strength for thousands of SLiM sequence variants, and a comprehensive interrogation of SLiM sequence features that control their recognition and potency. We show that multiple distinct classes of SLiM-binding domains can be analyzed by this method, and that the relative binding strength of peptides in vivo correlates with their biochemical affinities measured in vitro. Deep mutational scanning provides high-resolution definitions of motif recognition determinants and reveals how sequence variations at non-core positions can modulate binding strength. Furthermore, mutational scanning of multiple parent peptides that bind human tankyrase ARC or YAP WW domains identifies distinct binding modes and uncovers context effects in which the preferred residues at one position depend on residues elsewhere. The findings establish SIMBA as a fast and incisive approach for interrogating SLiM recognition via massively parallel quantification of protein-peptide binding strength in vivo.

CompariPSSM: a PSSM-PSSM comparison tool for motif-binding determinant analysis

MOTIVATION

Short linear motifs (SLiMs) are compact functional modules that mediate low-affinity protein-protein interactions. SLiMs direct the function of many dynamic signalling and regulatory complexes playing a central role in most biological processes of the cell. Motif-binding determinants describe the contribution of each residue in a motif-containing peptide to the affinity and specificity of binding to the motif-binding partner. Motif-binding determinants are generally defined as a motif consensus pattern or a position-specific scoring matrix (PSSM) encoding quantitative preferences. Motif-binding determinant comparison is an important motif analysis task and can be applied to motif annotation, classification, clustering, discovery and benchmarking. Currently, binding determinant comparison is generally performed by analysing consensus similarity; however, this ignores important quantitative information in both the consensus and non-consensus positions.

RESULTS

We have created a new tool, CompariPSSM, that quantifies the similarity between motif-binding determinants using sliding window PSSM-PSSM comparison and scores PSSM similarity using a randomisation-based probabilistic framework. The tool has been benchmarked on curated data from the eukaryotic linear motif database and experimental data from proteomic peptidephage display. CompariPSSM can be used for peptide classification to validate motif classes, peptide clustering to group functionally related conserved disordered regions, and benchmarking experimental motif discovery methods.

AVAILABILITY AND IMPLEMENTATION

CompariPSSM is available at https://slim.icr.ac.uk/projects/comparipssm.

Asp/ASPM phospho-regulation throughout the cell cycle

In mammals and Drosophila melanogaster, Asp/ASPM proteins contribute to cell proliferation and spindle formation. Recent evidence also suggests interphase roles for Asp/ASPM proteins, but little is known about the regulation allowing distinct roles in different cell cycle phases. In this review, we consider a cross-species comparison of Asp/ASPM protein sequences in light of cyclin-CDK literature, and suggest Asp/ASPM proteins to be prime candidates for cyclin-CDK regulation. Conserved regulatory features include an N-terminal proline directed serine/threonine (S/T-P) "supershift" phosphorylation domain common to proteins with bistable interphase and mitotic roles, as well as putative cyclin-binding sites positioned to allow multisite phosphorylation by cyclin-CDK complexes. Human, mouse, and Drosophila Asp/ASPM protein structural predictions show that multisite phosphorylation of the N-term supershift domain could alter the availability of CH-domains and HEAT-motifs, which can contribute to microtubule binding and protein aggregation likely required for spindle formation. Structural predictions of the smallest reported microcephaly patient truncation also emphasize the importance of the arrangement of these motifs. We position this in silico analysis within recent literature to build new hypotheses for Asp/ASPM regulation in interphase and mitosis, as well as de-regulation in microcephaly and cancer. We also highlight the utility of comparing structural/functional differences between human ASPM and Drosophila Asp to gain further insight.

AlphaFold-Multimer accurately captures interactions and dynamics of intrinsically disordered protein regions

Interactions mediated by intrinsically disordered protein regions (IDRs) pose formidable challenges in structural characterization. IDRs are highly versatile, capable of adopting diverse structures and engagement modes. Motivated by recent strides in protein structure prediction, we embarked on exploring the extent to which AlphaFold-Multimer can faithfully reproduce the intricacies of interactions involving IDRs. To this end, we gathered multiple datasets covering the versatile spectrum of IDR binding modes and used them to probe AlphaFold-Multimer's prediction of IDR interactions and their dynamics. Our analyses revealed that AlphaFold-Multimer is not only capable of predicting various types of bound IDR structures with high success rate, but that distinguishing true interactions from decoys, and unreliable predictions from accurate ones is achievable by appropriate use of AlphaFold-Multimer's intrinsic scores. We found that the quality of predictions drops for more heterogeneous, fuzzy interaction types, most likely due to lower interface hydrophobicity and higher coil content. Notably though, certain AlphaFold-Multimer scores, such as the Predicted Aligned Error and residue-ipTM, are highly correlated with structural heterogeneity of the bound IDR, enabling clear distinctions between predictions of fuzzy and more homogeneous binding modes. Finally, our benchmarking revealed that predictions of IDR interactions can also be successful when using full-length proteins, but not as accurate as with cognate IDRs. To facilitate identification of the cognate IDR of a given partner, we established "minD," which pinpoints potential interaction sites in a full-length protein. Our study demonstrates that AlphaFold-Multimer can correctly identify interacting IDRs and predict their mode of engagement with a given partner.

The transition state for coupled folding and binding of a disordered DNA binding domain resembles the unbound state

The basic zippers (bZIPs) are one of two large eukaryotic families of transcription factors whose DNA binding domains are disordered in isolation but fold into stable α-helices upon target DNA binding. Here, we systematically disrupt pre-existing helical propensity within the DNA binding region of the homodimeric bZIP domain of cAMP-response element binding protein (CREB) using Ala-Gly scanning and examine the impact on target binding kinetics. We find that the secondary structure of the transition state strongly resembles that of the unbound state. The residue closest to the dimerization domain is largely folded within both unbound and transition states; dimerization apparently propagates additional helical propensity into the basic region. The results are consistent with electrostatically-enhanced DNA binding, followed by rapid folding from the folded zipper outwards. Fly-casting theory suggests that protein disorder can accelerate binding. Interestingly however, we did not observe higher association rate constants for mutants with lower levels of residual structure in the unbound state.

SHARK enables sensitive detection of evolutionary homologs and functional analogs in unalignable and disordered sequences

Intrinsically disordered regions (IDRs) are structurally flexible protein segments with regulatory functions in multiple contexts, such as in the assembly of biomolecular condensates. Since IDRs undergo more rapid evolution than ordered regions, identifying homology of such poorly conserved regions remains challenging for state-of-the-art alignment-based methods that rely on position-specific conservation of residues. Thus, systematic functional annotation and evolutionary analysis of IDRs have been limited, despite them comprising ~21% of proteins. To accurately assess homology between unalignable sequences, we developed an alignment-free sequence comparison algorithm, SHARK (Similarity/Homology Assessment by Relating K-mers). We trained SHARK-dive, a machine learning homology classifier, which achieved superior performance to standard alignment-based approaches in assessing evolutionary homology in unalignable sequences. Furthermore, it correctly identified dissimilar but functionally analogous IDRs in IDR-replacement experiments reported in the literature, whereas alignment-based tools were incapable of detecting such functional relationships. SHARK-dive not only predicts functionally similar IDRs at a proteome-wide scale but also identifies cryptic sequence properties and motifs that drive remote homology and analogy, thereby providing interpretable and experimentally verifiable hypotheses of the sequence determinants that underlie such relationships. SHARK-dive acts as an alternative to alignment to facilitate systematic analysis and functional annotation of the unalignable protein universe.

Prediction and design of transcriptional repressor domains with large-scale mutational scans and deep learning

Regulatory proteins have evolved diverse repressor domains (RDs) to enable precise context-specific repression of transcription. However, our understanding of how sequence variation impacts the functional activity of RDs is limited. To address this gap, we generated a high-throughput mutational scanning dataset measuring the repressor activity of 115,000 variant sequences spanning more than 50 RDs in human cells. We identified thousands of clinical variants with loss or gain of repressor function, including TWIST1 HLH variants associated with Saethre-Chotzen syndrome and MECP2 domain variants associated with Rett syndrome. We also leveraged these data to annotate short linear interacting motifs (SLiMs) that are critical for repression in disordered RDs. Then, we designed a deep learning model called TENet ( T ranscriptional E ffector Net work) that integrates sequence, structure and biochemical representations of sequence variants to accurately predict repressor activity. We systematically tested generalization within and across domains with varying homology using the mutational scanning dataset. Finally, we employed TENet within a directed evolution sequence editing framework to tune the activity of both structured and disordered RDs and experimentally test thousands of designs. Our work highlights critical considerations for future dataset design and model training strategies to improve functional variant prioritization and precision design of synthetic regulatory proteins.

Structural Insights into the Role of the Proline Rich Region in Tau Function

Tau is a microtubule-associated protein that plays an important role in modulating axonal microtubules in neurons. Intracellular tau aggregates are found in a broad class of disorders, including Alzheimer's disease, termed tauopathies. Tau is an intrinsically disordered protein, and its structural disorder appears to be critical to its microtubule-related functions. Tubulin binding sites are found in tau's proline-rich region (PRR), microtubule binding repeats (MTBR: R1-R4), and pseudo-repeat, R'. While many post-translational modifications have been identified on tau, phosphorylation sites, which both regulate tubulin dimer and microtubule interactions and are correlated with disease, cluster with high frequency within the PRR. Here, we use fluorescence correlation spectroscopy and structural mass spectrometry techniques to characterize the impact of phosphomimic mutations in the PRR on tubulin dimer binding and probe the structure of the PRR-tubulin dimer complex. We find that phosphomimics cumulatively diminish tubulin dimer binding and slow microtubule polymerization. Additionally, we map two ∼15 residue regions of the PRR as primary tubulin dimer binding sites and propose a model in which PRR enhances lateral interactions between tubulin dimers, complementing the longitudinal interactions observed for MTBR. Together these measurements provide insight into the previously overlooked relevance of tau's PRR in functional interactions with tubulin.

GRAPHICAL ABSTRACT

Rational Tuning of the Concentration-independent Enrichment of Prion-like Domains in Stress Granules

Stress granules (SGs) are large ribonucleoprotein assemblies that form in response to acute stress in eukaryotes. SG formation is thought to be initiated by liquid-liquid phase separation (LLPS) of key proteins and RNA. These molecules serve as a scaffold for recruitment of client molecules. LLPS of scaffold proteins in vitro is highly concentration-dependent, yet biomolecular condensates in vivo contain hundreds of unique proteins, most of which are thought to be clients rather than scaffolds. Many proteins that localize to SGs contain low-complexity, prion-like domains (PrLDs) that have been implicated in LLPS and SG recruitment. The degree of enrichment of proteins in biomolecular condensates such as SGs can vary widely, but the underlying basis for these differences is not fully understood. Here, we develop a toolkit of model PrLDs to examine the factors that govern efficiency of PrLD recruitment to stress granules. Recruitment was highly sensitive to amino acid composition: enrichment in SGs could be tuned through subtle changes in hydrophobicity. By contrast, SG recruitment was largely insensitive to PrLD concentration at both a population level and single-cell level. These observations point to a model wherein PrLDs are enriched in SGs through either simple solvation effects or interactions that are effectively non-saturable even at high expression levels.

Proteome-scale characterisation of motif-based interactome rewiring by disease mutations

Whole genome and exome sequencing are reporting on hundreds of thousands of missense mutations. Taking a pan-disease approach, we explored how mutations in intrinsically disordered regions (IDRs) break or generate protein interactions mediated by short linear motifs. We created a peptide-phage display library tiling ~57,000 peptides from the IDRs of the human proteome overlapping 12,301 single nucleotide variants associated with diverse phenotypes including cancer, metabolic diseases and neurological diseases. By screening 80 human proteins, we identified 366 mutation-modulated interactions, with half of the mutations diminishing binding, and half enhancing binding or creating novel interaction interfaces. The effects of the mutations were confirmed by affinity measurements. In cellular assays, the effects of motif-disruptive mutations were validated, including loss of a nuclear localisation signal in the cell division control protein CDC45 by a mutation associated with Meier-Gorlin syndrome. The study provides insights into how disease-associated mutations may perturb and rewire the motif-based interactome.

Facile detection of peptide-protein interactions using an electrophoretic crosslinking shift assay

Protein-protein interactions with high specificity and low affinity are functionally important but are not comprehensively understood because they are difficult to identify. Particularly intriguing are the dynamic and specific interactions between folded protein domains and short unstructured peptides known as short linear motifs. Such domain-motif interactions (DMIs) are often difficult to identify and study because affinities are modest to weak. Here we describe "electrophoretic crosslinking shift assay" (ECSA), a simple in vitro approach that detects transient, low affinity interactions by covalently crosslinking a prey protein and a fluorescently labeled bait. We demonstrate this technique on the well characterized DMI between MAP kinases and unstructured D-motif peptide ligands. We show that ECSA detects sequence-specific micromolar interactions using less than a microgram of input prey protein per reaction, making it ideal for verifying candidate low-affinity DMIs of components that purify with low yield. We propose ECSA as an intermediate step in SLiM characterization that bridges the gap between high throughput techniques such as phage display and more resource-intensive biophysical and structural analysis.

Design of intrinsically disordered protein variants with diverse structural properties

Intrinsically disordered proteins (IDPs) perform a broad range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Computational design of IDPs with specific conformational properties has, however, been difficult because of their substantial dynamics and structural complexity. We describe a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs that differ in compaction, long-range contacts, and propensity to phase separate. We experimentally tested and validated our designs and analyzed the sequence features that determine conformations. We show how our results are captured by a machine learning model, enabling us to speed up the algorithm. Our work expands the toolbox for computational protein design and will facilitate the design of proteins whose functions exploit the many properties afforded by protein disorder.

α-Crystalline Domains and Intrinsically Disordered Regions Can Work in Parallel to Induce Accumulation of MBD6 at Chromocenters in Arabidopsis thaliana

Proteins are localized and concentrated at cellular and genomic locations for specific and efficient functions. Efforts to understand protein accumulation in eukaryotic organisms have primarily focused on multivalent interactions between intrinsically disordered regions (IDRs) as mediators of protein condensation. We previously showed that α-crystalline domain (ACD) proteins 15 (ACD15) and 21 (ACD21) were required for multimerization and the accumulation of gene-silencing methyl-CpG-binding domain protein 6 (MBD6) at chromocenters in Arabidopsis thaliana. Here, we demonstrate that ACDs and IDRs can act as parallel mechanisms, facilitating higher-order MBD6 assemblies. Using human IDRs known to be important for protein accumulation, we replicated and enhanced the accumulation of MBD6 at chromocenters. In addition, IDRs fused to MBD6 could substitute for ACD function and partially reconstitute the MBD6 gene-silencing function. However, the accumulation of MBD6 by IDRs still required ACD15 and ACD21 for full effect. These results establish that ACD-mediated protein accumulation is a mechanism that can function similarly to and together with IDR-mediated mechanisms.

High-throughput affinity measurements of direct interactions between activation domains and co-activators

Sequence-specific activation by transcription factors is essential for gene regulation1,2. Key to this are activation domains, which often fall within disordered regions of transcription factors3,4 and recruit co-activators to initiate transcription5. These interactions are difficult to characterize via most experimental techniques because they are typically weak and transient6,7. Consequently, we know very little about whether these interactions are promiscuous or specific, the mechanisms of binding, and how these interactions tune the strength of gene activation. To address these questions, we developed a microfluidic platform for expression and purification of hundreds of activation domains in parallel followed by direct measurement of co-activator binding affinities (STAMMPPING, for Simultaneous Trapping of Affinity Measurements via a Microfluidic Protein-Protein INteraction Generator). By applying STAMMPPING to quantify direct interactions between eight co-activators and 204 human activation domains (>1,500 K ds), we provide the first quantitative map of these interactions and reveal 334 novel binding pairs. We find that the metazoan-specific co-activator P300 directly binds >100 activation domains, potentially explaining its widespread recruitment across the genome to influence transcriptional activation. Despite sharing similar molecular properties (e.g. enrichment of negative and hydrophobic residues), activation domains utilize distinct biophysical properties to recruit certain co-activator domains. Co-activator domain affinity and occupancy are well-predicted by analytical models that account for multivalency, and in vitro affinities quantitatively predict activation in cells with an ultrasensitive response. Not only do our results demonstrate the ability to measure affinities between even weak protein-protein interactions in high throughput, but they also provide a necessary resource of over 1,500 activation domain/co-activator affinities which lays the foundation for understanding the molecular basis of transcriptional activation.

Evolution of Virus-like Features and Intrinsically Disordered Regions in Retrotransposon-derived Mammalian Genes

Several mammalian genes have originated from the domestication of retrotransposons, selfish mobile elements related to retroviruses. Some of the proteins encoded by these genes have maintained virus-like features; including self-processing, capsid structure formation, and the generation of different isoforms through -1 programmed ribosomal frameshifting. Using quantitative approaches in molecular evolution and biophysical analyses, we studied 28 retrotransposon-derived genes, with a focus on the evolution of virus-like features. By analyzing the rate of synonymous substitutions, we show that the -1 programmed ribosomal frameshifting mechanism in three of these genes (PEG10, PNMA3, and PNMA5) is conserved across mammals and originates alternative proteins. These genes were targets of positive selection in primates, and one of the positively selected sites affects a B-cell epitope on the spike domain of the PNMA5 capsid, a finding reminiscent of observations in infectious viruses. More generally, we found that retrotransposon-derived proteins vary in their intrinsically disordered region content and this is directly associated with their evolutionary rates. Most positively selected sites in these proteins are located in intrinsically disordered regions and some of them impact protein posttranslational modifications, such as autocleavage and phosphorylation. Detailed analyses of the biophysical properties of intrinsically disordered regions showed that positive selection preferentially targeted regions with lower conformational entropy. Furthermore, positive selection introduces variation in binary sequence patterns across orthologues, as well as in chain compaction. Our results shed light on the evolutionary trajectories of a unique class of mammalian genes and suggest a novel approach to study how intrinsically disordered region biophysical characteristics are affected by evolution.

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.

Chromosome compaction is triggered by an autonomous DNA-binding module within condensin

The compaction of chromatin into mitotic chromosomes is essential for faithful transmission of the genome during cell division. In eukaryotes, chromosome morphogenesis is regulated by the condensin complex, though the exact mechanism used to target condensin to chromatin and initiate condensation is not understood. Here, we reveal that condensin contains an intrinsically disordered region (IDR) that modulates its association with chromatin in early mitosis and exhibits phase separation. We describe DNA-binding motifs within the IDR that, upon deletion, inflict striking defects in chromosome condensation and segregation, ill-timed condensin turnover on chromatin, and cell death. Importantly, we demonstrate that the condensin IDR can impart cell cycle regulatory functions when transferred to other subunits within the complex, indicating its autonomous nature. Collectively, our study unveils the molecular basis for the initiation of chromosome condensation in early mitosis and how this process ultimately promotes genomic stability and faultless cell division.

The substrate quality of CK2 target sites has a determinant role on their function and evolution

Most biological processes are regulated by signaling modules that bind to short linear motifs. For protein kinases, substrates may have full or only partial matches to the kinase recognition motif, a property known as "substrate quality." However, it is not clear whether differences in substrate quality represent neutral variation or if they have functional consequences. We examine this question for the kinase CK2, which has many fundamental functions. We show that optimal CK2 sites are phosphorylated at maximal stoichiometries and found in many conditions, whereas minimal substrates are more weakly phosphorylated and have regulatory functions. Optimal CK2 sites tend to be more conserved, and substrate quality is often tuned by selection. For intermediate sites, increases or decreases in substrate quality may be deleterious, as we demonstrate for a CK2 substrate at the kinetochore. The results together suggest a strong role for substrate quality in phosphosite function and evolution. A record of this paper's transparent peer review process is included in the supplemental information.

Intrinsic Disorder in the Host Proteins Entrapped in Rabies Virus Particles

A proteomics analysis of purified rabies virus (RABV) revealed 47 entrapped host proteins within the viral particles. Out of these, 11 proteins were highly disordered. Our study was particularly focused on five of the RABV-entrapped mouse proteins with the highest levels of disorder: Neuromodulin, Chmp4b, DnaJB6, Vps37B, and Wasl. We extensively utilized bioinformatics tools, such as FuzDrop, D2P2, UniProt, RIDAO, STRING, AlphaFold, and ELM, for a comprehensive analysis of the intrinsic disorder propensity of these proteins. Our analysis suggested that these disordered host proteins might play a significant role in facilitating the rabies virus pathogenicity, immune system evasion, and the development of antiviral drug resistance. Our study highlighted the complex interaction of the virus with its host, with a focus on how the intrinsic disorder can play a crucial role in virus pathogenic processes, and suggested that these intrinsically disordered proteins (IDPs) and disorder-related host interactions can also be a potential target for therapeutic strategies.

The bacterial DNA sliding clamp, β-clamp: structure, interactions, dynamics and drug discovery

DNA replication is a tightly coordinated event carried out by a multiprotein replication complex. An essential factor in the bacterial replication complex is the ring-shaped DNA sliding clamp, β-clamp, ensuring processive DNA replication and DNA repair through tethering of polymerases and DNA repair proteins to DNA. β -clamp is a hub protein with multiple interaction partners all binding through a conserved clamp binding sequence motif. Due to its central role as a DNA scaffold protein, β-clamp is an interesting target for antimicrobial drugs, yet little effort has been put into understanding the functional interactions of β-clamp. In this review, we scrutinize the β-clamp structure and dynamics, examine how its interactions with a plethora of binding partners are regulated through short linear binding motifs and discuss how contexts play into selection. We describe the dynamic process of clamp loading onto DNA and cover the recent advances in drug development targeting β-clamp. Despite decades of research in β-clamps and recent landmark structural insight, much remains undisclosed fostering an increased focus on this very central protein.

Disordered clock protein interactions and charge blocks turn an hourglass into a persistent circadian oscillator

Organismal physiology is widely regulated by the molecular circadian clock, a feedback loop composed of protein complexes whose members are enriched in intrinsically disordered regions. These regions can mediate protein-protein interactions via SLiMs, but the contribution of these disordered regions to clock protein interactions had not been elucidated. To determine the functionality of these disordered regions, we applied a synthetic peptide microarray approach to the disordered clock protein FRQ in Neurospora crassa. We identified residues required for FRQ's interaction with its partner protein FRH, the mutation of which demonstrated FRH is necessary for persistent clock oscillations but not repression of transcriptional activity. Additionally, the microarray demonstrated an enrichment of FRH binding to FRQ peptides with a net positive charge. We found that positively charged residues occurred in significant "blocks" within the amino acid sequence of FRQ and that ablation of one of these blocks affected both core clock timing and physiological clock output. Finally, we found positive charge clusters were a commonly shared molecular feature in repressive circadian clock proteins. Overall, our study suggests a mechanistic purpose for positive charge blocks and yielded insights into repressive arm protein roles in clock function.

A complete map of specificity encoding for a partially fuzzy protein interaction

Thousands of human proteins function by binding short linear motifs embedded in intrinsically disordered regions. How affinity and specificity are encoded in these binding domains and the motifs themselves is not well understood. The evolvability of binding specificity - how rapidly and extensively it can change upon mutation - is also largely unexplored, as is the contribution of 'fuzzy' dynamic residues to affinity and specificity in protein-protein interactions. Here we report the first complete map of specificity encoding for a globular protein domain. Quantifying >200,000 energetic interactions between a PDZ domain and its ligand identifies 20 major energetically coupled pairs of sites that control specificity. These are organized into six modules, with most mutations in each module reprogramming specificity for a single position in the ligand. Nine of the major energetic couplings controlling specificity are between structural contacts and 11 have an allosteric mechanism of action. The dynamic tail of the ligand is more robust to mutation than the structured residues but contributes additively to binding affinity and communicates with structured residues to enable changes in specificity. Our results quantify the binding specificities of >1,800 globular proteins to reveal how specificity is encoded and provide a direct comparison of the encoding of affinity and specificity in structured and dynamic molecular recognition.

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.

Pathogenic mutations of human phosphorylation sites affect protein-protein interactions

Despite their lack of a defined 3D structure, intrinsically disordered regions (IDRs) of proteins play important biological roles. Many IDRs contain short linear motifs (SLiMs) that mediate protein-protein interactions (PPIs), which can be regulated by post-translational modifications like phosphorylation. 20% of pathogenic missense mutations are found in IDRs, and understanding how such mutations affect PPIs is essential for unraveling disease mechanisms. Here, we employ peptide-based interaction proteomics to investigate 36 disease-associated mutations affecting phosphorylation sites. Our results unveil significant differences in interactomes between phosphorylated and non-phosphorylated peptides, often due to disrupted phosphorylation-dependent SLiMs. We focused on a mutation of a serine phosphorylation site in the transcription factor GATAD1, which causes dilated cardiomyopathy. We find that this phosphorylation site mediates interaction with 14-3-3 family proteins. Follow-up experiments reveal the structural basis of this interaction and suggest that 14-3-3 binding affects GATAD1 nucleocytoplasmic transport by masking a nuclear localisation signal. Our results demonstrate that pathogenic mutations of human phosphorylation sites can significantly impact protein-protein interactions, offering insights into potential molecular mechanisms underlying pathogenesis.

Proteasome activation: A novel strategy for targeting undruggable intrinsically disordered proteins

Intrinsically disordered proteins (IDPs) are characterized by their inability to adopt well-defined tertiary structures under physiological conditions. Nonetheless, they often play pivotal roles in the progression of various diseases, including cancer, neurodegenerative disorders, and cardiovascular ailments. Owing to their inherent dynamism, conventional drug design approaches based on structural considerations encounter substantial challenges when applied to IDPs. Consequently, the pursuit of therapeutic interventions directed towards IDPs presents a complex endeavor. While there are indeed existing methodologies for targeting IDPs, they are encumbered by noteworthy constrains. Hence, there exists an imminent imperative to investigate more efficacious and universally applicable strategies for modulating IDPs. Here, we present an overview of the latest advancements in the research pertaining to IDPs, along with the indirect regulation approach involving the modulation of IDP degradation through proteasome. By comprehending these advancements in research, novel insights can be generated to facilitate the development of new drugs targeted at addressing the accumulation of IDPs in diverse pathological conditions.

Uncovering domain motif interactions using high-throughput protein-protein interaction detection methods

Protein-protein interactions (PPIs) are often mediated by short linear motifs (SLiMs) in one protein and domain in another, known as domain-motif interactions (DMIs). During the past decade, SLiMs have been studied to find their role in cellular functions such as post-translational modifications, regulatory processes, protein scaffolding, cell cycle progression, cell adhesion, cell signalling and substrate selection for proteasomal degradation. This review provides a comprehensive overview of the current PPI detection techniques and resources, focusing on their relevance to capturing interactions mediated by SLiMs. We also address the challenges associated with capturing DMIs. Moreover, a case study analysing the BioGrid database as a source of DMI prediction revealed significant known DMI enrichment in different PPI detection methods. Overall, it can be said that current high-throughput PPI detection methods can be a reliable source for predicting DMIs.

Evolutionary analyses of intrinsically disordered regions reveal widespread signals of conservation

Intrinsically disordered regions (IDRs) are segments of proteins without stable three-dimensional structures. As this flexibility allows them to interact with diverse binding partners, IDRs play key roles in cell signaling and gene expression. Despite the prevalence and importance of IDRs in eukaryotic proteomes and various biological processes, associating them with specific molecular functions remains a significant challenge due to their high rates of sequence evolution. However, by comparing the observed values of various IDR-associated properties against those generated under a simulated model of evolution, a recent study found most IDRs across the entire yeast proteome contain conserved features. Furthermore, it showed clusters of IDRs with common "evolutionary signatures," i.e. patterns of conserved features, were associated with specific biological functions. To determine if similar patterns of conservation are found in the IDRs of other systems, in this work we applied a series of phylogenetic models to over 7,500 orthologous IDRs identified in the Drosophila genome to dissect the forces driving their evolution. By comparing models of constrained and unconstrained continuous trait evolution using the Brownian motion and Ornstein-Uhlenbeck models, respectively, we identified signals of widespread constraint, indicating conservation of distributed features is mechanism of IDR evolution common to multiple biological systems. In contrast to the previous study in yeast, however, we observed limited evidence of IDR clusters with specific biological functions, which suggests a more complex relationship between evolutionary constraints and function in the IDRs of multicellular organisms.

The disordered C-terminal tail of fungal LPMOs from phytopathogens mediates protein dimerization and impacts plant penetration

Lytic polysaccharide monooxygenases (LPMOs) are monocopper enzymes that oxidatively degrade various polysaccharides, such as cellulose. Despite extensive research on this class of enzymes, the role played by their C-terminal regions predicted to be intrinsically disordered (dCTR) has been overlooked. Here, we investigated the function of the dCTR of an LPMO, called CoAA9A, up-regulated during plant infection by Colletotrichum orbiculare, the causative agent of anthracnose. After recombinant production of the full-length protein, we found that the dCTR mediates CoAA9A dimerization in vitro, via a disulfide bridge, a hitherto-never-reported property that positively affects both binding and activity on cellulose. Using SAXS experiments, we show that the homodimer is in an extended conformation. In vivo, we demonstrate that gene deletion impairs formation of the infection-specialized cell called appressorium and delays penetration of the plant. Using immunochemistry, we show that the protein is a dimer not only in vitro but also in vivo when secreted by the appressorium. As these peculiar LPMOs are also found in other plant pathogens, our findings open up broad avenues for crop protection.

Fast and scalable querying of eukaryotic linear motifs with gget elm

MOTIVATION

Eukaryotic linear motifs (ELMs), or Short Linear Motifs, are protein interaction modules that play an essential role in cellular processes and signaling networks and are often involved in diseases like cancer. The ELM database is a collection of manually curated motif knowledge from scientific papers. It has become a crucial resource for investigating motif biology and recognizing candidate ELMs in novel amino acid sequences. Users can search amino acid sequences or UniProt Accessions on the ELM resource web interface. However, as with many web services, there are limitations in the swift processing of large-scale queries through the ELM web interface or API calls, and, therefore, integration into protein function analysis pipelines is limited.

RESULTS

To allow swift, large-scale motif analyses on protein sequences using ELMs curated in the ELM database, we have extended the gget suite of Python and command line tools with a new module, gget elm, which does not rely on the ELM server for efficiently finding candidate ELMs in user-submitted amino acid sequences and UniProt Accessions. gget elm increases accessibility to the information stored in the ELM database and allows scalable searches for motif-mediated interaction sites in the amino acid sequences.

AVAILABILITY AND IMPLEMENTATION

The manual and source code are available at https://github.com/pachterlab/gget.

Proteome-wide assessment of human interactome as a source of capturing domain-motif and domain-domain interactions

Protein-protein interactions (PPIs) play a crucial role in various biological processes by establishing domain-motif (DMI) and domain-domain interactions (DDIs). While the existence of real DMIs/DDIs is generally assumed, it is rarely tested; therefore, this study extensively compared high-throughput methods and public PPI repositories as sources for DMI and DDI prediction based on the assumption that the human interactome provides sufficient data for the reliable identification of DMIs and DDIs. Different datasets from leading high-throughput methods (Yeast two-hybrid [Y2H], Affinity Purification coupled Mass Spectrometry [AP-MS], and Co-fractionation-coupled Mass Spectrometry) were assessed for their ability to capture DMIs and DDIs using known DMI/DDI information. High-throughput methods were not notably worse than PPI databases and, in some cases, appeared better. In conclusion, all PPI datasets demonstrated significant enrichment in DMIs and DDIs (p-value <0.001), establishing Y2H and AP-MS as reliable methods for predicting these interactions. This study provides valuable insights for biologists in selecting appropriate methods for predicting DMIs, ultimately aiding in SLiM discovery.

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 pathogen-encoded signalling receptor Tir exploits host-like intrinsic disorder for infection

The translocated intimin receptor (Tir) is an essential type III secretion system (T3SS) effector of attaching and effacing pathogens contributing to the global foodborne disease burden. Tir acts as a cell-surface receptor in host cells, rewiring intracellular processes by targeting multiple host proteins. We investigated the molecular basis for Tir's binding diversity in signalling, finding that Tir is a disordered protein with host-like binding motifs. Unexpectedly, also are several other T3SS effectors. By an integrative approach, we reveal that Tir dimerises via an antiparallel OB-fold within a highly disordered N-terminal cytosolic domain. Also, it has a long disordered C-terminal cytosolic domain partially structured at host-like motifs that bind lipids. Membrane affinity depends on lipid composition and phosphorylation, highlighting a previously unrecognised host interaction impacting Tir-induced actin polymerisation and cell death. Furthermore, multi-site tyrosine phosphorylation enables Tir to engage host SH2 domains in a multivalent fuzzy complex, consistent with Tir's scaffolding role and binding promiscuity. Our findings provide insights into the intracellular Tir domains, highlighting the ability of T3SS effectors to exploit host-like protein disorder as a strategy for host evasion.

Systematic discovery of protein interaction interfaces using AlphaFold and experimental validation

Structural resolution of protein interactions enables mechanistic and functional studies as well as interpretation of disease variants. However, structural data is still missing for most protein interactions because we lack computational and experimental tools at scale. This is particularly true for interactions mediated by short linear motifs occurring in disordered regions of proteins. We find that AlphaFold-Multimer predicts with high sensitivity but limited specificity structures of domain-motif interactions when using small protein fragments as input. Sensitivity decreased substantially when using long protein fragments or full length proteins. We delineated a protein fragmentation strategy particularly suited for the prediction of domain-motif interfaces and applied it to interactions between human proteins associated with neurodevelopmental disorders. This enabled the prediction of highly confident and likely disease-related novel interfaces, which we further experimentally corroborated for FBXO23-STX1B, STX1B-VAMP2, ESRRG-PSMC5, PEX3-PEX19, PEX3-PEX16, and SNRPB-GIGYF1 providing novel molecular insights for diverse biological processes. Our work highlights exciting perspectives, but also reveals clear limitations and the need for future developments to maximize the power of Alphafold-Multimer for interface predictions.

A cyclin D1 intrinsically disordered domain accesses modified histone motifs to govern gene transcription

The essential G1-cyclin, CCND1, is frequently overexpressed in cancer, contributing to tumorigenesis by driving cell-cycle progression. D-type cyclins are rate-limiting regulators of G1-S progression in mammalian cells via their ability to bind and activate CDK4 and CDK6. In addition, cyclin D1 conveys kinase-independent transcriptional functions of cyclin D1. Here we report that cyclin D1 associates with H2BS14 via an intrinsically disordered domain (IDD). The same region of cyclin D1 was necessary for the induction of aneuploidy, induction of the DNA damage response, cyclin D1-mediated recruitment into chromatin, and CIN gene transcription. In response to DNA damage H2BS14 phosphorylation occurs, resulting in co-localization with γH2AX in DNA damage foci. Cyclin D1 ChIP seq and γH2AX ChIP seq revealed ~14% overlap. As the cyclin D1 IDD functioned independently of the CDK activity to drive CIN, the IDD domain may provide a rationale new target to complement CDK-extinction strategies.

Protein intrinsically disordered regions have a non-random, modular architecture

MOTIVATION

Protein sequences can be broadly categorized into two classes: those which adopt stable secondary structure and fold into a domain (i.e. globular proteins), and those that do not. The sequences belonging to this latter class are conformationally heterogeneous and are described as being intrinsically disordered. Decades of investigation into the structure and function of globular proteins has resulted in a suite of computational tools that enable their sub-classification by domain type, an approach that has revolutionized how we understand and predict protein functionality. Conversely, it is unknown if sequences of disordered protein regions are subject to broadly generalizable organizational principles that would enable their sub-classification.

RESULTS

Here, we report the development of a statistical approach that quantifies linear variance in amino acid composition across a sequence. With multiple examples, we provide evidence that intrinsically disordered regions are organized into statistically non-random modules of unique compositional bias. Modularity is observed for both low and high-complexity sequences and, in some cases, we find that modules are organized in repetitive patterns. These data demonstrate that disordered sequences are non-randomly organized into modular architectures and motivate future experiments to comprehensively classify module types and to determine the degree to which modules constitute functionally separable units analogous to the domains of globular proteins.

AVAILABILITY AND IMPLEMENTATION

The source code, documentation, and data to reproduce all figures are freely available at https://github.com/MWPlabUTSW/Chi-Score-Analysis.git. The analysis is also available as a Google Colab Notebook (https://colab.research.google.com/github/MWPlabUTSW/Chi-Score-Analysis/blob/main/ChiScore_Analysis.ipynb).

The function of alternative splicing in the proteome: rewiring protein interactomes to put old functions into new contexts

Alternative splicing affects more than 95% of multi-exon genes in the human genome. These changes affect the proteome in a myriad of ways. Here, we review our understanding of the breadth of these changes from their effect on protein structure to their influence on interactions. These changes encompass effects on nucleic acid binding in the nucleus to protein-carbohydrate interactions in the extracellular milieu, altering interactions involving all major classes of biological molecules. Protein isoforms have profound influences on cellular and tissue physiology, for example, by shaping neuronal connections, enhancing insulin secretion by pancreatic beta cells and allowing for alternative viral defense strategies in stem cells. More broadly, alternative splicing enables repurposing proteins from one context to another and thereby contributes to both the evolution of new traits as well as the creation of disease-specific interactomes that drive pathological phenotypes. In this Review, we highlight this universal character of alternative splicing as a central regulator of protein function with implications for almost every biological process.

Glutamine-rich regions of the disordered CREB transactivation domain mediate dynamic intra- and intermolecular interactions

The cyclic AMP response element (CRE) binding protein (CREB) is a transcription factor that contains a 280-residue N-terminal transactivation domain and a basic leucine zipper that mediates interaction with DNA. The transactivation domain comprises three subdomains, the glutamine-rich domains Q1 and Q2 and the kinase inducible activation domain (KID). NMR chemical shifts show that the isolated subdomains are intrinsically disordered but have a propensity to populate local elements of secondary structure. The Q1 and Q2 domains exhibit a propensity for formation of short β-hairpin motifs that function as binding sites for glutamine-rich sequences. These motifs mediate intramolecular interactions between the CREB Q1 and Q2 domains as well as intermolecular interactions with the glutamine-rich Q1 domain of the TATA-box binding protein associated factor 4 (TAF4) subunit of transcription factor IID (TFIID). Using small-angle X-ray scattering, NMR, and single-molecule Förster resonance energy transfer, we show that the Q1, Q2, and KID regions remain dynamically disordered in a full-length CREB transactivation domain (CREBTAD) construct. The CREBTAD polypeptide chain is largely extended although some compaction is evident in the KID and Q2 domains. Paramagnetic relaxation enhancement reveals transient long-range contacts both within and between the Q1 and Q2 domains while the intervening KID domain is largely devoid of intramolecular interactions. Phosphorylation results in expansion of the KID domain, presumably making it more accessible for binding the CBP/p300 transcriptional coactivators. Our study reveals the complex nature of the interactions within the intrinsically disordered transactivation domain of CREB and provides molecular-level insights into dynamic and transient interactions mediated by the glutamine-rich domains.

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.

How different viruses perturb host cellular machinery via short linear motifs

The virus interacts with its hosts by developing protein-protein interactions. Most viruses employ protein interactions to imitate the host protein: A viral protein with the same amino acid sequence or structure as the host protein attaches to the host protein's binding partner and interferes with the host protein's pathways. Being opportunistic, viruses have evolved to manipulate host cellular mechanisms by mimicking short linear motifs. In this review, we shed light on the current understanding of mimicry via short linear motifs and focus on viral mimicry by genetically different viral subtypes by providing recent examples of mimicry evidence and how high-throughput methods can be a reliable source to study SLiM-mediated viral mimicry.

Design of intrinsically disordered protein variants with diverse structural properties

Intrinsically disordered proteins (IDPs) perform a wide range of functions in biology, suggesting that the ability to design IDPs could help expand the repertoire of proteins with novel functions. Designing IDPs with specific structural or functional properties has, however, been difficult, in part because determining accurate conformational ensembles of IDPs generally requires a combination of computational modelling and experiments. Motivated by recent advancements in efficient physics-based models for simulations of IDPs, we have developed a general algorithm for designing IDPs with specific structural properties. We demonstrate the power of the algorithm by generating variants of naturally occurring IDPs with different levels of compaction and that vary more than 100 fold in their propensity to undergo phase separation, even while keeping a fixed amino acid composition. We experimentally tested designs of variants of the low-complexity domain of hnRNPA1 and find high accuracy in our computational predictions, both in terms of single-chain compaction and propensity to undergo phase separation. We analyze the sequence features that determine changes in compaction and propensity to phase separate and find an overall good agreement with previous findings for naturally occurring sequences. Our general, physics-based method enables the design of disordered sequences with specified conformational properties. Our algorithm thus expands the toolbox for protein design to include also the most flexible proteins and will enable the design of proteins whose functions exploit the many properties afforded by protein disorder.

Recognition and coacervation of G-quadruplexes by a multifunctional disordered region in RECQ4 helicase

Biomolecular polyelectrolyte complexes can be formed between oppositely charged intrinsically disordered regions (IDRs) of proteins or between IDRs and nucleic acids. Highly charged IDRs are abundant in the nucleus, yet few have been functionally characterized. Here, we show that a positively charged IDR within the human ATP-dependent DNA helicase Q4 (RECQ4) forms coacervates with G-quadruplexes (G4s). We describe a three-step model of charge-driven coacervation by integrating equilibrium and kinetic binding data in a global numerical model. The oppositely charged IDR and G4 molecules form a complex in the solution that follows a rapid nucleation-growth mechanism leading to a dynamic equilibrium between dilute and condensed phases. We also discover a physical interaction with Replication Protein A (RPA) and demonstrate that the IDR can switch between the two extremes of the structural continuum of complexes. The structural, kinetic, and thermodynamic profile of its interactions revealed a dynamic disordered complex with nucleic acids and a static ordered complex with RPA protein. The two mutually exclusive binding modes suggest a regulatory role for the IDR in RECQ4 function by enabling molecular handoffs. Our study extends the functional repertoire of IDRs and demonstrates a role of polyelectrolyte complexes involved in G4 binding.

Intracellular peptides in SARS-CoV-2-infected patients

Intracellular peptides (InPeps) generated by the orchestrated action of the proteasome and intracellular peptidases have biological and pharmacological significance. Here, human plasma relative concentration of specific InPeps was compared between 175 patients infected with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and 45 SARS-CoV-2 non-infected patients; 2,466 unique peptides were identified, of which 67% were InPeps. The results revealed differences of a specific group of peptides in human plasma comparing non-infected individuals to patients infected by SARS-CoV-2, following the results of the semi-quantitative analyses by isotope-labeled electrospray mass spectrometry. The protein-protein interactions networks enriched pathways, drawn by genes encoding the proteins from which the peptides originated, revealed the presence of the coronavirus disease/COVID-19 network solely in the group of patients fatally infected by SARS-CoV-2. Thus, modulation of the relative plasma levels of specific InPeps could be employed as a predictive tool for disease outcome.

Annexin A1, A2, A5, and A6 involvement in human pathologies

The human genome codes for 12 annexins with highly homologous membrane-binding cores and unique amino termini, which endow each protein with its specific biological properties. Not unique to vertebrate biology, multiple annexin orthologs are present in almost all eukaryotes. Their ability to combine either dynamically or constitutively with membrane lipid bilayers is hypothetically the key property that has led to their retention and multiple adaptation in eukaryotic molecular cell biology. Annexin genes are differentially expressed in many cell types but their disparate functions are still being discovered after more than 40 years of international research. A picture is emerging from gene knock down and knock out studies of individual annexins that these are important supporters rather than critical players in organism development and normal cell and tissue function. However, they appear to be highly significant "early responders" toward challenges arising from cell and tissue abiotic or biotic stress. In humans, recent focus has been on involvement of the annexin family for its involvement in diverse pathologies, especially cancer. From what has become an exceedingly broad field of investigation, we have selected four annexins in particular: AnxA1, 2, 5, and 6. Present both within and external to cells, these annexins are currently under intensive investigation in translational research as biomarkers of cellular dysfunction and as potential therapeutic targets for inflammatory conditions, neoplasia, and tissue repair. Annexin expression and release in response to biotic stress appears to be a balancing act. Under- or over-expression in different circumstances appears to damage rather than restore a healthy homeostasis. This review reflects briefly on what is already known of the structures and molecular cell biology of these selected annexins and considers their actual and potential roles in human health and disease.

Conserved molecular recognition by an intrinsically disordered region in the absence of sequence conservation

Intrinsically disordered regions (IDRs) are critical for cellular function, yet often appear to lack sequence conservation when assessed by multiple sequence alignments. This raises the question of if and how function can be encoded and preserved in these regions despite massive sequence variation. To address this question, we have applied coarse-grained molecular dynamics simulations to investigate non-specific RNA binding of coronavirus nucleocapsid proteins. Coronavirus nucleocapsid proteins consist of multiple interspersed disordered and folded domains that bind RNA. We focussed here on the first two domains of coronavirus nucleocapsid proteins, the disordered N-terminal domain (NTD) followed by the folded RNA binding domain (RBD). While the NTD is highly variable across evolution, the RBD is structurally conserved. This combination makes the NTD-RBD a convenient model system to explore the interplay between an IDR adjacent to a folded domain, and how changes in IDR sequence can influence molecular recognition of a partner. Our results reveal a surprising degree of sequence-specificity encoded by both the composition and the precise order of the amino acids in the NTD. The presence of an NTD can - depending on the sequence - either suppress or enhance RNA binding. Despite this sensitivity, large-scale variation in NTD sequences is possible while certain sequence features are retained. Consequently, a conformationally-conserved fuzzy RNA:protein complex is found across nucleocapsid protein orthologs, despite large-scale changes in both NTD sequence and RBD surface chemistry. Taken together, these insights shed light on the ability of disordered regions to preserve functional characteristics despite their sequence variability.

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.