Intrinsic Disorder Content Research Articles

Disordered region encodes α-crystallin chaperone activity toward lens client γD-crystallin

Small heat-shock proteins (sHSPs) are a widely expressed family of ATP-independent molecular chaperones that are among the first responders to cellular stress. Mechanisms by which sHSPs delay aggregation of client proteins remain undefined. sHSPs have high intrinsic disorder content of up to ~60% and assemble into large, polydisperse homo- and hetero-oligomers, making them challenging structural and biochemical targets. Two sHSPs, HSPB4 and HSPB5, are present at millimolar concentrations in eye lens, where they are responsible for maintaining lens transparency over the lifetime of an organism. Together, HSPB4 and HSPB5 compose the hetero-oligomeric chaperone known as α-crystallin. To identify the determinants of sHSP function, we compared the effectiveness of HSPB4 and HSPB5 homo-oligomers and HSPB4/HSPB5 hetero-oligomers in delaying the aggregation of the lens protein γD-crystallin. In chimeric versions of HSPB4 and HSPB5, chaperone activity tracked with the identity of the 60-residue disordered N-terminal regions (NTR). A short 10-residue stretch in the middle of the NTR ("Critical sequence") contains three residues that are responsible for high HSPB5 chaperone activity toward γD-crystallin. These residues affect structure and dynamics throughout the NTR. Abundant interactions involving the NTR Critical sequence reveal it to be a hub for a network of interactions within oligomers. We propose a model whereby the NTR critical sequence influences local structure and NTR dynamics that modulate accessibility of the NTR, which in turn modulates chaperone activity.

A Thermodynamic Atlas of Proteomes Reveals Energetic Innovation across the Tree of Life.

Protein stability is a fundamental molecular property enabling organisms to adapt to their biological niches. How this is facilitated and whether there are kingdom specific or more general universal strategies are unknown. A principal obstacle to addressing this issue is that the vast majority of proteins lack annotation, specifically thermodynamic annotation, beyond the amino acid and chromosome information derived from genome sequencing. To address this gap and facilitate future investigation into large-scale patterns of protein stability and dynamics within and between organisms, we applied a unique ensemble-based thermodynamic characterization of protein folds to a substantial portion of extant sequenced genomes. Using this approach, we compiled a database resource focused on the position-specific variation in protein stability. Interrogation of the database reveals: 1) domains of life exhibit distinguishing thermodynamic features, with eukaryotes particularly different from both archaea and bacteria; 2) the optimal growth temperature of an organism is proportional to the average apolar enthalpy of its proteome; 3) intrinsic disorder content is also proportional to the apolar enthalpy (but unexpectedly not the predicted stability at 25 °C); and 4) secondary structure and global stability information of individual proteins is extractable. We hypothesize that wider access to residue-specific thermodynamic information of proteomes will result in deeper understanding of mechanisms driving functional adaptation and protein evolution. Our database is free for download at https://afc-science.github.io/thermo-env-atlas/ (last accessed January 18, 2022).

Genomic analysis of intrinsically disordered proteins in cereals: From mining to meaning

Intrinsically disordered proteins (IDPs) are highly abundant in eukaryotic proteomes and involved in key biological and cellular processes. Although some resources of disordered protein predictions are available from animal and plant proteomes, those related to cereals are largely unknown. Here, we present an overview of IDPomes from Oryza sativa, Zea mays, Sorghum bicolor and Brachypodium distachyon. The work includes a comparative analysis with the model plant Arabidopsis thaliana. The data show that the intrinsic disorder content increases with the proteome size. Gene Ontology analysis reveals that IDPs in the studied species are involved mainly in regulation of cellular and metabolic processes and responses to stimulus. Our findings strongly suggest that higher plants may use common cellular and regulatory mechanisms for adaptation to various environmental constraints.

A Suggestion of Converting Protein Intrinsic Disorder to Structural Entropy Using Shannon's Information Theory.

