Distribution and cluster analysis of predicted intrinsically disordered protein Pfam domains

Intrinsic Disorder in the Human Tear Proteome

Purpose

We aimed to characterize the proteome of human tears and assess for the presence of intrinsically disordered proteins (IDPs). IDPs, despite lacking a rigid three-dimensional structure, maintain biological functionality and could shed light on the molecular interactions within tears.

Methods

We analyzed a dataset of 1475 proteins identified in the tear film of three healthy subjects. We employed several computational tools, including the Compositional Profiler, Rapid Intrinsic Disorder Analysis Online, Search Tool for the Retrieval of Interacting Genes, and Database of Disordered Protein Predictors to evaluate the intrinsic disorder, protein interactions, and functional characterization of the disordered regions within this proteome.

Results

Our analysis showed a notable inclination toward intrinsic disorder. Two out of 10 order-promoting residues and five out of 10 disorder-promoting residues were found enriched. Using the Predictor of Natural Disordered Regions (PONDR) VSL2 output, 95% of these proteins were classified as highly or moderately disordered. We revealed an extensive protein-protein interaction network with significant interaction enrichment. The most disordered proteins exhibited higher disorder binding sites and diverse posttranslational modifications compared to the most ordered ones.

Conclusions

To the best of our knowledge, our study is the first comprehensive analysis of intrinsic disorder in the human tear film proteome, and it revealed an abundance of IDPs and their role in protein function and interaction networks. These findings suggest that variations in the intrinsic disorder of a tear film could be impacted by systemic and ocular conditions, offering promising avenues for disease biomarker identification and drug target development. Further research is needed to understand the implications of these findings in human health and disease.

Digging into the 3D Structure Predictions of AlphaFold2 with Low Confidence: Disorder and Beyond

AlphaFold2 (AF2) has created a breakthrough in biology by providing three-dimensional structure models for whole-proteome sequences, with unprecedented levels of accuracy. In addition, the AF2 pLDDT score, related to the model confidence, has been shown to provide a good measure of residue-wise disorder. Here, we combined AF2 predictions with pyHCA, a tool we previously developed to identify foldable segments and estimate their order/disorder ratio, from a single protein sequence. We focused our analysis on the AF2 predictions available for 21 reference proteomes (AFDB v1), in particular on their long foldable segments (>30 amino acids) that exhibit characteristics of soluble domains, as estimated by pyHCA. Among these segments, we provided a global analysis of those with very low pLDDT values along their entire length and compared their characteristics to those of segments with very high pLDDT values. We highlighted cases containing conditional order, as well as cases that could form well-folded structures but escape the AF2 prediction due to a shallow multiple sequence alignment and/or undocumented structure or fold. AF2 and pyHCA can therefore be advantageously combined to unravel cryptic structural features in whole proteomes and to refine predictions for different flavors of disorder.

Structural and Functional Characterization of the ABA-Water Deficit Stress Domain from Wheat and Barley: An Intrinsically Disordered Domain behind the Versatile Functions of the Plant Abscissic Acid, Stress and Ripening Protein Family

The ASR protein family has been discovered thirty years ago in many plant species and is involved in the tolerance of various abiotic stresses such as dehydration, salinity and heat. Despite its importance, nothing is known about the conserved ABA-Water Deficit Stress Domain (ABA-WDS) of the ASR gene family. In this study, we characterized two ABA-WDS domains, isolated from durum wheat (TtABA-WDS) and barley (HvABA-WDS). Bioinformatics analysis shows that they are both consistently predicted to be intrinsically disordered. Hydrodynamic and circular dichroism analysis indicate that both domains are largely disordered but belong to different structural classes, with HvABA-WDS and TtABA-WDS adopting a PreMolten Globule-like (PMG-like) and a Random Coil-like (RC-like) conformation, respectively. In the presence of the secondary structure stabilizer trifluoroethanol (TFE) or of increasing glycerol concentrations, which mimics dehydration, the two domains acquire an α-helical structure. Interestingly, both domains are able to prevent heat- and dehydration-induced inactivation of the enzyme lactate dehydrogenase (LDH). Furthermore, heterologous expression of TtABA-WDS and HvABA-WDS in the yeast Saccharomyces cerevisiae improves its tolerance to salt, heat and cold stresses. Taken together our results converge to show that the ABA-WDS domain is an intrinsically disordered functional domain whose conformational plasticity could be instrumental to support the versatile functions attributed to the ASR family, including its role in abiotic stress tolerance. Finally, and after validation in the plant system, this domain could be used to improve crop tolerance to abiotic stresses.

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.

