Linear motifs: evolutionary interaction switches

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.

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.

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.

Homeodomain complex formation and biomolecular condensates in Hox gene regulation

Hox genes are a family of homeodomain transcription factors that regulate specialized morphological structures along the anterior-posterior axis of metazoans. Over the past few decades, researchers have focused on defining how Hox factors with similar in vitro DNA binding activities achieve sufficient target specificity to regulate distinct cell fates in vivo. In this review, we highlight how protein interactions with other transcription factors, many of which are also homeodomain proteins, result in the formation of transcription factor complexes with enhanced DNA binding specificity. These findings suggest that Hox-regulated enhancers utilize distinct combinations of homeodomain binding sites, many of which are low-affinity, to recruit specific Hox complexes. However, low-affinity sites can only yield reproducible responses with high transcription factor concentrations. To overcome this limitation, recent studies revealed how transcription factors, including Hox factors, use intrinsically disordered domains (IDRs) to form biomolecular condensates that increase protein concentrations. Moreover, Hox factors with altered IDRs have been associated with altered transcriptional activity and human disease states, demonstrating the importance of IDRs in mediating essential Hox output. Collectively, these studies highlight how Hox factors use their DNA binding domains, protein-protein interaction domains, and IDRs to form specific transcription factor complexes that yield accurate gene expression.

The interaction network of the proteasome assembly chaperone PSMD9 regulates proteostasis

Functional networks in cells are created by physical, genetic, and regulatory interactions. Mapping them and annotating their functions by available methods remains a challenge. We use affinity purification mass spectrometry (AP-MS) coupled with SLiMFinder to discern such a network involving 26S proteasome non-ATPase regulatory subunit 9 (PSMD9), a chaperone of proteasome assembly. Approximately 20% of proteins within the PSMD9 interactome carry a short linear motif (SLiM) of the type 'EXKK'. The binding of purified PSMD9 with the peptide sequence ERKK, proteins heterogeneous nuclear ribonucleoproteins A2/B1 (hnRNPA2B1; containing ERKK), and peroxiredoxin-6 (PRDX6; containing EAKK) provided proof of principle for this motif-driven network. The EXKK motif in the peptide primarily interacts with the coiled-coil N domain of PSMD9, a unique interaction not reported for any coiled-coil domain. PSMD9 knockout (KO) HEK293 cells experience endoplasmic reticulum (ER) stress and respond by increasing the unfolded protein response (UPR) and reducing the formation of aggresomes and lipid droplets. Trans-expression of PSMD9 in the KO cells rescues lipid droplet formation. Overexpression of PSMD9 in HEK293 cells results in reduced UPR, and increased lipid droplet and aggresome formation. The outcome argues for the prominent role of PSMD9 in maintaining proteostasis. Probable mechanisms involve the binding of PSMD9 to binding immunoglobulin protein (BIP/GRP78; containing EDKK), an endoplasmic reticulum chaperone and key regulator of the UPR, and fatty acid synthase (FASN; containing ELKK), involved in fatty acid synthesis/lipid biogenesis. We propose that PSMD9 acts as a buffer in the cellular milieu by moderating the UPR and enhancing aggresome formation to reduce stress-induced proteotoxicity. Akin to waves created in ponds that perpetuate to a distance, perturbing the levels of PSMD9 would cause ripples down the networks, affecting final reactions in the pathway, one of which is altered proteostasis.

Is there a loophole in Dollo's law? A DevoEvo perspective on irreversibility (of felid dentition)

There is a longstanding interest in whether the loss of complex characters is reversible (so-called "Dollo's law"). Reevolution has been suggested for numerous traits but among the first was Kurtén, who proposed that the presence of the second lower molar (M2 ) of the Eurasian lynx (Lynx lynx) was a violation of Dollo's law because all other Felids lack M2 . While an early and often cited example for the reevolution of a complex trait, Kurtén and Werdelin used an ad hoc parsimony argument to support their. Here I revisit the evidence that M2 reevolved lynx using explicit parsimony and maximum likelihood models of character evolution and find strong evidence that Kurtén and Werdelin were correct-M2 reevolved in E. lynx. Next, I explore the developmental mechanisms which may explain this violation of Dollo's law and suggest that the reevolution of lost complex traits may arise from the reevolution of cis-regulatory elements and protein-protein interactions, which have a longer half-life after silencing that protein coding genes. Finally, I present a developmental model to explain the reevolution M2 in E. lynx, which suggest that the developmental programs required for the establishment of serially homologous characters may never really be lost so long as a single instance of the character remains-thus the gain and loss and regain of serially homologous characters, such mammalian molars, may be developmentally and evolutionarily "simple."

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.

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.

Conformational Manifold Sampled by Two Short Linear Motif Segments Probed by Circular Dichroism, Vibrational, and Nuclear Magnetic Resonance Spectroscopy

Disordered protein segments called short linear motifs (SLiM) serve as recognition sites for a variety of biological processes and act as targeting signals, modification, and ligand binding sites. While SLiMs do not adopt one of the known regular secondary structures, the conformational distribution might still reflect the structural propensities of their amino acid residues and possible interactions between them. In the past, conformational analyses of short peptides provided compelling evidence for the notion that individual residues are less conformationally flexible than locally expected for a random coil. Here, we combined various spectroscopies (NMR, IR, vibrational, and UV circular dichroism) to determine the Ramachandran plots of two SLiM motifs, i.e., GRRDSG and GRRTSG. They are two representatives of RxxS motifs that are capable of being phosphorylated by protein kinase A, an enzyme that plays a fundamental role in a variety of biological processes. Our results reveal that the nearest and non-nearest interactions between residues cause redistributions between polyproline II and β-strand basins while concomitantly stabilizing extended relative to turn-forming and helical structures. They also cause shifts in basin positions. With increasing temperature, β-strand populations become more populated at the expense of polyproline II. While molecular dynamics simulations with Amber ff14SB and CHARMM 36m force fields indicate residue-residue interactions, they do not account for the observed structural changes.

Diversity of Short Linear Interaction Motifs in SARS-CoV-2 Nucleocapsid Protein

Molecular mimicry of short linear interaction motifs has emerged as a key mechanism for viral proteins binding host domains and hijacking host cell processes. Here, we examine the role of RNA-virus sequence diversity in the dynamics of the virus-host interface, by analyzing the uniquely vast sequence record of viable SARS-CoV-2 species with focus on the multi-functional nucleocapsid protein. We observe the abundant presentation of motifs encoding several essential host protein interactions, alongside a majority of possibly non-functional and randomly occurring motif sequences absent in subsets of viable virus species. A large number of motifs emerge ex nihilo through transient mutations relative to the ancestral consensus sequence. The observed mutational landscape implies an accessible motif space that spans at least 25% of known eukaryotic motifs. This reveals motif mimicry as a highly dynamic process with the capacity to broadly explore host motifs, allowing the virus to rapidly evolve the virus-host interface.

