Structural disorder in eukaryotes

Ferroptosis - An iron- and disorder-dependent programmed cell death

Programmed cell death (PCD) is an integral component of both developmental and pathological features of an organism. Recently, ferroptosis, a new form of PCD that is dependent on reactive oxygen species and iron, has been described. As with apoptosis, necroptosis, and autophagy, ferroptosis is associated with a large set of proteins assembled in protein-protein interaction (PPI) networks, interactability of which is driven by the presence of intrinsically disordered proteins (IDPs) and IDP regions (IDPRs). Previous investigations have identified the prevalence and functionality of IDPs/IDPRs in non-ferroptosis PCD. The intrinsic disorder in protein structures is used to increase the regulatory powers of these processes. As uncontrolled PCD is associated with the onset of various pathological traits, uncovering the association between intrinsic disorder and ferroptosis-related proteins is crucial. To understand this association, 31 human ferroptosis-related proteins were analyzed via a multi-dimensional array of bioinformatics and computational techniques. In addition, a disorder meta-predictor (PONDR® FIT) was implored to look at the evolutionary conservation of intrinsic disorder in these proteins. This study presents evidence that IDPs and IDPRs are prevalent in ferroptosis. The intrinsic disorder found in ferroptosis has far-reaching functional implications related to ferroptosis-related PPIs and molecular interactions.

Intrinsic disorder associated with 14-3-3 proteins and their partners

Protein-protein interactions (PPIs) mediate a variety of cellular processes and form complex networks, where connectivity is achieved owing to the "hub" proteins whose interaction with multiple protein partners is facilitated by the intrinsically disordered protein regions (IDPRs) and posttranslational modifications (PTMs). Universal regulatory proteins of the eukaryotic 14-3-3 family nicely exemplify these concepts and are the focus of this chapter. The extremely wide interactome of 14-3-3 proteins is characterized by high levels of intrinsic disorder (ID) enabling protein phosphorylation and consequent specific binding to the well-structured 14-3-3 dimers, one of the first phosphoserine/phosphothreonine binding modules discovered. However, high ID enrichment also challenges structural studies, thereby limiting the progress in the development of small molecule modulators of the key 14-3-3 PPIs of increased medical importance. Besides the well-known structural flexibility of their variable C-terminal tails, recent studies revealed the strong and conserved ID propensity hidden in the N-terminal segment of 14-3-3 proteins (~40 residues), normally forming the α-helical dimerization region, that may have a potential role for the dimer/monomer dynamics and recently reported moonlighting chaperone-like activity of these proteins. We review the role of ID in the 14-3-3 structure, their interactome, and also in selected 14-3-3 complexes. In addition, we discuss approaches that, in the future, may help minimize the disproportion between the large amount of known 14-3-3 partners and the small number of 14-3-3 complexes characterized with atomic precision, to unleash the whole potential of 14-3-3 PPIs as drug targets.

Misprediction of Structural Disorder in Halophiles

Whereas the concept of intrinsic disorder derives from biophysical observations of the lack of structure of proteins or protein regions under native conditions, many of our respective concepts rest on proteome-scale bioinformatics predictions. It is established that most predictors work reliably on proteins commonly encountered, but it is often neglected that we know very little about their performance on proteins of microorganisms that thrive in environments of extreme temperature, pH, or salt concentration, which may cause adaptive sequence composition bias. To address this issue, we predicted structural disorder for the complete proteomes of different extremophile groups by popular prediction methods and compared them to those of the reference mesophilic group. While significant deviations from mesophiles could be explained by a lack or gain of disordered regions in hyperthermophiles and radiotolerants, respectively, we found systematic overprediction in the case of halophiles. Additionally, examples were collected from the Protein Data Bank (PDB) to demonstrate misprediction and to help understand the underlying biophysical principles, i.e., halophilic proteins maintain a highly acidic and hydrophilic surface to avoid aggregation in high salt conditions. Although sparseness of data on disordered proteins from extremophiles precludes the development of dedicated general predictors, we do formulate recommendations for how to address their disorder with current bioinformatics tools.

Large-Scale Analyses of Site-Specific Evolutionary Rates across Eukaryote Proteomes Reveal Confounding Interactions between Intrinsic Disorder, Secondary Structure, and Functional Domains

Various structural and functional constraints govern the evolution of protein sequences. As a result, the relative rates of amino acid replacement among sites within a protein can vary significantly. Previous large-scale work on Metazoan (Animal) protein sequence alignments indicated that amino acid replacement rates are partially driven by a complex interaction among three factors: intrinsic disorder propensity; secondary structure; and functional domain involvement. Here, we use sequence-based predictors to evaluate the effects of these factors on site-specific sequence evolutionary rates within four eukaryotic lineages: Metazoans; Plants; Saccharomycete Fungi; and Alveolate Protists. Our results show broad, consistent trends across all four Eukaryote groups. In all four lineages, there is a significant increase in amino acid replacement rates when comparing: (i) disordered vs. ordered sites; (ii) random coil sites vs. sites in secondary structures; and (iii) inter-domain linker sites vs. sites in functional domains. Additionally, within Metazoans, Plants, and Saccharomycetes, there is a strong confounding interaction between intrinsic disorder and secondary structure-alignment sites exhibiting both high disorder propensity and involvement in secondary structures have very low average rates of sequence evolution. Analysis of gene ontology (GO) terms revealed that in all four lineages, a high fraction of sequences containing these conserved, disordered-structured sites are involved in nucleic acid binding. We also observe notable differences in the statistical trends of Alveolates, where intrinsically disordered sites are more variable than in other Eukaryotes and the statistical interactions between disorder and other factors are less pronounced.

