Skip to main content
  • Methodology article
  • Open access
  • Published:

MFSPSSMpred: identifying short disorder-to-order binding regions in disordered proteins based on contextual local evolutionary conservation

Abstract

Background

Molecular recognition features (MoRFs) are short binding regions located in longer intrinsically disordered protein regions. Although these short regions lack a stable structure in the natural state, they readily undergo disorder-to-order transitions upon binding to their partner molecules. MoRFs play critical roles in the molecular interaction network of a cell, and are associated with many human genetic diseases. Therefore, identification of MoRFs is an important step in understanding functional aspects of these proteins and in finding applications in drug design.

Results

Here, we propose a novel method for identifying MoRFs, named as MFSPSSMpred (Masked, Filtered and Smoothed Position-Specific Scoring Matrix-based Predictor). Firstly, a masking method is used to calculate the average local conservation scores of residues within a masking-window length in the position-specific scoring matrix (PSSM). Then, the scores below the average are filtered out. Finally, a smoothing method is used to incorporate the features of flanking regions for each residue to prepare the feature sets for prediction. Our method employs no predicted results from other classifiers as input, i.e., all features used in this method are extracted from the PSSM of sequence only. Experimental results show that, comparing with other methods tested on the same datasets, our method achieves the best performance: achieving 0.004~0.079 higher AUC than other methods when tested on TEST419, and achieving 0.045~0.212 higher AUC than other methods when tested on TEST2012. In addition, when tested on an independent membrane proteins-related dataset, MFSPSSMpred significantly outperformed the existing predictor MoRFpred.

Conclusions

This study suggests that: 1) amino acid composition and physicochemical properties in the flanking regions of MoRFs are very different from those in the general non-MoRF regions; 2) MoRFs contain both highly conserved residues and highly variable residues and, on the whole, are highly locally conserved; and 3) combining contextual information with local conservation information of residues facilitates the prediction of MoRFs.

Background

With the breaking of the conventional protein paradigm of “sequence-structure-function”, the functional importance of intrinsically disordered proteins (IDPs) has become increasingly apparent. Although IDPs have no well-defined tertiary structures in their natural state, they possess essential biological functions. IDPs play critical roles in a variety of physiological processes such as signal transduction, translation regulation, and protein modification. Specifically, in interaction-mediated signaling events, IDPs possess unique advantages [1]: (a) they lack a stable three-dimensional structure and their conformations can fluctuate over time to time; (b) they possess a combination of high specificity and low affinity; (c) they can recognize multiple partners through adoption of different conformations; (d) multiple distinct partners can bind to a common binding site of IDPs, where these partners may assume different folds. Many proteins from higher organism have been found to be entirely disordered or to contain partly disordered regions [2]. Due to the particular properties of IDPs, conventional structure determination methods are simply inapplicable to IDPs in isolation, because no stable structure exists. However, when bound to their molecular partners, many IDPs undergo a disorder-to-order transition [1-3]. This characteristic makes it possible to obtain the structures of IDPs by crystallizing them in complexes with their molecular binding partners.

Molecular recognition features (MoRFs) are short binding regions (5-25 residues) located in IDPs regions, which easily undergo disorder-to-order transitions upon binding to partner proteins [4]. According to their structures in the bound state, MoRFs can be divided into at least three sub-types: α-MoRFs, β-MoRFs, and ι-MoRFs, which form α-helices, β-strands, and structures without a regular pattern of backbone hydrogen bonds [2]. MoRFs are enrichment in highly connected hub proteins, and their complexity reinforces the functional importance of the disordered regions [3]. They act as molecular switches in molecular-interaction networks of the cell, and are assumed to be implicated in the causes of many diseases [1]. Thus, identification of MoRFs is a key step in understanding the functions of these proteins and in finding applications in drug design.

Experimental methods for identifying MoRFs are expensive and time consuming, which makes computational methods indispensable for guiding experimental analysis. So far, because of the limited number of experimentally validated MoRFs, only four custom-built tools for predicting MoRFs are available: α-MoRF-PredI [5] and α-MoRF-PredII [6] are neural network-based predictors aimed at predicting α-MoRFs; ANCHOR [7] concentrates on the prediction of MoRFs which bind to globular proteins; and MoRFpred [4] is a comprehensive predictor which combines the annotations generated by sequence alignments with the prediction results generated by a support vector machine (SVM). These predictors use a variety of predicted features as their input, including predicted disorder probabilities [4-7], predicted solvent accessibility [4], predicted secondary structure propensities [4, 5], and predicted B-factors [4]. These predicted features not only easy to result in high-dimensional feature space, but also greatly increase the complexity of algorithms. Moreover, the performance of these predictors is also largely affected by other classifiers. Thus, a more simple and efficient method for identifying MoRFs is indispensable.

