Molecular Recognition Features in Zika Virus Proteome

Viruses have compact genomes that encode limited number of proteins in comparison to other biological entities. Interestingly, viral proteins have shown natural abundance of either completely disordered proteins that are recognized as intrinsically disorder proteins (IDPs) or partially disordered segments known as intrinsically disordered protein regions (IDPRs). IDPRs are involved in interactions with multiple binding partners to accomplish signaling, regulation, and control functions in cells. Tuning of IDPs and IDPRs are mediated through post-translational modification and alternative splicing. Often, the interactions of IDPRs with their binding protein partner(s) lead to transition from the state of disorder to ordered form. Such interaction-prone protein IDPRs are identified as molecular recognition features (MoRFs). Molecular recognition is an important initial step for the biomolecular interactions and their functional proceedings. Although previous studies have established occurrence of the IDPRs in Zika virus proteome, which provide the functional diversity and structural plasticity to viral proteins, the MoRF analysis has not been performed as of yet. Many computational methods have been developed for the identification of the MoRFs in protein sequences including ANCHOR, MoRFpred, DISOPRED3, and MoRFchibi_web server. In the current study, we have investigated the presence of MoRF regions in structural and non-structural proteins of Zika virus using an aforementioned set of computational techniques. Furthermore, we have experimentally validated the intrinsic disorderness of NS2B cofactor region of NS2B–NS3 protease. NS2B has one of the longest MoRF regions in Zika virus proteome. In future, this study may provide valuable information while investigating the virus host protein interaction networks.     

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.     

Exploring the binding diversity of intrinsically disordered proteins involved in one‐to‐many binding

Molecular recognition features (MoRFs) are intrinsically disordered protein regions that bind to partners via disorder‐to‐order transitions. In one‐to‐many binding, a single MoRF binds to two or more different partners individually. MoRF‐based one‐to‐many protein–protein interaction (PPI) examples were collected from the Protein Data Bank, yielding 23 MoRFs bound to 2–9 partners, with all pairs of same‐MoRF partners having less than 25% sequence identity. Of these, 8 MoRFs were bound to 2–9 partners having completely different folds, whereas 15 MoRFs were bound to 2–5 partners having the same folds but with low sequence identities. For both types of partner variation, backbone and side chain torsion angle rotations were used to bring about the conformational changes needed to enable close fits between a single MoRF and distinct partners. Alternative splicing events (ASEs) and posttranslational modifications (PTMs) were also found to contribute to distinct partner binding. Because ASEs and PTMs both commonly occur in disordered regions, and because both ASEs and PTMs are often tissue‐specific, these data suggest that MoRFs, ASEs, and PTMs may collaborate to alter PPI networks in different cell types. These data enlarge the set of carefully studied MoRFs that use inherent flexibility and that also use ASE‐based and/or PTM‐based surface modifications to enable the same disordered segment to selectively associate with two or more partners. The small number of residues involved in MoRFs and in their modifications by ASEs or PTMs may simplify the evolvability of signaling network diversity.     

Analysis of molecular recognition features (MoRFs)

Several proteomic studies in the last decade revealed that many proteins are either completely disordered or possess long structurally flexible regions. Many such regions were shown to be of functional importance, often allowing a protein to interact with a large number of diverse partners. Parallel to these findings, during the last five years structural bioinformatics has produced an explosion of results regarding protein-protein interactions and their importance for cell signaling. We studied the occurrence of relatively short (10-70 residues), loosely structured protein regions within longer, largely disordered sequences that were characterized as bound to larger proteins. We call these regions molecular recognition features (MoRFs, also known as molecular recognition elements, MoREs). Interestingly, upon binding to their partner(s), MoRFs undergo disorder-to-order transitions. Thus, in our interpretation, MoRFs represent a class of disordered region that exhibits molecular recognition and binding functions. This work extends previous research showing the importance of flexibility and disorder for molecular recognition. We describe the development of a database of MoRFs derived from the RCSB Protein Data Bank and present preliminary results of bioinformatics analyses of these sequences. Based on the structure adopted upon binding, at least three basic types of MoRFs are found: α-MoRFs, β-MoRFs, and ι-MoRFs, which form α-helices, β-strands, and irregular secondary structure when bound, respectively. Our data suggest that functionally significant residual structure can exist in MoRF regions prior to the actual binding event. The contribution of intrinsic protein disorder to the nature and function of MoRFs has also been addressed. The results of this study will advance the understanding of protein-protein interactions and help towards the future development of useful protein-protein binding site predictors.     