Co-Evolution of Intrinsically Disordered Proteins with Folded Partners Witnessed by Evolutionary Couplings

Although improved strategies for the detection and analysis of evolutionary couplings (ECs) between protein residues already enable the prediction of protein structures and interactions, they are mostly restricted to conserved and well-folded proteins. Whereas intrinsically disordered proteins (IDPs) are central to cellular interaction networks, due to the lack of strict structural constraints, they undergo faster evolutionary changes than folded domains. This makes the reliable identification and alignment of IDP homologs difficult, which led to IDPs being omitted in most large-scale residue co-variation analyses. By preforming a dedicated analysis of phylogenetically widespread bacterial IDP–partner interactions, here we demonstrate that partner binding imposes constraints on IDP sequences that manifest in detectable interprotein ECs. These ECs were not detected for interactions mediated by short motifs, rather for those with larger IDP–partner interfaces. Most identified coupled residue pairs reside close (<10 Å) to each other on the interface, with a third of them forming multiple direct atomic contacts. EC-carrying interfaces of IDPs are enriched in negatively charged residues, and the EC residues of both IDPs and partners preferentially reside in helices. Our analysis brings hope that IDP–partner interactions difficult to study could soon be successfully dissected through residue co-variation analysis.

The lineage-specific, intrinsically disordered N-terminal extension of monothiol glutaredoxin 1 from trypanosomes contains a regulatory region

Glutaredoxins (Grx) are small proteins conserved throughout all the kingdoms of life that are engaged in a wide variety of biological processes and share a common thioredoxin-fold. Among them, class II Grx are redox-inactive proteins involved in iron-sulfur (FeS) metabolism. They contain a single thiol group in their active site and use low molecular mass thiols such as glutathione as ligand for binding FeS-clusters. In this study, we investigated molecular aspects of 1CGrx1 from the pathogenic parasite Trypanosoma brucei brucei, a mitochondrial class II Grx that fulfills an indispensable role in vivo. Mitochondrial 1CGrx1 from trypanosomes differs from orthologues in several features including the presence of a parasite-specific N-terminal extension (NTE) whose role has yet to be elucidated. Previously we have solved the structure of a truncated form of 1CGrx1 containing only the conserved glutaredoxin domain but lacking the NTE. Our aim here is to investigate the effect of the NTE on the conformation of the protein. We therefore solved the NMR structure of the full-length protein, which reveals subtle but significant differences with the structure of the NTE-less form. By means of different experimental approaches, the NTE proved to be intrinsically disordered and not involved in the non-redox dependent protein dimerization, as previously suggested. Interestingly, the portion comprising residues 65–76 of the NTE modulates the conformational dynamics of the glutathione-binding pocket, which may play a role in iron-sulfur cluster assembly and delivery. Furthermore, we disclosed that the class II-strictly conserved loop that precedes the active site is critical for stabilizing the protein structure. So far, this represents the first communication of a Grx containing an intrinsically disordered region that defines a new protein subgroup within class II Grx.

Human long intrinsically disordered protein regions are frequent targets of positive selection

Intrinsically disordered regions occur frequently in proteins and are characterized by a lack of a well-defined three-dimensional structure. Although these regions do not show a higher order of structural organization, they are known to be functionally important. Disordered regions are rapidly evolving, largely attributed to relaxed purifying selection and an increased role of genetic drift. It has also been suggested that positive selection might contribute to their rapid diversification. However, for our own species, it is currently unknown whether positive selection has played a role during the evolution of these protein regions. Here, we address this question by investigating the evolutionary pattern of more than 6600 human proteins with intrinsically disordered regions and their ordered counterparts. Our comparative approach with data from more than 90 mammalian genomes uses a priori knowledge of disordered protein regions, and we show that this increases the power to detect positive selection by an order of magnitude. We can confirm that human intrinsically disordered regions evolve more rapidly, not only within humans but also across the entire mammalian phylogeny. They have, however, experienced substantial evolutionary constraint, hinting at their fundamental functional importance. We find compelling evidence that disordered protein regions are frequent targets of positive selection and estimate that the relative rate of adaptive substitutions differs fourfold between disordered and ordered protein regions in humans. Our results suggest that disordered protein regions are important targets of genetic innovation and that the contribution of positive selection in these regions is more pronounced than in other protein parts.

