CLIP: accurate prediction of disordered linear interacting peptides from protein sequences using co-evolutionary information.

One of key features of intrinsically disordered regions (IDRs) is facilitation of protein-protein and protein-nucleic acids interactions. These disordered binding regions include molecular recognition features (MoRFs), short linear motifs (SLiMs) and longer binding domains. Vast majority of current predictors of disordered binding regions target MoRFs, with a handful of methods that predict SLiMs and disordered protein-binding domains. A new and broader class of disordered binding regions, linear interacting peptides (LIPs), was introduced recently and applied in the MobiDB resource. LIPs are segments in protein sequences that undergo disorder-to-order transition upon binding to a protein or a nucleic acid, and they cover MoRFs, SLiMs and disordered protein-binding domains. Although current predictors of MoRFs and disordered protein-binding regions could be used to identify some LIPs, there are no dedicated sequence-based predictors of LIPs. To this end, we introduce CLIP, a new predictor of LIPs that utilizes robust logistic regression model to combine three complementary types of inputs: co-evolutionary information derived from multiple sequence alignments, physicochemical profiles and disorder predictions. Ablation analysis suggests that the co-evolutionary information is particularly useful for this prediction and that combining the three inputs provides substantial improvements when compared to using these inputs individually. Comparative empirical assessments using low-similarity test datasets reveal that CLIP secures area under receiver operating characteristic curve (AUC) of 0.8 and substantially improves over the results produced by the closest current tools that predict MoRFs and disordered protein-binding regions. The webserver of CLIP is freely available at http://biomine.cs.vcu.edu/servers/CLIP/ and the standalone code can be downloaded from http://yanglab.qd.sdu.edu.cn/download/CLIP/.

Intrinsic disorder in proteins associated with oxidative stress-induced JNK signaling.

The c-Jun N-terminal kinase (JNK) signaling cascade is a mitogen-activated protein kinase (MAPK) signaling pathway that can be activated in response to a wide range of environmental stimuli. Based on the type, degree, and duration of the stimulus, the JNK signaling cascade dictates the fate of the cell by influencing gene expression through its substrate transcription factors. Oxidative stress is a result of a disturbance in the pro-oxidant/antioxidant homeostasis of the cell and is associated with a large number of diseases, such as neurodegenerative disorders, cancer, diabetes, cardiovascular diseases, and disorders of the immune system, where it activates the JNK signaling pathway. Among different biological roles ascribed to the intrinsically disordered proteins (IDPs) and hybrid proteins containing ordered domains and intrinsically disordered protein regions (IDPRs) are signaling hub functions, as intrinsic disorder allows proteins to undertake multiple interactions, each with a different consequence. In order to ensure precise signaling, the cellular abundance of IDPs is highly regulated, and mutations or changes in abundance of IDPs/IDPRs are often associated with disease. In this study, we have used a combination of six disorder predictors to evaluate the presence of intrinsic disorder in proteins of the oxidative stress-induced JNK signaling cascade, and as per our findings, none of the 18 proteins involved in this pathway are ordered. The highest level of intrinsic disorder was observed in the scaffold proteins, JIP1, JIP2, JIP3; dual specificity phosphatases, MKP5, MKP7; 14-3-3ζ and transcription factor c-Jun. The MAP3Ks, MAP2Ks, MAPKs, TRAFs, and thioredoxin were the proteins that were predicted to be moderately disordered. Furthermore, to characterize the predicted IDPs/IDPRs in the proteins of the JNK signaling cascade, we identified the molecular recognition features (MoRFs), posttranslational modification (PTM) sites, and short linear motifs (SLiMs) associated with the disordered regions. These findings will serve as a foundation for experimental characterization of disordered regions in these proteins, which represents a crucial step for a better understanding of the roles of IDPRs in diseases associated with this important pathway.

Understanding COVID-19 via comparative analysis of dark proteomes of SARS-CoV-2, human SARS and bat SARS-like coronaviruses.