We propose a framework to convert the protein intrinsic disorder content to structural entropy (H) using Shannon’s information theory (IT). The structural capacity (C), which is the sum of H and structural information (I), is equal to the amino acid sequence length of the protein. The structural entropy of the residues expands a continuous spectrum, ranging from 0 (fully ordered) to 1 (fully disordered), consistent with Shannon’s IT, which scores the fully-determined state 0 and the fully-uncertain state 1. The intrinsically disordered proteins (IDPs) in a living cell may participate in maintaining the high-energy-low-entropy state. In addition, under this framework, the biological functions performed by proteins and associated with the order or disorder of their 3D structures could be explained in terms of information-gains or entropy-losses, or the reverse processes.

GENETIC OVERLAP AND CAUSALITY AMONG MAJOR DEPRESSIVE DISORDER, ALCOHOL DEPENDENCE, AND ALCOHOL CONSUMPTION: FINDINGS FROM THE PSYCHIATRIC GENOMICS CONSORTIUM

Intrinsically disordered proteins (IDPs) are highly abundant in eukaryotic proteomes and involved in key biological and cellular processes. Although some resources of disordered protein predictions are available from animal and plant proteomes, those related to cereals are largely unknown. Here, we present an overview of IDPomes from Oryza sativa, Zea mays, Sorghum bicolor and Brachypodium distachyon. The work includes a comparative analysis with the model plant Arabidopsis thaliana. The data show that the intrinsic disorder content increases with the proteome size. Gene Ontology analysis reveals that IDPs in the studied species are involved mainly in regulation of cellular and metabolic processes and responses to stimulus. Our findings strongly suggest that higher plants may use common cellular and regulatory mechanisms for adaptation to various environmental constraints.

The rheumatic disease-associated FAM167A-BLK locus encodes DIORA-1, a novel disordered protein expressed highly in bronchial epithelium and alveolar macrophages.

Triggering of autoimmunity that leads to rheumatic disease has been suggested to depend upon gene-environment interactions occurring in epithelial barriers and associated immune cells. Genetic studies have identified associations of the FAM167A-BLK locus with rheumatoid arthritis, systemic lupus erythematosus (SLE) and Sjögren's syndrome. While BLK (B lymphocyte kinase) has a well-established role in B cells, family with sequence similarity to 167 member A (FAM167A) and its gene family remain uncharacterized. To begin to understand the role of FAM167A in rheumatic disease pathogenesis, we explored this gene family and cloned and investigated the gene products. Expression of quantitative trait locus analysis was performed in immune cells. FAM167A and FAM167B were cloned from human peripheral blood mononuclear cells (PBMC). Gene conservation and protein properties were analysed by online tools, mRNA expression measured in mouse organs by quantitative polymerase chain reaction (qPCR) and protein expression investigated in human tissues by immunohistochemistry. We found that autoimmune risk genotypes within the FAM167A-BLK locus lead to increased expression of FAM167A. The FAM167 gene family includes two members, FAM167A and FAM167B, which are not homologous to any other annotated gene but are evolutionarily conserved. The encoded proteins, which we denote 'disordered autoimmunity' (DIORA)-1 and DIORA-2, respectively, are characterized by a high content of intrinsic disorder. Notably, DIORA-1 has its highest expression in the lung, detectable in both bronchial epithelium and alveolar macrophages with an endosomal localization pattern. In summary, the FAM167A gene is associated with several rheumatic diseases and encodes a novel disordered protein, DIORA-1, which is expressed highly in the lung, consistent with a potential role in disease pathogenesis.

First Experimental Assessment of Protein Intrinsic Disorder Involvement in an RNA Virus Natural Adaptive Process.