A number of studies have analyzed the attributes of MoRFs [1-9]. Norman et al. [3] found that several strong physicochemical preferences were shown in all MoRF types compared with in general disordered regions. Fuxreiter et al. [8] identified that amino acid composition and charge/hydropathy properties of MoRFs exhibited a mixture characteristic of folded and disordered structures. Chica et al. [9] found that the flanking regions of MoRFs were relevant to linear motif-mediated interactions at both the structural and sequence levels, and that the prediction of MoRFs can be facilitated by contextual information from the protein sequence.

Because the functional sites of proteins need to maintain a high degree of conservation to execute a given function, evolutionary information included in a position-specific scoring matrix (PSSM) has been considered as the most predictive feature for identifying the functional sites of ordered proteins. MoRFs are also found to be more conserved than their surrounding residues [9, 10]. However, disordered proteins usually evolve more rapidly than ordered proteins; therefore, standard PSSMs which incorporate the conservation information of proteins are ineffective when used directly. Fortunately, relative local conservation has been proven to be a good feature for motif discovery, and has been used in a number of studies [10-12]. In addition, Shimizu et al. [13] found that ignoring some redundant features in standard PSSMs can improve the prediction of protein disordered regions significantly.

Taking the above information into account, we developed a novel sequence-based method for identifying MoRFs in IDPs. In this method, firstly, a masking method is used to calculate the average local conservative scores. Then, a filtering method is used to drop the scores below the average. Finally, a smoothing method is used to incorporate features of the flanking regions for each residue in the PSSM. The masking and filtering steps strengthen the highly conserved information and filter out poorly conserved information for each residue, thereby ensuring that only highly local conserved features are considered in the prediction. Moreover, the smoothing method can incorporate contextual information of neighboring residues for any given residue. All the related features are extracted only from protein sequences themselves. We used a support vector machine (SVM) to build the classifier.

Methods

Benchmark datasets

In the present study, one training dataset and three test datasets were adopted. We employed the same training dataset as the one that was used by Fatemeh et al. [4] in their study of MoRFpred. This dataset, called as TRAINING421, includes 421 MoRF-containing chains, that contained 5,396 positive samples (MoRF residues) and 240,588 negative samples (non-MoRF residues). All the positive samples and an equal number of randomly selected negative samples were used for training our prediction model. After that, three test datasets were used for testing the developed model.

First, the two test dataset used in the study of MoRFpred [4] were adopted.

TEST419

It included 419 MoRF-containing chains, which were deposited in PDB from April 2008, named as TEST419. They shared up to 30% sequence identity with the training dataset.

TEST2012

It included 45 MoRF-containing chains, which were deposited in PDB from January 1 to March 11, 2012 and in UniProtKB from February 22, 2012, named as TEST2012. They shared up to 30% sequence identity with the training dataset.

