The Balancing Act of Intrinsically Disordered Proteins: Enabling Functional Diversity while Minimizing Promiscuity

Intrinsically disordered proteins (IDPs) or regions (IDRs) perform diverse cellular functions, but are also prone to forming promiscuous and potentially deleterious interactions. We investigate the extent to which the properties of, and content in, IDRs have adapted to enable functional diversity while limiting interference from promiscuous interactions in the crowded cellular environment. Information on protein sequences, their predicted intrinsic disorder, and 3D structure contents is related to data on protein cellular concentrations, gene co-expression, and protein–protein interactions in the well-studied yeast Saccharomyces cerevisiae. Results reveal that both the protein IDR content and the frequency of “sticky” amino acids in IDRs (those more frequently involved in protein interfaces) decrease with increasing protein cellular concentration. This implies that the IDR content and the amino acid composition of IDRs experience negative selection as the protein concentration increases. In the S. cerevisiae protein–protein interaction network, the higher a protein's IDR content, the more frequently it interacts with IDR-containing partners, and the more functionally diverse the partners are. Employing a clustering analysis of Gene Ontology terms, we newly identify ~ 600 putative multifunctional proteins in S. cerevisiae. Strikingly, these proteins are enriched in IDRs and contribute significantly to all the observed trends. In particular, IDRs of multi-functional proteins feature more sticky amino acids than IDRs of their non-multifunctional counterparts, or the surfaces of structured yeast proteins. This property likely affords sufficient binding affinity for the functional interactions, commonly mediated by short IDR segments, thereby counterbalancing the loss in overall IDR conformational entropy upon binding.     

Conformational flexibility of BECN1: Essential to its key role in autophagy and beyond

BECN1 (Beclin 1), a highly conserved eukaryotic protein, is a key regulator of autophagy, a cellular homeostasis pathway, and also participates in vacuolar protein sorting, endocytic trafficking, and apoptosis. BECN1 is important for embryonic development, the innate immune response, tumor suppression, and protection against neurodegenerative disorders, diabetes, and heart disease. BECN1 mediates autophagy as a core component of the class III phosphatidylinositol 3‐kinase complexes. However, the exact mechanism by which it regulates the activity of these complexes, or mediates its other diverse functions is unclear. BECN1 interacts with several diverse protein partners, perhaps serving as a scaffold or interaction hub for autophagy. Based on extensive structural, biophysical and bioinformatics analyses, BECN1 consists of an intrinsically disordered region (IDR), which includes a BH3 homology domain (BH3D); a flexible helical domain (FHD); a coiled‐coil domain (CCD); and a β‐α‐repeated autophagy‐specific domain (BARAD). Each of these BECN1 domains mediates multiple diverse interactions that involve concomitant conformational changes. Thus, BECN1 conformational flexibility likely plays a key role in facilitating diverse protein interactions. Further, BECN1 conformation and interactions are also modulated by numerous post‐translational modifications. A better structure‐based understanding of the interplay between different BECN1 conformational and binding states, and the impact of post‐translational modifications will be essential to elucidating the mechanism of its multiple biological roles.     

A Conserved Deubiquitinating Enzyme Uses Intrinsically Disordered Regions to Scaffold Multiple Protein Interaction Sites

In the canonical view of protein function, it is generally accepted that the three-dimensional structure of a protein determines its function. However, the past decade has seen a dramatic growth in the identification of proteins with extensive intrinsically disordered regions (IDRs), which are conformationally plastic and do not appear to adopt single three-dimensional structures. One current paradigm for IDR function is that disorder enables IDRs to adopt multiple conformations, expanding the ability of a protein to interact with a wide variety of disparate proteins. The capacity for many interactions is an important feature of proteins that occupy the hubs of protein networks, in particular protein-modifying enzymes that usually have a broad spectrum of substrates. One such protein modification is ubiquitination, where ubiquitin is attached to proteins through ubiquitin ligases (E3s) and removed through deubiquitinating enzymes. Numerous proteomic studies have found that thousands of proteins are dynamically regulated by cycles of ubiquitination and deubiquitination. Thus, how these enzymes target their wide array of substrates is of considerable importance for understanding the function of the cell''s diverse ubiquitination networks. Here, we characterize a yeast deubiquitinating enzyme, Ubp10, that possesses IDRs flanking its catalytic protease domain. We show that Ubp10 possesses multiple, distinct binding modules within its IDRs that are necessary and sufficient for directing protein interactions important for Ubp10''s known roles in gene silencing and ribosome biogenesis. The human homolog of Ubp10, USP36, also has IDRs flanking its catalytic domain, and these IDRs similarly contain binding modules important for protein interactions. This work highlights the significant protein interaction scaffolding abilities of IDRs in the regulation of dynamic protein ubiquitination.     

