ELM-the eukaryotic linear motif resource in 2020

Structure-based prediction of HDAC6 substrates validated by enzymatic assay reveals determinants of promiscuity and detects new potential substrates

Histone deacetylases play important biological roles well beyond the deacetylation of histone tails. In particular, HDAC6 is involved in multiple cellular processes such as apoptosis, cytoskeleton reorganization, and protein folding, affecting substrates such as ɑ-tubulin, Hsp90 and cortactin proteins. We have applied a biochemical enzymatic assay to measure the activity of HDAC6 on a set of candidate unlabeled peptides. These served for the calibration of a structure-based substrate prediction protocol, Rosetta FlexPepBind, previously used for the successful substrate prediction of HDAC8 and other enzymes. A proteome-wide screen of reported acetylation sites using our calibrated protocol together with the enzymatic assay provide new peptide substrates and avenues to novel potential functional regulatory roles of this promiscuous, multi-faceted enzyme. In particular, we propose novel regulatory roles of HDAC6 in tumorigenesis and cancer cell survival via the regulation of EGFR/Akt pathway activation. The calibration process and comparison of the results between HDAC6 and HDAC8 highlight structural differences that explain the established promiscuity of HDAC6.

Sensitive Identification of Known and Unknown Protease Activities by Unsupervised Linear Motif Deconvolution

The cleavage-site specificities for many proteases are not well understood, restricting the utility of supervised classification methods. We present an algorithm and web interface to overcome this limitation through the unsupervised detection of overrepresented patterns in protein sequence data, providing insight into the mixture of protease activities contributing to a complex system. Here, we apply the RObust LInear Motif Deconvolution (RoLiM) algorithm to confidently detect substrate cleavage patterns for SARS-CoV-2 MPro protease in the N-terminome data of an infected human cell line. Using mass spectrometry-based peptide data from a case-control comparison of 341 primary urothelial bladder cancer cases and 110 controls, we identified distinct sequence motifs indicative of increased matrix metallopeptidase activity in urine from cancer patients. The evaluation of N-terminal peptides from patient plasma post-chemotherapy detected novel granzyme B/corin activity. RoLiM will enhance the unbiased investigation of peptide sequences to establish the composition of known and uncharacterized protease activities in biological systems. RoLiM is available at http://langelab.org/rolim/.

Orchestrating serine/threonine phosphorylation and elucidating downstream effects by short linear motifs

Cellular function is based on protein-protein interactions. A large proportion of these interactions involves the binding of short linear motifs (SLiMs) by folded globular domains. These interactions are regulated by post-translational modifications, such as phosphorylation, that create and break motif binding sites or tune the affinity of the interactions. In addition, motif-based interactions are involved in targeting serine/threonine kinases and phosphatases to their substrate and contribute to the specificity of the enzymatic actions regulating which sites are phosphorylated. Here, we review how SLiM-based interactions assist in determining the specificity of serine/threonine kinases and phosphatases, and how phosphorylation, in turn, affects motif-based interactions. We provide examples of SLiM-based interactions that are turned on/off, or are tuned by serine/threonine phosphorylation and exemplify how this affects SLiM-based protein complex formation.

Linking post-translational modifications and protein turnover by site-resolved protein turnover profiling

Proteome-wide measurements of protein turnover have largely ignored the impact of post-translational modifications (PTMs). To address this gap, we employ stable isotope labeling and mass spectrometry to measure the turnover of >120,000 peptidoforms including >33,000 phosphorylated, acetylated, and ubiquitinated peptides for >9,000 native proteins. This site-resolved protein turnover (SPOT) profiling discloses global and site-specific differences in turnover associated with the presence or absence of PTMs. While causal relationships may not always be immediately apparent, we speculate that PTMs with diverging turnover may distinguish states of differential protein stability, structure, localization, enzymatic activity, or protein-protein interactions. We show examples of how the turnover data may give insights into unknown functions of PTMs and provide a freely accessible online tool that allows interrogation and visualisation of all turnover data. The SPOT methodology is applicable to many cell types and modifications, offering the potential to prioritize PTMs for future functional investigations.

Post-translational modifications (PTMs) can regulate cellular protein function but their global impact on protein turnover is largely unknown. Here, the authors develop proteomic workflows to profile PTM-resolved protein turnover and analyze the effects of phosphorylation, acetylation and ubiquitination.

Harnessing protein folding neural networks for peptide-protein docking

Highly accurate protein structure predictions by deep neural networks such as AlphaFold2 and RoseTTAFold have tremendous impact on structural biology and beyond. Here, we show that, although these deep learning approaches have originally been developed for the in silico folding of protein monomers, AlphaFold2 also enables quick and accurate modeling of peptide-protein interactions. Our simple implementation of AlphaFold2 generates peptide-protein complex models without requiring multiple sequence alignment information for the peptide partner, and can handle binding-induced conformational changes of the receptor. We explore what AlphaFold2 has memorized and learned, and describe specific examples that highlight differences compared to state-of-the-art peptide docking protocol PIPER-FlexPepDock. These results show that AlphaFold2 holds great promise for providing structural insight into a wide range of peptide-protein complexes, serving as a starting point for the detailed characterization and manipulation of these interactions.

An issue of concern: unique truncated ORF8 protein variants of SARS-CoV-2

Open reading frame 8 (ORF8) shows one of the highest levels of variability among accessory proteins in Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2), the causative agent of Coronavirus Disease 2019 (COVID-19). It was previously reported that the ORF8 protein inhibits the presentation of viral antigens by the major histocompatibility complex class I (MHC-I), which interacts with host factors involved in pulmonary inflammation. The ORF8 protein assists SARS-CoV-2 in evading immunity and plays a role in SARS-CoV-2 replication. Among many contributing mutations, Q27STOP, a mutation in the ORF8 protein, defines the B.1.1.7 lineage of SARS-CoV-2, engendering the second wave of COVID-19. In the present study, 47 unique truncated ORF8 proteins (T-ORF8) with the Q27STOP mutations were identified among 49,055 available B.1.1.7 SARS-CoV-2 sequences. The results show that only one of the 47 T-ORF8 variants spread to over 57 geo-locations in North America, and other continents, which include Africa, Asia, Europe and South America. Based on various quantitative features, such as amino acid homology, polar/non-polar sequence homology, Shannon entropy conservation, and other physicochemical properties of all specific 47 T-ORF8 protein variants, nine possible T-ORF8 unique variants were defined. The question as to whether T-ORF8 variants function similarly to the wild type ORF8 is yet to be investigated. A positive response to the question could exacerbate future COVID-19 waves, necessitating severe containment measures.

Computational Analysis of Alternative Splicing Using VAST-TOOLS and the VastDB Framework

Alternative splicing (AS) can vastly expand animal transcriptomes and proteomes. Two main open questions in the field are how AS is regulated across cell/tissue types and disease, and what roles different AS events play. To facilitate AS research, we have created the computational VastDB framework, which comprises a series of complementary software and resources that we describe in this chapter. The VastDB framework is especially designed to aid biomedical researchers without a strong computational background. It offers tools and resources to: (a) quantify AS and identify differentially spliced AS events using RNA-seq data (vast-tools), (b) perform multiple genomic and sequence analyses for investigating AS events (Matt), (c) identify AS events with genomic and regulatory conservation among species (ExOrthist), and (d) help with the biological interpretation of the results, and, ultimately, with the identification of interesting AS events to design wet-lab experiments (VastDB and PastDB).

Enhancing the Discovery of Functional Post-Translational Modification Sites with Machine Learning Models - Development, Validation, and Interpretation

Protein posttranslational modifications (PTMs) are a rapidly expanding feature class of significant importance in cell biology. Due to a high burden of experimental proof, the number of functionals PTMs in the eukaryotic proteome is currently underestimated. Furthermore, not all PTMs are functionally equivalent. Computational approaches that can confidently recommend PTMs of probable function can improve the heuristics of PTM investigation and alleviate these problems. To address this need, we developed SAPH-ire: a multifeature heuristic neural network model that takes community wisdom into account by recommending experimental PTMs similar to those which have previously been established as having regulatory impact. Here, we describe the principle behind the SAPH-ire model, how it is developed, how we evaluate its performance, and important caveats to consider when building and interpreting such models. Finally, we discus current limitations of functional PTM prediction models and highlight potential mechanisms for their improvement.