The training and the above two test datasets are available on the MoRFpred web server [http://biomine-ws.ece.ualberta.ca/MoRFpred/index.html] [4].

Because the TRAINING421 and TEST419 datasets were reported to contain a large number of immune response-related MoRFs (120 among the 840 MoRFs), we built another independent dataset to test whether the developed predictor was biased for some particular type of MoRFs.

TESTMem64

It included 64 non-redundant MoRF-containing membrane proteins (50 transmembrane proteins and 14 peripheral membrane proteins), which were extracted from the study of membrane proteins reported by Ioly et al. [14]. In their study, Ioly et al. collected 166 non-redundant MoRF-containing membrane protein sequences. We removed the sequences that either with sequence length >1000 residues, or with MoRF length >25 residues. Because sequences longer than 1000 residues cannot be processed by MoRFpred [4], in our study, we also aimed to predict short MoRFs (5-25 residues). After that, 64 sequences were retained, named as TESTMem64. Details of the TESTMem64 dataset are shown in Additional file 1: Table S1.

Feature analysis of MoRFs, their flanking regions, and general non-MoRF regions

Composition analysis

The composition of the 421 protein sequences related to TRAINING421dataset was analyzed. The sequences were characterized into three regions: MoRF region, regions flanking the MoRFs, and general non-MoRF regions.

When the length of the flanking region was 5-residue long, amino acid composition of the three regions is shown in Figure 1. The composition of the MoRFs and their flanking regions deviated significantly from that of the general non- MoRF regions. The MoRFs were enriched in hydrophobic residues Ile, Leu, Phe, Tyr and Lys, as well as in charged residues Arg, and Asp, while depleted in Ala and Gly. The flanking regions were enriched in variable (disorder-promoting) amino acids, such as Ala, Gly, Glu, Ser and Thr, and depleted in hydrophobic and rigid (order-promoting) residues. Both of MoRFs and the flanking regions exhibit a large excess of a small, polar residue--Serine. The general non-MoRF regions are overrepresented with Val, Cys, Ala, Gly, Asn, and Ser, most of which were small and tiny amino acids.

Figure 1
figure 1

Distribution of the amino acids in three regions: MoRF regions (blue), 5-residue long flanking regions (green), and non-MoRF regions (red).

We further analyzed the composition of flanking regions with different lengths. While the composition distribution of general non-MoRF regions was essentially unchanged over different lengths (see the Additional file 2: Figure S1), that of the flanking regions changed significantly (see the Additional file 3: Figure S2). The composition differences between the flanking regions and the general non-MoRF regions (% amino acids in flanking regions - % amino acids in general non-MoRF regions) are shown in Figure 2. We also calculated the Pearson product-moment correlation coefficients between the flanking length and the composition difference for each amino acid (Additional file 4: Table S2). The correlation coefficients related to Ile, Leu, Val, His, Arg, Glu and flanking regions Ser have absolute value > 0.75 at the p-value <0.001. It illustrates that, composition of these amino acids in the flanking regions and in the general non-MoRF regions are significantly different.

Figure 2
figure 2

Differences in amino acid composition between MoRF flanking and non-MoRF regions. W indicates the length of the flanking regions; for example, W5 indicates a 5-residue long MoRF flanking region.

Physicochemical properties

We next analyzed differences between flanking and general non-MoRF regions with respect to ten physicochemical properties, namely: hydrophobic, polar, small, proline, tiny, aliphatic, aromatic, positive, negative and charged. The difference of physicochemical propensity between the flanking regions and general non-MoRF regions (property percentage in flanking regions - property percentage in general non-MoRF regions) is shown in Figure 3. The corresponding Pearson correlation coefficients between the flanking length and each physicochemical properties difference are shown in Additional file 5: Table S3. Correlation coefficients related to hydrophobic, polar, aliphatic, aromatic, positive, negative and charged have absolute value > 0.75 at the p-value <0.001. This phenomenon illustrates that, physicochemical properties of flanking regions are very different from those of general non-MoRF regions, namely that, the arrangements and properties of neighboring residues significantly impact the plasticity of MoRFs [8].

Figure 3
figure 3

Differences in physicochemical properties between MoRF flanking and non-MoRF regions. W indicates the length of the flanking regions; for example, W5 indicates a 5-residue long MoRF flanking region.

The analyses of the composition and physicochemical properties for the three regions illustrated that flanking regions are highly relevant to the MoRFs. Therefore, we assumed that the MoRFs in protein sequences are highly contextual. Because MoRFs are found to be more conserved than surrounding residues [9, 10], we considered incorporating contextual information of residues with local evolutionary conservation to improve the prediction of MoRFs.

Evolutionary information from PSSM

Evolutionary information was obtained from PSSMs, which were generated by PSI-BLAST [15], searching against the NCBI non-redundant (nr) database [16] by three-time iteration with an e-value of 0.001. Evolutionary information for each amino acid was encapsulated in a vector of 20 dimensions, where the size of PSSM matrix of a protein with N residues is 20 × N. 20 dimensions were considered as a standard amino acid size, and N is the length of the protein sequence.

Masking and filtering the PSSM

Masked PSSM was used to describe the relative local evolutionary information of each residue in a protein. It was converted from a standard PSSM according to Formula (1).

Masking _ C i = C i - 1 2 n + 1 i - n i + n C j , i = 1 , , N , where N is the sequence length ,
(1)

where Masking_C i represents the relative local conservation score of residue i, Ci is the standard conservation score in the PSSM, and 2n+1 is the masking-window size.

In order to extract highly locally conserved information and filter out poorly locally conserved information from a sequence, the masked PSSM was further converted to a filtered PSSM according to Formula (2). All the scores below the mean score within a masking-window length were set to 0.

Filtering _ C i = Masking _ C i , Masking _ C i > 0 0 , Masking _ C i 0 , i = 1 , , N , where N is the sequence length
(2)

An example of a masked and filtered PSSM profile is shown in Figure 4. For amino acid 'A', the first column of the vector was masked by 4 - [(-1) + (-1) + 4 + (-1) + (-2)]/5 = 4.2, because 4.2 > 0, it was retained, otherwise, it would be set to 0. After modification, only the scores that represented strong local conservation were retained.

Figure 4
figure 4

Procedure of preparing feature sets for the predictor.

Smoothing the modified PSSM

Every value in a standard PSSM is calculated based on the assumption that the position of each value in the matrix is independent of the others. However, residues in MoRF regions appear continuously and the plasticity of MoRFs is also largely affected by neighboring residues [8]. Accordingly, in order to incorporate the dependency on surrounding neighbors of a central residue, we adopted a smoothing method [17] that was inspired by the smoothing skill used in image processing [18]. Firstly, to address the N-terminal and C-terminal of a protein sequence, m (m is an odd number) ZERO vectors were appended to the head and tail of a standard PSSM profile, where 2m+1 was the size of a smoothing sliding-window. The smoothing sliding-window was then used to incorporate the evolutionary information from the upstream and downstream residues. Each row vector of an amino acid residue C i was smoothed according to Formula (3). Finally, each value in the masked and smoothed PSSM matrix was scaled to a range of [-1, 1] according to a certain ratio.

Smoothing _ C i = i - m i + m Filtering _ C j , i = 1 , , N , N is the sequence length , 2 m + 1 is smoothing - window size
(3)

Prediction model

We adopt a modified PSSM, which incorporates the information of amino acid position, relative evolutionary information, and dependency on neighboring residues, to design our prediction model. Only highly conserved position scores calculated from PSSMs were considered in our prediction. Firstly, a masking method was used to calculate the average score of the local conservation within a masking-window length. Next, scores lower than the averages were filtered out. This procedure could strengthen the highly locally conserved information and filter out poorly locally conserved information of residues. Finally, a smoothing method was adopted to incorporate feature information from the neighbors flanking each residue. The prediction model is shown in Figure 5, and the detailed procedure for preparing the input feature set is shown in Figure 4.

Figure 5
figure 5

Prediction model design. 2h+1 is the length of the outside sliding-window, n is the sequence length, and v i represents the corresponding amino acids i in the feature vector.

SVM

Prediction of MoRFs can be addressed as a two-classification problem, namely, determining whether a given residue belongs to a MoRFs or not. Our prediction model was trained using the LIBSVM software package [19, 20]. Here, the Radial Basis Function (RBF kernel) was selected as the kernel function. The capacity parameter c and kernel width parameter g were then optimized using a grid search approach [19, 20].

Evaluation criteria

We adopted the evaluation criteria used in the 10th Critical Assessment of Techniques for Protein Structure Prediction (CASP10) experiment [21]. The area under ROC curve (AUC) and the accuracy (ACC) were adopted to evaluate the performance of the classifiers. Receiver operating characteristic (ROC) plots with AUC values were generated using the R statistical package [22]. The sensitivity, specificity, true positive rate (TPR), false negative rate (FPR), success rate, and ACC were defined as follows:

Specifity = TN TN + FP
(4)
TPR = Sensitivity = TP TP + FN
(5)
FPR = 1 - Specificity = FP TN + FP
(6)
Accuracy = TP + TN TP + TN + FP + FN
(7)
ACC = 1 2 Sensitivity + Specificity
(8)

Where TP, TN, FP and FN represents true positive, true negative, false positive and false negative respectively.

Results and discussion

Optimizing window size

To develop the MFSPSSMpred model, three window sizes were necessary: (i) the outside sliding-window size which ultimately determined the dimensions of feature vectors; (ii) the masking-window size which was used to calculate the average conservation scores in a local region; and (iii) the inside smoothing-window size which was used to strengthen locally conserved features. To make a fair comparison among the different predictors, we chose the same outside-sliding window size of '25’ with the research MoRFpred [4]. The masking-window size was also assigned to '25’, since it has the similar meaning with the outside sliding-window size, both of which indicated the flanking length that would be considered to affect a central residue. Here, because the average scores calculated from different lengths greater than 10 residues were highly similar, we found that variations in masking-window size had limited influence on the results (Results of MFSPSSMpred with different masking-window sizes are shown in Additional file 6: Figure S3). Next, MFSPSSMpred models with different smoothing-window sizes were tested by 5-cross-validation using the grid search approach [19]. The cross-validation accuracies according to different smoothing-window sizes are shown in Figure 6. The predictor performance stabilized with the smoothing-window sizes greater than 9. Therefore, we chose the relatively best size 13 as the inside smoothing-window size for our model.

Figure 6
figure 6

Cross-validation accuracy of MFSPSSMpred with different smoothing-window sizes. W indicates the length of windows, i.e., W3 means 3-residue long.

Effectiveness of the feature extracting methods

Effectiveness at the individual protein level

To confirm the effectiveness of our feature extracting method in distinguishing MoRF residues from non_MoRF residues and to determine how this might benefit the predication, we selected the yeast elongin C complex with a von Hippel-Lindau peptide [PDB: 1HV2] as an example. Yeast elongin C is a signaling protein, which is highly intrinsically disordered and has a MoRF located at the region of residues 100-212 of its B chain.

We extracted a column of position-specific scores randomly from the standard PSSM (there are a total of 20 columns, corresponding to the 20 standard amino acids). Then, the distribution of the scores in the standard PSSM, masked PSSM, and smoothed PSSM, was counted (Figure 7), in the standard PSSM, the distribution of scores between the MoRF and non_MoRF showed no distinct difference. The MoRFs regions included both highly conserved residues and highly variable residues, which is consistent with the property of MoRFs. i.e., to maintain a certain function, MoRFs must be highly conserved. On the other hand, MoRFs must be highly variable because of their intrinsically disordered nature. In the masked PSSM, scores below the average have been filtered out, thereby discarding the noise data (low conservative scores), that are undesirable for prediction. In the smoothed PSSM, scores of residues that are surrounded by highly conserved residues have been enhanced after smoothing, and scores of residues that are surrounded by poorly conserved residues have been weakened. Since MoRFs are usually composed of intensively conserved residues, they are easy to appear as highly conserved peaks compared with the non_MoRF regions. The results demonstrate that, after masking and smoothing the PSSM, the position-specific scores in MoRF regions become significantly higher than those in the non_MoRF regions. Here, we have shown the results for just one column of scores in the PSSM, similar results were obtained for the scores in the other 19 columns of the PSSM (data not shown).

Figure 7
figure 7

Distribution of the conservation scores for protein 1HV2 (PDB ID) related to the standard PSSM (a), masked PSSM (b), masked and smoothed PSSM (c). The red-dotted rectangles mark the position of MoRF (only residues 100~212 in the B chain are displayed).

Effectiveness at the whole dataset level

We also calculated the distributions of local conservation scores in each of the three PSSMs for all residues in the 840 chains (TRAINING421 and TEST419). We assumed that the summation of scores in each row of the corresponding PSSMs represented the local conservation of the related residues. The local conservation scores in each PSSM were calculated according to Formulas (9), (10) and (11).

Sum _ Original _ C i = j = 1 20 Original _ C j , i = 1 , , N , N is the sequence length
(9)
Sum _ Masked _ C i = j = 1 20 Masked _ C j , i = 1 , , N , N is the sequence length
(10)
Sum _ Smoothed _ C i = j = 1 20 Smoothed _ C j , i = 1 , , N , N is the sequence length
(11)

where '20’ is the number of standard amino acid residues, and also the number of columns in the PSSMs (PSSM is a 20 × N matrix), and N is the sequence length.

We compared the distribution of the amino acid residues in the MoRF and non_MoRF regions in each of the three PSSMs. In the standard PSSM, the MoRF residues are at a disadvantage in conservation compared with the non_MoRF residues (Figure 8a), because MoRFs contain highly variable residues. After masking and filtering, the differences between the residues in the MoRF and non_MoRF regions became less obvious, and some MoRF residues exhibited more strongly conserved than the non_MoRF residues (Figure 8b), showing that most of the filtered scores belonged to the non_MoRF residues. Finally, after smoothing, the residues in the two regions had distinctly different distributions; the MoRF residues were much more strongly conserved than most of non_MoRF residues (Figure 8c). The findings obtained using the whole of the TRAINING421 and TEST419 datasets show that our feature extracting method can effectively distinguish MoRF residues from non_MoRF residues as illustrated in Figure 8.

Figure 8
figure 8

Distribution of the summation scores of residues in MoRF and non-MoRF regions related to the standard PSSM (a), masked PSSM (b), masked and smoothed PSSM(c). Here the non-MoRF regions include the residues in the flanking regions.

Performance tested on TEST419, TEST2012 and comparison with other feature-based methods

First, TRAINING421 was used to train MSPSSMpred, which was then tested on the TEST419 and TEST2012 datasets. ROC plots of the results are shown in Figure 9(a ~b) respectively. For TEST419, the AUC was 0.677 and, for TEST2012, the AUC was 0.724.

Figure 9
figure 9

ROC plots of MFSPSSMpred tested on TEST419 (a) and TEST2012 (b).

Originally, the direct outputs of PSSMs from PSI-BLAST have provided conversation information by default and have been widely used to predict various protein functional sites. However, there is a room for improvement, because standard PSSMs contain redundant features.

Here, we have compared our method with four other PSSM_based methods: 1) the 'PSSM’ method, which uses the standard PSSM (the direct output of PSSMs) for prediction; 2) the 'Smooth_PSSM’ method, which uses smoothed PSSMs without masking and filtering; 3) the 'Mask_PSSM’ method, which uses masked and filtered PSSMs without smoothing; and 4) the 'MFS_Physi_PSSM’ method, which is similar to our MFSPSSMpred method but incorporates 10 physicochemical properties of residues as input. The performances of MFSPSSMpred and the other four methods based on TRAINING421 and TEST419 are shown in Table 1, and ROC plots for all the methods are shown in Figure 10. The results demonstrate that MFSPSSMpred achieves the best performance.