BAX and BAK1 are dispensable for ABT-737-induced dissociation of the BCL2-BECN1 complex and autophagy

Disruption of the complex of BECN1 with BCL2 or BCL2L1/BCL-X_L is an essential switch that turns on cellular autophagy in response to environmental stress or treatment with BH3 peptidomimetics. Recently, it has been proposed that BCL2 and BCL2L1/BCL-X_L may inhibit autophagy indirectly through a mechanism dependent on the proapoptotic BCL2 family members, BAX and BAK1. Here we report that the BH3 mimetic, ABT-737, induces autophagy in parallel with disruption of BCL2-BECN1 binding in 2 different apoptosis-deficient cell types lacking BAX and BAK1, namely in mouse embryonic fibroblasts cells and in human colon cancer HCT116 cells. We conclude that the BH3 mimetic ABT-737 induces autophagy through a BAX and BAK1-independent mechanism that likely involves disruption of BECN1 binding to antiapoptotic BCL2 family members.     

Features of molecular recognition of intrinsically disordered proteins via coupled folding and binding

The sequence–structure–function paradigm of proteins has been revolutionized by the discovery of intrinsically disordered proteins (IDPs) or intrinsically disordered regions (IDRs). In contrast to traditional ordered proteins, IDPs/IDRs are unstructured under physiological conditions. The absence of well‐defined three‐dimensional structures in the free state of IDPs/IDRs is fundamental to their function. Folding upon binding is an important mode of molecular recognition for IDPs/IDRs. While great efforts have been devoted to investigating the complex structures and binding kinetics and affinities, our knowledge on the binding mechanisms of IDPs/IDRs remains very limited. Here, we review recent advances on the binding mechanisms of IDPs/IDRs. The structures and kinetic parameters of IDPs/IDRs can vary greatly, and the binding mechanisms can be highly dependent on the structural properties of IDPs/IDRs. IDPs/IDRs can employ various combinations of conformational selection and induced fit in a binding process, which can be templated by the target and/or encoded by the IDP/IDR. Further studies should provide deeper insights into the molecular recognition of IDPs/IDRs and enable the rational design of IDP/IDR binding mechanisms in the future.     

Negatively Charged Disordered Regions are Prevalent and Functionally Important Across Proteomes

Intrinsically disordered regions (IDRs) of proteins are often characterized by a high fraction of charged residues, but differ in their overall net charge and in the organization of the charged residues. The function-encoding information stored via IDR charge composition and organization remains elusive. Here, we aim to decipher the sequence–function relationship in IDRs by presenting a comprehensive bioinformatic analysis of the charge properties of IDRs in the human, mouse, and yeast proteomes. About 50% of the proteins comprise at least a single IDR, which is either positively or negatively charged. Highly negatively charged IDRs are longer and possess greater net charge per residue compared with highly positively charged IDRs. A striking difference between positively and negatively charged IDRs is the characteristics of the repeated units, specifically, of consecutive Lys or Arg residues (K/R repeats) and Asp or Glu (D/E repeats) residues. D/E repeats are found to be about five times longer than K/R repeats, with the longest found containing 49 residues. Long stretches of consecutive D and E are found to be more prevalent in nucleic acid-related proteins. They are less common in prokaryotes, and in eukaryotes their abundance increases with genome size. The functional role of D/E repeats and the profound differences between them and K/R repeats are discussed.     

Intrinsic disorder mediates hepatitis C virus core–host cell protein interactions