Characterizing the Interaction between the HTLV-1 Transactivator Tax-1 with Transcription Elongation Factor ELL2 and Its Impact on Viral Transactivation

The human T-cell leukemia virus type 1 (HTLV-1)-encoded transactivator and oncoprotein Tax-1 is essential for HTLV-1 replication. We recently found that Tax-1 interacts with transcription elongation factor for RNA polymerase II 2, ELL2, which enhances Tax-1-mediated transactivation of the HTLV-1 promotor. Here, we characterize the Tax-1:ELL2 interaction and its impact on viral transactivation by confocal imaging, co-immunoprecipitation, and luciferase assays. We found that Tax-1 and ELL2 not only co-precipitate, but also co-localize in dot-like structures in the nucleus. Tax-1:ELL2 complex formation occurred independently of Tax-1 point mutations, which are crucial for post translational modifications (PTMs) of Tax-1, suggesting that these PTMs are irrelevant for Tax-1:ELL2 interaction. In contrast, Tax-1 deletion mutants lacking either N-terminal (aa 1-37) or C-terminal regions (aa 150-353) of Tax-1 were impaired in interacting with ELL2. Contrary to Tax-1, the related, non-oncogenic Tax-2B from HTLV-2B did not interact with ELL2. Finally, we found that ELL2-R1 (aa 1-353), which carries an RNA polymerase II binding domain, and ELL2-R3 (aa 515-640) are sufficient to interact with Tax-1; however, only ELL2-truncations expressing R1 could enhance Tax-1-mediated transactivation of the HTLV-1 promoter. Together, this study identifies domains in Tax-1 and ELL2 being required for Tax-1:ELL2 complex formation and for viral transactivation.

Functional enrichment of alternative splicing events with NEASE reveals insights into tissue identity and diseases

Alternative splicing (AS) is an important aspect of gene regulation. Nevertheless, its role in molecular processes and pathobiology is far from understood. A roadblock is that tools for the functional analysis of AS-set events are lacking. To mitigate this, we developed NEASE, a tool integrating pathways with structural annotations of protein-protein interactions to functionally characterize AS events. We show in four application cases how NEASE can identify pathways contributing to tissue identity and cell type development, and how it highlights splicing-related biomarkers. With a unique view on AS, NEASE generates unique and meaningful biological insights complementary to classical pathways analysis.

14-3-3ζ negatively regulates mitochondrial biogenesis in GBM residual cells

Glioblastoma (GBM) is the most lethal primary brain tumour with a median survival of only 15 months. We have previously demonstrated the generation of an in vitro therapy resistance model that captures the residual resistant (RR) disease cells of GBM post-radiation. We also reported the proteomic landscape of parent, residual, and relapse cells using iTRAQ based quantitative proteomics of glioma cells. The proteomics data revealed significant up-regulation (fold change >1.5) of 14-3-3ζ, specifically in GBM RR cells. This was further confirmed by western blots in residual cells generated from GBM cell lines and patient sample-derived short-term primary culture. ShRNA-mediated knockdown of 14-3-3ζ radio-sensitized GBM cells and further stimulated therapy-induced senescence (TIS) and multinucleated giant cells (MNGCs) phenotype in RR cells. Intriguingly, 14-3-3ζ knockdown residual cells also showed a significantly higher number of mitochondria and increased mtDNA content. Indeed, in vitro GST pull-down mass spectrometry analysis of GST tagged 14-3-3ζ from RR cells identified novel interacting partners of 14-3-3ζ involved in cellular metabolism. Taken together, here we identified novel interacting partners of 14-3-3ζ and proposed an unconventional function of 14-3-3ζ as a negative regulator of TIS and mitochondrial biogenesis in residual resistant cells and loss of which also radio-sensitize GBM cells.

•

14-3-3ζ is up-regulated in residual disease cells of GBM.

•

14-3-3ζ knockdown radiosensitizes GBM cells.

•

14-3-3ζ knockdown increases MNGCs formation and senescence in residual cells.

•

14-3-3ζ negatively regulates mitochondrial biogenesis of residual disease cells.

•

Novel interacting partners of 14-3-3ζ from residual cells are involved in cellular metabolism.

Dynamic, but Not Necessarily Disordered, Human-Virus Interactions Mediated through SLiMs in Viral Proteins

Most viruses have small genomes that encode proteins needed to perform essential enzymatic functions. Across virus families, primary enzyme functions are under functional constraint; however, secondary functions mediated by exposed protein surfaces that promote interactions with the host proteins may be less constrained. Viruses often form transient interactions with host proteins through conformationally flexible interfaces. Exposed flexible amino acid residues are known to evolve rapidly suggesting that secondary functions may generate diverse interaction potentials between viruses within the same viral family. One mechanism of interaction is viral mimicry through short linear motifs (SLiMs) that act as functional signatures in host proteins. Viral SLiMs display specific patterns of adjacent amino acids that resemble their host SLiMs and may occur by chance numerous times in viral proteins due to mutational and selective processes. Through mimicry of SLiMs in the host cell proteome, viruses can interfere with the protein interaction network of the host and utilize the host-cell machinery to their benefit. The overlap between rapidly evolving protein regions and the location of functionally critical SLiMs suggest that these motifs and their functional potential may be rapidly rewired causing variation in pathogenicity, infectivity, and virulence of related viruses. The following review provides an overview of known viral SLiMs with select examples of their role in the life cycle of a virus, and a discussion of the structural properties of experimentally validated SLiMs highlighting that a large portion of known viral SLiMs are devoid of predicted intrinsic disorder based on the viral SLiMs from the ELM database.

Subcellular Localization of Seed-Expressed LEA_4 Proteins Reveals Liquid-Liquid Phase Separation for LEA9 and for LEA48 Homo- and LEA42-LEA48 Heterodimers

LEA proteins are involved in plant stress tolerance. In Arabidopsis, the LEA_4 Pfam group is the biggest group with the majority of its members being expressed in dry seeds. To assess subcellular localization in vivo, we investigated 11 seed-expressed LEA_4 proteins in embryos dissected from dry seeds expressing LEA_4 fusion proteins under its native promoters with the Venus fluorescent protein (proLEA_4::LEA_4:Venus). LEA_4 proteins were shown to be localized in the endoplasmic reticulum, nucleus, mitochondria, and plastids. LEA9, in addition to the nucleus, was also found in cytoplasmic condensates in dry seeds dependent on cellular hydration level. Most investigated LEA_4 proteins were detected in 4-d-old seedlings. In addition, we assessed bioinformatic tools for predicting subcellular localization and promoter motifs of 11 seed-expressed LEA_4 proteins. Ratiometric bimolecular fluorescence complementation assays showed that LEA7, LEA29, and LEA48 form homodimers while heterodimers were formed between LEA7-LEA29 and LEA42-LEA48 in tobacco leaves. Interestingly, LEA48 homodimers and LEA42-LEA48 heterodimers formed droplets structures with liquid-like behavior. These structures, along with LEA9 cytoplasmic condensates, may have been formed through liquid-liquid phase separation. These findings suggest possible important roles of LLPS for LEA protein functions.

Genetic fusions favor tumorigenesis through degron loss in oncogenes

Chromosomal rearrangements can generate genetic fusions composed of two distinct gene sequences, many of which have been implicated in tumorigenesis and progression. Our study proposes a model whereby oncogenic gene fusions frequently alter the protein stability of the resulting fusion products, via exchanging protein degradation signal (degron) between gene sequences. Computational analyses of The Cancer Genome Atlas (TCGA) identify 2,406 cases of degron exchange events and reveal an enrichment of oncogene stabilization due to loss of degrons from fusion. Furthermore, we identify and experimentally validate that some recurrent fusions, such as BCR-ABL, CCDC6-RET and PML-RARA fusions, perturb protein stability by exchanging internal degrons. Likewise, we also validate that EGFR or RAF1 fusions can be stabilized by losing a computationally-predicted C-terminal degron. Thus, complementary to enhanced oncogene transcription via promoter swapping, our model of degron loss illustrates another general mechanism for recurrent fusion proteins in driving tumorigenesis.

Can the SARS-CoV-2 Spike Protein Bind Integrins Independent of the RGD Sequence?