The recently emerged coronavirus designated as SARS-CoV-2 (also known as 2019 novel coronavirus (2019-nCoV) or Wuhan coronavirus) is a causative agent of coronavirus disease 2019 (COVID-19), which is rapidly spreading throughout the world now. More than 1.21 million cases of SARS-CoV-2 infection and more than 67,000 COVID-19-associated mortalities have been reported worldwide till the writing of this article, and these numbers are increasing every passing hour. The World Health Organization (WHO) has declared the SARS-CoV-2 spread as a global public health emergency and admitted COVID-19 as a pandemic now. Multiple sequence alignment data correlated with the already published reports on SARS-CoV-2 evolution indicated that this virus is closely related to the bat severe acute respiratory syndrome-like coronavirus (bat SARS-like CoV) and the well-studied human SARS coronavirus (SARS-CoV). The disordered regions in viral proteins are associated with the viral infectivity and pathogenicity. Therefore, in this study, we have exploited a set of complementary computational approaches to examine the dark proteomes of SARS-CoV-2, bat SARS-like, and human SARS CoVs by analysing the prevalence of intrinsic disorder in their proteins. According to our findings, SARS-CoV-2 proteome contains very significant levels of structural order. In fact, except for nucleocapsid, Nsp8, and ORF6, the vast majority of SARS-CoV-2 proteins are mostly ordered proteins containing less intrinsically disordered protein regions (IDPRs). However, IDPRs found in SARS-CoV-2 proteins are functionally important. For example, cleavage sites in its replicase 1ab polyprotein are found to be highly disordered, and almost all SARS-CoV-2 proteins contains molecular recognition features (MoRFs), which are intrinsic disorder-based protein-protein interaction sites that are commonly utilized by proteins for interaction with specific partners. The results of our extensive investigation of the dark side of SARS-CoV-2 proteome will have important implications in understanding the structural and non-structural biology of SARS or SARS-like coronaviruses.

The dark side of Alzheimer's disease: unstructured biology of proteins from the amyloid cascade signaling pathway.

Alzheimer's disease (AD) is a leading cause of age-related dementia worldwide. Despite more than a century of intensive research, we are not anywhere near the discovery of a cure for this disease or a way to prevent its progression. Among the various molecular mechanisms proposed for the description of the pathogenesis and progression of AD, the amyloid cascade hypothesis, according to which accumulation of a product of amyloid precursor protein (APP) cleavage, amyloid ß (Aß) peptide, induces pathological changes in the brain observed in AD, occupies a unique niche. Although multiple proteins have been implicated in this amyloid cascade signaling pathway, their structure-function relationships are mostly unexplored. However, it is known that two major proteins related to AD pathology, Aß peptide, and microtubule-associated protein tau belong to the category of intrinsically disordered proteins (IDPs), which are the functionally important proteins characterized by a lack of fixed, ordered three-dimensional structure. IDPs and intrinsically disordered protein regions (IDPRs) play numerous vital roles in various cellular processes, such as signaling, cell cycle regulation, macromolecular recognition, and promiscuous binding. However, the deregulation and misfolding of IDPs may lead to disturbed signaling, interactions, and disease pathogenesis. Often, molecular recognition-related IDPs/IDPRs undergo disorder-to-order transition upon binding to their biological partners and contain specific disorder-based binding motifs, known as molecular recognition features (MoRFs). Knowing the intrinsic disorder status and disorder-based functionality of proteins associated with amyloid cascade signaling pathway may help to untangle the mechanisms of AD pathogenesis and help identify therapeutic targets. In this paper, we have used multiple computational tools to evaluate the presence of intrinsic disorder and MoRFs in 27 proteins potentially relevant to the amyloid cascade signaling pathway. Among these, BIN1, APP, APOE, PICALM, PSEN1 and CD33 were found to be highly disordered. Furthermore, their disorder-based binding regions and associated short linear motifs have also been identified. These findings represent important foundation for the future research, and experimental characterization of disordered regions in these proteins is required to better understand their roles in AD pathogenesis.

OPAL+: Length-Specific MoRF Prediction in Intrinsically Disordered Protein Sequences.

Intrinsically disordered proteins (IDPs) contain long unstructured regions, which play an important role in their function. These intrinsically disordered regions (IDRs) participate in binding events through regions called molecular recognition features (MoRFs). Computational prediction of MoRFs helps identify the potentially functional regions in IDRs. In this study, OPAL+, a novel MoRF predictor, is presented. OPAL+ uses separate models to predict MoRFs of varying lengths along with incorporating the hidden Markov model (HMM) profiles and physicochemical properties of MoRFs and their flanking regions. Together, these features help OPAL+ achieve a marginal performance improvement of 0.4-0.7% over its predecessor for diverse MoRF test sets. This performance improvement comes at the expense of increased run time as a result of the requirement of HMM profiles. OPAL+ is available for download at https://github.com/roneshsharma/OPAL-plus/wiki/OPAL-plus-Download.