Table 1 Performance comparison with four other PSSM-based methods
Figure 10
figure 10

ROC plots for MFSPSSMpred and four other PSSM-based methods.

Performance comparison with existing predictors

Some existing tools that are publicly available for MoRFs prediction have been tested on the TEST419 and TEST2012 datasets [4]. Here, we list them out for a comparison. Results are shown in Table 2, and the results of other classifiers are quoted from the study of Fatemeh et al. [4]. Details of ROC and ACC for all the predictors are shown in Figure 11. The results demonstrate that, MFSPSSMpred outperformed the other predictors with respect to ACC and AUC on both the TEST419 and TEST2012 datasets.

Table 2 Performance comparisons tested on the TEST419 and TEST2012 datasets
Figure 11
figure 11

ROC and ACC for all the predictors tested on TEST419 (a) and TEST2012 (b). The AUC cannot be computed for α-MoRF-predI and α-MoRF-predII because they generate only binary predictions.

Performance on unbalanced training samples

The TRAINING421 dataset contains 5,601positive samples and 262,732 negative samples, and the ratio between them is 1:46.9. In order to analyze whether this imbalance biased the prediction method, we developed another training model with a 1:2 ratio between the MoRF and non-MoRF residues, that is, 5,601 MoRFs with 112,02 non-MoRF residues, and tested it on the TEST419 and TEST2012 datasets (Figure 12(a ~b)). Our results demonstrate that there was no significant difference in the performance between the 2:1 and 1:1 ratios.