Disordered epitopes as peptide vaccines

The development of clinically useful peptide-based vaccines remains a long-standing goal. This review highlights that intrinsically disordered protein antigens, which lack an ordered three-dimensional structure, represent excellent starting points for the development of such vaccines. Disordered proteins represent an important class of antigen in a wide range of human pathogens, and, contrary to widespread belief, they are frequently targets of protective antibody responses. Importantly, disordered epitopes appear invariably to be linear epitopes, rendering them ideally suited to incorporation into a peptide vaccine. Nonetheless, the conformational properties of disordered antigens, and hence their recognition by antibodies, frequently depend on the interactions they make and the context in which they are presented to the immune system. These effects must be considered in the design of an effective vaccine. Here we discuss these issues and propose design principles that may facilitate the development of peptide vaccines targeting disordered antigens.

The evolutionary origins of cell type diversification and the role of intrinsically disordered proteins

The evolution of complex multicellular life forms occurred multiple times and was attended by cell type specialization. We review seven lines of evidence indicating that intrinsically disordered/ductile proteins (IDPs) played a significant role in the evolution of multicellularity and cell type specification: (i) most eukaryotic transcription factors (TFs) and multifunctional enzymes contain disproportionately long IDP sequences (≥30 residues in length), whereas highly conserved enzymes are normally IDP region poor; (ii) ~80% of the proteome involved in development are IDPs; (iii) the majority of proteins undergoing alternative splicing (AS) of pre-mRNA contain significant IDP regions; (iv) proteins encoded by DNA regions flanking crossing-over 'hot spots' are significantly enriched in IDP regions; (v) IDP regions are disproportionately subject to combinatorial post-translational modifications (PTMs) as well as AS; (vi) proteins involved in transcription and RNA processing are enriched in IDP regions; and (vii) a strong positive correlation exists between the number of different cell types and the IDP proteome fraction across a broad spectrum of uni- and multicellular algae, plants, and animals. We argue that the multifunctionalities conferred by IDPs and the disproportionate involvement of IDPs with AS and PTMs provided a IDP-AS-PTM 'motif' that significantly contributed to the evolution of multicellularity in all major eukaryotic lineages.

Structural Analysis of Anesthetics in Complex with Soluble Proteins

Anesthetics interact with a broad range of different targets, including both soluble and membrane-bound proteins. Understanding these interactions at the molecular level requires detailed structural knowledge of anesthetic-protein complexes, and one of the most productive routes to such knowledge is X-ray crystallography. In this chapter we discuss the application of this technique to the analysis of complexes of anesthetics with soluble proteins. The model protein apoferritin is highlighted, and protocols are presented for obtaining diffraction-quality crystals of this protein in complex with different general anesthetics.

pKID Binds to KIX via an Unstructured Transition State with Nonnative Interactions

Understanding the detailed mechanism of interaction of intrinsically disordered proteins with their partners is crucial to comprehend their functions in signaling and transcription. Through its interaction with KIX, the disordered pKID region of CREB protein is central in the transcription of cAMP responsive genes, including those involved in long-term memory. Numerous simulation studies have investigated these interactions. Combined with experimental results, these can provide valuable and comprehensive understanding of the mechanisms involved. Here, we probe the transition state of this interaction experimentally through analyzing the kinetic effect of mutating both interface and solvent exposed residues in pKID. We show that very few specific interactions between pKID and KIX are required in the initial binding process. Only a small number of weak interactions are formed at the transition state, including nonnative interactions, and most of the folding occurs after the initial binding event. These properties are consistent with computational results and also the majority of experimental studies of intrinsically disordered protein coupled folding and binding in other protein systems, suggesting that these may be common features.

Not an exception to the rule: the functional significance of intrinsically disordered protein regions in enzymes

Intrinsically disordered protein regions (IDPRs) are remarkably common and have unique and important biological functions. Enzymes have long been considered an exception to the rule of protein intrinsic disorder due to the structural requirements for catalysis. Although functionally significant IDPRs have been described in several enzymes, there has been no study quantifying the extent of this phenomenon. We have conducted a multilevel computational analysis of missing regions in X-ray crystal structures in the PDB and predicted disorder in 66 representative proteomes. We found that the fraction of predicted disorder was higher in non-enzymes than enzymes, because non-enzymes were more likely to be fully disordered. However, we also found that transferases, hydrolases and enzymes with multiple assigned functional classifications were similar to non-enzymes in terms of the length of the longest continuous stretch of predicted disorder. Both eukaryotic enzymes and non-enzymes had a greater disorder content than was seen in bacteria. Disorder at the proteome level appears to emerge in response to organismic and functional complexity, and enzymes are not an exception to this rule.

CIDER: Resources to Analyze Sequence-Ensemble Relationships of Intrinsically Disordered Proteins