Viral proteins bind to numerous cellular and viral proteins throughout the infection cycle. However, the mechanisms by which viral proteins interact with such large numbers of factors remain unknown. Cellular proteins that interact with multiple, distinct partners often do so through short sequences known as molecular recognition features (MoRFs) embedded within intrinsically disordered regions (IDRs). In this study, we report the first evidence that MoRFs in viral proteins play a similar role in targeting the host cell. Using a combination of evolutionary modeling, protein–protein interaction analyses and forward genetic screening, we systematically investigated two computationally predicted MoRFs within the N‐terminal IDR of the hepatitis C virus (HCV) Core protein. Sequence analysis of the MoRFs showed their conservation across all HCV genotypes and the canine and equine Hepaciviruses. Phylogenetic modeling indicated that the Core MoRFs are under stronger purifying selection than the surrounding sequence, suggesting that these modules have a biological function. Using the yeast two‐hybrid assay, we identified three cellular binding partners for each HCV Core MoRF, including two previously characterized cellular targets of HCV Core (DDX3X and NPM1). Random and site‐directed mutagenesis demonstrated that the predicted MoRF regions were required for binding to the cellular proteins, but that different residues within each MoRF were critical for binding to different partners. This study demonstrated that viruses may use intrinsic disorder to target multiple cellular proteins with the same amino acid sequence and provides a framework for characterizing the binding partners of other disordered regions in viral and cellular proteomes.     

Antagonism of Beclin 1‐dependent autophagy by BCL‐2 at the endoplasmic reticulum requires NAF‐1

In addition to mitochondria, BCL‐2 is located at the endoplasmic reticulum (ER) where it is a constituent of several distinct complexes. Here, we identify the BCL‐2‐interacting protein at the ER, nutrient‐deprivation autophagy factor‐1 (NAF‐1)—a bitopic integral membrane protein whose defective expression underlies the aetiology of the neurodegenerative disorder Wolfram syndrome 2 (WFS2). NAF‐1 contains a two iron–two sulphur coordinating domain within its cytosolic region, which is necessary, but not sufficient for interaction with BCL‐2. NAF‐1 is displaced from BCL‐2 by the ER‐restricted BH3‐only protein BIK and contributes to regulation of BIK‐initiated autophagy, but not BIK‐dependent activation of caspases. Similar to BCL‐2, NAF‐1 is found in association with the inositol 1,4,5‐triphosphate receptor and is required for BCL‐2‐mediated depression of ER Ca²⁺ stores. During nutrient deprivation as a physiological stimulus of autophagy, BCL‐2 is known to function through inhibition of the autophagy effector and tumour suppressor Beclin 1. NAF‐1 is required in this pathway for BCL‐2 at the ER to functionally antagonize Beclin 1‐dependent autophagy. Thus, NAF‐1 is a BCL‐2‐associated co‐factor that targets BCL‐2 for antagonism of the autophagy pathway at the ER.     

Protein kinases phosphorylate long disordered regions in intrinsically disordered proteins

Phosphorylation is a major post‐translational modification that plays a central role in signaling pathways. Protein kinases phosphorylate substrates (phosphoproteins) by adding phosphate at Ser/Thr or Tyr residues (phosphosites). A large amount of data identifying and describing phosphosites in phosphoproteins has been reported but the specificity of phosphorylation is not fully resolved. In this report, data of kinase‐substrate pairs identified by the Kinase‐Interacting Substrate Screening (KISS) method were used to analyze phosphosites in intrinsically disordered regions (IDRs) of intrinsically disordered proteins. We compared phosphorylated and nonphosphorylated IDRs and found that the phosphorylated IDRs were significantly longer than nonphosphorylated IDRs. The phosphorylated IDR is often the longest IDR (71%) in a phosphoprotein when only a single phosphosite exists in the IDR, and when the phosphoprotein has multiple phosphosites in an IDR(s), the phosphosites are primarily localized in a single IDR (78%) and this IDR is usually the longest one (81%). We constructed a stochastic model of phosphorylation to estimate the effect of IDR length. The model that accounted for IDR length produced more realistic results when compared with a model that excluded the IDR length. We propose that the IDR length is a significant determinant for locating kinase phosphorylation sites in phosphoproteins.     