Figure 12
figure 12

Performance of MFSPSSMpred on different ratios of training samples, tested on TEST419 (a) and TEST2012 (b). Red plots represent the result based on a training dataset with 1:2 ratio of MoRFs to non_MoRFs; the black plots represent the result based on a training dataset with 1:1 ratio of MoRFs to non_MoRFs.

Performance tested on TESTMem64 and comparison with MoRFpred

The TRAINING421 and TEST419 datasets contain a large proportion (120 among the 840 MoRFs) of immune response-related MoRFs [4]. In order to test whether MFSPSSMpred was biased for some particular type of MoRFs, we built an independent test dataset --TESTMem64, which was extracted from another independent study of membrane proteins by Ioly et al. [14]. MoRFpred [4] was also tested on TESTMem64 for a comparison. Results are shown in Figure 13 and a detailed comparison of ACC, TPR, FPR and AUC is shown in Table 3. MFSPSSMpred performed much better than MoRFpred, achieving significantly higher ACC and AUC than MoRFpred.

Figure 13
figure 13

ROC plots of MFSPSSMpred (black) and MoRFpred (red) tested on TESTMem64.

Table 3 Performance comparison between MFSPSSMPred and MoRFpred tested on TESTMem64

We speculate that the reasons for the better performance of MFSPSSMpred include: (1) the MoRFpred method incorporated many predicted results, such as predicted disorder probabilities, predicted B-factor and predicted relative solvent accessibility derived from other predictors, as input for the prediction. These predicted features themselves are largely affected by the other classifiers those were used. Moreover, incorporating many predicted features can easy to result in a high-dimensional feature space; (2) MoRFpred [4] merges the result generated by an SVM and the result generated by sequence alignment with the MoRFs database into their final prediction result. Since there are so many immune response-related MoRFs in their database, MoRFpred is inevitably biased towards this type of MoRFs; (3) MFSPSSMpred used only the PSSM as input for prediction. It caught the point that, MoRF regions in a sequence are mingled with highly conserved residues and highly variable residues. Therefore, our approach is independent of the type or binding partners of MoRFs, and can be applied to the prediction of all MoRFs.