Intrinsically disordered proteins and regions (IDPs) represent a large class of proteins that are defined by conformational heterogeneity and lack of persistent tertiary/secondary structure. IDPs play important roles in a range of biological functions, and their dysregulation is central to numerous diseases, including neurodegeneration and cancer. The conformational ensembles of IDPs are encoded by their amino acid sequences. Here, we present two computational tools that are designed to enable rapid and high-throughput analyses of a wide range of physicochemical properties encoded by IDP sequences. The first, CIDER, is a user-friendly webserver that enables rapid analysis of IDP sequences. The second, localCIDER, is a high-performance software package that enables a wide range of analyses relevant to IDP sequences. In addition to introducing the two packages, we demonstrate the utility of these resources using examples where sequence analysis offers biophysical insights.

Transient helicity in intrinsically disordered Axin-1 studied by NMR spectroscopy and molecular dynamics simulations

Many natural proteins are, as a whole or in part, intrinsically disordered. Frequently, such intrinsically disordered regions (IDRs) undergo a transition to a defined and often helical conformation upon binding to partner molecules. The intrinsic propensity of an IDR sequence to fold into a helical conformation already in the absence of a binding partner can have a decisive influence on the binding process and affinity. Using a combination of NMR spectroscopy and molecular dynamics (MD) simulations we have investigated the tendency of regions of Axin-1, an intrinsically disordered scaffolding protein of the WNT signaling pathway, to form helices in segments interacting with binding partners. Secondary chemical shifts from NMR measurements show an increased helical population in these regions. Systematic application of MD advanced sampling approaches on peptide segments of Axin-1 reproduces the experimentally observed tendency and allows insights into the distribution of segment conformations and free energies of helix formation. The results, however, were found to dependent on the force field water model. Recent water models specifically designed for IDRs significantly reduce the predicted helical content and do not improve the agreement with experiment.

Unfoldomics of prostate cancer: on the abundance and roles of intrinsically disordered proteins in prostate cancer

Prostatic diseases such as prostate cancer and benign prostatic hyperplasia are highly prevalent among men. The number of studies focused on the abundance and roles of intrinsically disordered proteins in prostate cancer is rather limited. The goal of this study is to analyze the prevalence and degree of disorder in proteins that were previously associated with the prostate cancer pathogenesis and to compare these proteins to the entire human proteome. The analysis of these datasets provides means for drawing conclusions on the roles of disordered proteins in this common male disease. We also hope that the results of our analysis can potentially lead to future experimental studies of these proteins to find novel pathways associated with this disease.

Evaluation of the Coupled Two-Dimensional Main Chain Torsional Potential in Modeling Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) carry out crucial biological functions in essential biological processes of life. Because of the highly dynamic and conformationally heterogeneous nature of the disordered states of IDPs, molecular dynamics simulations are becoming an indispensable tool for the investigation of the conformational ensembles and dynamic properties of IDPs. Nevertheless, there is still no consensus on the most reliable force field in molecular dynamics simulations for IDPs hitherto. In this work, the recently proposed AMBER99SB^2D force field is evaluated in modeling some disordered polypeptides and proteins by checking its ability to reproduce experimental NMR data. The results highlight that when the ildn side-chain corrections are included, AMBER99SB^2D-ildn exhibits reliable results that agree with experiments compared with its predecessors, the AMBER14SB, AMBER99SB, AMBER99SB-ildn, and AMBER99SB^2D force fields, and that decreasing the overall magnitude of protein-protein interactions in favor of protein-water interactions is a key ingredient behind the improvement.

Are Charge-State Distributions a Reliable Tool Describing Molecular Ensembles of Intrinsically Disordered Proteins by Native MS?

Native mass spectrometry (MS) has become a central tool of structural proteomics, but its applicability to the peculiar class of intrinsically disordered proteins (IDPs) is still object of debate. IDPs lack an ordered tridimensional structure and are characterized by high conformational plasticity. Since they represent valuable targets for cancer and neurodegeneration research, there is an urgent need of methodological advances for description of the conformational ensembles populated by these proteins in solution. However, structural rearrangements during electrospray-ionization (ESI) or after the transfer to the gas phase could affect data obtained by native ESI-MS. In particular, charge-state distributions (CSDs) are affected by protein conformation inside ESI droplets, while ion mobility (IM) reflects protein conformation in the gas phase. This review focuses on the available evidence relating IDP solution ensembles with CSDs, trying to summarize cases of apparent consistency or discrepancy. The protein-specificity of ionization patterns and their responses to ligands and buffer conditions suggests that CSDs are imprinted to protein structural features also in the case of IDPs. Nevertheless, it seems that these proteins are more easily affected by electrospray conditions, leading in some cases to rearrangements of the conformational ensembles. Graphical Abstract ᅟ.

Sampling conformational space of intrinsically disordered proteins in explicit solvent: Comparison between well-tempered ensemble approach and solute tempering method