Computational prediction of MoRFs based on protein sequences and minimax probability machine.

BACKGROUND: Molecular recognition features (MoRFs) are one important type of disordered segments that can promote specific protein-protein interactions. They are located within longer intrinsically disordered regions (IDRs), and undergo disorder-to-order transitions upon binding to their interaction partners. The functional importance of MoRFs and the limitation of experimental identification make it necessary to predict MoRFs accurately with computational methods. RESULTS: In this study, a new sequence-based method, named as MoRFMPM, is proposed for predicting MoRFs. MoRFMPM uses minimax probability machine (MPM) to predict MoRFs based on 16 features and 3 different windows, which neither relying on other predictors nor calculating the properties of the surrounding regions of MoRFs separately. Comparing with ANCHOR, MoRFpred and MoRFCHiBi on the same test sets, MoRFMPM not only obtains higher AUC, but also obtains higher TPR at low FPR. CONCLUSIONS: The features used in MoRFMPM can effectively predict MoRFs, especially after preprocessing. Besides, MoRFMPM uses a linear classification algorithm and does not rely on results of other predictors which makes it accessible and repeatable.

Computational Disorder Analysis in Ethylene Response Factors Uncovers Binding Motifs Critical to Their Diverse Functions.

APETALA2/ETHYLENE RESPONSE FACTOR transcription factors (AP2/ERFs) play crucial roles in adaptation to stresses such as those caused by pathogens, wounding and cold. Although their name suggests a specific role in ethylene signalling, some ERF members also co-ordinate signals regulated by other key plant stress hormones such as jasmonate, abscisic acid and salicylate. We analysed a set of ERF proteins from three divergent plant species for intrinsically disorder regions containing conserved segments involved in protein-protein interaction known as Molecular Recognition Features (MoRFs). Then we correlated the MoRFs identified with a number of known functional features where these could be identified. Our analyses suggest that MoRFs, with plasticity in their disordered surroundings, are highly functional and may have been shuffled between related protein families driven by selection. A particularly important role may be played by the alpha helical component of the structured DNA binding domain to permit specificity. We also present examples of computationally identified MoRFs that have no known function and provide a valuable conceptual framework to link both disordered and ordered structural features within this family to diverse function.

Prediction of MoRFs in Protein Sequences with MLPs Based on Sequence Properties and Evolution Information.

Molecular recognition features (MoRFs) are one important type of intrinsically disordered proteins functional regions that can undergo a disorder-to-order transition through binding to their interaction partners. Prediction of MoRFs is crucial, as the functions of MoRFs are associated with many diseases and can therefore become the potential drug targets. In this paper, a method of predicting MoRFs is developed based on the sequence properties and evolutionary information. To this end, we design two distinct multi-layer perceptron (MLP) neural networks and present a procedure to train them. We develop a preprocessing process which exploits different sizes of sliding windows to capture various properties related to MoRFs. We then use the Bayes rule together with the outputs of two trained MLP neural networks to predict MoRFs. In comparison to several state-of-the-art methods, the simulation results show that our method is competitive.

The Disordered C-Terminus of Yeast Hsf1 Contains a Cryptic Low-Complexity Amyloidogenic Region.

Response mechanisms to external stress rely on networks of proteins able to activate specific signaling pathways to ensure the maintenance of cell proteostasis. Many of the proteins mediating this kind of response contain intrinsically disordered regions, which lack a defined structure, but still are able to interact with a wide range of clients that modulate the protein function. Some of these interactions are mediated by specific short sequences embedded in the longer disordered regions. Because the physicochemical properties that promote functional and abnormal interactions are similar, it has been shown that, in globular proteins, aggregation-prone and binding regions tend to overlap. It could be that the same principle applies for disordered protein regions. In this context, we show here that a predicted low-complexity interacting region in the disordered C-terminus of the stress response master regulator heat shock factor 1 (Hsf1) protein corresponds to a cryptic amyloid region able to self-assemble into fibrillary structures resembling those found in neurodegenerative disorders.