A conserved motif in the disordered linker of human MLH1 is vital for DNA mismatch repair and its function is diminished by a cancer family mutation

DNA mismatch repair (MMR) is essential for correction of DNA replication errors. Germline mutations of the human MMR gene MLH1 are the major cause of Lynch syndrome, a heritable cancer predisposition. In the MLH1 protein, a non-conserved, intrinsically disordered region connects two conserved, catalytically active structured domains of MLH1. This region has as yet been regarded as a flexible spacer, and missense alterations in this region have been considered non-pathogenic. However, we have identified and investigated a small motif (ConMot) in this linker which is conserved in eukaryotes. Deletion of the ConMot or scrambling of the motif abolished mismatch repair activity. A mutation from a cancer family within the motif (p.Arg385Pro) also inactivated MMR, suggesting that ConMot alterations can be causative for Lynch syndrome. Intriguingly, the mismatch repair defect of the ConMot variants could be restored by addition of a ConMot peptide containing the deleted sequence. This is the first instance of a DNA mismatch repair defect conferred by a mutation that can be overcome by addition of a small molecule. Based on the experimental data and AlphaFold2 predictions, we suggest that the ConMot may bind close to the C-terminal MLH1-PMS2 endonuclease and modulate its activation during the MMR process.

The relevance of short peptides for an understanding of unfolded and intrinsically disordered proteins

Over the last thirty years the unfolded state of proteins has attracted considerable interest owing to the discovery of intrinsically disordered proteins which perform a plethora of functions despite resembling unfolded proteins to a significant extent. Research on both, unfolded and disordered proteins has revealed that their conformational properties can deviate locally from random coil behavior. In this context results from work on short oligopeptides suggest that individual amino acid residues sample the sterically allowed fraction of the Ramachandran plot to a different extent. Alanine has been found to exhibit a peculiarity in that it has a very high propensity for adopting polyproline II like conformations. This Perspectives article reviews work on short peptides aimed at exploring the Ramachandran distributions of amino acid residues in different contexts with experimental and computational means. Based on the thus provided overview the article discussed to what extent short peptides can serve as tools for exploring unfolded and disordered proteins and as benchmarks for the development of a molecular dynamics force field.

Computational prediction of disordered binding regions

One of the key features of intrinsically disordered regions (IDRs) is their ability to interact with a broad range of partner molecules. Multiple types of interacting IDRs were identified including molecular recognition fragments (MoRFs), short linear sequence motifs (SLiMs), and protein-, nucleic acids- and lipid-binding regions. Prediction of binding IDRs in protein sequences is gaining momentum in recent years. We survey 38 predictors of binding IDRs that target interactions with a diverse set of partners, such as peptides, proteins, RNA, DNA and lipids. We offer a historical perspective and highlight key events that fueled efforts to develop these methods. These tools rely on a diverse range of predictive architectures that include scoring functions, regular expressions, traditional and deep machine learning and meta-models. Recent efforts focus on the development of deep neural network-based architectures and extending coverage to RNA, DNA and lipid-binding IDRs. We analyze availability of these methods and show that providing implementations and webservers results in much higher rates of citations/use. We also make several recommendations to take advantage of modern deep network architectures, develop tools that bundle predictions of multiple and different types of binding IDRs, and work on algorithms that model structures of the resulting complexes.

Problem of Domain/Building Block Preservation in the Evolution of Biological Macromolecules and Evolutionary Computation

Structurally and functionally isolated domains in biological macromolecular evolution, both natural and artificial, are largely similar to "schemata", building blocks (BBs), in evolutionary computation (EC). The problem of preserving in subsequent evolutionary searches the already found domains / BBs is well known and quite relevant in biology as well as in EC. Both biology and EC are seeing parallel and independent development of several approaches to identifying and preserving previously identified domains / BBs. First, we notice the similarity of DNA shuffling methods in synthetic biology and multi-parent recombination algorithms in EC. Furthermore, approaches to computer identification of domains in proteins that are being developed in biology can be aligned with BB identification methods in EC. Finally, approaches to chimeric protein libraries optimization in biology can be compared to evolutionary search methods based on probabilistic models in EC. We propose to validate the prospects of mutual exchange of ideas and transfer of algorithms and approaches between evolutionary systems biology and EC in these three principal directions. A crucial aim of this transfer is the design of new advanced experimental techniques capable of solving more complex problems of in vitro evolution.

Targeting peptide-mediated interactions in omics

Peptide-mediated interactions (PMIs) play a crucial role in cell signaling network, which are responsible for about half of cellular protein-protein associations in the human interactome and have recently been recognized as a new kind of promising druggable target for drug development and disease therapy. In this article, we give a systematic review regarding the proteome-wide discovery of PMIs and targeting druggable PMIs (dPMIs) with chemical drugs, self-inhibitory peptides (SIPs) and protein agents, particularly focusing on their implications and applications for therapeutic purpose in omics. We also introduce computational peptidology strategies used to model, analyze, and design PMI-targeted molecular entities and further extend the concepts of protein context, direct/indirect readout, and enthalpy/entropy effect involved in PMIs. Current issues and future perspective on this topic are discussed. There is still a long way to go before establishment of efficient therapeutic strategies to target PMIs on the omics scale.

CLIP: accurate prediction of disordered linear interacting peptides from protein sequences using co-evolutionary information

One of key features of intrinsically disordered regions (IDRs) is facilitation of protein-protein and protein-nucleic acids interactions. These disordered binding regions include molecular recognition features (MoRFs), short linear motifs (SLiMs) and longer binding domains. Vast majority of current predictors of disordered binding regions target MoRFs, with a handful of methods that predict SLiMs and disordered protein-binding domains. A new and broader class of disordered binding regions, linear interacting peptides (LIPs), was introduced recently and applied in the MobiDB resource. LIPs are segments in protein sequences that undergo disorder-to-order transition upon binding to a protein or a nucleic acid, and they cover MoRFs, SLiMs and disordered protein-binding domains. Although current predictors of MoRFs and disordered protein-binding regions could be used to identify some LIPs, there are no dedicated sequence-based predictors of LIPs. To this end, we introduce CLIP, a new predictor of LIPs that utilizes robust logistic regression model to combine three complementary types of inputs: co-evolutionary information derived from multiple sequence alignments, physicochemical profiles and disorder predictions. Ablation analysis suggests that the co-evolutionary information is particularly useful for this prediction and that combining the three inputs provides substantial improvements when compared to using these inputs individually. Comparative empirical assessments using low-similarity test datasets reveal that CLIP secures area under receiver operating characteristic curve (AUC) of 0.8 and substantially improves over the results produced by the closest current tools that predict MoRFs and disordered protein-binding regions. The webserver of CLIP is freely available at http://biomine.cs.vcu.edu/servers/CLIP/ and the standalone code can be downloaded from http://yanglab.qd.sdu.edu.cn/download/CLIP/.

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.