Synthesis and Biological Activity of Scyllatoxin-Based BH3 Domain Mimetics Containing Two Disulfide Linkages

The B cell lymphoma 2 (BCL2) proteins are a family of evolutionarily related proteins that act as positive or negative regulators of the intrinsic apoptosis pathway. Overexpression of anti-apoptotic BCL2 proteins in cells is associated with apoptotic resistance, which can result in cancerous phenotypes and pathogenic cell survival. Consequently, anti-apoptotic BCL2 proteins have attracted considerable interest as therapeutic targets. We recently reported the development of a novel class of synthetic protein based on scyllatoxin (ScTx) designed to mimic the helical BH3 interaction domain of the pro-apoptotic BCL2 protein Bax. These studies showed that the number and position of native disulfide linkages contained within the ScTx-Bax structure significantly influences the ability for these constructs to target anti-apoptotic BCL2 proteins in vitro. The goal of the present study is to investigate the contribution of two disulfide linkages in the folding and biological activity of ScTx-Bax proteins. Here, we report the full chemical synthesis of three ScTx-Bax sequence variants, each presenting two native disulfide linkages at different positions within the folded structure. It was observed that two disulfide linkages were sufficient to fold ScTx-Bax proteins into native-like architectures reminiscent of wild-type ScTx. Furthermore, we show that select (bis)disulfide ScTx-Bax variants can target Bcl-2 (proper) in vitro and that the position of the disulfide bonds significantly influences binding affinity. Despite exhibiting only modest binding to Bcl-2, the successful synthesis of ScTx-Bax proteins containing two disulfide linkages represents a viable route to ScTx-based BH3 domain mimetics that preserve native-like conformations. Finally, structural models of ScTx-Bax proteins in complex with Bcl-2 indicate that these helical mimetics bind in similar configurations as wild-type Bax BH3 domains. Taken together, these results suggest that ScTx-Bax proteins may serve as potent lead compounds that expand the repertoire of “druggable” protein–protein interactions.     

Membraneless organelles: P granules in Caenorhabditis elegans

Membraneless organelles are distinct compartments within a cell that are not enclosed by a traditional lipid membrane and instead form through a process called liquid‐liquid phase separation. Examples of these non‐membrane‐bound organelles include nucleoli, stress granules, P bodies, pericentriolar material and germ granules. Many recent studies have used Caenorhabditis elegans germ granules, known as P granules, to expand our understanding of the formation of these unique cellular compartments. From this work, we know that proteins with intrinsically disordered regions (IDRs) play a critical role in the process of phase separation. IDR phase separation is further tuned through their interactions with RNA and through protein modifications such as phosphorylation and methylation. These findings from C elegans, combined with work done in other model organisms, continue to provide insight into the formation of membraneless organelles and the important role they play in compartmentalizing cellular processes.     

Intrinsically disordered regions of nucleophosmin/B23 regulate its RNA binding activity through their inter- and intra-molecular association

Nucleophosmin (NPM1/B23) is a nucleolar protein implicated in growth-associated functions, in which the RNA binding activity of B23 plays essential roles in ribosome biogenesis. The C-terminal globular domain (CTD) of B23 has been believed to be the RNA binding domain because the splicing variant B23.2 lacking the CTD binds considerably less efficiently to RNA. However, the recognition of target RNAs by B23 remains poorly understood. Herein, we report a novel mechanism by which B23 recognizes specific RNA targets. We observed that the nucleolar retention of B23.3 lacking the basic region of B23.1 was lower than that of B23.1 because of its low RNA binding activity. Circular dichroism measurements indicated that the basic region and adjacent acidic regions of B23 are intrinsically disordered regions (IDRs). Biochemical analyses revealed that the basic IDR alone strongly binds to RNA with low specificity. The excessive RNA binding activity of the basic IDR was restrained by intra-molecular interaction with the acidic IDR of B23. Chemical cross-linking experiments and fluorescent labeling of bipartite tetracysteine-tagged proteins suggested that the inter- and intra-molecular interactions between the two IDRs contribute to the regulation of the RNA binding activity of CTD to control the cellular localization and functions of B23.     

<Emphasis Type="Italic">Escherichia coli</Emphasis> AlkB interacts with single-stranded DNA binding protein SSB by an intrinsically disordered region of SSB