Conclusions

In this study, we propose a novel method which adopts a modified PSSM encoding scheme for MoRFs prediction. Our method employs no predicted results as input, and all input features are extracted only from the PSSMs of sequences. By means of masking, filtering and smoothing, the modified PSSMs combine predictive features, which can effectively distinguish MoRF from non_MoRF residues. When comparing with other existing methods on the same datasets, MFSPSSMpred outperformed them all, achieving 3.3%~9.3% higher ACC and 0.004~0.079 higher AUC than other methods when tested on TEST419, and 10.3%~17.2% higher ACC and 0.045~0.212 higher AUC than other methods when tested on TEST2012. Moreover, despite the training dataset being biased toward immune response-related proteins, when MFSPSSMpred was tested on an independent membrane proteins-related dataset---TESTMem64, it showed good adaptability and significantly outperformed the existing MoRFpred predictor [4].

In summary, this study shows that combining contextual information with local conservation information of residues is predictive for identifying MoRFs. In addition, our study revealed some hallmarks of MoRFs; namely, MoRFs are mingled with highly conserved residues and highly variable residues, and MoRFs, on the whole, are highly locally conserved and are flanked by less conserved residues. A free Web server has been developed, which allows users to identify MoRFs in a given sequence using the model trained on our dataset. It is available from URL http://webapp.yama.info.waseda.ac.jp/fang/MoRFs.php.

Abbreviations

MoRFs:

Molecular recognition features

PSSM:

Position-specific scoring matrix

SVM:

Support vector machine

IDPs:

Intrinsically disordered proteins

ROC:

Receiver operating characteristic

AUC:

Area under the corresponding ROC curve