Intrinsic disorder (ID) in proteins is defined as a lack of stable structure in physiological conditions. Intrinsically disordered regions (IDRs) are highly abundant in some RNA virus proteomes. Low topological constraints exerted on IDRs are expected to buffer the effect of numerous deleterious mutations and could be related to the remarkable adaptive potential of RNA viruses to overcome resistance of their host. To experimentally test this hypothesis in a natural pathosystem, a set of four variants of Potato virus Y (PVY; Potyvirus genus) containing various ID degrees in the Viral genome-linked (VPg) protein, a key determinant of potyvirus adaptation, was designed. To estimate the ID contribution to the VPg-based PVY adaptation, the adaptive ability of the four PVY variants was monitored in the pepper host (Capsicum annuum) carrying a recessive resistance gene. Intriguingly, the two mutants with the highest ID content displayed a significantly higher ability to restore infection in the resistant host, whereas the less intrinsically disordered mutant was unable to restore infection. The role of ID on virus adaptation may be due either to a larger exploration of evolutionary pathways or the minimization of fitness penalty caused by resistance-breaking mutations. This pioneering study strongly suggests the positive impact of ID in an RNA virus adaptive capacity.

Content of intrinsic disorder influences the outcome of cell-free protein synthesis.

Cell-free protein synthesis is used to produce proteins with various structural traits. Recent bioinformatics analyses indicate that more than half of eukaryotic proteins possess long intrinsically disordered regions. However, no systematic study concerning the connection between intrinsic disorder and expression success of cell-free protein synthesis has been presented until now. To address this issue, we examined correlations of the experimentally observed cell-free protein expression yields with the contents of intrinsic disorder bioinformatically predicted in the expressed sequences. This analysis revealed strong relationships between intrinsic disorder and protein amenability to heterologous cell-free expression. On the one hand, elevated disorder content was associated with the increased ratio of soluble expression. On the other hand, overall propensity for detectable protein expression decreased with disorder content. We further demonstrated that these tendencies are rooted in some distinct features of intrinsically disordered regions, such as low hydrophobicity, elevated surface accessibility and high abundance of sequence motifs for proteolytic degradation, including sites of ubiquitination and PEST sequences. Our findings suggest that identification of intrinsically disordered regions in the expressed amino acid sequences can be of practical use for predicting expression success and optimizing cell-free protein synthesis.

Elucidating evolutionary features and functional implications of orphan genes in Leishmania major

Orphan genes are protein coding genes that lack recognizable homologs in other organisms. These genes were reported to comprise a considerable fraction of coding regions in all sequenced genomes and thought to be allied with organism’s lineage-specific traits. However, their evolutionary persistence and functional significance still remain elusive. Due to lack of homologs with the host genome and for their probable lineage-specific functional roles, orphan gene product of pathogenic protozoan might be considered as the possible therapeutic targets. Leishmania major is an important parasitic protozoan of the genus Leishmania that is associated with the disease cutaneous leishmaniasis. Therefore, evolutionary and functional characterization of orphan genes in this organism may help in understanding the factors prevailing pathogen evolution and parasitic adaptation. In this study, we systematically identified orphan genes of L. major and employed several in silico analyses for understanding their evolutionary and functional attributes. To trace the signatures of molecular evolution, we compared their evolutionary rate with non-orphan genes. In agreement with prior observations, here we noticed that orphan genes evolve at a higher rate as compared to non-orphan genes. Lower sequence conservation of orphan genes was previously attributed solely due to their younger gene age. However, here we observed that together with gene age, a number of genomic (like expression level, GC content, variation in codon usage) and proteomic factors (like protein length, intrinsic disorder content, hydropathicity) could independently modulate their evolutionary rate. We considered the interplay of all these factors and analyzed their relative contribution on protein evolutionary rate by regression analysis. On the functional level, we observed that orphan genes are associated with regulatory, growth factor and transport related processes. Moreover, these genes were found to be enriched with various types of interaction and trafficking motifs, implying their possible involvement in host–parasite interactions. Thus, our comprehensive analysis of L. major orphan genes provided evidence for their extensive roles in host–pathogen interactions and virulence.