Intrinsically disordered proteins (IDPs) are a class of proteins that expected to be largely unstructured under physiological conditions. Due to their heterogeneous nature, experimental characterization of IDP is challenging. Temperature replica exchange molecular dynamics (T-REMD) is a widely used enhanced sampling method to probe structural characteristics of these proteins. However, its application has been hindered due to its tremendous computational cost, especially when simulating large systems in explicit solvent. Two methods, parallel tempering well-tempered ensemble (PT-WTE) and replica exchange with solute tempering (REST), have been proposed to alleviate the computational expense of T-REMD. In this work, we select three different IDP systems to compare the sampling characteristics and efficiencies of the two methods Both the two methods could efficiently sample the conformational space of IDP and yield highly consistent results for all the three IDPs. The efficiencies of the two methods: are compatible, with about 5-6 times better than the plain T-REMD. Besides, the advantages and disadvantages of each method are also discussed. Specially, the PT-WTE method could provide temperature dependent data of the system which could not be achieved by REST, while the REST method could readily be used to a part of the system, which is quite efficient to simulate some biological processes.

Chicken cathelicidins as potent intrinsically disordered biocides with antimicrobial activity against infectious pathogens

This study was performed to identify the expression patterns of the cathelicidin genes in a local chicken breed and to evaluate the antimicrobial activities of the cathelicidin peptides against pathogenic bacteria. This analysis revealed that the coding regions of CATH-1, -2, and -3 genes contain 447 bp, 465 bp, and 456 bp, respectively, and encode proteins of 148, 154, 151 amino acids, respectively. The complete amino acid sequences of the cathelicidin peptides are similar to those found in Meleagris gallopavo, Phasianus colchicus, and Coturnix coturnix, and show high sequence identity to their Columba livia and Anas platyrhynchos counterparts. In contrast, these avian peptides shared a very low sequence identity with the mammalian cathelicidins. The analysis further revealed that the cathelicidin genes are expressed in various organ and tissues. We also show that the CATH peptides 1, 2, 3 and their amide-modified structures possess potent antimicrobial activities against both Gram-positive and Gram-negative pathogens, with these bacteria being affected to different extents. The antimicrobial activities of the peptides are slightly lower than those of their amide analogs. Computational analysis revealed that pre-pro-cathelicidins are hybrid proteins that contain ordered domains and functional intrinsically disordered regions. Furthermore, high structural and sequence variability of mature cathelicidins is a strong indication of their rather disordered nature. It is likely that intrinsic disorder is needed for the multifarious functionality of these antimicrobial peptides. Our analyses indicated that cathelicidin peptides require further study to better understand their full potentials in the treatment of diseases in both humans and animals. The data obtained for synthetic avian peptides will help elucidating of their potential applications in the pharmaceutical industry.

An Integrated Mass Spectrometry Based Approach to Probe the Structure of the Full-Length Wild-Type Tetrameric p53 Tumor Suppressor

We present an integrated approach for investigating the topology of proteins through native mass spectrometry (MS) and cross-linking/MS, which we applied to the full-length wild-type p53 tetramer. For the first time, the two techniques were combined in one workflow to obtain not only structural insight in the p53 tetramer, but also information on the cross-linking efficiency and the impact of cross-linker modification on the conformation of an intrinsically disordered protein (IDP). P53 cross-linking was monitored by native MS and as such, our strategy serves as a quality control for different cross-linking reagents. Our approach can be applied to the structural investigation of various protein systems, including IDPs and large protein assemblies, which are challenging to study by the conventional methods used for protein structure characterization.

Large-scale analysis of intrinsic disorder flavors and associated functions in the protein sequence universe

Intrinsic disorder (ID) in proteins has been extensively described for the last decade; a large-scale classification of ID in proteins is mostly missing. Here, we provide an extensive analysis of ID in the protein universe on the UniProt database derived from sequence-based predictions in MobiDB. Almost half the sequences contain an ID region of at least five residues. About 9% of proteins have a long ID region of over 20 residues which are more abundant in Eukaryotic organisms and most frequently cover less than 20% of the sequence. A small subset of about 67,000 (out of over 80 million) proteins is fully disordered and mostly found in Viruses. Most proteins have only one ID, with short ID evenly distributed along the sequence and long ID overrepresented in the center. The charged residue composition of Das and Pappu was used to classify ID proteins by structural propensities and corresponding functional enrichment. Swollen Coils seem to be used mainly as structural components and in biosynthesis in both Prokaryotes and Eukaryotes. In Bacteria, they are confined in the nucleoid and in Viruses provide DNA binding function. Coils & Hairpins seem to be specialized in ribosome binding and methylation activities. Globules & Tadpoles bind antigens in Eukaryotes but are involved in killing other organisms and cytolysis in Bacteria. The Undefined class is used by Bacteria to bind toxic substances and mediate transport and movement between and within organisms in Viruses. Fully disordered proteins behave similarly, but are enriched for glycine residues and extracellular structures.

Protein intrinsic disorder within the Potyvirus genus: from proteome-wide analysis to functional annotation