TPR:

True positive rate

FPR:

False negative rate.

References

  1. Mohan A, Oldfield CJ, Radivojac P, Vacic V, Cortese MS, Dunker AK, Uversky VN: Analysis of molecular recognition features (MoRFs). J Mol Biol. 2006, 362: 1043-1059. 10.1016/j.jmb.2006.07.087.

    Article  CAS  PubMed  Google Scholar 

  2. Vacic V, Christopher JO, Amrita M, Predrag R, Marc SC, Vladimir NU, Dunker AK: Characterization of molecular recognition features, MoRFs, and their binding partners. J Proteome Res. 2007, 6 (6): 2351-2366. 10.1021/pr0701411.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Norman ED, Kim VR, Robert JW: Attributes of short linear motifs. Mol Biosyst. 2012, 8: 268-281. 10.1039/c1mb05231d.

    Article  Google Scholar 

  4. Fatemeh MD, Wei-Lun H, Marcin JM, Christopher JO, Bin X, Dunker AK, Vladimir NU, Lukasz K: MoRFpred, a computational tool for sequence-based prediction and characterization of short disorder-to-order transitioning binding regions in proteins. Bioinformatics. 2012, 28 (12): i75-i83. 10.1093/bioinformatics/bts209.

    Article  Google Scholar 

  5. Oldfield CJ, Cheng Y, Cortese MS, Romero P, Uversky VN, Dunker AK: Coupled folding and binding with alpha-helix-forming molecular recognition elements. Biochemistry. 2005, 44: 12454-12470. 10.1021/bi050736e.

    Article  CAS  PubMed  Google Scholar 

  6. Yugong C, Christopher JO, Jingwei M, Pedro R, Vladimir NU, Dunker AK: Mining α-helix-forming molecular recognition features with cross species sequence alignments. Biochemistry. 2007, 46 (47): 13468-13477. 10.1021/bi7012273.

    Article  Google Scholar 

  7. Dosztanyi Z, Mészáros SI: ANCHOR: web server for predicting protein binding regions in disordered proteins. Bioinformatics. 2009, 25 (20): 2745-2746. 10.1093/bioinformatics/btp518.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  8. Fuxreiter M, Peter T, Istvan S: Local structural disorder imparts plasticity on linear motifs. Bioinformatics. 2007, 23 (8): 950-956. 10.1093/bioinformatics/btm035.

    Article  CAS  PubMed  Google Scholar 

  9. Chica C, Diella F, Gibson TJ: Evidence for the concerted evolution between short linear protein motifs and their flanking regions. PLoS ONE. 2009, 4 (7): e6052-10.1371/journal.pone.0006052.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Norman ED, Denis CS, Richard JE: Masking residues using context-specific evolutionary conservation significantly improves short linear motif discovery. Bioinformatics. 2009, 25 (4): 443-450. 10.1093/bioinformatics/btn664.

    Article  Google Scholar 

  11. Norman ED, Joanne LC, Denis CS, Toby JG, Mark JC, Richard JE: SLiMPrints: conservation-based discovery of functional motif fingerprints in intrinsically disordered protein regions. Nucleic Acids Res. 2012, September, 12: 1-14.

    Google Scholar 

  12. Niall JH, Denis CS: Profile-based short linear protein motif discovery. BMC Bioinformatics. 2012, 13: 104-10.1186/1471-2105-13-104.

    Article  Google Scholar 

  13. Shimizu K, Hirose S, Tamotsu N: POODLE-S: web application for predicting protein disorder by using physicochemical features and reduced amino acid set of a position-specific scoring matrix. Bioinformatics. 2007, 23 (17): 2337-2338. 10.1093/bioinformatics/btm330.

    Article  CAS  PubMed  Google Scholar 

  14. Ioly KL, Georgios NT, Stavros JH: Analysis of molecular recognition features (MoRFs) in membrane proteins. Biochimica et Biophysica Acta (BBA) - Proteins and Proteomics. 2013, 1834 (4): 798-807. 10.1016/j.bbapap.2013.01.006.

    Article  Google Scholar 

  15. Stephen FA, Thomas LM, Alejandro AS, Jinghui Z, Zheng Z: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25 (17): 3389-3402. 10.1093/nar/25.17.3389.

    Article  Google Scholar 

  16. NR. ftp://ftp.ncbi.nih.gov/blast/db/fasta/nr.gz

  17. Cheng-Wei C, Emily CYS, Jenn-Kang H, Ting-Yi S, Wen-Lian H: Predicting RNA-binding sites of proteins using support vector machines and evolutionary information. BMC Bioinformatics. 2008, 9 (Suppl 12): S6-10.1186/1471-2105-9-S12-S6.

    Article  Google Scholar 

  18. Gonzalez RC, Woods RE: Digital Image Processing. The Second Edition. 2002, Prentice Hall

    Google Scholar 

  19. Chih-Chung C, Chih-Jen L: LIBSVM: a library for support vector machines. ACM Transactions on Intelligent Systems and Technology. 2011, 2 (3): 27-

    Google Scholar 

  20. A library for support vector machines. [http://www.csie.ntu.edu.tw/~cjlin/libsvm/]

  21. CASP10. [http://predictioncenter.org/casp10/index.cgi]

  22. R statistical package. [http://www.r-project.org/]

  23. Avner S, Marco P, Guy Y, Laszlo K, Burkhard R: Improved disorder prediction by combination of orthogonal approaches. PLoS One. 2009, 4 (2): e4433-10.1371/journal.pone.0004433.

    Article  Google Scholar 

  24. Zsuzsanna D, Veronika C, Peter T, István S: IUPred: web server for the prediction of intrinsically unstructured regions of proteins based on estimated energy content. Bioinformatics. 2005, 21 (16): 3433-3434. 10.1093/bioinformatics/bti541.

    Article  Google Scholar 

  25. Marcin JM, Wojciech S, Ke C, Kanaka DK, Fatemeh MD, Lukasz K: Improved sequence-based prediction of disordered regions with multilayer fusion of multiple information sources. Bioinformatics. 2010, 26 (18): i489-i496. 10.1093/bioinformatics/btq373.

    Article  Google Scholar 

  26. Tuo Z, Eshel F, Bin X, Dunker AK, Vladimir NU, Yaoqi Z: SPINE-D: accurate prediction of short and long disordered regions by a single neural-network based method. Journal of Biomolecular Structure and Dynamics. 2012, 4 (29): 799-813.

    Google Scholar 

  27. Jonathan JW, Liam JM, Kevin B, Bernard FB, David TJ: The DISOPRED server for the prediction of protein disorder. Bioinformatics. 2004, 20 (13): 2138-2139. 10.1093/bioinformatics/bth195.

    Article  Google Scholar 

  28. McGuffin L: Intrinsic disorder prediction from the analysis of multiple protein fold recognition models. Bioinformatics. 2008, 24: 1798-1804. 10.1093/bioinformatics/btn326.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgments

We would like to thank Prof. Toh of Computational Biological Research Center (CBRC) providing a lot of convenience for our study, and Fatemeh’s laboratory providing their prediction results of the MoRFpred predictor. We also thank an anonymous reviewer for his/her helpful comments, which improved the manuscript.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Chun Fang.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

C.F carried out the implementation and drafted the manuscript, T.N and H.Y read and revised the final manuscript, all authors read and approved the final manuscript.

Electronic supplementary material

12859_2013_6120_MOESM1_ESM.pdf

Additional file 1: Table S1: Information about the TESTMem64 dataset. Information of binding partners is included. (PDF 128 KB)

12859_2013_6120_MOESM2_ESM.tif

Additional file 2: Figure S1: Composition distribution of general non-MoRF regions with different flanking-window lengths (excluding flanking regions, W indicates the flanking window length). (TIF 73 KB)

12859_2013_6120_MOESM3_ESM.tif

Additional file 3: Figure S2: Composition distribution of the flanking regions. W indicates the flanking window length. (TIF 81 KB)

12859_2013_6120_MOESM4_ESM.tif

Additional file 4: Table S2: Correlation coefficients between flanking length and composition difference (Flanking regions - general non-MoRF regions) for each amino acid. (TIF 32 KB)

12859_2013_6120_MOESM5_ESM.tif

Additional file 5: Table S3: Correlation coefficients between flanking length and physicochemical properties difference (Flanking regions - general non-MoRF regions) for each property. (TIF 21 KB)

12859_2013_6120_MOESM6_ESM.tif

Additional file 6: Figure S3: ROC plots of MFSPSSMpred tested with different masking-window sizes (from 3 to 19). (TIF 31 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Authors’ original file for figure 4

Authors’ original file for figure 5

Authors’ original file for figure 6

Authors’ original file for figure 7

Authors’ original file for figure 8

Authors’ original file for figure 9

Authors’ original file for figure 10

Authors’ original file for figure 11

Authors’ original file for figure 12

Authors’ original file for figure 13

Authors’ original file for figure 14

Authors’ original file for figure 16

Authors’ original file for figure 17

Authors’ original file for figure 18

Authors’ original file for figure 19

Authors’ original file for figure 20

Authors’ original file for figure 21

Authors’ original file for figure 22

Authors’ original file for figure 23

Authors’ original file for figure 24

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Fang, C., Noguchi, T., Tominaga, D. et al. MFSPSSMpred: identifying short disorder-to-order binding regions in disordered proteins based on contextual local evolutionary conservation. BMC Bioinformatics 14, 300 (2013). https://doi.org/10.1186/1471-2105-14-300

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-14-300

Keywords