The RGD motif in the Severe Acute Syndrome Coronavirus 2 (SARS-CoV-2) spike protein has been predicted to bind RGD-recognizing integrins. Recent studies have shown that the spike protein does, indeed, interact with αVβ3 and α5β1 integrins, both of which bind to RGD-containing ligands. However, computational studies have suggested that binding between the spike RGD motif and integrins is not favourable, even when unfolding occurs after conformational changes induced by binding to the canonical host entry receptor, angiotensin-converting enzyme 2 (ACE2). Furthermore, non-RGD-binding integrins, such as αx, have been suggested to interact with the SARS-CoV-2 spike protein. Other viral pathogens, such as rotaviruses, have been recorded to bind integrins in an RGD-independent manner to initiate host cell entry. Thus, in order to consider the potential for the SARS-CoV-2 spike protein to bind integrins independent of the RGD sequence, we investigate several factors related to the involvement of integrins in SARS-CoV-2 infection. First, we review changes in integrin expression during SARS-CoV-2 infection to identify which integrins might be of interest. Then, all known non-RGD integrin-binding motifs are collected and mapped to the spike protein receptor-binding domain and analyzed for their 3D availability. Several integrin-binding motifs are shown to exhibit high sequence similarity with solvent accessible regions of the spike receptor-binding domain. Comparisons of these motifs with other betacoronavirus spike proteins, such as SARS-CoV and RaTG13, reveal that some have recently evolved while others are more conserved throughout phylogenetically similar betacoronaviruses. Interestingly, all of the potential integrin-binding motifs, including the RGD sequence, are conserved in one of the known pangolin coronavirus strains. Of note, the most recently recorded mutations in the spike protein receptor-binding domain were found outside of the putative integrin-binding sequences, although several mutations formed inside and close to one motif, in particular, may potentially enhance binding. These data suggest that the SARS-CoV-2 spike protein may interact with integrins independent of the RGD sequence and may help further explain how SARS-CoV-2 and other viruses can evolve to bind to integrins.

Expanding the Disorder-Function Paradigm in the C-Terminal Tails of Erbbs

ErbBs are receptor tyrosine kinases involved not only in development, but also in a wide variety of diseases, particularly cancer. Their extracellular, transmembrane, juxtamembrane, and kinase folded domains were described extensively over the past 20 years, structurally and functionally. However, their whole C-terminal tails (CTs) following the kinase domain were only described at atomic resolution in the last 4 years. They were shown to be intrinsically disordered. The CTs are known to be tyrosine-phosphorylated when the activated homo- or hetero-dimers of ErbBs are formed. Their phosphorylation triggers interaction with phosphotyrosine binding (PTB) or Src Homology 2 (SH2) domains and activates several signaling pathways controling cellular motility, proliferation, adhesion, and apoptosis. Beyond this passive role of phosphorylated domain and site display for partners, recent structural and function studies unveiled active roles in regulation of phosphorylation and interaction: the CT regulates activity of the kinase domain; different phosphorylation states have different compaction levels, potentially modulating the succession of phosphorylation events; and prolines have an important role in structure, dynamics, and possibly regulatory interactions. Here, we review both the canonical role of the disordered CT domains of ErbBs as phosphotyrosine display domains and the recent findings that expand the known range of their regulation functions linked to specific structural and dynamic features.

CYP11B1 variants influence skeletal maturation via alternative splicing

We performed genome-wide association study meta-analysis to identify genetic determinants of skeletal age (SA) deviating in multiple growth disorders. The joint meta-analysis (N = 4557) in two multiethnic cohorts of school-aged children identified one locus, CYP11B1 (expression confined to the adrenal gland), robustly associated with SA (rs6471570-A; β = 0.14; P = 6.2 × 10^-12). rs6410 (a synonymous variant in the first exon of CYP11B1 in high LD with rs6471570), was prioritized for functional follow-up being second most significant and the one closest to the first intron-exon boundary. In 208 adrenal RNA-seq samples from GTEx, C-allele of rs6410 was associated with intron 3 retention (P = 8.11 × 10^-40), exon 4 inclusion (P = 4.29 × 10^-34), and decreased exon 3 and 5 splicing (P = 7.85 × 10^-43), replicated using RT-PCR in 15 adrenal samples. As CYP11B1 encodes 11-β-hydroxylase, involved in adrenal glucocorticoid and mineralocorticoid biosynthesis, our findings highlight the role of adrenal steroidogenesis in SA in healthy children, suggesting alternative splicing as a likely underlying mechanism.

In-Depth Annotation of the Drosophila Bithorax-Complex Reveals the Presence of Several Alternative ORFs That Could Encode for Motif-Rich Peptides

It is recognized that a large proportion of eukaryotic RNAs and proteins is not produced from conventional genes but from short and alternative (alt) open reading frames (ORFs) that are not captured by gene prediction programs. Here we present an in silico prediction of altORFs by applying several selecting filters based on evolutionary conservation and annotations of previously characterized altORF peptides. Our work was performed in the Bithorax-complex (BX-C), which was one of the first genomic regions described to contain long non-coding RNAs in Drosophila. We showed that several altORFs could be predicted from coding and non-coding sequences of BX-C. In addition, the selected altORFs encode for proteins that contain several interesting molecular features, such as the presence of transmembrane helices or a general propensity to be rich in short interaction motifs. Of particular interest, one altORF encodes for a protein that contains a peptide sequence found in specific isoforms of two Drosophila Hox proteins. Our work thus suggests that several altORF proteins could be produced from a particular genomic region known for its critical role during Drosophila embryonic development. The molecular signatures of these altORF proteins further suggests that several of them could make numerous protein–protein interactions and be of functional importance in vivo.

Investigation of Evolutionary History and Origin of the Tre1 Family Suggests a Role in Regulating Hemocytes Cells Infiltration of the Blood-Brain Barrier

Simple Summary

Understanding the evolutionary association between immune cells and the blood–brain barrier (BBB) is vital to develop therapeutic approaches. In Drosophila, glial cells form the BBB that regulates the access of hemocytes to the brain. It is still not known which diapedesis route hemocytes cells follow. In vertebrates, paracellular migration is dependent on PECAM1, while transcellular migration is dependent on the expression of CAV1. The drosophila genome lacks both genes. The Tre1 family (Tre1, moody, and Dmel_CG4313) contribute to regulating transepithelial migration in Drosophila. However, its evolutionary history is not known. We performed phylogenetic analysis to reconstruct the evolutionary history of the Tre1 family. We found Dmel_CG4313 only in insects. Tre1 exists only in invertebrates and is highly conserved. moody evolutionary history is more spread as it appears from Cnidaria up to mammals and is less conserved. The Tre1 family origin seems to be related to opsins. We have identified an SH3 motif in Tre1, moody, and Dmel_CG4313. SH3 regulates actin movement in a Rho-dependent manner in PECAM1. Our results suggest that the Tre1 family could be playing an important role in paracellular diapedesis in Drosophila. Thus, targeting the Tre1 family could help us regulate access to the brain.

Abstract

Understanding the evolutionary relationship between immune cells and the blood–brain barrier (BBB) is important to devise therapeutic strategies. In vertebrates, immune cells follow either a paracellular or a transcellular pathway to infiltrate the BBB. In Drosophila, glial cells form the BBB that regulates the access of hemocytes to the brain. However, it is still not known which diapedesis route hemocytes cells follow. In vertebrates, paracellular migration is dependent on PECAM1, while transcellular migration is dependent on the expression of CAV1. Interestingly Drosophila genome lacks both genes. Tre1 family (Tre1, moody, and Dmel_CG4313) play a diverse role in regulating transepithelial migration in Drosophila. However, its evolutionary history and origin are not yet known. We performed phylogenetic analysis, together with HH search, positive selection, and ancestral reconstruction to investigate the Tre1 family. We found that Tre1 exists in Mollusca, Arthropoda, Ambulacraria, and Scalidophora. moody is shown to be a more ancient protein and it has existed since Cnidaria emergence and has a homolog (e.g., GPCR84) in mammals. The third family member (Dmel_CG4313) seems to only exist in insects. The origin of the family seems to be related to the rhodopsin-like family and in particular family α. We found that opsin is the nearest receptor to have a common ancestor with the Tre1 family that has diverged in sponges. We investigated the positive selection of the Tre1 family using PAML. Tre1 seems to have evolved under negative selection, whereas moody has evolved during positive selection. The sites that we found under positive selection are likely to play a role in the speciation of function in the case of moody. We have identified an SH3 motif, in Tre1 and, moody and Dmel_CG4313. SH3 is known to play a fundamental role in regulating actin movement in a Rho-dependent manner in PECAM1. Our results suggest that the Tre1 family could be playing an important role in paracellular diapedesis in Drosophila.