Portability of a Small-Molecule Binding Site between Disordered Proteins

Intrinsically disordered proteins (IDPs) are important in both normal and disease states. Small molecules can be targeted to disordered regions, but we currently have only a limited understanding of the nature of small-molecule binding sites in IDPs. Here, we show that a minimal small-molecule binding sequence of eight contiguous residues derived from the Myc protein can be ported into a different disordered protein and recapitulate small-molecule binding activity in the new context. We also find that the residue immediately flanking the binding site can have opposing effects on small-molecule binding in the different disordered protein contexts. The results demonstrate that small-molecule binding sites can act modularly and are portable between disordered protein contexts but that residues outside of the minimal binding site can modulate binding affinity.

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.

Soft disorder modulates the assembly path of protein complexes

The relationship between interactions, flexibility and disorder in proteins has been explored from many angles over the years: folding upon binding, flexibility of the core relative to the periphery, entropy changes, etc. In this work, we provide statistical evidence for the involvement of highly mobile and disordered regions in complex assembly. We ordered the entire set of X-ray crystallographic structures in the Protein Data Bank into hierarchies of progressive interactions involving identical or very similar protein chains, yielding 40205 hierarchies of protein complexes with increasing numbers of partners. We then examine them as proxies for the assembly pathways. Using this database, we show that upon oligomerisation, the new interfaces tend to be observed at residues that were characterised as softly disordered (flexible, amorphous or missing residues) in the complexes preceding them in the hierarchy. We also rule out the possibility that this correlation is just a surface effect by restricting the analysis to residues on the surface of the complexes. Interestingly, we find that the location of soft disordered residues in the sequence changes as the number of partners increases. Our results show that there is a general mechanism for protein assembly that involves soft disorder and modulates the way protein complexes are assembled. This work highlights the difficulty of predicting the structure of large protein complexes from sequence and emphasises the importance of linking predictors of soft disorder to the next generation of predictors of complex structure. Finally, we investigate the relationship between the Alphafold2’s confidence metric pLDDT for structure prediction in unbound versus bound structures, and soft disorder. We show a strong correlation between Alphafold2 low confidence residues and the union of all regions of soft disorder observed in the hierarchy. This paves the way for using the pLDDT metric as a proxy for predicting interfaces and assembly paths.

Intrinsic protein disorder uncouples affinity from binding specificity

Intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs) of proteins often function by molecular recognition, in which they undergo induced folding. Based on prior generalizations, the idea prevails in the IDP field that due to the entropic penalty of induced folding, the major functional advantage associated with this binding mode is "uncoupling" specificity from binding strength. Nevertheless, both weaker binding and high specificity of IDPs/IDRs rest on limited experimental observations, making these assumptions more speculations than evidence-supported facts. The issue is also complicated by the rather vague concept of specificity that lacks an exact measure, such as the K_d for binding strength. We addressed these issues by creating and analyzing a comprehensive dataset of well-characterized ID/globular protein complexes, for which both the atomic structure of the complex and free energy (ΔG, K_d ) of interaction is known. Through this analysis, we provide evidence that the affinity distributions of IDP/globular and globular/globular complexes show different trends, whereas specificity does not connote to weaker binding strength of IDPs/IDRs. Furthermore, protein disorder extends the spectrum in the direction of very weak interactions, which may have important regulatory consequences and suggest that, in a biological sense, strict correlation of specificity and binding strength are uncoupled by structural disorder.

Short Linear Motifs Orchestrate Functioning of Human Proteins during Embryonic Development, Redox Regulation, and Cancer

Short linear motifs (SLiMs) are evolutionarily conserved functional modules of proteins that represent amino acid stretches composed of 3 to 10 residues. The biological activities of two short peptide segments of human alpha-fetoprotein (AFP), a major embryo-specific and cancer-related protein, have been confirmed experimentally. This is a heptapeptide segment LDSYQCT in domain I designated as AFP14–20 and a nonapeptide segment EMTPVNPGV in domain III designated as GIP-9. In our work, we searched the UniprotKB database for human proteins that contain SLiMs with sequence similarity to the both segments of human AFP and undertook gene ontology (GO)-based functional categorization of retrieved proteins. Gene set enrichment analysis included GO terms for biological process, molecular function, metabolic pathway, KEGG pathway, and protein–protein interaction (PPI) categories. We identified the SLiMs of interest in a variety of non-homologous proteins involved in multiple cellular processes underlying embryonic development, cancer progression, and, unexpectedly, the regulation of redox homeostasis. These included transcription factors, cell adhesion proteins, ubiquitin-activating and conjugating enzymes, cell signaling proteins, and oxidoreductase enzymes. They function by regulating cell proliferation and differentiation, cell cycle, DNA replication/repair/recombination, metabolism, immune/inflammatory response, and apoptosis. In addition to the retrieved genes, new interacting genes were identified. Our data support the hypothesis that conserved SLiMs are incorporated into non-homologous proteins to serve as functional blocks for their orchestrated functioning.

Insights into Membrane Curvature Sensing and Membrane Remodeling by Intrinsically Disordered Proteins and Protein Regions

Cellular membranes are highly dynamic in shape. They can rapidly and precisely regulate their shape to perform various cellular functions. The protein’s ability to sense membrane curvature is essential in various biological events such as cell signaling and membrane trafficking. As they are bound, these curvature-sensing proteins may also change the local membrane shape by one or more curvature driving mechanisms. Established curvature-sensing/driving mechanisms rely on proteins with specific structural features such as amphipathic helices and intrinsically curved shapes. However, the recent discovery and characterization of many proteins have shattered the protein structure–function paradigm, believing that the protein functions require a unique structural feature. Typically, such structure-independent functions are carried either entirely by intrinsically disordered proteins or hybrid proteins containing disordered regions and structured domains. It is becoming more apparent that disordered proteins and regions can be potent sensors/inducers of membrane curvatures. In this article, we outline the basic features of disordered proteins and regions, the motifs in such proteins that encode the function, membrane remodeling by disordered proteins and regions, and assays that may be employed to investigate curvature sensing and generation by ordered/disordered proteins.

Short Linear Motifs (SLiMs) in “Core” RxLR Effectors of Phytophthora parasitica var. nicotianae : a Case of PpRxLR1 Effector