Protein intrinsic disorder and network connectivity. The case of 14-3-3 proteins

The understanding of networks is a common goal of an unprecedented array of traditional disciplines. One of the protein network properties most influenced by the structural contents of its nodes is the inter-connectivity. Recent studies in which structural information was included into the topological analysis of protein networks revealed that the content of intrinsic disorder in the nodes could modulate the network topology, rewire networks, and change their inter-connectivity, which is defined by its clustering coefficient. Here, we review the role of intrinsic disorder present in the partners of the highly conserved 14-3-3 protein family on its interaction networks. The 14-3-3s are phospho-serine/threonine binding proteins that have strong influence in the regulation of metabolism and signal transduction networks. Intrinsic disorder increases the clustering coefficients, namely the inter-connectivity of the nodes within each 14-3-3 paralog networks. We also review two new ideas to measure intrinsic disorder independently of the primary sequence of proteins, a thermodynamic model and a method that uses protein structures and their solvent environment. This new methods could be useful to explain unsolved questions about versatility and fixation of intrinsic disorder through evolution. The relation between the intrinsic disorder and network topologies could be an interesting model to investigate new implicitness of the graph theory into biology.

RAPID: Fast and accurate sequence-based prediction of intrinsic disorder content on proteomic scale

Recent research in the protein intrinsic disorder was stimulated by the availability of accurate computational predictors. However, most of these methods are relatively slow, especially considering proteome-scale applications, and were shown to produce relatively large errors when estimating disorder at the protein- (in contrast to residue-) level, which is defined by the fraction/content of disordered residues. To this end, we propose a novel support vector Regression-based Accurate Predictor of Intrinsic Disorder (RAPID). Key advantages of RAPID are speed (prediction of an average-size eukaryotic proteome takes <1h on a modern desktop computer); sophisticated design (multiple, complementary information sources that are aggregated over an input chain are combined using feature selection); and high-quality and robust predictive performance. Empirical tests on two diverse benchmark datasets reveal that RAPID's predictive performance compares favorably to a comprehensive set of state-of-the-art disorder and disorder content predictors. Drawing on high speed and good predictive quality, RAPID was used to perform large-scale characterization of disorder in 200+ fully sequenced eukaryotic proteomes. Our analysis reveals interesting relations of disorder with structural coverage and chain length, and unusual distribution of fully disordered chains. We also performed a comprehensive (using 56000+ annotated chains, which doubles the scope of previous studies) investigation of cellular functions and localizations that are enriched in the disorder in the human proteome. RAPID, which allows for batch (proteome-wide) predictions, is available as a web server at http://biomine.ece.ualberta.ca/RAPID/.

Insights into the Evolutionary Features of Human Neurodegenerative Diseases

Comparative analyses between human disease and non-disease genes are of great interest in understanding human disease gene evolution. However, the progression of neurodegenerative diseases (NDD) involving amyloid formation in specific brain regions is still unknown. Therefore, in this study, we mainly focused our analysis on the evolutionary features of human NDD genes with respect to non-disease genes. Here, we observed that human NDD genes are evolutionarily conserved relative to non-disease genes. To elucidate the conserved nature of NDD genes, we incorporated the evolutionary attributes like gene expression level, number of regulatory miRNAs, protein connectivity, intrinsic disorder content and relative aggregation propensity in our analysis. Our studies demonstrate that NDD genes have higher gene expression levels in favor of their lower evolutionary rates. Additionally, we observed that NDD genes have higher number of different regulatory miRNAs target sites and also have higher interaction partners than the non-disease genes. Moreover, miRNA targeted genes are known to have higher disorder content. In contrast, our analysis exclusively established that NDD genes have lower disorder content. In favor of our analysis, we found that NDD gene encoded proteins are enriched with multi interface hubs (party hubs) with lower disorder contents. Since, proteins with higher disorder content need to adapt special structure to reduce their aggregation propensity, NDD proteins found to have elevated relative aggregation propensity (RAP) in support of their lower disorder content. Finally, our categorical regression analysis confirmed the underlined relative dominance of protein connectivity, 3′UTR length, RAP, nature of hubs (singlish/multi interface) and disorder content for such evolutionary rates variation between human NDD genes and non-disease genes.