Intrinsically disordered regions (IDRs) of proteins often regulate function through interactions with folded domains. Escherichia coli single-stranded DNA binding protein SSB binds and stabilizes single-stranded DNA (ssDNA). The N-terminal of SSB contains characteristic OB (oligonucleotide/oligosaccharide-binding) fold which binds ssDNA tightly but non-specifically. SSB also forms complexes with a large number proteins via the C-terminal interaction domain consisting mostly of acidic amino acid residues. The amino acid residues located between the OB-fold and C-terminal acidic domain are known to constitute an IDR and no functional significance has been attributed to this region. Although SSB is known to bind many DNA repair protein, it is not known whether it binds to DNA dealkylation repair protein AlkB. Here, we characterize AlkB SSB interaction and demonstrate that SSB binds to AlkB via the IDR. We have established that AlkB-SSB interaction by in vitro pull-down and yeast two-hybrid analysis. We mapped the site of contact to be the residues 152–169 of SSB. Unlike most of the SSB-binding proteins which utilize C-terminal acidic domain for interaction, IDR of SSB is necessary and sufficient for AlkB interaction.     

Characterization of intrinsically disordered regions in proteins informed by human genetic diversity

All proteomes contain both proteins and polypeptide segments that don’t form a defined three-dimensional structure yet are biologically active—called intrinsically disordered proteins and regions (IDPs and IDRs). Most of these IDPs/IDRs lack useful functional annotation limiting our understanding of their importance for organism fitness. Here we characterized IDRs using protein sequence annotations of functional sites and regions available in the UniProt knowledgebase (“UniProt features”: active site, ligand-binding pocket, regions mediating protein-protein interactions, etc.). By measuring the statistical enrichment of twenty-five UniProt features in 981 IDRs of 561 human proteins, we identified eight features that are commonly located in IDRs. We then collected the genetic variant data from the general population and patient-based databases and evaluated the prevalence of population and pathogenic variations in IDPs/IDRs. We observed that some IDRs tolerate 2 to 12-times more single amino acid-substituting missense mutations than synonymous changes in the general population. However, we also found that 37% of all germline pathogenic mutations are located in disordered regions of 96 proteins. Based on the observed-to-expected frequency of mutations, we categorized 34 IDRs in 20 proteins (DDX3X, KIT, RB1, etc.) as intolerant to mutation. Finally, using statistical analysis and a machine learning approach, we demonstrate that mutation-intolerant IDRs carry a distinct signature of functional features. Our study presents a novel approach to assign functional importance to IDRs by leveraging the wealth of available genetic data, which will aid in a deeper understating of the role of IDRs in biological processes and disease mechanisms.     

High-throughput prediction of RNA,DNA and protein binding regions mediated by intrinsic disorder

Intrinsically disordered proteins and regions (IDPs and IDRs) lack stable 3D structure under physiological conditions in-vitro, are common in eukaryotes, and facilitate interactions with RNA, DNA and proteins. Current methods for prediction of IDPs and IDRs do not provide insights into their functions, except for a handful of methods that address predictions of protein-binding regions. We report first-of-its-kind computational method DisoRDPbind for high-throughput prediction of RNA, DNA and protein binding residues located in IDRs from protein sequences. DisoRDPbind is implemented using a runtime-efficient multi-layered design that utilizes information extracted from physiochemical properties of amino acids, sequence complexity, putative secondary structure and disorder and sequence alignment. Empirical tests demonstrate that it provides accurate predictions that are competitive with other predictors of disorder-mediated protein binding regions and complementary to the methods that predict RNA- and DNA-binding residues annotated based on crystal structures. Application in Homo sapiens, Mus musculus, Caenorhabditis elegans and Drosophila melanogaster proteomes reveals that RNA- and DNA-binding proteins predicted by DisoRDPbind complement and overlap with the corresponding known binding proteins collected from several sources. Also, the number of the putative protein-binding regions predicted with DisoRDPbind correlates with the promiscuity of proteins in the corresponding protein–protein interaction networks. Webserver: http://biomine.ece.ualberta.ca/DisoRDPbind/