Predicting MoRFs in protein sequences using HMM profiles.

BACKGROUND: Intrinsically Disordered Proteins (IDPs) lack an ordered three-dimensional structure and are enriched in various biological processes. The Molecular Recognition Features (MoRFs) are functional regions within IDPs that undergo a disorder-to-order transition on binding to a partner protein. Identifying MoRFs in IDPs using computational methods is a challenging task. METHODS: In this study, we introduce hidden Markov model (HMM) profiles to accurately identify the location of MoRFs in disordered protein sequences. Using windowing technique, HMM profiles are utilised to extract features from protein sequences and support vector machines (SVM) are used to calculate a propensity score for each residue. Two different SVM kernels with high noise tolerance are evaluated with a varying window size and the scores of the SVM models are combined to generate the final propensity score to predict MoRF residues. The SVM models are designed to extract maximal information between MoRF residues, its neighboring regions (Flanks) and the remainder of the sequence (Others). RESULTS: To evaluate the proposed method, its performance was compared to that of other MoRF predictors; MoRFpred and ANCHOR. The results show that the proposed method outperforms these two predictors. CONCLUSIONS: Using HMM profile as a source of feature extraction, the proposed method indicates improvement in predicting MoRFs in disordered protein sequences.

MoRF_ESM: Prediction of MoRFs in disordered proteins based on a deep transformer protein language model.

Molecular recognition features (MoRFs) are particular functional segments of disordered proteins, which play crucial roles in regulating the phase transition of membrane-less organelles and frequently serve as central sites in cellular interaction networks. As the association between disordered proteins and severe diseases continues to be discovered, identifying MoRFs has gained growing significance. Due to the limited number of experimentally validated MoRFs, the performance of existing MoRF's prediction algorithms is not good enough and still needs to be improved. In this research, we present a model named MoRF_ESM, which utilizes deep-learning protein representations to predict MoRFs in disordered proteins. This approach employs a pretrained ESM-2 protein language model to generate embedding representations of residues in the form of attention map matrices. These representations are combined with a self-learned TextCNN model for feature extraction and prediction. In addition, an averaging step was incorporated at the end of the MoRF_ESM model to refine the output and generate final prediction results. In comparison to other impressive methods on benchmark datasets, the MoRF_ESM approach demonstrates state-of-the-art performance, achieving [Formula: see text] higher AUC than other methods when tested on TEST1 and achieving [Formula: see text] higher AUC than other methods when tested on TEST2. These results imply that the combination of ESM-2 and TextCNN can effectively extract deep evolutionary features related to protein structure and function, along with capturing shallow pattern features located in protein sequences, and is well qualified for the prediction task of MoRFs. Given that ESM-2 is a highly versatile protein language model, the methodology proposed in this study can be readily applied to other tasks involving the classification of protein sequences.

CoMemMoRFPred: Sequence-based Prediction of MemMoRFs by Combining Predictors of Intrinsic Disorder, MoRFs and Disordered Lipid-binding Regions.

Molecular recognition features (MoRFs) are a commonly occurring type of intrinsically disordered regions (IDRs) that undergo disorder-to-order transition upon binding to partner molecules. We focus on recently characterized and functionally important membrane-binding MoRFs (MemMoRFs). Motivated by the lack of computational tools that predict MemMoRFs, we use a dataset of experimentally annotated MemMoRFs to conceptualize, design, evaluate and release an accurate sequence-based predictor. We rely on state-of-the-art tools that predict residues that possess key characteristics of MemMoRFs, such as intrinsic disorder, disorder-to-order transition and lipid-binding. We identify and combine results from three tools that include flDPnn for the disorder prediction, DisoLipPred for the prediction of disordered lipid-binding regions, and MoRFCHiBiLight for the prediction of disorder-to-order transitioning protein binding regions. Our empirical analysis demonstrates that combining results produced by these three methods generates accurate predictions of MemMoRFs. We also show that use of a smoothing operator produces predictions that closely mimic the number and sizes of the native MemMoRF regions. The resulting CoMemMoRFPred method is available as an easy-to-use webserver at http://biomine.cs.vcu.edu/servers/CoMemMoRFPred. This tool will aid future studies of MemMoRFs in the context of exploring their abundance, cellular functions, and roles in pathologic phenomena.