Enhancer recruitment of transcription repressors RUNX1 and TLE3 by mis-expressed FOXC1 blocks differentiation in acute myeloid leukemia

Despite absent expression in normal hematopoiesis, the Forkhead factor FOXC1, a critical mesenchymal differentiation regulator, is highly expressed in ∼30% of HOXAhigh acute myeloid leukemia (AML) cases to confer blocked monocyte/macrophage differentiation. Through integrated proteomics and bioinformatics, we find that FOXC1 and RUNX1 interact through Forkhead and Runt domains, respectively, and co-occupy primed and active enhancers distributed close to differentiation genes. FOXC1 stabilizes association of RUNX1, HDAC1, and Groucho repressor TLE3 to limit enhancer activity: FOXC1 knockdown induces loss of repressor proteins, gain of CEBPA binding, enhancer acetylation, and upregulation of nearby genes, including KLF2. Furthermore, it triggers genome-wide redistribution of RUNX1, TLE3, and HDAC1 from enhancers to promoters, leading to repression of self-renewal genes, including MYC and MYB. Our studies highlight RUNX1 and CEBPA transcription factor swapping as a feature of leukemia cell differentiation and reveal that FOXC1 prevents this by stabilizing enhancer binding of a RUNX1/HDAC1/TLE3 transcription repressor complex to oncogenic effect.

Detecting Rewiring Events in Protein-Protein Interaction Networks Based on Transcriptomic Data

Proteins rarely carry out their cellular functions in isolation. Instead, eukaryotic proteins engage in about six interactions with other proteins on average. The aggregated protein interactome of an organism forms a “hairy ball”-type protein-protein interaction (PPI) network. Yet, in a typical human cell, only about half of all proteins are expressed at a particular time. Hence, it has become common practice to prune the full PPI network to the subset of expressed proteins. If RNAseq data is available, one can further resolve the specific protein isoforms present in a cell or tissue. Here, we review various approaches, software tools and webservices that enable users to construct context-specific or tissue-specific PPI networks and how these are rewired between two cellular conditions. We illustrate their different functionalities on the example of the interactions involving the human TNR6 protein. In an outlook, we describe how PPI networks may be integrated with epigenetic data or with data on the activity of splicing factors.

On the roles of intrinsically disordered proteins and regions in cell communication and signaling

For proteins, the sequence → structure → function paradigm applies primarily to enzymes, transmembrane proteins, and signaling domains. This paradigm is not universal, but rather, in addition to structured proteins, intrinsically disordered proteins and regions (IDPs and IDRs) also carry out crucial biological functions. For these proteins, the sequence → IDP/IDR ensemble → function paradigm applies primarily to signaling and regulatory proteins and regions. Often, in order to carry out function, IDPs or IDRs cooperatively interact, either intra- or inter-molecularly, with structured proteins or other IDPs or intermolecularly with nucleic acids. In this IDP/IDR thematic collection published in Cell Communication and Signaling, thirteen articles are presented that describe IDP/IDR signaling molecules from a variety of organisms from humans to fruit flies and tardigrades ("water bears") and that describe how these proteins and regions contribute to the function and regulation of cell signaling. Collectively, these papers exhibit the diverse roles of disorder in responding to a wide range of signals as to orchestrate an array of organismal processes. They also show that disorder contributes to signaling in a broad spectrum of species, ranging from micro-organisms to plants and animals.

CEACAM Gene Family Mutations Associated With Inherited Breast Cancer Risk - A Comparative Oncology Approach to Discovery

Introduction

Recent studies comparing canine mammary tumors (CMTs) and human breast cancers have revealed remarkable tumor similarities, identifying shared expression profiles and acquired mutations. CMTs can also provide a model of inherited breast cancer susceptibility in humans; thus, we investigated breed-specific whole genome sequencing (WGS) data in search for novel CMT risk factors that could subsequently explain inherited breast cancer risk in humans.

Methods

WGS was carried out on five CMT-affected Gold Retrievers from a large pedigree of 18 CMT-affected dogs. Protein truncating variants (PTVs) detected in all five samples (within human orthlogs) were validated and then genotyped in the 13 remaining CMT-affected Golden Retrievers. Allele frequencies were compared to canine controls. Subsequently, human blood-derived exomes from The Cancer Genome Atlas breast cancer cases were analyzed and allele frequencies were compared to Exome Variant Server ethnic-matched controls.

Results

Carcinoembryonic Antigen-related Cell Adhesion Molecule 24 (CEACAM24) c.247dupG;p.(Val83Glyfs^∗48) was the only validated variant and had a frequency of 66.7% amongst the 18 Golden Retrievers with CMT. This was significant compared to the European Variation Archive (p-value 1.52 × 10^–8) and non-Golden Retriever American Kennel Club breeds (p-value 2.48 × 10^–5). With no direct ortholog of CEACAM24 in humans but high homology to all CEACAM gene family proteins, all human CEACAM genes were investigated for PTVs. A total of six and sixteen rare PTVs were identified in African and European American breast cancer cases, respectively. Single variant assessment revealed five PTVs associated with breast cancer risk. Gene-based aggregation analyses revealed that rare PTVs in CEACAM6, CEACAM7, and CEACAM8 are associated with European American breast cancer risk, and rare PTVs in CEACAM7 are associated with breast cancer risk in African Americans. Ultimately, rare PTVs in the entire CEACAM gene family are associated with breast cancer risk in both European and African Americans with respective p-values of 1.75 × 10^–13 and 1.87 × 10^–04.

Conclusion

This study reports the first association of inherited CEACAM mutations and breast cancer risk, and potentially implicates the whole gene family in genetic risk. Precisely how these mutations contribute to breast cancer needs to be determined; especially considering our current knowledge on the role that the CEACAM gene family plays in tumor development, progression, and metastasis.

Diversification of CD1 Molecules Shapes Lipid Antigen Selectivity

Molecular studies of host-pathogen evolution have largely focused on the consequences of variation at protein-protein interaction surfaces. The potential for other microbe-associated macromolecules to promote arms race dynamics with host factors remains unclear. The cluster of differentiation 1 (CD1) family of vertebrate cell surface receptors plays a crucial role in adaptive immunity through binding and presentation of lipid antigens to T-cells. Although CD1 proteins present a variety of endogenous and microbial lipids to various T-cell types, they are less diverse within vertebrate populations than the related major histocompatibility complex (MHC) molecules. We discovered that CD1 genes exhibit a high level of divergence between simian primate species, altering predicted lipid-binding properties and T-cell receptor interactions. These findings suggest that lipid-protein conflicts have shaped CD1 genetic variation during primate evolution. Consistent with this hypothesis, multiple primate CD1 family proteins exhibit signatures of repeated positive selection at surfaces impacting antigen presentation, binding pocket morphology, and T-cell receptor accessibility. Using a molecular modeling approach, we observe that interspecies variation as well as single mutations at rapidly-evolving sites in CD1a drastically alter predicted lipid binding and structural features of the T-cell recognition surface. We further show that alterations in both endogenous and microbial lipid-binding affinities influence the ability of CD1a to undergo antigen swapping required for T-cell activation. Together these findings establish lipid-protein interactions as a critical force of host-pathogen conflict and inform potential strategies for lipid-based vaccine development.

Delineating the molecular and phenotypic spectrum of the SETD1B-related syndrome

Purpose

Pathogenic variants in SETD1B have been associated with a syndromic neurodevelopmental disorder including intellectual disability, language delay, and seizures. To date, clinical features have been described for 11 patients with (likely) pathogenic SETD1B sequence variants. This study aims to further delineate the spectrum of the SETD1B-related syndrome based on characterizing an expanded patient cohort.

Methods

We perform an in-depth clinical characterization of a cohort of 36 unpublished individuals with SETD1B sequence variants, describing their molecular and phenotypic spectrum. Selected variants were functionally tested using in vitro and genome-wide methylation assays.

Results