Within proteins, intrinsically disordered regions (IDRs) are devoid of stable secondary and tertiary structures under physiological conditions and rather exist as dynamic ensembles of inter-converting conformers. Although ubiquitous in all domains of life, the intrinsic disorder content is highly variable in viral genomes. Over the years, functional annotations of disordered regions at the scale of the whole proteome have been conducted for several animal viruses. But to date, similar studies applied to plant viruses are still missing. Based on disorder prediction tools combined with annotation programs and evolutionary studies, we analyzed the intrinsic disorder content in Potyvirus, using a 10-species dataset representative of this genus diversity. In this paper, we revealed that: (i) the Potyvirus proteome displays high disorder content, (ii) disorder is conserved during Potyvirus evolution, suggesting a functional advantage of IDRs, (iii) IDRs evolve faster than ordered regions, and (iv) IDRs may be associated with major biological functions required for the Potyvirus cycle. Notably, the proteins P1, Coat protein (CP) and Viral genome-linked protein (VPg) display a high content of conserved disorder, enriched in specific motifs mimicking eukaryotic functional modules and suggesting strategies of host machinery hijacking. In these three proteins, IDRs are particularly conserved despite their high amino acid polymorphism, indicating a link to adaptive processes. Through this comprehensive study, we further investigate the biological relevance of intrinsic disorder in Potyvirus biology and we propose a functional annotation of potyviral proteome IDRs.

Essential functions linked with structural disorder in organisms of minimal genome

Abstract

Intrinsically disordered regions (IDRs) of proteins fulfill important regulatory roles in most organisms. However, the proteins of certain endosymbiont and intracellular pathogenic bacteria with extremely reduced genomes contain disproportionately small amounts of IDRs, consisting almost entirely of folded domains. As their genomes co-evolving with their hosts have been reduced in unrelated lineages, the proteomes of these bacteria represent independently evolved minimal protein sets. We systematically analyzed structural disorder in a representative set of such minimal organisms to see which types of functionally relevant longer IDRs are invariably retained in them. We found that a few characteristic functions are consistently linked with conformational disorder: ribosomal proteins, key components of the protein production machinery, a central coordinator of DNA metabolism and certain housekeeping chaperones seem to strictly rely on structural disorder even in genome-reduced organisms. We propose that these functions correspond to the most essential and probably also the most ancient ones fulfilled by structural disorder in cellular organisms.

Reviewers

This article was reviewed by Michael Gromiha, Zoltan Gaspari and Sandor Pongor.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-016-0149-y) contains supplementary material, which is available to authorized users.

Computational analysis of translational readthrough proteins in Drosophila and yeast reveals parallels to alternative splicing

In translational readthrough (TR) the ribosome continues extending the nascent protein beyond the first in-frame termination codon. Due to the lack of dedicated analyses of eukaryotic TR cases, the associated functional-evolutionary advantages are still unclear. Here, based on a variety of computational methods, we describe the structural and functional properties of previously proposed D. melanogaster and S. cerevisiae TR proteins and extensions. We found that in D. melanogaster TR affects long proteins in mainly regulatory roles. Their TR-extensions are structurally disordered and rich in binding motifs, which, together with their cell-type- and developmental stage-dependent inclusion, suggest that similarly to alternatively spliced exons they rewire cellular interaction networks in a temporally and spatially controlled manner. In contrast, yeast TR proteins are rather short and fulfil mainly housekeeping functions, like translation. Yeast extensions usually lack disorder and linear motifs, which precludes elucidating their functional relevance with sufficient confidence. Therefore we propose that by being much more restricted and by lacking clear functional hallmarks in yeast as opposed to fruit fly, TR shows remarkable parallels with alternative splicing. Additionally, the lack of conservation of TR extensions among orthologous TR proteins suggests that TR-mediated functions may be generally specific to lower taxonomic levels.

Order, Disorder, and Everything in Between

In addition to the “traditional” proteins characterized by the unique crystal-like structures needed for unique functions, it is increasingly recognized that many proteins or protein regions (collectively known as intrinsically disordered proteins (IDPs) and intrinsically disordered protein regions (IDPRs)), being biologically active, do not have a specific 3D-structure in their unbound states under physiological conditions. There are also subtler categories of disorder, such as conditional (or dormant) disorder and partial disorder. Both the ability of a protein/region to fold into a well-ordered functional unit or to stay intrinsically disordered but functional are encoded in the amino acid sequence. Structurally, IDPs/IDPRs are characterized by high spatiotemporal heterogeneity and exist as dynamic structural ensembles. It is important to remember, however, that although structure and disorder are often treated as binary states, they actually sit on a structural continuum.

Simulation of coupled folding and binding of an intrinsically disordered protein in explicit solvent with metadynamics