Oomycetes of the genus Phytophthora encompass several of the most successful plant pathogens described to date. The success of infection by Phytophthora species is attributed to the pathogens’ ability to secrete effector proteins that alter the host’s physiological processes. Structural analyses of effector proteins mainly from bacterial and viral pathogens have revealed the presence of intrinsically disordered regions that host short linear motifs (SLiMs). These motifs play important biological roles by facilitating protein-protein interactions as well as protein translocation. Nonetheless, SLiMs in Phytophthora species RxLR effectors have not been investigated previously and their roles remain unknown. Using a bioinformatics pipeline, we identified 333 candidate RxLR effectors in the strain INRA 310 of Phytophthora parasitica. Of these, 71 (21%) were also found to be present in 10 other genomes of P. parasitica, and hence, these were designated core RxLR effectors (CREs). Within the CRE sequences, the N terminus exhibited enrichment in intrinsically disordered regions compared to the C terminus, suggesting a potential role of disorder in effector translocation. Although the disorder content was reduced in the C-terminal regions, it is important to mention that most SLiMs were in this terminus. PpRxLR1 is one of the 71 CREs identified in this study, and its genes encode a 6-amino acid (aa)-long SLiM at the C terminus. We showed that PpRxLR1 interacts with several host proteins that are implicated in defense. Structural analysis of this effector using homology modeling revealed the presence of potential ligand-binding sites. Among key residues that were predicted to be crucial for ligand binding, L¹⁰² and Y¹⁰⁶ were of interest since they form part of the 6-aa-long PpRxLR1 SLiM. In silico substitution of these two residues to alanine was predicted to have a significant effect on both the function and the structure of PpRxLR1 effector. Molecular docking simulations revealed possible interactions between PpRxLR1 effector and ubiquitin-associated proteins. The ubiquitin-like SLiM carried in this effector was shown to be a potential mediator of these interactions. Further studies are required to validate and elucidate the underlying molecular mechanism of action.

IMPORTANCE The continuous gain and loss of RxLR effectors makes the control of Phytophthora spp. difficult. Therefore, in this study, we endeavored to identify RxLR effectors that are highly conserved among species, also known as “core” RxLR effectors (CREs). We reason that these highly conserved effectors target conserved proteins or processes; thus, they can be harnessed in breeding for durable resistance in plants. To further understand the mechanisms of action of CREs, structural dissection of these proteins is crucial. Intrinsically disordered regions (IDRs) that do not adopt a fixed, three-dimensional fold carry short linear motifs (SLiMs) that mediate biological functions of proteins. The presence and potential role of these SLiMs in CREs of Phytophthora spp. have been overlooked. To our knowledge, we have effectively identified CREs as well as SLiMs with the potential of promoting effector virulence. Together, this work has advanced our comprehension of Phytophthora RxLR effector function and may facilitate the development of innovative and effective control strategies.

Host diversification is concurrent with linear motif evolution in a Mastadenovirus hub protein

Over one hundred Mastadenovirus types infect seven orders of mammals. Virus-host coevolution may involve cospeciation, duplication, host switch and partial extinction events. We reconstruct Mastadenovirus diversification, finding that while cospeciation is dominant, the other three events are also common in Mastadenovirus evolution. Linear motifs are fast-evolving protein functional elements and key mediators of virus-host interactions, thus likely to partake in adaptive viral evolution. We study the evolution of eleven linear motifs in the Mastadenovirus E1A protein, a hub of virus-host protein-protein interactions, in the context of host diversification. The reconstruction of linear motif gain and loss events shows fast linear motif turnover, corresponding a virus-host protein-protein interaction turnover orders of magnitude faster than in model host proteomes. Evolution of E1A linear motifs is coupled, indicating functional coordination at the protein scale, yet presents motif-specific patterns suggestive of convergent evolution. We report a pervasive association between Mastadenovirus host diversification events and the evolution of E1A linear motifs. Eight of 17 host switches associate with the gain of one linear motif and the loss of four different linear motifs, while five of nine partial extinctions associate with the loss of one linear motif. The specific changes in E1A linear motifs during a host switch or a partial extinction suggest that changes in the host molecular environment lead to modulation of the interactions with the retinoblastoma protein and host transcriptional regulators. Altogether, changes in the linear motif repertoire of a viral hub protein are associated with adaptive evolution events during Mastadenovirus evolution.

Repeating Aspartic Acid Residues Prefer Turn-like Conformations in the Unfolded State: Implications for Early Protein Folding

Protein folding can be described as a motion of the polypeptide chain in a potential energy funnel, where the conformational manifold is narrowed as the chain traverses from a completely unfolded state until it reaches the folded (native) state. The initial folding stages set the tone for this process by substantially narrowing the manifold of accessible conformations. In an ideally unfolded state with no long-range stabilizing forces, local conformations (i.e., residual structures) are likely to drive the folding process. While most amino acid residues tend to predominantly adopt extended structures in unfolded proteins and peptides, aspartic acid exhibits a relatively high intrinsic preference for turn-forming conformations. Regions in an unfolded polypeptide or protein that are rich in aspartic acid residues may therefore be crucial sites for protein folding steps. By combining NMR and vibrational spectroscopies, we observed that the conformational sampling of multiple sequentially neighbored aspartic acid residues in the model peptides GDDG and GDDDG even show an on average higher propensity for turn-forming structures than the intrinsic reference system D in GDG, which suggests that nearest neighbor interactions between adjacent aspartic acid residues stabilize local turn-forming structures. In the presence of the unlike neighbor phenylalanine, nearest neighbor interactions are of a totally different nature in that it they decrease the turn-forming propensities and mutually increase the sampling of polyproline II (pPII) conformations. We hypothesize the structural role of aspartic residues in intrinsically disordered proteins in general, and particularly in small linear motifs, that are very much determined by their respective neighbors.

The HGR motif is the antiangiogenic determinant of vasoinhibin: implications for a therapeutic orally active oligopeptide

The hormone prolactin acquires antiangiogenic and antivasopermeability properties after undergoing proteolytic cleavage to vasoinhibin, an endogenous prolactin fragment of 123 or more amino acids that inhibits the action of multiple proangiogenic factors. Preclinical and clinical evidence supports the therapeutic potential of vasoinhibin against angiogenesis-related diseases including diabetic retinopathy, peripartum cardiomyopathy, rheumatoid arthritis, and cancer. However, the use of vasoinhibin in the clinic has been limited by difficulties in its production. Here, we removed this barrier to using vasoinhibin as a therapeutic agent by showing that a short linear motif of just three residues (His46-Gly47-Arg48) (HGR) is the functional determinant of vasoinhibin. The HGR motif is conserved throughout evolution, its mutation led to vasoinhibin loss of function, and oligopeptides containing this sequence inhibited angiogenesis and vasopermeability with the same potency as whole vasoinhibin. Furthermore, the oral administration of an optimized cyclic retro-inverse vasoinhibin heptapeptide containing HGR inhibited melanoma tumor growth and vascularization in mice and exhibited equal or higher antiangiogenic potency than other antiangiogenic molecules currently used as anti-cancer drugs in the clinic. Finally, by unveiling the mechanism that obscures the HGR motif in prolactin, we anticipate the development of vasoinhibin-specific antibodies to solve the on-going challenge of measuring endogenous vasoinhibin levels for diagnostic and interventional purposes, the design of vasoinhibin antagonists for managing insufficient angiogenesis, and the identification of putative therapeutic proteins containing HGR.