Our data present evidence for a loss-of-function mechanism of SETD1B variants, resulting in a core clinical phenotype of global developmental delay, language delay including regression, intellectual disability, autism and other behavioral issues, and variable epilepsy phenotypes. Developmental delay appeared to precede seizure onset, suggesting SETD1B dysfunction impacts physiological neurodevelopment even in the absence of epileptic activity. Males are significantly overrepresented and more severely affected, and we speculate that sex-linked traits could affect susceptibility to penetrance and the clinical spectrum of SETD1B variants.

Conclusion

Insights from this extensive cohort will facilitate the counseling regarding the molecular and phenotypic landscape of newly diagnosed patients with the SETD1B-related syndrome.

Matrix-screening reveals a vast potential for direct protein-protein interactions among RNA binding proteins

RNA-binding proteins (RBPs) are crucial factors of post-transcriptional gene regulation and their modes of action are intensely investigated. At the center of attention are RNA motifs that guide where RBPs bind. However, sequence motifs are often poor predictors of RBP-RNA interactions in vivo. It is hence believed that many RBPs recognize RNAs as complexes, to increase specificity and regulatory possibilities. To probe the potential for complex formation among RBPs, we assembled a library of 978 mammalian RBPs and used rec-Y2H matrix screening to detect direct interactions between RBPs, sampling > 600 K interactions. We discovered 1994 new interactions and demonstrate that interacting RBPs bind RNAs adjacently in vivo. We further find that the mRNA binding region and motif preferences of RBPs deviate, depending on their adjacently binding interaction partners. Finally, we reveal novel RBP interaction networks among major RNA processing steps and show that splicing impairing RBP mutations observed in cancer rewire spliceosomal interaction networks. The dataset we provide will be a valuable resource for understanding the combinatorial interactions of RBPs with RNAs and the resulting regulatory outcomes.

Light-Controllable PROTACs for Temporospatial Control of Protein Degradation

PROteolysis-TArgeting Chimeras (PROTACs) is an emerging and promising approach to target intracellular proteins for ubiquitination-mediated degradation, including those so-called undruggable protein targets, such as transcriptional factors and scaffold proteins. To date, plenty of PROTACs have been developed to degrade various disease-relevant proteins, such as estrogen receptor (ER), androgen receptor (AR), RTK, and CDKs. However, the on-target off-tissue and off-target effect is one of the major limitation that prevents the usage of PROTACs in clinic. To this end, we and several other groups have recently developed light-controllable PROTACs, as the representative for the third generation controllable PROTACs, by using either photo-caging or photo-switch approaches. In this review, we summarize the emerging light-controllable PROTACs and the prospective for other potential ways to achieve temporospatial control of PROTACs.

Computational Phosphorylation Network Reconstruction: An Update on Methods and Resources

Most proteins undergo some form of modification after translation, and phosphorylation is one of the most relevant and ubiquitous post-translational modifications. The succession of protein phosphorylation and dephosphorylation catalyzed by protein kinase and phosphatase, respectively, constitutes a key mechanism of molecular information flow in cellular systems. The protein interactions of kinases, phosphatases, and their regulatory subunits and substrates are the main part of phosphorylation networks. To elucidate the landscape of phosphorylation events has been a central goal pursued by both experimental and computational approaches. Substrate specificity (e.g., sequence, structure) or the phosphoproteome has been utilized in an array of different statistical learning methods to infer phosphorylation networks. In this chapter, different computational phosphorylation network inference-related methods and resources are summarized and discussed.

Predicted Cellular Interactors of the Endogenous Retrovirus-K Integrase Enzyme

Integrase (IN) enzymes are found in all retroviruses and are crucial in the retroviral integration process. Many studies have revealed how exogenous IN enzymes, such as the human immunodeficiency virus (HIV) IN, contribute to altered cellular function. However, the same consideration has not been given to viral IN originating from symbionts within our own DNA. Endogenous retrovirus-K (ERVK) is pathologically associated with neurological and inflammatory diseases along with several cancers. The ERVK IN interactome is unknown, and the question of how conserved the ERVK IN protein-protein interaction motifs are as compared to other retroviral integrases is addressed in this paper. The ERVK IN protein sequence was analyzed using the Eukaryotic Linear Motif (ELM) database, and the results are compared to ELMs of other betaretroviral INs and similar eukaryotic INs. A list of putative ERVK IN cellular protein interactors was curated from the ELM list and submitted for STRING analysis to generate an ERVK IN interactome. KEGG analysis was used to identify key pathways potentially influenced by ERVK IN. It was determined that the ERVK IN potentially interacts with cellular proteins involved in the DNA damage response (DDR), cell cycle, immunity, inflammation, cell signaling, selective autophagy, and intracellular trafficking. The most prominent pathway identified was viral carcinogenesis, in addition to select cancers, neurological diseases, and diabetic complications. This potentiates the role of ERVK IN in these pathologies via protein-protein interactions facilitating alterations in key disease pathways.

IUPred3: prediction of protein disorder enhanced with unambiguous experimental annotation and visualization of evolutionary conservation

Intrinsically disordered proteins and protein regions (IDPs/IDRs) exist without a single well-defined conformation. They carry out important biological functions with multifaceted roles which is also reflected in their evolutionary behavior. Computational methods play important roles in the characterization of IDRs. One of the commonly used disorder prediction methods is IUPred, which relies on an energy estimation approach. The IUPred web server takes an amino acid sequence or a Uniprot ID/accession as an input and predicts the tendency for each amino acid to be in a disordered region with an option to also predict context-dependent disordered regions. In this new iteration of IUPred, we added multiple novel features to enhance the prediction capabilities of the server. First, learning from the latest evaluation of disorder prediction methods we introduced multiple new smoothing functions to the prediction that decreases noise and increases the performance of the predictions. We constructed a dataset consisting of experimentally verified ordered/disordered regions with unambiguous annotations which were added to the prediction. We also introduced a novel tool that enables the exploration of the evolutionary conservation of protein disorder coupled to sequence conservation in model organisms. The web server is freely available to users and accessible at https://iupred3.elte.hu.

Mechnetor: a web server for exploring protein mechanism and the functional context of genetic variants

Advances in DNA sequencing and proteomics mean that researchers must now regularly interrogate thousands of positional gene/protein changes in order to find those relevant for potential clinical application or biological insights. The abundance of already known information on protein interactions, mechanism, and tertiary structure provides the possible means to understand these changes rapidly, though a careful and systematic integration of these diverse datasets is first needed. For this purpose, we developed Mechnetor, a tool that allows users to quickly explore and visualize integrated mechanistic data for proteins or interactions of interest. Central to the system is a careful cataloguing of diverse sources of protein interaction mechanism, and an efficient means to visualize interactions between relevant and/or known protein regions. The result is a finer resolution interaction network that provides more immediate clues as to points of intervention or mechanistic understanding. Users can import protein, interactions, genetic variants or post-translational modifications and see these data in the best known mechanistic context. We demonstrate the tool with topical examples in human genetic diseases and cancer genomics. The tool is freely available at: mechnetor.russelllab.org.

P2T2: Protein Panoramic annoTation Tool for the interpretation of protein coding genetic variants

MOTIVATION

Genomic data are prevalent, leading to frequent encounters with uninterpreted variants or mutations with unknown mechanisms of effect. Researchers must manually aggregate data from multiple sources and across related proteins, mentally translating effects between the genome and proteome, to attempt to understand mechanisms.

MATERIALS AND METHODS

P2T2 presents diverse data and annotation types in a unified protein-centric view, facilitating the interpretation of coding variants and hypothesis generation. Information from primary sequence, domain, motif, and structural levels are presented and also organized into the first Paralog Annotation Analysis across the human proteome.

RESULTS

Our tool assists research efforts to interpret genomic variation by aggregating diverse, relevant, and proteome-wide information into a unified interactive web-based interface. Additionally, we provide a REST API enabling automated data queries, or repurposing data for other studies.

CONCLUSION

The unified protein-centric interface presented in P2T2 will help researchers interpret novel variants identified through next-generation sequencing. Code and server link available at github.com/GenomicInterpretation/p2t2.

On the Need to Tell Apart Fraternal Twins eEF1A1 and eEF1A2, and Their Respective Outfits