"Order from disordered": Potential role of intrinsically disordered regions in phytopathogenic oomycete intracellular effector proteins.

There is a continuous arms race between pathogens and their host plants. However, successful pathogens, such as phytopathogenic oomycetes, secrete effector proteins to manipulate host defense responses for disease development. Structural analyses of these effector proteins reveal the existence of regions that fail to fold into three-dimensional structures, intrinsically disordered regions (IDRs). Because of their flexibility, these regions are involved in important biological functions of effector proteins, such as effector-host protein interactions that perturb host immune responses. Despite their significance, the role of IDRs in phytopathogenic oomycete effector-host protein interactions is not clear. This review, therefore, searched the literature for functionally characterized oomycete intracellular effectors with known host interactors. We further classify regions that mediate effector-host protein interactions into globular or disordered binding sites in these proteins. To fully appreciate the potential role of IDRs, five effector proteins encoding potential disordered binding sites were used as case studies. We also propose a pipeline that can be used to identify, classify as well as characterize potential binding regions in effector proteins. Understanding the role of IDRs in these effector proteins can aid in the development of new disease-control strategies.

Coronaviruses spike glycoprotein endodomains: The sequence and structure-based comprehensive study.

Any protein's flexibility or region makes it available to interact with many biomolecules in the cell. Specifically, such interactions in viruses help them to perform more functions despite having a smaller genome. Therefore, these flexible regions can be exciting and essential targets to be explored for their role in pathogenicity and therapeutic developments as they achieve essential interactions. In the continuation with our previous study on disordered analysis of SARS-CoV-2 spike protein's cytoplasmic tail (CTR), or endodomain, here we have explored the endodomain's disordered potential of six other coronaviruses using multiple bioinformatics approaches and molecular dynamics simulations. Based on the comprehensive analysis of its sequence and structural composition, we report the varying disorder propensity in endodomains of spike proteins of coronaviruses. The observations of this study may help to understand the importance of spike glycoprotein endodomain and creating therapeutic interventions against them.

Computational Methods to Predict Intrinsically Disordered Regions and Functional Regions in Them.

Intrinsically disordered regions (IDRs) are protein regions that do not adopt fixed tertiary structures. Since these regions lack ordered three-dimensional structures, they should be excluded from the target portions of homology modeling. IDRs can be predicted from the amino acid sequences, because their amino acid compositions are different from that of the structured domains. This chapter provides a review of the prediction methods of IDRs and a case study of IDR prediction.

Computational prediction of disordered binding regions.

One of the key features of intrinsically disordered regions (IDRs) is their ability to interact with a broad range of partner molecules. Multiple types of interacting IDRs were identified including molecular recognition fragments (MoRFs), short linear sequence motifs (SLiMs), and protein-, nucleic acids- and lipid-binding regions. Prediction of binding IDRs in protein sequences is gaining momentum in recent years. We survey 38 predictors of binding IDRs that target interactions with a diverse set of partners, such as peptides, proteins, RNA, DNA and lipids. We offer a historical perspective and highlight key events that fueled efforts to develop these methods. These tools rely on a diverse range of predictive architectures that include scoring functions, regular expressions, traditional and deep machine learning and meta-models. Recent efforts focus on the development of deep neural network-based architectures and extending coverage to RNA, DNA and lipid-binding IDRs. We analyze availability of these methods and show that providing implementations and webservers results in much higher rates of citations/use. We also make several recommendations to take advantage of modern deep network architectures, develop tools that bundle predictions of multiple and different types of binding IDRs, and work on algorithms that model structures of the resulting complexes.

MoRF-FUNCpred: Molecular Recognition Feature Function Prediction Based on Multi-Label Learning and Ensemble Learning.