The C-terminal domain of measles virus nucleoprotein is an intrinsically disordered protein that could bind to the X domain (XD) of phosphoprotein P to exert its physiological function. Experiments reveal that the minimal binding unit is a 21-residue α-helical molecular recognition element (α-MoRE-MeV), which adopts a fully helical conformation upon binding to XD. Due to currently limited computing power, direct simulation of this coupled folding and binding process with atomic force field in explicit solvent cannot be achieved. In this work, two advanced sampling methods, metadynamics and parallel tempering, are combined to characterize the free energy surface of this process and investigate the underlying mechanism. Starting from an unbound and partially folded state of α-MoRE-MeV, multiple folding and binding events are observed during the simulation and the energy landscape was well estimated. The results demonstrate that the isolated α-MoRE-MeV resembles a molten globule and rapidly interconverts between random coil and multiple partially helical states in solution. The coupled folding and binding process occurs through the induced fit mechanism, with the residual helical conformations providing the initial binding sites. Upon binding, α-MoRE-MeV can easily fold into helical conformation without obvious energy barriers. Two mechanisms, namely, the system tending to adopt the structure in which the free energy of isolated α-MoRE-MeV is the minimum, and the binding energy of α-MoRE-MeV to its partner protein XD tending to the minimum, jointly dominate the coupled folding and binding process. With the advanced sampling approach, more IDP systems could be simulated and common mechanisms concerning the coupled folding and binding process could be investigated in the future.

Three reasons protein disorder analysis makes more sense in the light of collagen

We have identified that the collagen helix has the potential to be disruptive to analyses of intrinsically disordered proteins. The collagen helix is an extended fibrous structure that is both promiscuous and repetitive. Whilst its sequence is predicted to be disordered, this type of protein structure is not typically considered as intrinsic disorder. Here, we show that collagen-encoding proteins skew the distribution of exon lengths in genes. We find that previous results, demonstrating that exons encoding disordered regions are more likely to be symmetric, are due to the abundance of the collagen helix. Other related results, showing increased levels of alternative splicing in disorder-encoding exons, still hold after considering collagen-containing proteins. Aside from analyses of exons, we find that the set of proteins that contain collagen significantly alters the amino acid composition of regions predicted as disordered. We conclude that research in this area should be conducted in the light of the collagen helix.

The BECN1 N-terminal domain is intrinsically disordered

BECN1/Beclin 1 has a critical role in the early stages of autophagosome formation. Recently, structures of its central and C-terminal domains were reported, however, little structural information is available on the N-terminal domain, comprising a third of the protein. This lack of structural information largely stems from the inability to produce this region in a purified form. Here, we describe the expression and purification of the N-terminal domain of BECN1 (residues 1 to 150) and detailed biophysical characterization, including NMR spectroscopy. Combined, our studies demonstrated at the atomic level that the BECN1 N-terminal domain is intrinsically disordered, and apart from the BH3 subdomain, remains disordered following interaction with a binding partner, BCL2L1/BCL-XL. In addition, the BH3 domain α-helix induced upon interaction with BCL2L1 reverts to a disordered state when the complex is dissociated by exposure to a competitive inhibitor. No significant interactions between N- and C-terminal domains were detected.

Accelerating the Conformational Sampling of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) are a class of proteins lacking a well-defined secondary structure. Instead, they are able to attain multiple conformations, bind to multiple targets, and respond to changes in their surroundings. Functionally, IDPs have been associated with molecular recognition, cell regulation, and signal transduction. The dynamic conformational ensemble of IDPs is highly environmental and binding partner dependent, rendering the characterization of IDPs extremely challenging. Here, we compare the sampling efficiencies of conventional molecular dynamics (MD), well-tempered metadynamics (WT-META), and bias-exchange metadynamics (BE-META). The total simulation time was over 10 μs, and a 20-mer peptide derived from the Neh2 domain of the Nuclear factor erythroid 2-related factor 2 (Nrf2) protein was simulated. BE-META, with a neutral replica and seven biased replicas employing a set of seven relevant collective variables (CVs), provided the most reliable and efficient sampling. Finally, we propose a free-energy reconstruction method based on the probability distribution of the secondary structure contents. This postprocessing analysis confirms the presence of not only the β-hairpin conformation of the free Neh2 peptide but also its rare bound-state-like conformation, both of that have been experimentally observed. In addition, our simulations also predict other possible conformations to be verified with future experiments.

Biophysical Studies on BEX3, the p75NTR-Associated Cell Death Executor, Reveal a High-Order Oligomer with Partially Folded Regions

BEX3 (Brain Expressed X-linked protein 3) is a member of a mammal-specific placental protein family. Several studies have found the BEX proteins to be associated with neurodegeneration, the cell cycle and cancer. BEX3 has been predicted to be intrinsically disordered and also to represent an intracellular hub for cell signaling. The pro-apoptotic activity of BEX3 in association with a number of additional proteins has been widely supported; however, to the best of our knowledge, very limited data are available on the conformation of any of the members of the BEX family. In this study, we structurally characterized BEX3 using biophysical experimental data. Small angle X-ray scattering and atomic force microscopy revealed that BEX3 forms a specific higher-order oligomer that is consistent with a globular molecule. Solution nuclear magnetic resonance, partial proteinase K digestion, circular dichroism spectroscopy, and fluorescence techniques that were performed on the recombinant protein indicated that the structure of BEX3 is composed of approximately 31% α-helix and 20% β-strand, contains partially folded regions near the N- and C-termini, and a core which is proteolysis-resistant around residues 55-120. The self-oligomerization of BEX3 has been previously reported in cell culture and is consistent with our in vitro data.