eEF1A1 and eEF1A2 are paralogous proteins whose presence in most normal eukaryotic cells is mutually exclusive and developmentally regulated. Often described in the scientific literature under the collective name eEF1A, which stands for eukaryotic elongation factor 1A, their best known activity (in a monomeric, GTP-bound conformation) is to bind aminoacyl-tRNAs and deliver them to the A-site of the 80S ribosome. However, both eEF1A1 and eEF1A2 are endowed with multitasking abilities (sometimes performed by homo- and heterodimers) and can be located in different subcellular compartments, from the plasma membrane to the nucleus. Given the high sequence identity of these two sister proteins and the large number of post-translational modifications they can undergo, we are often confronted with the dilemma of discerning which is the particular proteoform that is actually responsible for the ascribed biochemical or cellular effects. We argue in this review that acquiring this knowledge is essential to help clarify, in molecular and structural terms, the mechanistic involvement of these two ancestral and abundant G proteins in a variety of fundamental cellular processes other than translation elongation. Of particular importance for this special issue is the fact that several de novo heterozygous missense mutations in the human EEF1A2 gene are associated with a subset of rare but severe neurological syndromes and cardiomyopathies.

Ubiquitin Interacting Motifs: Duality Between Structured and Disordered Motifs

Ubiquitin is a small protein at the heart of many cellular processes, and several different protein domains are known to recognize and bind ubiquitin. A common motif for interaction with ubiquitin is the Ubiquitin Interacting Motif (UIM), characterized by a conserved sequence signature and often found in multi-domain proteins. Multi-domain proteins with intrinsically disordered regions mediate interactions with multiple partners, orchestrating diverse pathways. Short linear motifs for binding are often embedded in these disordered regions and play crucial roles in modulating protein function. In this work, we investigated the structural propensities of UIMs using molecular dynamics simulations and NMR chemical shifts. Despite the structural portrait depicted by X-crystallography of stable helical structures, we show that UIMs feature both helical and intrinsically disordered conformations. Our results shed light on a new class of disordered UIMs. This group is here exemplified by the C-terminal domain of one isoform of ataxin-3 and a group of ubiquitin-specific proteases. Intriguingly, UIMs not only bind ubiquitin. They can be a recruitment point for other interactors, such as parkin and the heat shock protein Hsc70-4. Disordered UIMs can provide versatility and new functions to the client proteins, opening new directions for research on their interactome.

The diversity of molecular interactions involving intrinsically disordered proteins: A molecular modeling perspective

Intrinsically Disordered Proteins and Regions (IDPs/IDRs) are key components of a multitude of biological processes. Conformational malleability enables IDPs/IDRs to perform very specialized functions that cannot be accomplished by globular proteins. The functional role for most of these proteins is related to the recognition of other biomolecules to regulate biological processes or as a part of signaling pathways. Depending on the extent of disorder, the number of interacting sites and the type of partner, very different architectures for the resulting assemblies are possible. More recently, molecular condensates with liquid-like properties composed of multiple copies of IDPs and nucleic acids have been proven to regulate key processes in eukaryotic cells. The structural and kinetic details of disordered biomolecular complexes are difficult to unveil experimentally due to their inherent conformational heterogeneity. Computational approaches, alone or in combination with experimental data, have emerged as unavoidable tools to understand the functional mechanisms of this elusive type of assemblies. The level of description used, all-atom or coarse-grained, strongly depends on the size of the molecular systems and on the timescale of the investigated mechanism. In this mini-review, we describe the most relevant architectures found for molecular interactions involving IDPs/IDRs and the computational strategies applied for their investigation.

Cellular Chaperone Function of Intrinsically Disordered Dehydrin ERD14

Disordered plant chaperones play key roles in helping plants survive in harsh conditions, and they are indispensable for seeds to remain viable. Aside from well-known and thoroughly characterized globular chaperone proteins, there are a number of intrinsically disordered proteins (IDPs) that can also serve as highly effective protecting agents in the cells. One of the largest groups of disordered chaperones is the group of dehydrins, proteins that are expressed at high levels under different abiotic stress conditions, such as drought, high temperature, or osmotic stress. Dehydrins are characterized by the presence of different conserved sequence motifs that also serve as the basis for their categorization. Despite their accepted importance, the exact role and relevance of the conserved regions have not yet been formally addressed. Here, we explored the involvement of each conserved segment in the protective function of the intrinsically disordered stress protein (IDSP) A. thaliana's Early Response to Dehydration (ERD14). We show that segments that are directly involved in partner binding, and others that are not, are equally necessary for proper function and that cellular protection emerges from the balanced interplay of different regions of ERD14.

Model-based analysis uncovers mutations altering autophagy selectivity in human cancer

Autophagy can selectively target protein aggregates, pathogens, and dysfunctional organelles for the lysosomal degradation. Aberrant regulation of autophagy promotes tumorigenesis, while it is far less clear whether and how tumor-specific alterations result in autophagic aberrance. To form a link between aberrant autophagy selectivity and human cancer, we establish a computational pipeline and prioritize 222 potential LIR (LC3-interacting region) motif-associated mutations (LAMs) in 148 proteins. We validate LAMs in multiple proteins including ATG4B, STBD1, EHMT2 and BRAF that impair their interactions with LC3 and autophagy activities. Using a combination of transcriptomic, metabolomic and additional experimental assays, we show that STBD1, a poorly-characterized protein, inhibits tumor growth via modulating glycogen autophagy, while a patient-derived W203C mutation on LIR abolishes its cancer inhibitory function. This work suggests that altered autophagy selectivity is a frequently-used mechanism by cancer cells to survive during various stresses, and provides a framework to discover additional autophagy-related pathways that influence carcinogenesis.

Although autophagy has been linked to tumourigenesis, it is unclear how genomic alterations affect autophagy selectivity in tumours. Here, the authors establish a pipeline that integrates computational and experimental approaches to show that altered autophagy selectivity is frequent in cancer cells and link glycogen autophagy with tumourigenesis.

Characterization of a Novel SARS-CoV-2 Genetic Variant with Distinct Spike Protein Mutations

The COVID-19 pandemic, which began in Wuhan (Hubei, China), has been ongoing for about a year and a half. An unprecedented number of people around the world have been infected with SARS-CoV-2, the etiological agent of COVID-19. Despite the fact that the mortality rate for COVID-19 is relatively low, the total number of deaths has currently already reached more than three million and continues to increase due to high incidence. Since the beginning of the pandemic, a large number of sequences have been obtained and many genetic variants have been identified. Some of them bear significant mutations that affect biological properties of the virus. These genetic variants, currently Variants of Concern (VoC), include the so-called United Kingdom variant (20I/501Y), the Brazilian variant (20J/501Y.V3), and the South African variant (20H/501Y.V2). We describe here a novel SARS-CoV-2 variant with distinct spike protein mutations, first obtained at the end of January 2021 in northwest Russia. Therefore, it is necessary to pay attention to the dynamics of its spread among patients with COVID-19, as well as to study in detail its biological properties.

Pepscan Approach for the Identification of Protein-Protein Interfaces: Lessons from Experiment

PEPscan is an old approach that has recently gained renewed interest for the identification of interfering peptides (IPs), i.e., peptides able to interfere with protein-protein interactions (PPIs). Its principle is to slice a protein sequence as a series of short overlapping peptides that are synthesized on a peptide array and tested for their ability to bind a partner, with positive spots corresponding to candidate IPs. PEPscan has been applied with a rather large success in various contexts, but the structural determinants underlying this success remain obscure. Here, we analyze the results of 14 PEPscan experiments, and confront the in vitro results with the available structural information. PEPscan identifies candidate IPs in limited numbers that in all cases correspond to solvent-accessible regions of the structures, their location at the protein-protein interface remaining to be further demonstrated. A strong point of PEPscan seems to be its ability to identify specific IPs. IPs identified from the same protein differ depending on the target PPI, and correspond to patches not frequently involved in the interactions seen in the 3D structures available. Overall, PEPscan seems to provide a cheap and rapid manner to identify candidate IPs, that also comes with room for improvement.

Structural Profiling of Bacterial Effectors Reveals Enrichment of Host-Interacting Domains and Motifs