Genomic analysis uncovers functional variation in the C-terminus of anthocyanin-activating MYB transcription factors

MYB transcription factors regulate diverse aspects of plant development and secondary metabolism, often by partnering in transcriptional regulatory complexes. Here, we harness genomic resources to identify novel MYBs, thereby producing an updated eudicot MYB phylogeny with revised relationships among subgroups as well as new information on sequence variation in the disordered C-terminus of anthocyanin-activating MYBs. BLAST® and hidden Markov model scans of gene annotations identified a total of 714 MYB transcription factors across the genomes of four crops that span the eudicots: apple, grape, kiwifruit and tomato. Codon model-based phylogenetic inference identified novel members of previously defined subgroups, and the function of specific anthocyanin-activating subgroup 6 members was assayed transiently in tobacco leaves. Sequence conservation within subgroup 6 highlighted one previously described and two novel short linear motifs in the disordered C-terminal region. The novel motifs have a mix of hydrophobic and acidic residues and are predicted to be relatively ordered compared with flanking protein sequences. Comparison of motifs with the Eukaryotic Linear Motif database suggests roles in protein-protein interaction. Engineering of motifs and their flanking regions from strong anthocyanin activators into weak activators, and vice versa, affected function. We conclude that, although the MYB C-terminal sequence diverges greatly even within MYB clades, variation within the C-terminus at and near relatively ordered regions offers opportunities for exploring MYB function and developing superior alleles for plant breeding.

Short peptides as predictors for the structure of polyarginine sequences in disordered proteins

Intrinsically disordered proteins and intrinsically disordered regions are frequently enriched in charged amino acids. Intrinsically disordered regions are regularly involved in important biological processes in which one or more charged residues is the driving force behind a protein-biomolecule interaction. Several lines of experimental and computational evidence suggest that polypeptides and proteins that carry high net charges have a high preference for extended conformations with average end-to-end distances exceeding expectations for self-avoiding random coils. Here, we show that charged arginine residues even in short glycine-capped model peptides (GRRG and GRRRG) significantly affect the conformational propensities of each other when compared with the intrinsic propensities of a mostly unperturbed arginine in the tripeptide GRG. A conformational analysis based on experimentally determined J-coupling constants from heteronuclear NMR spectroscopy and amide I' band profiles from vibrational spectroscopy reveals that nearest-neighbor interactions stabilize extended β-strand conformations at the expense of polyproline II and turn conformations. The results from molecular dynamics simulations with a CHARMM36m force field and TIP3P water reproduce our results only to a limited extent. The use of the Ramachandran distribution of the central residue of GRRRG in a calculation of end-to-end distances of polyarginines of different length yielded the expected power law behavior. The scaling coefficient of 0.66 suggests that such peptides would be more extended than predicted by a self-avoiding random walk. Our findings thus support in principle theoretical predictions.

The Phytophthora infestans AVR2 Effector Escapes R2 Recognition Through Effector Disordering

Intrinsic disorder is a common structural characteristic of proteins and a central player in the biochemical processes of species. However, the role of intrinsic disorder in the evolution of plant-pathogen interactions is rarely investigated. Here, we explored the role of intrinsic disorder in the development of the pathogenicity in the RXLR AVR2 effector of Phytophthora infestans. We found AVR2 exhibited high nucleotide diversity generated by point mutation, early-termination, altered start codon, deletion/insertion, and intragenic recombination and is predicted to be an intrinsically disordered protein. AVR2 amino acid sequences conferring a virulent phenotype had a higher disorder tendency in both the N- and C-terminal regions compared with sequences conferring an avirulent phenotype. In addition, we also found virulent AVR2 mutants gained one or two short linear interaction motifs, the critical components of disordered proteins required for protein-protein interactions. Furthermore, virulent AVR2 mutants were predicted to be unstable and have a short protein half-life. Taken together, these results support the notion that intrinsic disorder is important for the effector function of pathogens and demonstrate that SLiM-mediated protein-protein interaction in the C-terminal effector domain might contribute greatly to the evasion of resistance-protein detection in P. infestans.

Intrinsically disordered proteins of viruses: Involvement in the mechanism of cell regulation and pathogenesis

Intrinsically disordered proteins (IDPs) possess the property of inherent flexibility and can be distinguished from other proteins in terms of lack of any fixed structure. Such dynamic behavior of IDPs earned the name "Dancing Proteins." The exploration of these dancing proteins in viruses has just started and crucial details such as correlation of rapid evolution, high rate of mutation and accumulation of disordered contents in viral proteome at least understood partially. In order to gain a complete understanding of this correlation, there is a need to decipher the complexity of viral mediated cell hijacking and pathogenesis in the host organism. Further there is necessity to identify the specific patterns within viral and host IDPs such as aggregation; Molecular recognition features (MoRFs) and their association to virulence, host range and rate of evolution of viruses in order to tackle the viral-mediated diseases. The current book chapter summarizes the aforementioned details and suggests the novel opportunities for further research of IDPs senses in viruses.

One Actor, Multiple Roles: The Performances of Cryptochrome in Drosophila

Cryptochromes (CRYs) are flavoproteins that are sensitive to blue light, first identified in Arabidopsis and then in Drosophila and mice. They are evolutionarily conserved and play fundamental roles in the circadian clock of living organisms, enabling them to adapt to the daily 24-h cycles. The role of CRYs in circadian clocks differs among different species: in plants, they have a blue light-sensing activity whereas in mammals they act as light-independent transcriptional repressors within the circadian clock. These two different functions are accomplished by two principal types of CRYs, the light-sensitive plant/insect type 1 CRY and the mammalian type 2 CRY acting as a negative autoregulator in the molecular circadian clockwork. Drosophila melanogaster possesses just one CRY, belonging to type 1 CRYs. Nevertheless, this single CRY appears to have different functions, specific to different organs, tissues, and even subset of cells in which it is expressed. In this review, we will dissect the multiple roles of this single CRY in Drosophila, focusing on the regulatory mechanisms that make its pleiotropy possible.

Understanding protein multifunctionality: from short linear motifs to cellular functions