Orderly order in protein intrinsic disorder distribution: disorder in 3500 proteomes from viruses and the three domains of life

Intrinsically disordered proteins and intrinsically disordered protein regions are highly abundant in nature. However, the quantitative and qualitative measures of protein intrinsic disorder in species with known genomes are still not available. Furthermore, although the correlation between high fraction of disordered residues and advanced species has been reported, the details of this correlation and the connection between the disorder content and proteome complexity have not been reported as of yet. To fill this gap, we analysed entire proteomes of 3484 species from three domains of life (archaea, bacteria and eukaryotes) and from viruses. Our analysis revealed that the evolution process is characterized by distinctive patterns of changes in the protein intrinsic disorder content. We are showing here that viruses are characterized by the widest spread of the proteome disorder content (the percentage of disordered residues ranges from 7.3% in human coronavirus NL63 to 77.3% in Avian carcinoma virus). For several organisms, a clear correlation is seen between their disorder contents and habitats. In multicellular eukaryotes, there is a weak correlation between the complexity of an organism (evaluated as a number of different cell types) and its overall disorder content. For both the prokaryotes and eukaryotes, the disorder content is generally independent of the proteome size. However, disorder shows a sharp increase associated with the transition from prokaryotic to eukaryotic cells. This suggests that the increased disorder content in eukaryotic proteomes might be used by nature to deal with the increased cell complexity due to the appearance of the various cellular compartments.

Abundance of Intrinsic Disorder in Protein Associated with Cardiovascular Disease

Evidence that many protein regions and even entire proteins lacking stable tertiary and/or secondary structure in solution (i.e., intrinsically disordered proteins) might be involved in protein-protein interactions, regulation, recognition, and signal transduction is rapidly accumulating. These signaling proteins play a crucial role in the development of several pathological conditions, including cancer. To test a hypothesis that intrinsic disorder is also abundant in cardiovascular disease (CVD), a data set of 487 CVD-related proteins was extracted from SWISS-PROT. CVD-related proteins are depleted in major order-promoting residues (Trp, Phe, Tyr, Ile, and Val) and enriched in some disorder-promoting residues (Arg, Gln, Ser, Pro, and Glu). The application of a neural network predictor of natural disordered regions (PONDR VL-XT) together with cumulative distribution function (CDF) analysis, charge-hydropathy plot (CH plot) analysis, and alpha-helical molecular recognition feature (alpha-MoRF) indicator revealed that CVD-related proteins are enriched in intrinsic disorder. In fact, the percentage of proteins with 30 or more consecutive residues predicted by PONDR VL-XT to be disordered was 57 +/- 4% for CVD-associated proteins. This value is close that described earlier for signaling proteins (66 +/- 6%) and is significantly larger than the content of intrinsic disorder in eukaryotic proteins from SWISS-PROT (47 +/- 4%) and in nonhomologous protein segments with a well-defined three-dimensional structure (13 +/- 4%). Furthermore, CDF and CH-plot analyses revealed that 120 and 36 CVD-related proteins, respectively, are wholly disordered. This high level of intrinsic disorder could be important for the function of CVD-related proteins and for the control and regulation of processes associated with cardiovascular disease. In agreement with this hypothesis, 198 alpha-MoRFs were predicted in 101 proteins from the CVD data set. A comparison of disorder predictions with the experimental structural and functional data for a subset of the CVD-associated proteins indicated good agreement between predictions and observations. Thus, our data suggest that intrinsically disordered proteins might play key roles in cardiovascular disease.