Flexibility and Disorder in Gene Regulation: LacI/GalR and Hox Proteins

To modulate transcription, a variety of input signals must be sensed by genetic regulatory proteins. In these proteins, flexibility and disorder are emerging as common themes. Prokaryotic regulators generally have short, flexible segments, whereas eukaryotic regulators have extended regions that lack predicted secondary structure (intrinsic disorder). Two examples illustrate the impact of flexibility and disorder on gene regulation: the prokaryotic LacI/GalR family, with detailed information from studies on LacI, and the eukaryotic family of Hox proteins, with specific insights from investigations of Ultrabithorax (Ubx). The widespread importance of structural disorder in gene regulatory proteins may derive from the need for flexibility in signal response and, particularly in eukaryotes, in protein partner selection.

Unfoldome variation upon plant-pathogen interactions: strawberry infection by Colletotrichum acutatum

Intrinsically disordered proteins (IDPs) are proteins that lack secondary and/or tertiary structure under physiological conditions. These proteins are very abundant in eukaryotic proteomes and play crucial roles in all molecular mechanisms underlying the response to environmental challenges. In plants, different IDPs involved in stress response have been identified and characterized. Nevertheless, a comprehensive evaluation of protein disorder in plant proteomes under abiotic or biotic stresses is not available so far. In the present work the transcriptome dataset of strawberry (Fragaria X ananassa) fruits interacting with the fungal pathogen Colletotrichum acutatum was actualized onto the woodland strawberry (Fragaria vesca) genome. The obtained cDNA sequences were translated into protein sequences, which were subsequently subjected to disorder analysis. The results, providing the first estimation of disorder abundance associated to plant infection, showed that the proteome activated in the strawberry red fruit during the active fungal propagation is remarkably depleted in disorder. On the other hand, in the resistant white fruit, no significant disorder reduction is observed in the proteins expressed in response to fungal infection. Four representative proteins, FvSMP, FvPRKRIP, FvPCD-4 and FvFAM32A-like, predicted as mainly disordered and never experimentally characterized before, were isolated, and the absence of structure was validated at the secondary and tertiary level using circular dichroism and differential scanning fluorimetry. Their quaternary structure was also established using light scattering. The results are discussed considering the role of protein disorder in plant defense.

In Silico Analysis of Correlations between Protein Disorder and Post-Translational Modifications in Algae

Recent proteome analyses have reported that intrinsically disordered regions (IDRs) of proteins play important roles in biological processes. In higher plants whose genomes have been sequenced, the correlation between IDRs and post-translational modifications (PTMs) has been reported. The genomes of various eukaryotic algae as common ancestors of plants have also been sequenced. However, no analysis of the relationship to protein properties such as structure and PTMs in algae has been reported. Here, we describe correlations between IDR content and the number of PTM sites for phosphorylation, glycosylation, and ubiquitination, and between IDR content and regions rich in proline, glutamic acid, serine, and threonine (PEST) and transmembrane helices in the sequences of 20 algae proteomes. Phosphorylation, O-glycosylation, ubiquitination, and PEST preferentially occurred in disordered regions. In contrast, transmembrane helices were favored in ordered regions. N-glycosylation tended to occur in ordered regions in most of the studied algae; however, it correlated positively with disordered protein content in diatoms. Additionally, we observed that disordered protein content and the number of PTM sites were significantly increased in the species-specific protein clusters compared to common protein clusters among the algae. Moreover, there were specific relationships between IDRs and PTMs among the algae from different groups.

Disordered regions in transmembrane proteins

The functions of transmembrane proteins in living cells are widespread; they range from various transport processes to energy production, from cell-cell adhesion to communication. Structurally, they are highly ordered in their membrane-spanning regions, but may contain disordered regions in the cytosolic and extra-cytosolic parts. In this study, we have investigated the disordered regions in transmembrane proteins by a stringent definition of disordered residues on the currently available largest experimental dataset, and show a significant correlation between the spatial distributions of positively charged residues and disordered regions. This finding suggests a new role of disordered regions in transmembrane proteins by providing structural flexibility for stabilizing interactions with negatively charged head groups of the lipid molecules. We also find a preference of structural disorder in the terminal--as opposed to loop--regions in transmembrane proteins, and survey the respective functions involved in recruiting other proteins or mediating allosteric signaling effects. Finally, we critically compare disorder prediction methods on our transmembrane protein set. While there are no major differences between these methods using the usual statistics, such as per residue accuracies, Matthew's correlation coefficients, etc.; substantial differences can be found regarding the spatial distribution of the predicted disordered regions. We conclude that a predictor optimized for transmembrane proteins would be of high value to the field of structural disorder.