Moonlighting proteins perform multiple unrelated functions without any change in polypeptide sequence. They can coordinate cellular activities, serving as switches between pathways and helping to respond to changes in the cellular environment. Therefore, regulation of the multiple protein activities, in space and time, is likely to be important for the homeostasis of biological systems. Some moonlighting proteins may perform their multiple functions simultaneously while others alternate between functions due to certain triggers. The switch of the moonlighting protein's functions can be regulated by several distinct factors, including the binding of other molecules such as proteins. We here review the approaches used to identify moonlighting proteins and existing repositories. We particularly emphasise the role played by short linear motifs and PTMs as regulatory switches of moonlighting functions.

Evolution of a Landscape Phage Library in a Mouse Xenograft Model of Human Breast Cancer

Peptide-displayed phage libraries are billion-clone collections of diverse chimeric bacteriophage particles, decorated by genetically fused peptides built from a random combination of natural amino acids. Studying the molecular evolution of peptide-displayed libraries in mammalian model systems, using in vivo phage display techniques, can provide invaluable knowledge about the underlying physiology of the vasculature system, allow recognition of organ- and tissue-specific networks of protein-protein interactions, and provide ligands for targeted diagnostics and therapeutics. Recently, we discovered that landscape phage libraries, a specific type of multivalent peptide phage display library, expose on their surface comprehensive collections of elementary binding units (EBUs), which can form short linear motifs (SLiMs) that interact with functional domains of physiologically relevant proteins. Because of their unique structural and functional features, landscape phages can use an alternative mechanism of directed molecular evolution, i.e., combinatorial avidity selection. These discoveries fueled our interest in revisiting the in vivo evolution of phage displayed libraries using another format of display, i.e., landscape phages. In this study, we monitored the evolution of a landscape phage library in a mouse model with and without an implanted human breast cancer tumor xenograft. As expected, the multivalent architecture of landscape phage displayed proteins provided strong tissue selectivity and resulted in a huge diversity of tissue penetrating, chimeric phage particles. We identified several types of EBU interactions that evolved during the course of tissue distribution, which included interactions of EBUs with all tissue types, those EBUs that interacted selectively with specific organs or tissues with shared gene expression profiles or functionalities, and other EBUs that interacted in a tissue-selective manner. We demonstrated that landscape phage libraries are a rich collection of unique nanobioparticles that can be used to identify functional organ and tissue-binding elements after the evolution of a phage display library in vivo.

Combinatorial Avidity Selection of Mosaic Landscape Phages Targeted at Breast Cancer Cells-An Alternative Mechanism of Directed Molecular Evolution

Low performance of actively targeted nanomedicines required revision of the traditional drug targeting paradigm and stimulated the development of novel phage-programmed, self-navigating drug delivery vehicles. In the proposed smart vehicles, targeting peptides, selected from phage libraries using traditional principles of affinity selection, are substituted for phage proteins discovered through combinatorial avidity selection. Here, we substantiate the potential of combinatorial avidity selection using landscape phage in the discovery of Short Linear Motifs (SLiMs) and their partner domains. We proved an algorithm for analysis of phage populations evolved through multistage screening of landscape phage libraries against the MDA-MB-231 breast cancer cell line. The suggested combinatorial avidity selection model proposes a multistage accumulation of Elementary Binding Units (EBU), or Core Motifs (CorMs), in landscape phage fusion peptides, serving as evolutionary initiators for formation of SLiMs. Combinatorial selection has the potential to harness directed molecular evolution to create novel smart materials with diverse novel, emergent properties.

RNA-protein interactions: disorder, moonlighting and junk contribute to eukaryotic complexity

RNA-protein interactions are crucial for most biological processes in all organisms. However, it appears that the complexity of RNA-based regulation increases with the complexity of the organism, creating additional regulatory circuits, the scope of which is only now being revealed. It is becoming apparent that previously unappreciated features, such as disordered structural regions in proteins or non-coding regions in DNA leading to higher plasticity and pliability in RNA-protein complexes, are in fact essential for complex, precise and fine-tuned regulation. This review addresses the issue of the role of RNA-protein interactions in generating eukaryotic complexity, focusing on the newly characterized disordered RNA-binding motifs, moonlighting of metabolic enzymes, RNA-binding proteins interactions with different RNA species and their participation in regulatory networks of higher order.

The carboxy-terminus, a key regulator of protein function

All proteins end with a carboxyl terminus that has unique biophysical properties and is often disordered. Although there are examples of important C-termini functions, a more global role for the C-terminus is not yet established. In this review, we summarize research on C-termini, a unique region in proteins that cells exploit. Alternative splicing and proteolysis increase the diversity of proteins and peptides in cells with unique C-termini. The C-termini of proteins contain minimotifs, short peptides with an encoded function generally characterized as binding, posttranslational modifications, and trafficking. Many of these activities are specific to minimotifs on the C-terminus. Approximately 13% of C-termini in the human proteome have a known minimotif, and the majority, if not all of the remaining termini have conserved motifs inferring a function that remains to be discovered. C-termini, their predictions, and their functions are collated in the C-terminome, Proteus, and Terminus Oriented Protein Function INferred Database (TopFIND) database/web systems. Many C-termini are well conserved, and some have a known role in health and disease. We envision that this summary of C-termini will guide future investigation of their biochemical and physiological significance.

Autism risk genes are evolutionarily ancient and maintain a unique feature landscape that echoes their function

Previous research on autism risk (ASD), developmental regulatory (DevReg), and central nervous system (CNS) genes suggests they tend to be large in size, enriched in nested repeats, and mutation intolerant. The relevance of these genomic features is intriguing yet poorly understood. In this study, we investigated the feature landscape of these gene groups to discover structural themes useful in interpreting their function, developmental patterns, and evolutionary history. ASD, DevReg, CNS, housekeeping, and whole genome control (WGC) groups were compiled using various resources. Multiple gene features of interest were extracted from NCBI/UCSC Bioinformatics. Residual variation intolerance scores, Exome Aggregation Consortium pLI scores, and copy number variation data from Decipher were used to estimate variation intolerance. Gene age and protein-protein interactions (PPI) were estimated using Ensembl and EBI Intact databases, respectively. Compared to WGC: ASD, DevReg, and CNS genes are longer, produce larger proteins, maintain greater numbers/density of conserved noncoding elements and transposable elements, produce more transcript variants, and are comparatively variation intolerant. After controlling for gene size, mutation tolerance, and clinical association, ASD genes still retain many of these same features. In addition, we also found that ASD genes that are extremely mutation intolerant have larger PPI networks. These data support many of the recent findings within the field of autism genetics but also expand our understanding of the evolution of these broad gene groups, their potential regulatory complexity, and the extent to which they interact with the cellular network. Autism Res 2019, 12: 860-869. © 2019 International Society for Autism Research, Wiley Periodicals, Inc. LAY SUMMARY: Autism risk genes are more ancient compared to other genes in the genome. As such, they exhibit physical features related to their age, including long gene and protein size and regulatory sequences that help to control gene expression. They share many of these same features with other genes that are expressed in the brain and/or are associated with prenatal development.