Abundance and functional roles of intrinsic disorder in allergenic proteins and allergen representative peptides

The pathological process of allergies generally involves an initial activation of certain immune cells, tied to an ensuing inflammatory reaction on renewed contact with the allergen. In IgE-mediated hypersensitivity, this typically occurs in response to otherwise harmless food- or air-borne proteins. As some members of certain protein families carry special properties that make them allergenic, exploring protein allergens at the molecular level is instrumental to an improved understanding of the disease mechanisms, including the identification of relevant antigen features. For this purpose, we inspected a previously identified set of allergen representative peptides (ARPs) to scrutinize protein intrinsic disorder. The resulting study presented here focused on the association between these ARPs and protein intrinsic disorder. In addition, the connection between the disorder-enriched ARPs and UniProt functional keywords was considered. Our analysis revealed that ～ 20% of the allergen peptides are highly disordered, and that ～ 77% of ARPs are either located within disordered regions of corresponding allergenic proteins or show more disorder/flexibility than their neighbor regions. Furthermore, among the subset of allergenic proteins, ～ 70% of the predicted molecular recognition features (MoRFs that consist of short interactive disordered regions undergoing disorder-to-order transitions at interaction with binding partners) were identified as ARPs. These results suggest that intrinsic disorder and MoRFs may play functional roles in IgE-mediated allergy.     

Characterization of molecular recognition features, MoRFs, and their binding partners

Molecular Recognition Features (MoRFs) are short, interaction-prone segments of protein disorder that undergo disorder-to-order transitions upon specific binding, representing a specific class of intrinsically disordered regions that exhibit molecular recognition and binding functions. MoRFs are common in various proteomes and occupy a unique structural and functional niche in which function is a direct consequence of intrinsic disorder. Example MoRFs collected from the Protein Data Bank (PDB) have been divided into three subtypes according to their structures in the bound state: alpha-MoRFs form alpha-helices, beta-MoRFs form beta-strands, and iota-MoRFs form structures without a regular pattern of backbone hydrogen bonds. These example MoRFs were indicated to be intrinsically disordered in the absence of their binding partners by several criteria. In this study, we used several geometric and physiochemical criteria to examine the properties of 62 alpha-, 20 beta-, and 176 iota-MoRF complex structures. Interface residues were examined by calculating differences in accessible surface area between the complex and isolated monomers. The compositions and physiochemical properties of MoRF and MoRF partner interface residues were compared to the interface residues of homodimers, heterodimers, and antigen-antibody complexes. Our analysis indicates that there are significant differences in residue composition and several geometric and physicochemical properties that can be used to discriminate, with a high degree of accuracy, between various interfaces in protein interaction data sets. Implications of these findings for the development of MoRF-partner interaction predictors are discussed. In addition, structural changes upon MoRF-to-partner complex formation were examined for several illustrative examples.     

Defining the structural basis for assembly of a transmembrane cytochrome

To define the structural basis for cofactor binding to membrane proteins, we introduce a manageable model system, which allows us, for the first time, to study the influence of individual transmembrane helices and of single amino acid residues on the assembly of a transmembrane cytochrome. In vivo as well as in vitro analyses indicate central roles of single amino acid residues for either interaction of the transmembrane helices or for binding of the cofactor. The results clearly show that interaction of the PsbF transmembrane helix is independent from binding of the heme cofactor. On the other hand, binding of the cofactor highly depends on helix-helix interactions. By site-directed mutagenesis critical amino acid residues were identified, which are involved in the assembly of a functional transmembrane cytochrome. Especially, a highly conserved glycine residue is critical for interaction of the transmembrane helices and assembly of the cytochrome. Based on the two-stage-model of alpha-helical membrane protein folding, the presented results clearly indicate a third stage of membrane protein folding, in which a cofactor binds to a pre-assembled transmembrane protein.     

Functional competition within a membrane: Lipid recognition vs. transmembrane helix oligomerization