Effector proteins are bacterial virulence factors secreted directly into host cells and, through extensive interactions with host proteins, rewire host signaling pathways to the advantage of the pathogen. Despite the crucial role of globular domains as mediators of protein-protein interactions (PPIs), previous structural studies of bacterial effectors are primarily focused on individual domains, rather than domain-mediated PPIs, which limits their ability to uncover systems-level molecular recognition principles governing host-bacteria interactions. Here, we took an interaction-centric approach and systematically examined the potential of structural components within bacterial proteins to engage in or target eukaryote-specific domain-domain interactions (DDIs). Our results indicate that: 1) effectors are about six times as likely as non-effectors to contain host-like domains that mediate DDIs exclusively in eukaryotes; 2) the average domain in effectors is about seven times as likely as that in non-effectors to co-occur with DDI partners in eukaryotes rather than in bacteria; and 3) effectors are about nine times as likely as non-effectors to contain bacteria-exclusive domains that target host domains mediating DDIs exclusively in eukaryotes. Moreover, in the absence of host-like domains or among pathogen proteins without domain assignment, effectors harbor a higher variety and density of short linear motifs targeting host domains that mediate DDIs exclusively in eukaryotes. Our study lends novel quantitative insight into the structural basis of effector-induced perturbation of host-endogenous PPIs and may aid in the design of selective inhibitors of host-pathogen interactions.

The Arp2/3 Inhibitory Protein Arpin Is Required for Intestinal Epithelial Barrier Integrity

The intestinal epithelial barrier (IEB) depends on stable interepithelial protein complexes such as tight junctions (TJ), adherens junctions (AJ), and the actin cytoskeleton. During inflammation, the IEB is compromised due to TJ protein internalization and actin remodeling. An important actin regulator is the actin-related protein 2/3 (Arp2/3) complex, which induces actin branching. Activation of Arp2/3 by nucleation-promoting factors is required for the formation of epithelial monolayers, but little is known about the relevance of Arp2/3 inhibition and endogenous Arp2/3 inhibitory proteins for IEB regulation. We found that the recently identified Arp2/3 inhibitory protein arpin was strongly expressed in intestinal epithelial cells. Arpin expression decreased in response to tumor necrosis factor (TNF)α and interferon (IFN)γ treatment, whereas the expression of gadkin and protein interacting with protein C-kinase α-subunit 1 (PICK1), other Arp2/3 inhibitors, remained unchanged. Of note, arpin coprecipitated with the TJ proteins occludin and claudin-1 and the AJ protein E-cadherin. Arpin depletion altered the architecture of both AJ and TJ, increased actin filament content and actomyosin contractility, and significantly increased epithelial permeability, demonstrating that arpin is indeed required for maintaining IEB integrity. During experimental colitis in mice, arpin expression was also decreased. Analyzing colon tissues from ulcerative colitis patients by Western blot, we found different arpin levels with overall no significant changes. However, in acutely inflamed areas, arpin was significantly reduced compared to non-inflamed areas. Importantly, patients receiving mesalazine had significantly higher arpin levels than untreated patients. As arpin depletion (theoretically meaning more active Arp2/3) increased permeability, we wanted to know whether Arp2/3 inhibition would show the opposite. Indeed, the specific Arp2/3 inhibitor CK666 ameliorated TNFα/IFNγ-induced permeability in established Caco-2 monolayers by preventing TJ disruption. CK666 treatment also attenuated colitis development, colon tissue damage, TJ disruption, and permeability in dextran sulphate sodium (DSS)-treated mice. Our results demonstrate that loss of arpin triggers IEB dysfunction during inflammation and that low arpin levels can be considered a novel hallmark of acute inflammation.

Thousands of protein linear motif classes may still be undiscovered

Linear motifs are short protein subsequences that mediate protein interactions. Hundreds of motif classes including thousands of motif instances are known. Our theory estimates how many motif classes remain undiscovered. As commonly done, we describe motif classes as regular expressions specifying motif length and the allowed amino acids at each motif position. We measure motif specificity for a pair of motif classes by quantifying how many motif-discriminating positions prevent a protein subsequence from matching the two classes at once. We derive theorems for the maximal number of motif classes that can simultaneously maintain a certain number of motif-discriminating positions between all pairs of classes in the motif universe, for a given amino acid alphabet. We also calculate the fraction of all protein subsequences that would belong to a motif class if all potential motif classes came into existence. Naturally occurring pairs of motif classes present most often a single motif-discriminating position. This mild specificity maximizes the potential number of coexisting motif classes, the expansion of the motif universe due to amino acid modifications and the fraction of amino acid sequences that code for a motif instance. As a result, thousands of linear motif classes may remain undiscovered.

Order from disorder in the sarcomere: FATZ forms a fuzzy but tight complex and phase-separated condensates with α-actinin

In sarcomeres, α-actinin cross-links actin filaments and anchors them to the Z-disk. FATZ (filamin-, α-actinin-, and telethonin-binding protein of the Z-disk) proteins interact with α-actinin and other core Z-disk proteins, contributing to myofibril assembly and maintenance. Here, we report the first structure and its cellular validation of α-actinin-2 in complex with a Z-disk partner, FATZ-1, which is best described as a conformational ensemble. We show that FATZ-1 forms a tight fuzzy complex with α-actinin-2 and propose an interaction mechanism via main molecular recognition elements and secondary binding sites. The obtained integrative model reveals a polar architecture of the complex which, in combination with FATZ-1 multivalent scaffold function, might organize interaction partners and stabilize α-actinin-2 preferential orientation in Z-disk. Last, we uncover FATZ-1 ability to phase-separate and form biomolecular condensates with α-actinin-2, raising the question whether FATZ proteins can create an interaction hub for Z-disk proteins through membraneless compartmentalization during myofibrillogenesis.

Synergy and allostery in ligand binding by HIV-1 Nef

The Nef protein of human and simian immunodeficiency viruses boosts viral pathogenicity through its interactions with host cell proteins. By combining the polyvalency of its large unstructured regions with the binding selectivity and strength of its folded core domain, Nef can associate with many different host cell proteins, thereby disrupting their functions. For example, the combination of a linear proline-rich motif and hydrophobic core domain surface allows Nef to bind tightly and specifically to SH3 domains of Src family kinases. We investigated whether the interplay between Nef's flexible regions and its core domain could allosterically influence ligand selection. We found that the flexible regions can associate with the core domain in different ways, producing distinct conformational states that alter the way in which Nef selects for SH3 domains and exposes some of its binding motifs. The ensuing crosstalk between ligands might promote functionally coherent Nef-bound protein ensembles by synergizing certain subsets of ligands while excluding others. We also combined proteomic and bioinformatics analyses to identify human proteins that select SH3 domains in the same way as Nef. We found that only 3% of clones from a whole-human fetal library displayed Nef-like SH3 selectivity. However, in most cases, this selectivity appears to be achieved by a canonical linear interaction rather than by a Nef-like ‘tertiary' interaction. Our analysis supports the contention that Nef's mode of hijacking SH3 domains is a virus-specific adaptation with no or very few cellular counterparts. Thus, the Nef tertiary binding surface is a promising virus-specific drug target.

Ehrlichia chaffeensis TRP120 Is a Wnt Ligand Mimetic That Interacts with Wnt Receptors and Contains a Novel Repetitive Short Linear Motif That Activates Wnt Signaling