Histidine, the less interactive cousin of arginine

Electrostatic interactions are one of the main factors influencing biomolecular conformation. The formation of noncovalent complexes by electrostatic interactions is governed by certain amino acid residues and post-translational modifications. It has been demonstrated that adjacent arginine forms noncovalent complex with phosphate; however, histidine noncovalent complexes have rarely been investigated. In the present work, we compare the interaction between basic epitopes (NLRRITRVN, SHHGLHSTPD) and diverse acidic and aromatic-rich peptides using both MALDI and ESI Mass spectrometry. We show that adjacent histidines can also form stable noncovalent bonds and that those bonds are probably formed by a salt bridge between the phosphate or the acid residues and the histidines. However, noncovalent complexes with the arginine epitopes form more readily and are stronger than those with histidine-containing epitopes.

Rewiring of RSK-PDZ Interactome by Linear Motif Phosphorylation

Phosphorylation of short linear peptide motifs is a widespread process for the dynamic regulation of protein-protein interactions. However, the global impact of phosphorylation events on the protein-protein interactome is rarely addressed. The disordered C-terminal tail of ribosomal S6 kinase 1 (RSK1) binds to PDZ domain-containing scaffold proteins, and it harbors a phosphorylatable PDZ-binding motif (PBM) responsive to epidermal growth factor stimulation. Here, we examined binding of two versions of the RSK1 PBM, either phosphorylated or unphosphorylated at position -3, to almost all (95%) of the 266 PDZ domains of the human proteome. PBM phosphorylation dramatically altered the PDZ domain-binding landscape of RSK1, by strengthening or weakening numerous interactions to various degrees. The RSK-PDZome interactome analyzed in this study reveals how linear motif-based phospho-switches convey stimulus-dependent changes in the context of related network components.

Screening and Identification of PLK1-Polo Box Binding Peptides by High-Throughput Sequencing of Phage-Selected Libraries

BACKGROUND

Human proteome contains a plethora of short linear peptide motifs that is crucial for signaling and other cellular processes. These motifs are difficult to identify due to lack of systematic approach for their detection.

OBJECTIVES

Here we demonstrate the use of peptide phage display in combination with high throughput next generation sequencing to identify enriched peptide sequences through biopanning process against polo box domain (PBD) of mitotic polo like kinase 1 (Plk1).

METHODS

Purified recombinant Plk1 and two unrelated controls namely B-lymphocyte antigen (CD20) and fluorescent protein (mCherry) were subjected to peptide phage display analysis. Bacterially-propagated phage DNA was amplified by PCR using triplet bar coded primers to tag the pool from each amplicon.

RESULTS

Proteomic peptide phage display along with next generation sequencing and Bioinformatics analysis demonstrated several known and putative novel interactions which were potentially related to Plk1-PBD. With our strategy, we were able to identify and characterize several Plk1-PBD binding peptides, as well as define more precisely, consensus sequences.

CONCLUSION

We believe that this information could provide valuable tools for exploring novel interaction involved in Plk1 signaling as well as to choose peptides for Plk1 specific drug development.

Evidence for Stabilizing Selection Driving Mutational Turnover of Short Motifs in the Eukaryotic Complementary Sex Determiner (Csd) Protein

Short linear motifs (SLiMs) can play pivotal functional roles in proteins, such as targeting proteins to specific subcellular localizations, modulating the efficiency of translation and tagging proteins for degradation. Until recently we had little knowledge about SLiM evolution. Only a few amino acids in these motifs are functionally important, making them likely to evolve ex nihilo and suggesting that they can play key roles in protein evolution. Several reports now suggest that these motifs can appear and disappear while their function in the protein is preserved, a process sometimes referred to as “turnover”. However, there has been a lack of specific experiments to determine whether independently evolved motifs do indeed have the same function, which would conclusively determine whether the process of turnover actually occurs. In this study, we experimentally detected evidence for such a mutational turnover process for nuclear localization signals (NLS) during the post-duplication divergence of the Complementary sex determiner (Csd) and Feminizer (Fem) proteins in the honeybee (Apis mellifera) lineage. Experiments on the nuclear transport activity of protein segments and those of the most recent common ancestor (MRCA) sequences revealed that three new NLS motifs evolved in the Csd protein during the post-duplication divergence while other NLS motifs were lost that existed before duplication. A screen for essential and newly evolved amino acids revealed that new motifs in the Csd protein evolved by one or two missense mutations coding for lysine. Amino acids that were predating the duplication were also essential in the acquisition of the C1 motif suggesting that the ex nihilo origin was constrained by preexisting amino acids in the physical proximity. Our data support a model in which stabilizing selection maintains the constancy of nuclear transport function but allowed mutational turnover of the encoding NLS motifs.

The Transcriptional Repressor BS69 is a Conserved Target of the E1A Proteins from Several Human Adenovirus Species

Early region 1A (E1A) is the first viral protein produced upon human adenovirus (HAdV) infection. This multifunctional protein transcriptionally activates other HAdV early genes and reprograms gene expression in host cells to support productive infection. E1A functions by interacting with key cellular regulatory proteins through short linear motifs (SLiMs). In this study, the molecular determinants of interaction between E1A and BS69, a cellular repressor that negatively regulates E1A transactivation, were systematically defined by mutagenesis experiments. We found that a minimal sequence comprised of MPNLVPEV, which contains a conserved PXLXP motif and spans residues 112–119 in HAdV-C5 E1A, was necessary and sufficient in binding to the myeloid, Nervy, and DEAF-1 (MYND) domain of BS69. Our study also identified residues P113 and L115 as critical for this interaction. Furthermore, the HAdV-C5 and -A12 E1A proteins from species C and A bound BS69, but those of HAdV-B3, -E4, -D9, -F40, and -G52 from species B, E, D, F, and G, respectively, did not. In addition, BS69 functioned as a repressor of E1A-mediated transactivation, but only for HAdV-C5 and HAdV-A12 E1A. Thus, the PXLXP motif present in a subset of HAdV E1A proteins confers interaction with BS69, which serves as a negative regulator of E1A mediated transcriptional activation.

Modulation of the aggregation of an amyloidogenic sequence by flanking-disordered region in the intrinsically disordered antigen merozoite surface protein 2