Intrinsically disordered regions (IDRs) without stable structure are important for protein structures and functions. Some IDRs can be combined with molecular fragments to make itself completed the transition from disordered to ordered, which are called molecular recognition features (MoRFs). There are five main functions of MoRFs: molecular recognition assembler (MoR_assembler), molecular recognition chaperone (MoR_chaperone), molecular recognition display sites (MoR_display_sites), molecular recognition effector (MoR_effector), and molecular recognition scavenger (MoR_scavenger). Researches on functions of molecular recognition features are important for pharmaceutical and disease pathogenesis. However, the existing computational methods can only predict the MoRFs in proteins, failing to distinguish their different functions. In this paper, we treat MoRF function prediction as a multi-label learning task and solve it with the Binary Relevance (BR) strategy. Finally, we use Support Vector Machine (SVM), Logistic Regression (LR), Decision Tree (DT), and Random Forest (RF) as basic models to construct MoRF-FUNCpred through ensemble learning. Experimental results show that MoRF-FUNCpred performs well for MoRF function prediction. To the best knowledge of ours, MoRF-FUNCpred is the first predictor for predicting the functions of MoRFs. Availability and Implementation: The stand alone package of MoRF-FUNCpred can be accessed from https://github.com/LiangYu-Xidian/MoRF-FUNCpred.

Anchoring intrinsically disordered proteins to multiple targets: lessons from N-terminus of the p53 protein.

Anchor residues, which are deeply buried upon binding, play an important role in protein-protein interactions by providing recognition specificity and facilitating the binding kinetics. Up to now, studies on anchor residues have been focused mainly on ordered proteins. In this study, we investigated anchor residues in intrinsically disordered proteins (IDPs) which are flexible in the free state. We identified the anchor residues of the N-terminus of the p53 protein (Glu17-Asn29, abbreviated as p53N) which are involved in binding with two different targets (MDM2 and Taz2), and analyzed their side chain conformations in the unbound states. The anchor residues in the unbound p53N were found to frequently sample conformations similar to those observed in the bound complexes (i.e., Phe19, Trp23, and Leu26 in the p53N-MDM2 complex, and Leu22 in the p53N-Taz2 complex). We argue that the bound-like conformations of the anchor residues in the unbound state are important for controlling the specific interactions between IDPs and their targets. Further, we propose a mechanism to account for the binding promiscuity of IDPs in terms of anchor residues and molecular recognition features (MoRFs).

Prediction of MoRFs based on sequence properties and convolutional neural networks.

BACKGROUND: Intrinsically disordered proteins possess flexible 3-D structures, which makes them play an important role in a variety of biological functions. Molecular recognition features (MoRFs) act as an important type of functional regions, which are located within longer intrinsically disordered regions and undergo disorder-to-order transitions upon binding their interaction partners. RESULTS: We develop a method, MoRFCNN, to predict MoRFs based on sequence properties and convolutional neural networks (CNNs). The sequence properties contain structural and physicochemical properties which are used to describe the differences between MoRFs and non-MoRFs. Especially, to highlight the correlation between the target residue and adjacent residues, three windows are selected to preprocess the selected properties. After that, these calculated properties are combined into the feature matrix to predict MoRFs through the constructed CNN. Comparing with other existing methods, MoRFCNN obtains better performance. CONCLUSIONS: MoRFCNN is a new individual MoRFs prediction method which just uses protein sequence properties without evolutionary information. The simulation results show that MoRFCNN is effective and competitive.

MedProDB: A database of Mediator proteins.

In the last three decades, the multi-subunit Mediator complex has emerged as the key component of transcriptional regulation of eukaryotic gene expression. Although there were initial hiccups, recent advancements in bioinformatics tools contributed significantly to in-silico prediction and characterization of Mediator subunits from several organisms belonging to different eukaryotic kingdoms. In this study, we have developed the first database of Mediator proteins named MedProDB with 33,971 Mediator protein entries. Out of those, 12531, 11545, and 9895 sequences belong to metazoans, plants, and fungi, respectively. Apart from the core information consisting of sequence, length, position, organism, molecular weight, and taxonomic lineage, additional information of each Mediator sequence like aromaticity, hydropathy, instability index, isoelectric point, functions, interactions, repeat regions, diseases, sequence alignment to Mediator subunit family, Intrinsically Disordered Regions (IDRs), Post-translation modifications (PTMs), and Molecular Recognition Features (MoRFs) may be of high utility to the users. Furthermore, different types of search and browse options with four different tools namely BLAST, Smith-Waterman Align, IUPred, and MoRF-Chibi_Light are provided at MedProDB to perform different types of analysis. Being a critical component of the transcriptional machinery and regulating almost all the aspects of transcription, it generated lots of interest in structural and functional studies of Mediator functioning. So, we think that the MedProDB database will be very useful for researchers studying the process of transcription. This database is freely available at www.nipgr.ac.in/MedProDB.