B cell epitopes in the intrinsically disordered regions of neuraminidase and hemagglutinin proteins of H5N1 and H9N2 avian influenza viruses for peptide-based vaccine development

Avian influenza viruses (AIV) are very active in several parts of the globe and are the cause of huge economic loss for the poultry industry and also human fatalities. Three dimensional modeling was carried out for neuraminidase (NA) and hemagglutinin (HA) proteins of AIV. The C-score, estimated TM-Score, and estimated root-mean-square deviation (RMSD) score for NA of H5N1 were -1.18, 0.57 ± 0.15, and 9.8 ± 7.6, respectively. The C-score, estimated TM-Score, and estimated RMSD score for NA of H9N2 were -1.43, 0.54 ± 0.15, and 10.5 ± 4.6, respectively. The C-score, estimated TM-Score, and estimated RMSD score for HA of H5N1 were -0.03, 0.71 ± 0.12, and 7.7 ± 4.3, respectively. The C-score, estimated TM-Score, and estimated RMSD score for HA of H9N2 were -0.57, 0.64 ± 0.13, and 8.9 ± 4.6, respectively. Intrinsically disordered regions were identified for the NA and HA proteins of H5N1 and H9N2 with the use of PONDR program. Linear B cell epitope was predicted using BepiPred 2 program for NA and HA of H5N1 and H9N2 avian influenza strains. Discontinuous epitopes were predicted by Discotope 2 program. The linear epitopes that were considered likely to be immunogenic and within the intrinsically disordered region for the NA of H5N1 was TKSTNSRSGFEMIWDPNGWTGTDSSFSVK, and for H9N2 it was VGDTPRNDDSSSSSNCRDPNNERGAP. In the case of HA of H5N1, it was QRLVPKIATRSKVNGQSG and ATGLRNSPQRERRRKK; for H9N2 it was INRTFKPLIGPRPLVNGLQG and SLKLAVGLRNVPARSSR. The discontinuous epitopes of NA of H5N1 and H9N2 were identified at various regions of the protein structure spanning from amino acid residue positions 90 to 449 and 107 to 469, respectively. Similarly, the discontinuous epitopes of HA of H5N1 and H9N2 were identified in the amino acid residue positions 27 to 517 and 136 to 521, respectively. This study has identified potential and highly immunogenic linear and conformational B-cell epitopes towards developing a vaccine against AIV both for human and poultry use.

Sequence fingerprints distinguish erroneous from correct predictions of intrinsically disordered protein regions

More than 60 prediction methods for intrinsically disordered proteins (IDPs) have been developed over the years, many of which are accessible on the World Wide Web. Nearly, all of these predictors give balanced accuracies in the ~65%-~80% range. Since predictors are not perfect, further studies are required to uncover the role of amino acid residues in native IDP as compared to predicted IDP regions. In the present work, we make use of sequences of 100% predicted IDP regions, false positive disorder predictions, and experimentally determined IDP regions to distinguish the characteristics of native versus predicted IDP regions. A higher occurrence of asparagine is observed in sequences of native IDP regions but not in sequences of false positive predictions of IDP regions. The occurrences of certain combinations of amino acids at the pentapeptide level provide a distinguishing feature in the IDPs with respect to globular proteins. The distinguishing features presented in this paper provide insights into the sequence fingerprints of amino acid residues in experimentally determined as compared to predicted IDP regions. These observations and additional work along these lines should enable the development of improvements in the accuracy of disorder prediction algorithm.

Polymorphism Analysis Reveals Reduced Negative Selection and Elevated Rate of Insertions and Deletions in Intrinsically Disordered Protein Regions