Ehrlichia chaffeensis expresses the TRP120 multifunctional effector, which is known to play a role in phagocytic entry, on the surface of infectious dense-cored ehrlichiae, but a cognate host receptor has not been identified. We recently reported that E. chaffeensis activates canonical Wnt signaling in monocytes to promote bacterial uptake and intracellular survival and that TRP120 was involved in this activation event. To identify the specific mechanism of pathway activation, we hypothesized that TRP120 is a Wnt signaling ligand mimetic that initiates Wnt pathway activity through direct interaction with the Wnt pathway Frizzled family of receptors. In this study, we used confocal immunofluorescence microscopy to demonstrate very strong colocalization between E. chaffeensis and Fzd2, 4, 5, 7, and 9 as well as coreceptor LRP5 at 1 to 3 h postinfection. Direct binding between TRP120 and multiple Fzd receptors was further confirmed by enzyme-linked immunosorbent assay (ELISA) and surface plasmon resonance (SPR). Interfering RNA knockdown of Wnt receptors, coreceptors, and signaling pathway components significantly reduced E. chaffeensis infection, demonstrating that complex and redundant interactions are involved in Wnt pathway exploitation. We utilized in silico approaches to identify a repetitive short linear motif (SLiM) in TRP120 that is homologous to Wnt ligands and used mutant SLiM peptides and an α-TRP120-Wnt-SLiM antibody to demonstrate that the TRP120 Wnt SLiM activates the canonical Wnt pathway and promotes E. chaffeensis infection. This study reports the first example of bacterial mimicry of Wnt pathway ligands and highlights a pathogenic mechanism with potential for targeting by antimicrobial therapeutics.IMPORTANCE Upon infecting mammalian hosts, Ehrlichia chaffeensis establishes a replicative niche in microbe-eating immune system cells where it expertly orchestrates infection and spread. One of the ways Ehrlichia survives within these phagocytes is by activating evolutionarily conserved signaling pathways including the Wnt pathway; however, the molecular details of pathway hijacking have not been defined. This study is significant because it identifies an ehrlichial protein that directly interacts with components of the Wnt receptor complex, influencing pathway activity and promoting infection. Consequentially, Ehrlichia serves as a unique tool to investigate the intricacies of how pathogens repurpose human immune cell signaling and provides an opportunity to better understand many cellular processes in health and disease. Furthermore, understanding how this bacterium utilizes its small genome to survive within cells that evolved to destroy pathogens will facilitate the development of antibacterial therapeutics that could target Ehrlichia as well as other intracellular agents of human disease.

Short Linear Motifs Characterizing Snake Venom and Mammalian Phospholipases A2

Snake venom phospholipases A2 (PLA2s) have sequences and structures very similar to those of mammalian group I and II secretory PLA2s, but they possess many toxic properties, ranging from the inhibition of coagulation to the blockage of nerve transmission, and the induction of muscle necrosis. The biological properties of these proteins are not only due to their enzymatic activity, but also to protein–protein interactions which are still unidentified. Here, we compare sequence alignments of snake venom and mammalian PLA2s, grouped according to their structure and biological activity, looking for differences that can justify their different behavior. This bioinformatics analysis has evidenced three distinct regions, two central and one C-terminal, having amino acid compositions that distinguish the different categories of PLA2s. In these regions, we identified short linear motifs (SLiMs), peptide modules involved in protein–protein interactions, conserved in mammalian and not in snake venom PLA2s, or vice versa. The different content in the SLiMs of snake venom with respect to mammalian PLA2s may result in the formation of protein membrane complexes having a toxic activity, or in the formation of complexes whose activity cannot be blocked due to the lack of switches in the toxic PLA2s, as the motif recognized by the prolyl isomerase Pin1.

UbiNet 2.0: a verified, classified, annotated and updated database of E3 ubiquitin ligase-substrate interactions

Ubiquitination is an important post-translational modification, which controls protein turnover by labeling malfunctional and redundant proteins for proteasomal degradation, and also serves intriguing non-proteolytic regulatory functions. E3 ubiquitin ligases, whose substrate specificity determines the recognition of target proteins of ubiquitination, play crucial roles in ubiquitin-proteasome system. UbiNet 2.0 is an updated version of the database UbiNet. It contains 3332 experimentally verified E3-substrate interactions (ESIs) in 54 organisms and rich annotations useful for investigating the regulation of ubiquitination and the substrate specificity of E3 ligases. Based on the accumulated ESIs data, the recognition motifs in substrates for each E3 were also identified and a functional enrichment analysis was conducted on the collected substrates. To facilitate the research on ESIs with different categories of E3 ligases, UbiNet 2.0 performed strictly evidence-based classification of the E3 ligases in the database based on their mechanisms of ubiquitin transfer and substrate specificity. The platform also provides users with an interactive tool that can visualize the ubiquitination network of a group of self-defined proteins, displaying ESIs and protein-protein interactions in a graphical manner. The tool can facilitate the exploration of inner regulatory relationships mediated by ubiquitination among proteins of interest. In summary, UbiNet 2.0 is a user-friendly web-based platform that provides comprehensive as well as updated information about experimentally validated ESIs and a visualized tool for the construction of ubiquitination regulatory networks available at http://awi.cuhk.edu.cn/~ubinet/index.php.

TNF receptor-associated factor 3 restrains B-cell receptor signaling in normal and malignant B cells

TRAF3 has diverse signaling functions, which vary by cell type. Uniquely in B lymphocytes, TRAF3 inhibits homeostatic survival. Highlighting the role of TRAF3 as a tumor suppressor, loss-of-function TRAF3 mutations are associated with human B-cell malignancies, while B-cell-specific deletion of TRAF3 in mice leads to autoimmunity and lymphoma development. The role of TRAF3 in inhibiting noncanonical NF-κB activation, CD40 and BAFF-R signaling to B cells is well documented. In contrast, TRAF3 enhances many T-cell effector functions, through associating with and enhancing signaling by the T-cell receptor (TCR)-CD28 complex. The present study was designed to determine the role of TRAF3 in signaling via the B-cell antigen receptor (BCR). The BCR is crucial for antigen recognition, survival, proliferation, and antibody production, and defects in BCR signaling can promote abnormal survival of malignant B cells. Here, we show that TRAF3 is associated with both CD79B and the BCR-activated kinases Syk and Btk following BCR stimulation. BCR-induced phosphorylation of Syk and additional downstream kinases was increased in TRAF3^−/− B cells, with regulation observed in both follicular and marginal zone B-cell subsets. BCR stimulation of TRAF3^−/− B cells resulted in increased surface expression of MHC-II, CD80, and CD86 molecules. Interestingly, increased survival of TRAF3^−/− primary B cells was resistant to inhibition of Btk, while TRAF3-deficient malignant B-cell lines showed enhanced sensitivity. TRAF3 serves to restrain normal and malignant BCR signaling, with important implications for its role in normal B-cell biology and abnormal survival of malignant B cells.

Interesting Biochemistries in the Structure and Function of Bacterial Effectors

Bacterial effector proteins, delivered into host cells by specialized multiprotein secretion systems, are a key mediator of bacterial pathogenesis. Following delivery, they modulate a range of host cellular processes and functions. Strong selective pressures have resulted in bacterial effectors evolving unique structures that can mimic host protein biochemical activity or enable novel and distinct biochemistries. Despite the protein structure-function paradigm, effectors from different bacterial species that share biochemical activities, such as the conjugation of ubiquitin to a substrate, do not necessarily share structural or sequence homology to each other or the eukaryotic proteins that carry out the same function. Furthermore, some bacterial effectors have evolved structural variations to known protein folds which enable different or additional biochemical and physiological functions. Despite the overall low occurrence of intrinsically disordered proteins or regions in prokaryotic proteomes compared to eukaryotes proteomes, bacterial effectors appear to have adopted intrinsically disordered regions that mimic the disordered regions of eukaryotic signaling proteins. In this review, we explore examples of the diverse biochemical properties found in bacterial effectors that enable effector-mediated interference of eukaryotic signaling pathways and ultimately support pathogenesis. Despite challenges in the structural and functional characterisation of effectors, recent progress has been made in understanding the often unusual and fascinating ways in which these virulence factors promote pathogenesis. Nevertheless, continued work is essential to reveal the array of remarkable activities displayed by effectors.

On the Prevalence and Potential Functionality of an Intrinsic Disorder in the MERS-CoV Proteome

Middle East respiratory syndrome is a severe respiratory illness caused by an infectious coronavirus. This virus is associated with a high mortality rate, but there is as of yet no effective vaccine or antibody available for human immunity/treatment. Drug design relies on understanding the 3D structures of viral proteins; however, arriving at such understanding is difficult for intrinsically disordered proteins, whose disorder-dependent functions are key to the virus’s biology. Disorder is suggested to provide viral proteins with highly flexible structures and diverse functions that are utilized when invading host organisms and adjusting to new habitats. To date, the functional roles of intrinsically disordered proteins in the mechanisms of MERS-CoV pathogenesis, transmission, and treatment remain unclear. In this study, we performed structural analysis to evaluate the abundance of intrinsic disorder in the MERS-CoV proteome and in individual proteins derived from the MERS-CoV genome. Moreover, we detected disordered protein binding regions, namely, molecular recognition features and short linear motifs. Studying disordered proteins/regions in MERS-CoV could contribute to unlocking the complex riddles of viral infection, exploitation strategies, and drug development approaches in the near future by making it possible to target these important (yet challenging) unstructured regions.