The abundant Plasmodium falciparum merozoite surface protein MSP2, a potential malaria vaccine candidate, is an intrinsically disordered protein with some nascent secondary structure present in its conserved N-terminal region. This relatively ordered region has been implicated in both membrane interactions and amyloid-like aggregation of the protein, while the significance of the flanking-disordered region is unclear. In this study, we show that aggregation of the N-terminal conserved region of MSP2 is influenced in a length- and sequence-dependent fashion by the disordered central variable sequences. Intriguingly, MSP2 peptides containing the conserved region and the first five residues of the variable disordered regions aggregated more rapidly than a peptide corresponding to the conserved region alone. In contrast, MSP2 peptides extending 8 or 12 residues into the disordered region aggregated more slowly, consistent with the expected inhibitory effect of flanking-disordered sequences on the aggregation of amyloidogenic ordered sequences. Computational analyses indicated that the helical propensity of the ordered region of MSP2 was modulated by the adjacent disordered five residues in a sequence-dependent manner. Nuclear magnetic resonance and circular dichroism spectroscopic studies with synthetic peptides confirmed the computational predictions, emphasizing the correlation between aggregation propensity and conformation of the ordered region and the effects thereon of the adjacent disordered region. These results show that the effects of flanking-disordered sequences on a more ordered sequence may include enhancement of aggregation through modulation of the conformational properties of the more ordered sequence.

Order in Disorder as Observed by the "Hydrophobic Cluster Analysis" of Protein Sequences

Hydrophobic cluster analysis (HCA) is an original approach for protein sequence analysis, which provides access to the foldable repertoire of the protein universe, including yet unannotated protein segments ("dark proteome"). Foldable segments correspond to ordered regions, as well as to intrinsically disordered regions (IDRs) undergoing disorder to order transitions. In this review, how HCA can be used to give insight into this last category of foldable segments is illustrated, with examples matching known 3D structures. After reviewing the HCA principles, examples of short foldable segments are given, which often contain short linear motifs, typically matching hydrophobic clusters. These segments become ordered upon contact with partners, with secondary structure preferences generally corresponding to those observed in the 3D structures within the complexes. Such small foldable segments are sometimes larger than the segments of known 3D structures, including flanking hydrophobic clusters that may be critical for interaction specificity or regulation, as well as intervening sequences allowing fuzziness. Cases of larger conditionally disordered domains are also presented, with lower density in hydrophobic clusters than well-folded globular domains or with exposed hydrophobic patches, which are stabilized by interaction with partners.

SLiM-Enrich: computational assessment of protein-protein interaction data as a source of domain-motif interactions

Many important cellular processes involve protein–protein interactions (PPIs) mediated by a Short Linear Motif (SLiM) in one protein interacting with a globular domain in another. Despite their significance, these domain-motif interactions (DMIs) are typically low affinity, which makes them challenging to identify by classical experimental approaches, such as affinity pulldown mass spectrometry (AP-MS) and yeast two-hybrid (Y2H). DMIs are generally underrepresented in PPI networks as a result. A number of computational methods now exist to predict SLiMs and/or DMIs from experimental interaction data but it is yet to be established how effective different PPI detection methods are for capturing these low affinity SLiM-mediated interactions. Here, we introduce a new computational pipeline (SLiM-Enrich) to assess how well a given source of PPI data captures DMIs and thus, by inference, how useful that data should be for SLiM discovery. SLiM-Enrich interrogates a PPI network for pairs of interacting proteins in which the first protein is known or predicted to interact with the second protein via a DMI. Permutation tests compare the number of known/predicted DMIs to the expected distribution if the two sets of proteins are randomly associated. This provides an estimate of DMI enrichment within the data and the false positive rate for individual DMIs. As a case study, we detect significant DMI enrichment in a high-throughput Y2H human PPI study. SLiM-Enrich analysis supports Y2H data as a source of DMIs and highlights the high false positive rates associated with naïve DMI prediction. SLiM-Enrich is available as an R Shiny app. The code is open source and available via a GNU GPL v3 license at: https://github.com/slimsuite/SLiMEnrich. A web server is available at: http://shiny.slimsuite.unsw.edu.au/SLiMEnrich/.

Large-scale docking predicts that sORF-encoded peptides may function through protein-peptide interactions in Arabidopsis thaliana

Several recent studies indicate that small Open Reading Frames (sORFs) embedded within multiple eukaryotic non-coding RNAs can be translated into bioactive peptides of up to 100 amino acids in size. However, the functional roles of the 607 Stress Induced Peptides (SIPs) previously identified from 189 Transcriptionally Active Regions (TARs) in Arabidopsis thaliana remain unclear. To provide a starting point for functional annotation of these plant-derived peptides, we performed a large-scale prediction of peptide binding sites on protein surfaces using coarse-grained peptide docking. The docked models were subjected to further atomistic refinement and binding energy calculations. A total of 530 peptide-protein pairs were successfully docked. In cases where a peptide encoded by a TAR is predicted to bind at a known ligand or cofactor-binding site within the protein, it can be assumed that the peptide modulates the ligand or cofactor-binding. Moreover, we predict that several peptides bind at protein-protein interfaces, which could therefore regulate the formation of the respective complexes. Protein-peptide binding analysis further revealed that peptides employ both their backbone and side chain atoms when binding to the protein, forming predominantly hydrophobic interactions and hydrogen bonds. In this study, we have generated novel predictions on the potential protein-peptide interactions in A. thaliana, which will help in further experimental validation.

Identification and Analysis of OVATE Family Members from Genome of the Early Land Plants Provide Insights into Evolutionary History of OFP Family and Function

Mosses, liverworts, hornworts and lycophytes represent transition stages between the aquatic to terrestrial/land plants. Several morphological and adaptive novelties driven by genomic components including emergence and expansion of new or existing gene families have played a critical role during and after the transition, and contributed towards successful colonization of terrestrial ecosystems. It is crucial to decipher the evolutionary transitions and natural selection on the gene structure and function to understand the emergence of phenotypic and adaptive diversity. Plants at the "transition zone", between aquatic and terrestrial ecosystem, are also the most vulnerable because of climate change and may contain clues for successful mitigation of the challenges of climate change. Identification and comparative analyses of such genetic elements and gene families are few in mosses, liverworts, hornworts and lycophytes. Ovate family proteins (OFPs) are plant-specific transcriptional repressors and are acknowledged for their roles in important growth and developmental processes in land plants, and information about the functional aspects of OFPs in early land plants is fragmentary. As a first step towards addressing this gap, a comprehensive in silico analysis was carried out utilizing publicly available genome sequences of Marchantia polymorpha (Mp), Physcomitrella patens (Pp), Selaginella moellendorffii (Sm) and Sphagnum fallax (Sf). Our analysis led to the identification of 4 MpOFPs, 19 PpOFPs, 6 SmOFPs and 3 SfOFPs. Cross-genera analysis revealed a drastic change in the structure and physiochemical properties in OFPs suggesting functional diversification and genomic plasticity during the evolutionary course. Knowledge gained from this comparative analysis will form the framework towards deciphering and dissection of their developmental and adaptive role/s in early land plants and could provide insights into evolutionary strategies adapted by land plants.