Binding of specific lipids to large, polytopic membrane proteins is well described, and it is clear that such lipids are crucial for protein stability and activity. In contrast, binding of defined lipid species to individual transmembrane helices and regulation of transmembrane helix monomer–oligomer equilibria by binding of distinct lipids is a concept, which has emerged only lately. Lipids bind to single-span membrane proteins, both in the juxta-membrane region as well as in the hydrophobic membrane core. While some interactions counteract transmembrane helix oligomerization, in other cases lipid binding appears to enhance oligomerization. As reversible oligomerization is involved in activation of many membrane proteins, binding of defined lipids to single-span transmembrane proteins might be a mechanism to regulate and/or fine-tune the protein activity. But how could lipid binding trigger the activity of a protein? How can binding of a single lipid molecule to a transmembrane helix affect the structure of a transmembrane helix oligomer, and consequently its signaling state? These questions are discussed in the present article based on recent results obtained with simple, single-span transmembrane proteins. This article is part of a Special Issue entitled: Lipid–protein interactions.     

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.     

Identification of Novel Death-Associated Protein Kinase 2 Interaction Partners by Proteomic Screening Coupled with Bimolecular Fluorescence Complementation

Death-associated protein kinase 2 (DAPK2) is a Ca²⁺/calmodulin-dependent Ser/Thr kinase that possesses tumor-suppressive functions and regulates programmed cell death, autophagy, oxidative stress, hematopoiesis, and motility. As only few binding partners of DAPK2 have been determined, the molecular mechanisms governing these biological functions are largely unknown. We report the identification of 180 potential DAPK2 interaction partners by affinity purification-coupled mass spectrometry, 12 of which are known DAPK binding proteins. A small subset of established and potential binding proteins detected in this screen was further investigated by bimolecular fluorescence complementation (BiFC) assays, a method to visualize protein interactions in living cells. These experiments revealed that α-actinin-1 and 14-3-3-β are novel DAPK2 binding partners. The interaction of DAPK2 with α-actinin-1 was localized at the plasma membrane, resulting in massive membrane blebbing and reduced cellular motility, whereas the interaction of DAPK2 with 14-3-3-β was localized to the cytoplasm, with no impact on blebbing, motility, or viability. Our results therefore suggest that DAPK2 effector functions are influenced by the protein''s subcellular localization and highlight the utility of combining mass spectrometry screening with bimolecular fluorescence complementation to identify and characterize novel protein-protein interactions.     

In vitro analysis of the complete CIDE domain interactions of the Drep system in fly

CIDE domain containing proteins are involved in apoptosis and lipid metabolism, and four CIDE containing proteins, Drep1, Drep2, Drep3, and Drep4, have been identified in fly. In this study, we found that Drep3 interacts with Drep4 via the CIDE domain specifically, which completes the interaction map of Drep system in fly, cyclic interactions: Drep1–Drep2–Drep3–Drep4–Drep1. In addition, we analyzed the dynamic stoichiometry changes of Drep proteins upon binding to their binding partners. Our current studies will help us to understand Drep system in fly as well as CIDE domain for protein–protein interactions.     

Evolutionary Capacitance and Control of Protein Stability in Protein-Protein Interaction Networks

In addition to their biological function, protein complexes reduce the exposure of the constituent proteins to the risk of undesired oligomerization by reducing the concentration of the free monomeric state. We interpret this reduced risk as a stabilization of the functional state of the protein. We estimate that protein-protein interactions can account for of additional stabilization; a substantial contribution to intrinsic stability. We hypothesize that proteins in the interaction network act as evolutionary capacitors which allows their binding partners to explore regions of the sequence space which correspond to less stable proteins. In the interaction network of baker''s yeast, we find that statistically proteins that receive higher energetic benefits from the interaction network are more likely to misfold. A simplified fitness landscape wherein the fitness of an organism is inversely proportional to the total concentration of unfolded proteins provides an evolutionary justification for the proposed trends. We conclude by outlining clear biophysical experiments to test our predictions.     

Membrane interactions of intrinsically disordered proteins: The example of alpha-synuclein