Intrinsically disordered protein regions are abundant in eukaryotic proteins and lack stable tertiary structures and enzymatic functions. Previous studies of disordered region evolution based on interspecific alignments have revealed an increased propensity for indels and rapid rates of amino acid substitution. How disordered regions are maintained at high abundance in the proteome and across taxa, despite apparently weak evolutionary constraints, remains unclear. Here, we use single nucleotide and indel polymorphism data in yeast and human populations to survey the population variation within disordered regions. First, we show that single nucleotide polymorphisms in disordered regions are under weaker negative selection compared with more structured protein regions and have a higher proportion of neutral non-synonymous sites. We also confirm previous findings that nonframeshifting indels are much more abundant in disordered regions relative to structured regions. We find that the rate of nonframeshifting indel polymorphism in intrinsically disordered regions resembles that of noncoding DNA and pseudogenes, and that large indels segregate in disordered regions in the human population. Our survey of polymorphism confirms patterns of evolution in disordered regions inferred based on longer evolutionary comparisons.

Prediction of protein structural features from sequence data based on Shannon entropy and Kolmogorov complexity

While the genome for a given organism stores the information necessary for the organism to function and flourish it is the proteins that are encoded by the genome that perhaps more than anything else characterize the phenotype for that organism. It is therefore not surprising that one of the many approaches to understanding and predicting protein folding and properties has come from genomics and more specifically from multiple sequence alignments. In this work I explore ways in which data derived from sequence alignment data can be used to investigate in a predictive way three different aspects of protein structure: secondary structures, inter-residue contacts and the dynamics of switching between different states of the protein. In particular the use of Kolmogorov complexity has identified a novel pathway towards achieving these goals.

Identifying novel cell cycle proteins in Apicomplexa parasites through co-expression decision analysis

Hypothetical proteins comprise roughly half of the predicted gene complement of Toxoplasma gondii and Plasmodium falciparum and represent the largest class of uniquely functioning proteins in these parasites. Following the idea that functional relationships can be informed by the timing of gene expression, we devised a strategy to identify the core set of apicomplexan cell division cycling genes with important roles in parasite division, which includes many uncharacterized proteins. We assembled an expanded list of orthologs from the T. gondii and P. falciparum genome sequences (2781 putative orthologs), compared their mRNA profiles during synchronous replication, and sorted the resulting set of dual cell cycle regulated orthologs (744 total) into protein pairs conserved across many eukaryotic families versus those unique to the Apicomplexa. The analysis identified more than 100 ortholog gene pairs with unknown function in T. gondii and P. falciparum that displayed co-conserved mRNA abundance, dynamics of cyclical expression and similar peak timing that spanned the complete division cycle in each parasite. The unknown cyclical mRNAs encoded a diverse set of proteins with a wide range of mass and showed a remarkable conservation in the internal organization of ordered versus disordered structural domains. A representative sample of cyclical unknown genes (16 total) was epitope tagged in T. gondii tachyzoites yielding the discovery of new protein constituents of the parasite inner membrane complex, key mitotic structures and invasion organelles. These results demonstrate the utility of using gene expression timing and dynamic profile to identify proteins with unique roles in Apicomplexa biology.

Intrinsically disordered proteins and intrinsically disordered protein regions

Intrinsically disordered proteins (IDPs) and IDP regions fail to form a stable structure, yet they exhibit biological activities. Their mobile flexibility and structural instability are encoded by their amino acid sequences. They recognize proteins, nucleic acids, and other types of partners; they accelerate interactions and chemical reactions between bound partners; and they help accommodate posttranslational modifications, alternative splicing, protein fusions, and insertions or deletions. Overall, IDP-associated biological activities complement those of structured proteins. Recently, there has been an explosion of studies on IDP regions and their functions, yet the discovery and investigation of these proteins have a long, mostly ignored history. Along with recent discoveries, we present several early examples and the mechanisms by which IDPs contribute to function, which we hope will encourage comprehensive discussion of IDPs and IDP regions in biochemistry textbooks. Finally, we propose future directions for IDP research.