Peripheral membrane proteins associate reversibly with biological membranes that, compared to protein binding partners, are structurally labile and devoid of specific binding pockets. Membranes in different subcellular compartments vary primarily in their chemical composition and physical properties, and recognition of these features is therefore critical for allowing such proteins to engage their proper membrane targets. Intrinsically disordered proteins (IDPs) are well-suited to accomplish this task using highly specific and low- to moderate-affinity interactions governed by recognition principles that are both similar to and different from those that mediate the membrane interactions of rigid proteins. IDPs have also evolved multiple mechanisms to regulate membrane (and other) interactions and achieve their impressive functional diversity. Moreover, IDP-membrane interactions may have a kinetic advantage in fast processes requiring rapid control of such interactions, such as synaptic transmission or signaling. Herein we review the biophysics, regulation and functional implications of IDP-membrane interactions and include a brief overview of some of the methods that can be used to study such interactions. At each step, we use the example of alpha-synuclein, a protein involved in the pathogenesis of Parkinson’s disease and one of the best characterized membrane-binding IDP, to illustrate some of the principles discussed.     

Great interactions: How binding incorrect partners can teach us about protein recognition and function

Protein–protein interactions play a key part in most biological processes and understanding their mechanism is a fundamental problem leading to numerous practical applications. The prediction of protein binding sites in particular is of paramount importance since proteins now represent a major class of therapeutic targets. Amongst others methods, docking simulations between two proteins known to interact can be a useful tool for the prediction of likely binding patches on a protein surface. From the analysis of the protein interfaces generated by a massive cross‐docking experiment using the 168 proteins of the Docking Benchmark 2.0, where all possible protein pairs, and not only experimental ones, have been docked together, we show that it is also possible to predict a protein's binding residues without having any prior knowledge regarding its potential interaction partners. Evaluating the performance of cross‐docking predictions using the area under the specificity‐sensitivity ROC curve (AUC) leads to an AUC value of 0.77 for the complete benchmark (compared to the 0.5 AUC value obtained for random predictions). Furthermore, a new clustering analysis performed on the binding patches that are scattered on the protein surface show that their distribution and growth will depend on the protein's functional group. Finally, in several cases, the binding‐site predictions resulting from the cross‐docking simulations will lead to the identification of an alternate interface, which corresponds to the interaction with a biomolecular partner that is not included in the original benchmark. Proteins 2016; 84:1408–1421. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Biochemical and computational analysis of LNX1 interacting proteins

PDZ (Post-synaptic density, 95 kDa, Discs large, Zona Occludens-1) domains are protein interaction domains that bind to the carboxy-terminal amino acids of binding partners, heterodimerize with other PDZ domains, and also bind phosphoinositides. PDZ domain containing proteins are frequently involved in the assembly of multi-protein complexes and clustering of transmembrane proteins. LNX1 (Ligand of Numb, protein X 1) is a RING (Really Interesting New Gene) domain-containing E3 ubiquitin ligase that also includes four PDZ domains suggesting it functions as a scaffold for a multi-protein complex. Here we use a human protein array to identify direct LNX1 PDZ domain binding partners. Screening of 8,000 human proteins with isolated PDZ domains identified 53 potential LNX1 binding partners. We combined this set with LNX1 interacting proteins identified by other methods to assemble a list of 220 LNX1 interacting proteins. Bioinformatic analysis of this protein list was used to select interactions of interest for future studies. Using this approach we identify and confirm six novel LNX1 binding partners: KCNA4, PAK6, PLEKHG5, PKC-alpha1, TYK2 and PBK, and suggest that LNX1 functions as a signalling scaffold.     

Intrinsically disordered domains: Sequence ➔ disorder ➔ function relationships

Disordered domains are long regions of intrinsic disorder that ideally have conserved sequences, conserved disorder, and conserved functions. These domains were first noticed in protein–protein interactions that are distinct from the interactions between two structured domains and the interactions between structured domains and linear motifs or molecular recognition features (MoRFs). So far, disordered domains have not been systematically characterized. Here, we present a bioinformatics investigation of the sequence–disorder–function relationships for a set of probable disordered domains (PDDs) identified from the Pfam database. All the Pfam seed proteins from those domains with at least one PDD sequence were collected. Most often, if a set contains one PDD sequence, then all members of the set are PDDs or nearly so. However, many seed sets have sequence collections that exhibit diverse proportions of predicted disorder and structure, thus giving the completely unexpected result that conserved sequences can vary substantially in predicted disorder and structure. In addition to the induction of structure by binding to protein partners, disordered domains are also induced to form structure by disulfide bond formation, by ion binding, and by complex formation with RNA or DNA. The two new findings, (a) that conserved sequences can vary substantially in their predicted disorder content and (b) that homologues from a single domain can evolve from structure to disorder (or vice versa), enrich our understanding of the sequence ? disorder ensemble ? function paradigm.