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

In-silico prediction of disorder content using hybrid sequence representation

Abstract

Background

Intrinsically disordered proteins play important roles in various cellular activities and their prevalence was implicated in a number of human diseases. The knowledge of the content of the intrinsic disorder in proteins is useful for a variety of studies including estimation of the abundance of disorder in protein families, classes, and complete proteomes, and for the analysis of disorder-related protein functions. The above investigations currently utilize the disorder content derived from the per-residue disorder predictions. We show that these predictions may over-or under-predict the overall amount of disorder, which motivates development of novel tools for direct and accurate sequence-based prediction of the disorder content.

Results

We hypothesize that sequence-level aggregation of input information may provide more accurate content prediction when compared with the content extracted from the local window-based residue-level disorder predictors. We propose a novel predictor, DisCon, that takes advantage of a small set of 29 custom-designed descriptors that aggregate and hybridize information concerning sequence, evolutionary profiles, and predicted secondary structure, solvent accessibility, flexibility, and annotation of globular domains. Using these descriptors and a ridge regression model, DisCon predicts the content with low, 0.05, mean squared error and high, 0.68, Pearson correlation. This is a statistically significant improvement over the content computed from outputs of ten modern disorder predictors on a test dataset with proteins that share low sequence identity with the training sequences. The proposed predictive model is analyzed to discuss factors related to the prediction of the disorder content.

Conclusions

DisCon is a high-quality alternative for high-throughput annotation of the disorder content. We also empirically demonstrate that the DisCon's predictions can be used to improve binary annotations of the disordered residues from the real-value disorder propensities generated by current residue-level disorder predictors. The web server that implements the DisCon is available at http://biomine.ece.ualberta.ca/DisCon/.

Background

The intrinsically disordered proteins (IDPs), also referred to as natively unfolded or intrinsically unstructured proteins, lack stable tertiary structure in vitro. These proteins are implicated in numerous processes including cellular signal transduction, transcriptional regulation, and translation [1], and their prevalence was demonstrated in several human diseases [2, 3], including cancer [4], cardiovascular disease [5], neurodegenerative diseases [6, 7], genetic diseases [8], and amyloidoses [9]. At the same time, the annotations of the IDPs are accumulated at a relatively low pace when compared with the growth of the number of known, non-redundant protein sequences. Over the last decade numerous sequence-derived characteristics, including low complexity [10], which was proposed in [11], high net charge and low content of hydrophobic amino acids [12, 13], lack of regular secondary structure [14], to name just a few, were found to differentiate between disordered and ordered regions. The abovementioned results suggest that disorder can be predicted from the sequence and they motivate the development of computational models for the prediction of the disordered regions. Several such predictors were already developed; see [15] for a recent review. Majority of the existing predictors generate the disorder predictions for each residue in the input protein chain. These per-residue predictors can be divided into 4 categories: i) methods that utilize the relative propensity of amino acids to form disorder/ordered regions which include GlobPlot [16], IUPred [17], FoldIndex [18], and Ucon [19]; ii) methods that are based on classifiers generated with the help of machine learning algorithms, such as DISpro [20], DISOPRED [21], DISOPRED2 [22], PrDOS [23], POODLE predictors [24, 25], PONDR predictors [10, 26, 27], Spritz [28], PROFbval [29], DisPSSMP [30], DisPSSMP2 [31], IUP [32], NORSnet [33] and OnD-CRFs [34]; iii) meta-approach methods that are based on a consensus of multiple base predictors including MULTICOM (also called PreDisorder) [35, 36], metaPrDOS [37] and recent MD [38], MFDp [39], and PONDR-FIT [40] predictors; and iv) approaches that find disordered residues through an analysis of the predicted 3D structural models such as PrDOS [23] and DISOclust [41]. There are also methods that predict the propensity of the entire protein chain to be unstructured [13, 42–44]. One of these approaches is based on the charge-hydropathy plots [13] and another utilizes distributions of the predicted per-residue disorder scores [42–44]. The abovementioned per-residue and per-chain methods perform the predictions in a high-throughput manner and consequently they can be used as a possible solution to close the annotation gap.

Although the per-residue methods are successful in the disorder prediction at the residue level, i.e., they achieve AUC (area under the ROC curve) of about 0.8 [38–40] and MCC (Matthews Correlation Coefficient) of about 0.45 [39] when tested on large benchmark datasets, we observe that they typically make relatively substantial mistakes at the sequence-level. More specifically, these methods may over-or under-predict the overall amount of disorder in the sequence. Tests of 10 recent disorder predictors that include methods from all four groups on a benchmark dataset of 200 chains, see Table 1, show that the mean average (over the dataset) squared errors between the native and the predicted amount of disorder vary between 0.07 and 0.18 (see the Results and Discussion sections). One of the potential reasons for these errors is the fact that virtually all of the most accurate recent predictors, such as NORSnet, DISOPRED2, MD, PONDR-FIT, and MFDp, use a local sequence window to predict the disorder while the information encoded in the entire chain may reveal an overall sequence-level disorder bias. The disorder predictors use identical (for all chains) cut-off values to annotate disordered residues based on the predicted real-value propensities, and we show that these annotations can be improved if the cut-off is adjusted to match the native amount of the disorder in the entire chain, which suggests that the knowledge of this sequence bias could be useful.

Table 1 Comparison of predictive quality of the DisCon and the disorder content extracted from the predictions of the 10 considered modern disorder predictors on the test dataset.

The overall disorder content was used in the past to estimate the abundance of intrinsic disorder in several protein databases [45, 46], in various protein families and classes [47–58], and in complete proteomes [59–62]. The high values of the disorder content were reported for several disease-related proteins [2–9]. The content was also used for the analysis of intrinsic disorder-related protein functions [63–65]. Importantly, in all these and similar cases, the disorder content was evaluated based on the results of either binary classifiers or was derived from the per-residue disorder predictions. As mentioned above, these per-residue disorder prediction methods may over-or under-predict the overall amount of disorder in the sequence. This observation and the fact that the knowledge of the disorder content in a given protein or in a set of proteins of interest or in an entire proteome can be utilized to investigate numerous important hypotheses motivate the development of new computational tools for the accurate prediction of the disorder content.

We propose a novel method, named DisCon (Disorder Content predictor), that aims to provide accurate sequence-based predictions of the disorder content. Our approach is based on the premise that sequence-level aggregation of information may provide more accurate content prediction when compared with the content extracted from the local window-based residue-level disorder predictors. DisCon extends the capabilities of the binary predictors from [13, 42–44] as it provides a real-value, instead of binary, estimates of the amount of the disorder. Our solution has two key characteristics. Firstly, we use a comprehensive selection of the input information sources including sequence, evolutionary profiles generated with PSI-BLAST, and predicted secondary structure, solvent accessibility, B-factors, signal peptides and globular domains. The main reason to use such a diverse set of inputs is to capture different aspects/flavours of disorder [66]. The selection of the first five sources is motivated by their successful use for the residue-level disorder predictions, see Table 2. Similarly as for the ordered proteins, for which the correct folding into biologically active conformations is determined by their amino acid chain, the absence of rigid structure in the intrinsically disordered proteins or regions is also encoded in their amino acid sequences [67, 68]. The disordered regions are usually depleted in so-called order-promoting residues (Trp, Tyr, Phe, Ile, Leu, Val, Cys, and Asn) while they include larger numbers of the disorder-promoting residues (Ala, Arg, Gly, Gln, Ser, Glu, Pro, and Lys) [10, 42, 67, 69, 70]. Moreover, the disorder is often observed in parts of the sequence that are characterized by low complexity, higher number of Pro and charged residues, and lower amount of hydrophobic and bulky amino acids [10, 12, 13, 27, 71], which motivates the use of the input protein chain. The predicted secondary structure is useful since many of the disordered regions are characterized by lack of the secondary structure [14, 33, 66, 69, 71]. These unstructured regions usually have a large solvent-accessible area [38] which motivates the application of the predicted solvent accessibility. High values of B-factors are often associated with disordered regions [72], which is why we use their predictions to implement inputs to DisCon. The signal peptides were previously found relevant in the context of analysis of differences between disorder predictors [73] and residues in domains are less likely to be disordered and this information was used in the NORSnet predictor [33]. Secondly, we perform a careful design of the input features that are computed based on the aggregation of values for each information source and also by combining information sources. We also designed several customized features that quantify the size and relative location of the predicted secondary structure segments.

Table 2 List of input information sources used by the disorder predictors considered in this work.

We empirically demonstrate that the DisCon's predictions are more accurate than the content extracted from the residue-level annotations generated by modern disorder predictors, including methods listed in Table 2. One of the potential applications of the predicted disorder content is to adjust the cut-offs used by the disorder predictors to annotate the disordered residues. We show that these annotations can be improved when the threshold values is adjusted for each chain such that the amount of the predicted disordered amino acids matches not only the native but also the predicted content.

Methods

Definition of Disorder

In the past CASP experiments the disordered residues were defined as the amino acids that lack coordinates in their crystal structures and, in the case of the structures solved by NMR, as the amino acids that exhibit high variability within an ensemble or that were annotated by experimentalists as disordered in the REMARK 465 [74, 75]. Another commonly used source for the disorder annotations is based on the experimentally-validated and biologically relevant disordered segments from the DisProt database [76]. We note that the assignment of the disordered regions using different experimental methods was previously shown to be potentially inconsistent [66]. Consequently, the disorder predictors that were developed using annotations provided with one method could lead to larger errors when used to predict annotations generated with the help of other methods [19, 33]. Therefore, we created a dataset that combines the CASP-defined annotations with the DisProt annotations.

Datasets

The proposed method was designed and tested using a dataset that was developed to validate a recent meta-predictor of disordered residues, the MFDp [39]. The protein chains were collected from the Protein Data Bank (PDB) [77] and the DisProt [76] databases. The culled PDB list [78] was used to derive a high-quality and low sequence identity subset of the PDB protein. More specifically, only the proteins for which the structure is characterized by R-factor < 0.2 and resolution < 2.0Ã…, and that are characterized by sequence identity < 25% were kept. We randomly selected 20% of the fully structured proteins among the resulting chains. This is motivated by the fact that many of chains selected using the culled PDB list are annotated as ordered while a recent study shows that completely ordered proteins are not highly abundant in PDB [46]. The PDB chains were combined with all 523 proteins from the release 4.9 of the DisProt. The resulting dataset was filtered to reduce the pairwise sequence identity to below 25% by removing similar sequence with fewer disordered residues. Among the remaining 514 chains we removed four for which MD failed to produce predictions; this also resulted in lack of predictions from Ucon, PROFbval and NORSnet that are bundled with the MD predictions. Moreover, we improved the annotations of the DisProt chains using the procedure described in [79]. We applied the approach based on the SL dataset [79] that combines the disorder annotations from the DisProt with the annotations of disorder and order based on the corresponding structural domains that can be found in PDB. We note that in contrast to the SL dataset that is based on the release 4.5 of DisProt, our annotations are based on the newer release 4.9. Finally, we also removed the HIS-tags that are introduced to ease the crystallization. The resulting dataset includes 305 chains from DisProt and 205 from PDB. This dataset was divided at random into two subsets, the training dataset with 310 chains and the test dataset with the remaining 200 chains.

We note that although there is some overlap between the training and test sequences (depending on the alignment tool used), they are mostly independent at the 25% similarity level. The training dataset was used to develop the predictor including selection of the input features and the parameterization of the prediction model, which were performed based on the 5-fold cross validation protocol. Next, our predictor that was computed using the training dataset was compared with the existing per-residue prediction methods using the test dataset. The training and test datasets are available at http://biomine.ece.ualberta.ca/DisCon/.

Evaluation criteria

The disorder content predictions are evaluated using three measures:

where y i ∈ Y is the native and x i ∈ X is predicted disorder content for the ith protein chain, avg X and avg Y are the sample means of X and Y, s x and s y are the sample standard deviations of X and Y.

Following [19, 22, 33, 39] the binary, per-residue disorder predictions and the per-chain predictions of disorder content (protein categorized based on a given disorder content cut-off) were assessed using two measures:

where TP is the number of true positives (correctly predicted disordered residues), FP denotes false positives (structured residues that were predicted as disordered), TN denotes true negatives (correctly predicted structured residues), and FN denotes false negatives (disordered residues that were predicted as structured). Accuracy quantifies the overall success rate, i.e., fraction of correct predictions among all prediction, but since it may lead to misleading results when the dataset is unbalanced (which is the case here since majority of residues are structured) we also use MCC. The MCC values range between -1 and 1 and they are equal zero when all residues are predicted to be structured or to be disordered. Higher values of PCC, accuracy and MCC and lower values of MSE and MAE correspond to better predictions. We also evaluated the real-value, per-residue disorder predictions based on the area under the ROC curve (AUC) measure.

Overview of the proposed predictor

The prediction of the disorder content is performed in three steps, see Figure 1. First, the input protein chain is processed through PSI-BLAST [80] to generate Position Specific Scoring Matrix (PSSM) and Weighted Observed Percentage (WOP) profiles using the NCBI's nr database downloaded on Nov 19th 2009, which was filtered using PFILT [81] to remove low-complexity regions, transmembrane regions, and coiled-coil segments. We use PSIPRED [82] to obtain the 3-state predicted secondary structure (SS), Real-SPINE3 [83] for the prediction of the relative solvent accessibility (RSA), PROFbval [29] for the B-factor and residue flexibility predictions, and IUPred [17] to predict globular domains. The former three methods use the PSSM profiles as their inputs. Second, the above predictions, profiles and sequence are used to generate a set of numerical descriptors/features. These features quantify the information encoded in each of the predictions/profiles and also between multiple predictions/profiles. We performed feature selection to select a small subset of 29 features that are relevant to the prediction of the disorder content. Third, the selected descriptors are fed into a ridge regression model to generate the predictions.

Figure 1
figure 1

Architecture of the DisCon predictor.

Feature-based encoding of the input protein sequence

The input sequence is processed to generate predictions of the 3-state SS, RSA, normalized [72] real-value B-factors, binary annotation of the residue flexibility as provided by PROFbval in two modes, the strict and the non-strict [29], binary annotation of residues that form globular domains, and sequence profiles encoded using PSSM and WOP. We normalized the ASA values predicted by Real-SPINE3 using the maximal ASA values provided in [83] and we preprocessed the 3-state SS by converting the predicted helices that had < 3 residues into coils. We also binarized the real-values RSA to annotate the residues as either solvent exposed when RSA > 0.25 or buried when RSA ≤ 0.25; this cut-off value was used in past studies [72, 84]. We also attempted to use signal peptide prediction provided by RPSP [85], but these features were removed during the feature selection. Detailed description of features is provided in Table A1 in the Additional File 1. We generated total of 614 features that are based on

  • composition of amino acids

  • length and relative location of predicted helix, strand and coil segments

  • composition of solvent exposed residues

  • composition of flexible residues and sequence segments composed of flexible residues

  • number and size of sequence segments that correspond to predicted globular domains

  • composition of residue predicted as signal peptides

  • fusion of the information coming from multiple predictions, including SS states, solvent exposure, flexibility, and domain annotations. We consider all combinations of two, three and four of the above predictions.

  • aggregations of the sequence profiles using entropy and relative (using background probability) entropy by both rows and columns of the PSSM and WOP

  • entropy-based aggregations of the sequence profiles encoded with PSSM and WOP which is performed for specific amino acid types, and for residues characterized by specific SS state, solvent exposure, flexibility, and domain annotations.

We emphasize that most of the features, in particular the features that are based on the secondary structure segments, flexible sequence segments, and that combine multiple predicted properties, are novel and unique to this work.

Design of the predictive model

The features were generated to comprehensively cover information that can be extracted from each predicted property, sequence and sequence profile, and their combinations. Consequently, some of these inputs may not be relevant to the prediction of the disorder content and some could be redundant with each other. We performed two-step feature selection to find a small set of non-redundant and relevant features; the second step also includes computation and parameterization of the predictive model. First, we remove the irrelevant and redundant features using a coarse-grained evaluation based on correlation, and next we perform a wrapper-based selection using the remaining features.

In the first step, for each feature we compute its average PCC with the disorder content (the PCC values are based on 5-fold cross validation on the training dataset and they are averages of the coefficients computed in the five training folds) and we remove the features with average absolute PCC value < 0.2. We selected the 0.2 cut-off as this value corresponds to a visible dip in the distribution of the correlation values, see Figure A1 in the Additional File 1. Next, we filtered the remaining 322 features to remove redundancy by assuring that the maximal average absolute PCC value between any pair of these features is < 0.9. Starting with the feature that has the highest average absolute PCC with the native content, we added another feature into the set of filtered features if the average (over the five training folds) absolute PCC between this feature and each feature which is already in the set of filtered features was < 0.9.

In the second step, we use the remaining 152 features to perform wrapper-based selection in which a subset of features that results in favourable performance in prediction of the disorder content is retained. We consider two types of predictors, ridge regression and Support Vector Regression (SVR) [86]. The selection of the regression model is motivated by its successful application in several related areas, including evaluation of peptide identification [87] and prediction of folding rates [88, 89], solvent accessibility [90], secondary structure content [91], and affinity of protein-ligand complexes [92], to name a few. Similarly, the SVR also enjoys a wide range of relevant applications including prediction of B-factors [93], solvent accessibility [94], and half-sphere exposure [95]. The values of the regression coefficients and the SVR models were estimated from the data in the training folds using WEKA workbench [96]. We consider three types of kernel functions to build SVR models, polynomial, Radial Basis Function (RBF), and Pearson VII function-based Universal Kernel (PUK) [97]. We parameterized each kernel and the complexity constant C by performing grid search. We use linear and quadratic polynomials, and C equal 2 x where x = -8, -7,..., 2; the RBF kernel with gamma (spread) equal 2 y where y = -11, -10,..., 2, and C values where x = -3, -2,..., 6; and the PUK kernel with omega equal 2 z where z = -4, -3,..., 1, and C values where x = -4, -3,..., 5. We also parameterized the ridge parameter in the ridge regression; we considered ridge values equal 10 w where w = -11, -9,..., 2. We first parameterized these 4 predictors (3 SVM types + 1 ridge regression) using a representative subset of the 152 features. We selected one features with the highest average absolute PCC from each of the feature groups defined in Table A1 in the Additional File 1. The representative subset includes 23 features since that number of groups was covered among the 152 features. Next, these parameterized predictors were used to perform feature selection in which we searched for a subset of features that results in the best MSE value. We performed forward and backward best first searches. The forward/backward best first search starts with the empty/entire (152 features) set of features, and it adds/removes one feature at the time if it decreases/increases the MSE value. The search stops when the entire list of features is scanned. As a result, we obtained 8 configurations of 4 predictors with 2 search types. The predictors in each configuration were parameterized using the grid search as described above. The parameterizations and all steps of the feature selection were executed based on multiple repetitions of 5-fold cross validations on the training dataset, and they aimed to minimize the average MSE score between the predicted and the native disorder content. We repeated the cross validations for up to five times using randomized division into the 5 folds for as long as the coefficient of variation (the ratio of the standard deviation to the mean) was below 0.02; this approach should assure a robust estimate of the MSE values. The parameters of the four predictors and the corresponding number of the selected features are given in Table A2 in the Additional File 1. The predictive performance, which was evaluated based on 5-fold cross validation on the training dataset, for the eight configurations is summarized in Table A3 in the Additional File 1. The best performance, in terms of the MSE and PCC values, is achieved with the ridge regression that uses 29 features selected using the forward best first search, and this configuration is used to implement the proposed DisCon predictor.

Results and Discussion

Disorder content prediction

We compare the performance of the DisCon with the results obtained using the disorder content computed from the disorder predictions generated by DISOPRED2, IUPred (both versions, IUPredL and IUPredS), PROFbval, NORSnet, Ucon, DISOclust, MD, PONDR-FIT, and MFDp methods. For the per-residue predictors we used the web servers or standalone implementations provided by the authors, and we calculated the content by counting the number of residues predicted as disordered and dividing it by the length of the corresponding protein chain. The results are computed on the test dataset with 200 chains which shares low identity to chains in our training dataset. We note that the methods we compare with use training datasets that may share higher similarity with the chains in our dataset, which could inflate their predictive quality. We also analyze statistical significance of the differences between the content predicted by DisCon and the other methods. We compare the per-chain values of the absolute errors and the squared errors over the 200 chains in the test dataset and the Pearson correlation coefficients computed for 200 randomly selected sets of 100 proteins from the test dataset. Since the measurements follow normal distribution (evaluated using Shapiro-Wilk test at 0.05 significance) we apply the paired t-test and we measure significance of the differences at 0.05 and 0.01 levels. We evaluate the extent of the over-and under-prediction of the disorder content by quantifying the number of the over-and under-predicted chains and the corresponding MAE values and we also provide the AUC, accuracy, and MCC values for the 10 considered per-residue predictors. The results are summarized in Table 1.

The DisCon is shown to provide favourable predictive performance. It obtains MSE equal 0.05 and PCC equal 0.68 on the test dataset. We note that these results are consistent with the results obtained based on the 5-fold cross validation on the training dataset (PCC = 0.70, MSE = 0.05; see Table A3 in the Additional File 1). On the test set, the best performing per-residue disorder predictors are worse than DisCon by 0.016 MSE and 0.07 PCC for the disorder content prediction. The average absolute error of DisCon equals 0.156 when compared with value at or over 0.167 obtained with the current disorder predictors, except for IUPredS for which MAE = 0.155. The improvements in MSE and PCC offered by DisCon are shown to be statistically significant when compared with all considered competitors. The MAE values computed from our predictions are significantly better than the errors based on the predictions with four existing methods and are equivalent with the remaining six predictors. Further analysis reveals that the quality of the DisCon predictions is better for longer chains, while some other methods may produce favorable predictions for short chains. Figure A2 in the Additional File 1, which shows the relation between the chain length and the absolute errors generated by the DisCon and the top-three methods from Table 1, i.e., IUPRedS, Ucon, and PONDR-FIT, demonstrates that the proposed predictor is characterized by smaller absolute errors for longer chains, while the other three predictors on average provide more accurate predictions for short chains. DisCon provides relatively balanced predictions with similar number of over-and under-predicted chains and low MAE values for these two types of errors. We observe that PROFbval and DISOclust are characterized by substantial levels of the over-prediction of the disorder content that are expressed by the large number of the over-predicted chains and/or high MAE for the over-predicted chains. The under-prediction of the disorder content is characteristic for the NORSnet method. Table 1 also shows that the Ucon which obtains relatively low MSE and mid-range PCC is characterized by lower quality of the per-residue predictions with MCC = 0.28.

Binary prediction of the disorder amount

We apply the predicted disorder content on the test dataset to perform binary prediction of chains that are characterized by the amount of the disorder below/above a specific cut-off value. The cut-offs at 1 and 0 corresponds to detection of fully-disordered and fully-ordered proteins, respectively, while the intermediate cut-off values could be used to find partially structured chains. We measure the MCC of these binary predictions defined using the cut-off values between 0 and 1 with step of 0.05, and the content predicted by DisCon and extracted from the outputs of the 10 considered disorder predictors, see Figure 2. We observe that the usage of the content predicted by DisCon results in the highest MCC values, which range 0.53 and 0.61, for cut-offs between 0.35 and 0.65 inclusive. The best results for the cut-off values above 0.65 are achieved with the content predicted with MD while several predictors do comparably well for the cut-offs below 0.35. We note that the prediction of the fully-structured and fully-disordered proteins could be also accomplished using specialized predictors, such as NN-CDF [44].

Figure 2
figure 2

The MCC values ( y -axis) for the binary prediction of chains that are characterized by the amount of the disorder below/above a cut-off value shown on the x -axis. The binary predictions are computed by thresholding the predicted disorder content generated by DisCon and the 10 considered disorder predictors on the test dataset.

Content guided thresholding of the real-value disorder prediction

The disorder predictors usually provide both the real-value propensity of the disorder and the binary order/disorder assignment for each residue. The binary assignment is usually based on thresholding of the real-value propensities with a fixed cut-off. We investigate whether the predicted disorder content could be used to guide the selection of the threshold value. This means that instead of using the fixed cut-off we adjust the threshold such that the amount of the residues annotated as disordered equals to the amount of the disorder content predicted with DisCon (using the predictions on the test dataset). We evaluate the binary residue-level disorder prediction of the original and the content-adjusted disorder predictors using MCC, see Figure 3, and accuracy, see Figure A3 in the Additional File 1. We note that the 10 methods that are considered in these figures provide continuous (real-value) prediction values that represent the propensity for a residue to be disordered. We compare the original binary disorder predictions, predictions that are based on a fixed cut-off that maximizes the MCC of a given predictor on our test dataset (to remove a potential bias due to the usage of a different or less complete, i.e., older, disorder annotation to select the original cut-off), and predictions where the cut-off is selected to match the content predicted by DisCon. Based on the Figure 2 that shows that MD outperforms DisCon when the native content < 0.1 or > 0.65, we also consider combining content predicted by the MD and DisCon to adjust the threshold. If the content predicted by MD is > 0.65 or < 0.1 for a given chain then we use the MD predicted content; otherwise we use the content predicted by the DisCon method.

Figure 3
figure 3

The MCC values for the residue-level disorder prediction adjusted using content predicted by DisCon. The bar chart includes the original predictions (densely dotted red bars), predictions with a fixed cut-off that is optimized to maximize MCC on the entire test dataset (sparsely dotted blue bars), predictions where the content predicted by DisCon is used to adjust the cut-off (solid black bars) and where the content predicted by MD if its values are > 0.65 or < 0.1 and otherwise content predicted by DisCon are used to adjust the cut-off (solid green bars). The results were computed on the test dataset and the methods on the x-axis are sorted by their original MCC values.

The thresholding of the predicted real-value disorder using the content predicted by DisCon leads to improvements in both MCC and accuracy for all predictors except for the MD, in which case the accuracy is slightly improved but the MCC is lower. The average (across all methods) improvement in MCC and accuracy equal 0.03 and 0.05, respectively. When we use the combination of the content predicted with DisCon and MD the improvements are more substantial and they range between 0.01 and 0.14 for the MCC (on average 0.06) and 0.01 and 0.24 for the accuracy (on average 0.05); the best MCC is obtained using the predictions from the MFDp and it equals 0.45 when compared with 0.425 that was obtained without the content-based adjustment. Interestingly, using this cut-off adjustment the MCC values obtained by seven out of ten considered predictors are > 0.4 while originally (with the default cut-offs), see Table 1, only two methods have MCC > 0.4. This suggests that majority of the considered disorder predictors differentiate between structured and disordered residues based on their real-value propensities in a given chain with relatively similar quality, but only a few of them can accurately scale the range of the real-value propensities between sequences. The content-guided selection of the cut-offs alleviates the prediction bias, i.e., the tendency to under-or over-predict the amount of disorder. The binary predictions of PROFbval, DISOclust, and NORSnet that are originally characterized by relatively low MCC values and a bias towards either over-or under-prediction, see Table 1, are shown to improve by a wide margin when using the disorder content predicted by DisCon or by the combination of DisCon and MD. We observe that the relatively poor performance of the Ucon method does not stem from the prediction bias but rather from its overall problems with the quality of the residue-level annotations, as evidenced by the relatively low AUC and MCC in Table 1 which is in contrast to the sequence-level amount of disorder that is predicted quite accurately by this method.

We visualize the improvements that result from the cut-off adjustments using two case studies, one where the original predictions over-estimate the native amount of disorder and another where the predictions are under-estimated. In both cases, we compare the original binary annotations of disordered residues with the annotations that are adjusted using the content predicted with DiscCon, and we include predictions from the top six methods from Table 1, i.e., Ucon, MD, MFDp, PONDR-FIT, IUPredS, and DISOPRED2. The first example is the apoptosis-inducing ligand 2 (Apo2L) protein (PDB ID 1DG6 chain A), see Figure 4. This protein was solved using high-resolution (at 1.3 Ã…) X-ray crystallography and the structure includes two relatively short disordered segments in the vicinity of the N-terminus [98]. Figure 4 reveals that all six predictors annotate disorder at the N-terminus and that PONDR-FIT and IUPredS also predict a short disordered segment at the C-terminus. However, the disorder at the N-terminus is over-predicted; the residues between positions 30 and 40 and between 52 to about 60 are predicted as disordered, while the X-ray structure shows them as structured. These over-predictions were minimized when the cut-off was adjusted to match the content predicted by the DisCon. Importantly, the adjusted predictions identify the two disordered segments, with particularly good results for the MD and MFDp predictors that quite accurately identify both of the disordered segments and no other disordered residues. After the adjustment, the predictions from PONDR-FIT show the two disordered segments at the N-terminus, although the first segment is predicted to be 11 residues too short, and the C-terminus is predicted as structured. Similarly, the adjusted predictions from IUPredS and DISOPRED2 show disorder in the vicinity of the N-terminus, while the disorder predictions in other parts of this chain are removed.

Figure 4
figure 4

Prediction of disordered residues in the apoptosis-inducing ligand 2 (Apo2L) protein (PDB ID 1DG6 chain A) by Ucon (thin blue line), PONDR-FIT (thin red line), MD (thin green line), MFDp (thin gray line), IUPredS (thin pink line), and DISOPRED2 (thin cyan line) predictors. The original cut-offs are shown using dashed lines. The native disordered regions are annotated using black horizontal line. The original binary predictions from Ucon, PONDR-FIT, MD, MFDp, IUPredS, and DISOPRED2 are denoted using blue (at the -0.1 point on the y-axis), red (at the -0.2), green (at the -0.3), gray (at the -0.4), pink (at the -0.5), and cyan (at the -0.6) horizontal lines. The binary predictions that were adjusted to match content predicted with DisCon are shown using horizontal bright green lines located immediately under the lines that show the original predictions.

The second case is the inosine-5'-monophosphate dehydrogenase (DisProt ID DP00399) which was also solved using X-ray crystallography, see Figure 5. This protein includes four disordered segments, one longer between positions 102 and 221, and three shorter towards the C-terminus [99]. Overall, the six predictors under-predicted the disorder levels in this protein. They predicted only a few disordered residues at both termini, with the exception of Ucon that predicted about a dozen of short disordered segments throughout the entire chain and MFDp that predicted three disordered segments, including both termini and a segment between positions 421 and 434. The DisCon predicted 28.6% of residues to be disordered, when compared with 1%, 1.2%, 2.8%, 3%, 6.6%, and 12% that were predicted by the DISOPRED2, MD, IUPredS, PONDR-FIT, MFDp, and Ucon, respectively; the native amount of disorder is 33.7%. The content-adjusted annotation of the disordered residues captures a large number of disordered residues in the long segment between positions 102 and 221, as well as the two disordered segments nearest to the C-terminus. These improvements come as a trade-off for an over-prediction of the disorder at the N-terminus, particularly for the PONDR-FIT and DISOPRED2 predictors. Overall, we observe that the adjusted predictions show a denser concentration of the disordered residues around the natively disordered regions.

Figure 5
figure 5

Prediction of disordered residues in the inosine-5'-monophosphate dehydrogenase protein (DisProt ID DP00399) by Ucon (thin blue line), PONDR-FIT (thin red line), MD (thin green line), MFDp (thin gray line), IUPredS (thin pink line), and DISOPRED2 (thin cyan line) predictors. The original cut-offs are shown using dashed lines. The native disordered regions are annotated using black horizontal line. The original binary predictions from Ucon, PONDR-FIT, MD, MFDp, IUPredS, and DISOPRED2 are denoted using blue (at the -0.1 point on the y-axis), red (at the -0.2), green (at the -0.3), gray (at the -0.4), pink (at the -0.5), and cyan (at the -0.6) horizontal lines. The binary predictions that were adjusted to match content predicted with DisCon are shown using horizontal bright green lines located immediately under the lines that show the original predictions.

We conclude that although predictions shown in the two case studies should not be assumed typical, they demonstrate that the content predicted with DisCon offers valuable assistance in selection of the cut-offs to annotate the disordered residues based on the real-values predictions from modern disorder predictors.

Factors related to the amount of disorder/order

We convert the input protein sequence into a custom-designed set of selected 29 numerical descriptors which utilize information related to the evolutionary profiles, sequence itself, and predicted secondary structure (SS), solvent accessible residues (RSA), B-factors and globular domains; see Table A4 in the Additional File 1. Majority of the selected features combine multiple input sources. For instance, the largest group of 5 similar features is based on counting the residues in certain predicted SS states, with certain levels of predicted solvent exposure and B-factors which are located within the predicted domains. For instance, the SSHEDOMinBFNSlowRSAB feature counts the predicted helix and strand residues that are in globular domains, have low B-factors (are structurally rigid) and are buried. As expected, these values are negatively correlated with the content of the native disorder (with PCC = -0.54 in our dataset), see Figure 6, since this descriptor highlights hallmarks of well-structured proteins, i.e., they usually include conserved domains with buried strands and/or helices that are usually structurally rigid. Another feature that attains negative, -0.34, correlation with the native content is BFNSlowSeg10, which quantifies the number of predicted rigid residues (with low B-factor) that are grouped together in the sequence in segments of size at least 10. Figure 6 shows a cluster of BFNSlowSeg10 values between 0.4 and 0.7 for chains with low amount of native disorder, which suggests that well-structured proteins include significant amount of rigid residues that are grouped together in the sequence, while disordered chains contain fewer of such rigid residues. We also discuss two features that have relatively high positive correlation with the disorder content. The SSCHBFNShighDOMnotin feature counts the number of coil and helix residues with high B-factor that are not in the globular domains. This feature that has PCC = 0.54 in our dataset, see Figure 6, agrees with characteristic properties of the disorder, which often concerns flexible residues that are outside of the domains. The CHC...CHSeg feature (PCC = 0.45 in our dataset) computes the number of residues in the longest segment in the predicted SS that does not include strands, i.e., stretch of the sequence that is composed only of a coils or of alternating helix-coil-helix-coil... segments. This feature demonstrates that disordered proteins are often depleted of beta-sheets.

Figure 6
figure 6

Scatter plots of the relations between the values of selected four input features, SS_HE-DOM_in-BFNS_low-RSA_B (green circle markers), BFNS_low-Seg_10 (red triangle markers), SS_CH-BFNS_high-DOM_notin (black × markers), and CHC ... CHSeg (blue hollow circle markers) shown on y -axi s and the native disorder content given on x -axis. The lines correspond to linear regressions with the corresponding R2 values.

The input features also highlight the importance of the relation between sequence conservation and the amount of the disorder, i.e., 13 out of the 29 features utilize entropy computed from the PSSM or WOP profiles. For instance, the EntAvePSSM feature, which computes the entropy of the average PSSM scores for each column (amino acid type) in the matrix along the sequence, has PCC = -0.5. This means that well-structured proteins are characterized on average by a stronger degree of sequence conservation when compared with the disordered proteins. Our observation is in agreement with the results of previous study, where the evolution rates of ordered and intrinsically disordered regions were compared using the pairwise genetic distances between the ordered and the disordered regions of 26 protein families having at least one member with a structurally characterized region of disorder of 30 or more consecutive residues [100]. This study demonstrated that the disordered regions evolved significantly more rapidly than the ordered regions in 19 of the 26 families studied [100].

Conclusions

In spite of the fact that the quality of the high-throughput disorder prediction continues to improve [75], researchers recognize that new and more accurate predictors are still needed [38, 39]. We address the shortage of accurate methods that predict the overall amount of disorder in a given protein chain, which is motivated by the fact that current disorder predictors tend to provide relatively inaccurate estimates of the disorder content. We propose a novel approach, called DisCon, which combines information derived from sequence, sequence profiles, and predicted secondary structure, solvent accessibility, flexibility, and annotation of globular domains. We custom designed feature-based representation of the input protein chain that aggregates and combines these inputs and we performed feature selection that found a small set of 29 complementary features that are well correlated with the native disorder. Using these features and a ridge regression-based model, the DisCon predicts the disorder content with low, 0.05, mean squared error and high, 0.68, correlation, as evaluated on an independent test dataset. These predictions are empirically shown to be significantly better than the disorder content estimates derived from outputs of ten modern disorder predictors. The DisCon's predictions provide a high-quality alternative for high-throughput annotation of the disorder content. They are also shown to provide useful input to improve binary annotations of the disordered residues from the real-value disorder propensities generated by current disorder prediction methods.

Abbreviations

AUC:

Area Under the ROC Curve

ASA:

Absolute Solvent Accessibility

CASP:

Critical Assessment of Techniques for Protein Structure Prediction

CDF:

Cumulative Distribution Function

CH-plot:

Charge-Hydropathy plot

DisCon:

Disorder Content predictor

IDP:

Intrinsically Disordered Protein

MAE:

Mean Absolute Error

MCC:

Matthews Correlation Coefficient

MSE:

Mean Squared Error

NMR:

Nuclear Magnetic Resonance

PCC:

Pearson Correlation Coefficient

PDB:

Protein Data Bank

PSSM:

Position Specific Scoring Matrix

PUK:

Pearson VII function-based Universal Kernel

RBF:

Radial Basis Function

RSA:

Relative Solvent Accessibility

SS:

Secondary Structure

SVR:

Support Vector Regression

WOP:

Weighted Observed Percentage.

References

  1. Dunker AK, Oldfield CJ, Meng J, Romero P, Yang JY, Chen JW, Vacic V, Obradovic Z, Uversky V: The unfoldomics decade: an update on intrinsically disordered proteins. BMC Genomics 2008, 9(Suppl 2):S1. 10.1186/1471-2164-9-S2-S1

    Article  Google Scholar 

  2. Uversky VN, Oldfield CJ, Midic U, Xie H, Vucetic S, Xue B, Iakoucheva LM, Obradovic Z, Dunker AK: Unfoldomics of human diseases: Linking protein intrinsic disorder with diseases. BMC Genomics 2009, 10(Suppl 1):S7. 10.1186/1471-2164-10-S1-S7

    Article  PubMed Central  PubMed  Google Scholar 

  3. Uversky VN, Oldfield CJ, Dunker AK: Intrinsically disordered proteins in human diseases: introducing the D2 concept. Annu Rev Biophys 2008, 37: 215–246. 10.1146/annurev.biophys.37.032807.125924

    Article  CAS  PubMed  Google Scholar 

  4. Iakoucheva LM, Brown CJ, Lawson JD, Obradovic Z, Dunker AK: Intrinsic disorder in cell-signaling and cancer-associated proteins. J Mol Biol 2002, 323: 573–584. 10.1016/S0022-2836(02)00969-5

    Article  CAS  PubMed  Google Scholar 

  5. Cheng Y, LeGall T, Oldfield CJ, Dunker AK, Uversky VN: Abundance of intrinsic disorder in protein associated with cardiovascular disease. Biochemistry 2006, 45: 10448–10460. 10.1021/bi060981d

    Article  CAS  PubMed  Google Scholar 

  6. Raychaudhuri S, Dey S, Bhattacharyya NP, Mukhopadhyay D: The role of intrinsically unstructured proteins in neurodegenerative diseases. PLoS One 2009, 4(5):e5566. 10.1371/journal.pone.0005566

    Article  PubMed Central  PubMed  Google Scholar 

  7. Uversky VN: Intrinsic disorder in proteins associated with neurodegenerative diseases. Front Biosci 2009, 14: 5188–5238. 10.2741/3594

    Article  CAS  Google Scholar 

  8. Midic U, Oldfield CJ, Dunker AK, Obradovic Z, Uversky VN: Protein disorder in the human diseasome: Unfoldomics of human genetic diseases. BMC Genomics 2009, 10(Suppl 1):S12. 10.1186/1471-2164-10-S1-S12

    Article  PubMed Central  PubMed  Google Scholar 

  9. Uversky VN: Amyloidogenesis of natively unfolded proteins. Curr. Alzheimer Res 2008, 5(3):260–287. 10.2174/156720508784533312

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  10. Romero P, Obradovic Z, Li X, Garner EC, Brown CJ, Dunker AK: Sequence complexity of disordered protein. Proteins 2001, 42: 38–48. 10.1002/1097-0134(20010101)42:1<38::AID-PROT50>3.0.CO;2-3

    Article  CAS  PubMed  Google Scholar 

  11. Wootton JC, Federhen S: Statistics of local complexity in amino acid sequences and sequence databases. Comput Chem 1993, 17: 149–163. 10.1016/0097-8485(93)85006-X

    Article  CAS  Google Scholar 

  12. Dyson HJ, Wright PE: Intrinsically unstructured proteins and their functions. Nat Rev Mol Cell Biol 2005, 6: 197–208. 10.1038/nrm1589

    Article  CAS  PubMed  Google Scholar 

  13. Uversky VN, Gillespie JR, Fink AL: Why are "natively unfolded" proteins unstructured under physiologic conditions? Proteins 2000, 41: 415–427. 10.1002/1097-0134(20001115)41:3<415::AID-PROT130>3.0.CO;2-7

    Article  CAS  PubMed  Google Scholar 

  14. Liu J, Tan H, Rost B: Loopy proteins appear conserved in evolution. J Mol Biol 2002, 322: 53–64. 10.1016/S0022-2836(02)00736-2

    Article  CAS  PubMed  Google Scholar 

  15. He B, Wang K, Liu YL, Xue B, Uversky VN, Dunker AK: Predicting intrinsic disorder in proteins: An overview. Cell Research 2009, 19(8):929–949. 10.1038/cr.2009.87

    Article  CAS  PubMed  Google Scholar 

  16. Linding R, Russell RB, Neduva V, Gibson TJ: GlobPlot: Exploring protein sequences for globularity and disorder. Nucleic Acids Res 2003, 31: 3701–3708. 10.1093/nar/gkg519

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Dosztányi Z, Csizmok V, Tompa P, Simon I: IUPred: web server for the pre-diction of intrinsically unstructured regions of proteins based on estimated energy content. Bioinformatics 2005, 21: 3433–3434. 10.1093/bioinformatics/bti541

    Article  PubMed  Google Scholar 

  18. Prilusky J, Felder CE, Zeev-Ben-Mordehai T, Rydberg EH, Man O, Beckmann JS, Silman I, Sussman JL: FoldIndex: a simple tool to predict whether a given protein sequence is intrinsically unfolded. Bioinformatics 2005, 21: 3435–3438. 10.1093/bioinformatics/bti537

    Article  CAS  PubMed  Google Scholar 

  19. Schlessinger A, Punta M, Rost B: Natively unstructured regions in proteins identified from contact predictions. Bioinformatics 2007, 23: 2376–2384. 10.1093/bioinformatics/btm349

    Article  CAS  PubMed  Google Scholar 

  20. Hecker J, Yang JY, Cheng J: Protein disorder prediction at multiple levels of sensitivity and specificity. BMC Genomics 2008, 9(Suppl 1):S9. 10.1186/1471-2164-9-S1-S9

    Article  PubMed Central  PubMed  Google Scholar 

  21. Jones DT, Ward JJ: Prediction of disordered regions in proteins from position specific score matrices. Proteins 2003, 53(Suppl 6):573–578.

    Article  CAS  PubMed  Google Scholar 

  22. Ward JJ, McGuffin LJ, Bryson K, Buxton BF, Jones DT: The DISOPRED server for the prediction of protein disorder. Bioinformatics 2004, 20: 2138–2139. 10.1093/bioinformatics/bth195

    Article  CAS  PubMed  Google Scholar 

  23. Ishida T, Kinoshita K: PrDOS: prediction of disordered protein regions from amino acid sequence. Nucleic Acids Res 2007, 35: W460–464. 10.1093/nar/gkm363

    Article  PubMed Central  PubMed  Google Scholar 

  24. Shimizu K, Muraoka Y, Hirose S, Tomii K, Noguchi T: Predicting mostly disordered proteins by using structure-unknown protein data. BMC Bioinformatics 2007, 8: 78. 10.1186/1471-2105-8-78

    Article  PubMed Central  PubMed  Google Scholar 

  25. Hirose S, Shimizu K, Kanai S, Kuroda Y, Noguchi T: POODLE-L: a two-level SVM prediction system for reliably predicting long disordered regions. Bioinformatics 2007, 23: 2046–2053. 10.1093/bioinformatics/btm302

    Article  CAS  PubMed  Google Scholar 

  26. Peng K, Vucetic S, Radivojac P, Brown CJ, Dunker AK, Obradovic Z: Optimizing long intrinsic disorder predictors with protein evolutionary information. J. Bioinform. Comput. Biol 2005, 3: 35–60. 10.1142/S0219720005000886

    Article  CAS  PubMed  Google Scholar 

  27. Peng K, Radivojac P, Vucetic S, Dunker AK, Obradovic Z: Length-dependent prediction of protein intrinsic disorder. BMC Bioinformatics 2006, 7: 208. 10.1186/1471-2105-7-208

    Article  PubMed Central  PubMed  Google Scholar 

  28. Vullo A, Bortolami O, Pollastri G, Tosatto SC: Spritz: a server for the predic-tion of intrinsically disordered regions in protein sequences using kernel machines. Nucleic Acids Res 2006, 34: W164–168. 10.1093/nar/gkl166

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  29. Schlessinger A, Yachdav G, Rost B: PROFbval: predict flexible and rigid residues in proteins. Bioinformatics 2006, 22: 891–893. 10.1093/bioinformatics/btl032

    Article  CAS  PubMed  Google Scholar 

  30. Su CT, Chen CY, Ou YY: Protein disorder prediction by condensed PSSM considering propensity for order or disorder. BMC Bioinformatics 2006, 7: 319. 10.1186/1471-2105-7-319

    Article  PubMed Central  PubMed  Google Scholar 

  31. Su CT, Chen CY, Hsu CM: iPDA: integrated protein disorder analyzer. Nucleic Acids Res 2007, 35: 465–472. 10.1093/nar/gkm353

    Article  Google Scholar 

  32. Yang MQ, Yang JY: IUP: intrinsically unstructured protein predictor-a software tool for analyzing polypeptide sequences. Sixth IEEE Symposium on BioInformatics and BioEngineering: 16–18 October 2006; Arlington, Virginia, USA 2006, 3–11.

    Chapter  Google Scholar 

  33. Schlessinger A, Liu J, Rost B: Natively unstructured loops differ from other loops. PLoS Comput Biol 2007, 3: e140. 10.1371/journal.pcbi.0030140

    Article  PubMed Central  PubMed  Google Scholar 

  34. Wang L, Sauer UH: OnD-CRF: predicting order and disorder in proteins using conditional random fields. Bioinformatics 2008, 24: 1401–1402. 10.1093/bioinformatics/btn132

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Cheng J, Sweredoski M, Baldi P: Accurate prediction of protein disordered regions by mining protein structure data. Data Mining Knowl Disc 2005, 11: 213–222. 10.1007/s10618-005-0001-y

    Article  Google Scholar 

  36. Deng X, Eickholt J, Cheng J: PreDisorder: Ab initio sequence-based prediction of protein disordered regions. BMC Bioinformatics 2009, 10: 436. 10.1186/1471-2105-10-436

    Article  PubMed Central  PubMed  Google Scholar 

  37. Ishida T, Kinoshita K: Prediction of disordered regions in proteins based on the meta approach. Bioinformatics 2008, 24: 1344–1348. 10.1093/bioinformatics/btn195

    Article  CAS  PubMed  Google Scholar 

  38. Schlessinger A, Punta M, Yachdav G, et al.: Improved disorder prediction by combination of orthogonal approaches. PLoS One 2009, 4: e4433. 10.1371/journal.pone.0004433

    Article  PubMed Central  PubMed  Google Scholar 

  39. Mizianty MJ, Stach W, Chen K, Kedarisetti KD, Disfani F, Kurgan L: 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  PubMed Central  CAS  PubMed  Google Scholar 

  40. Xue B, Dunbrack RL, Williams RW, Dunker AK, Uversky VN: PONDR-FIT: a meta-predictor of intrinsically disordered amino acids. Biochim Biophys Acta 2010, 1804(4):996–1010.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  41. McGuffin LJ: 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 

  42. Dunker AK, Lawson JD, Brown CJ, Williams RM, Romero P, Oh JS, Oldfield CJ, Campen AM, Ratliff CM, Hipps KW, et al.: Intrinsically disordered protein. J Mol Graph Model 2001, 19: 26–59. 10.1016/S1093-3263(00)00138-8

    Article  CAS  PubMed  Google Scholar 

  43. Oldfield CJ, Cheng Y, Cortese MS, Brown CJ, Uversky VN, Dunker AK: Comparing and combining predictors of mostly disordered proteins. Biochemistry 2005, 44: 1989–2000. 10.1021/bi047993o

    Article  CAS  PubMed  Google Scholar 

  44. Xue B, Oldfield CJ, Dunker AK, Uversky VN: CDF it all: consensus prediction of intrinsically disordered proteins based on various cumulative distribution functions. FEBS Lett 2009, 583(9):1469–1474. 10.1016/j.febslet.2009.03.070

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Romero P, Obradovic Z, Kissinger CR, Villafranca JE, Garner E, Guilliot S, Dunker AK: Thousands of proteins likely to have long disordered regions. Proceedings of the Pac Symp Biocomput.: 4–9 January 1998; Hawaii 1998, 437–448.

    Google Scholar 

  46. Le Gall T, Romero P, Cortese MS, Uversky VN, Dunker AK: Intrinsic disorder in the Protein Data Bank. J. Biomol. Struct. Dyn 2007, 24(4):303–428.

    Article  Google Scholar 

  47. Haynes C, Ji F, Oldfield CJ, Klitgord N, Cusick ME, Radivojac P, Uversky VN, Vidal M, Iakoucheva LM: Intrinsic disorder is a common feature of hub proteins from four eukaryotic interactomes. PLoS Comput Biol 2006, 2(8):e100. 10.1371/journal.pcbi.0020100

    Article  PubMed Central  PubMed  Google Scholar 

  48. Liu J, Perumal NB, Oldfield CJ, Su EW, Uversky VN, Dunker AK: Intrinsic disorder in transcription factors. Biochemistry 2006, 45(22):6773–6888. 10.1021/bi0523815

    Article  Google Scholar 

  49. Uversky VN, Roman A, Oldfield CJ, Dunker AK: Protein intrinsic disorder and human papillomaviruses: Increased amount of disorder in E6 and E7 oncoproteins from high risk HPVs. J Proteome Res 2006, 5(8):1829–1842. 10.1021/pr0602388

    Article  CAS  PubMed  Google Scholar 

  50. Dosztányi Z, Chen J, Dunker AK, Simon I, Tompa P: Disorder and sequence repeats in hub proteins and their implications for network evolution. J Proteome Res 2006, 5(11):2985–2995. 10.1021/pr060171o

    Article  PubMed  Google Scholar 

  51. Goh GK-M, Dunker AK, Uversky VN: A comparative analysis of viral matrix proteins using disorder predictors. Virology J 2008, 5: 126. 10.1186/1743-422X-5-126

    Article  Google Scholar 

  52. Cortese MS, Uversky VN, Dunker AK: Intrinsic disorder in scaffold proteins: Getting more from less. Progress Bioph Mol Biol 2008, 98(1):85–106. 10.1016/j.pbiomolbio.2008.05.007

    Article  CAS  Google Scholar 

  53. De Biasio A, Guarnaccia C, Popovic M, Uversky VN, Pintar P, Pongor S: Prevalence of intrinsic disorder in the intracellular region of human single-pass type I proteins: The case of the Notch ligand Delta-4. J Proteome Res 2008, 7(6):2496–2506. 10.1021/pr800063u

    Article  CAS  PubMed  Google Scholar 

  54. Hébrard E, Bessin Y, Michon T, Longhi S, Uversky VN, Delalande F, Van Dorsselaer A, Romero P, Walter J, Declerk N, et al.: Intrinsic disorder in viral proteins genome-linked: Experimental and predictive analyses. Virology J 2009, 6: 23. 10.1186/1743-422X-6-23

    Article  Google Scholar 

  55. Balázs A, Csizmok V, Buday L, Rakács M, Kiss R, Bokor M, Udupa R, Tompa K, Tompa P: High levels of structural disorder in scaffold proteins as exemplified by a novel neuronal protein, CASK-interactive protein1. FEBS J 2009, 276(14):3744–3756. 10.1111/j.1742-4658.2009.07090.x

    Article  PubMed  Google Scholar 

  56. Hegyi H, Buday L, Tompa P: Intrinsic structural disorder confers cellular viability on oncogenic fusion proteins. PLoS Comput Biol 2009, 5(10):e1000552. 10.1371/journal.pcbi.1000552

    Article  PubMed Central  PubMed  Google Scholar 

  57. Tompa P, Kovacs D: Intrinsically disordered chaperones in plants and animals. Biochem Cell Biol 2010, 88(2):167–174. 10.1139/O09-163

    Article  CAS  PubMed  Google Scholar 

  58. Xue B, Williams RW, Oldfield CJ, Goh GK-M, Dunker AK, Uversky VN: Viral disorder or disordered viruses: Do viral proteins possess unique features? Prot. Pept. Lett 2010, 17(8):932–951. 10.2174/092986610791498984

    Article  CAS  Google Scholar 

  59. Dunker AK, Obradovic Z, Romero P, Garner EC, Brown CJ: Intrinsic protein disorder in complete genomes. Genome Inform 2000, 11: 161–171.

    CAS  Google Scholar 

  60. Ward JJ, Sodhi JS, McGuffin LJ, Buxton BF, Jones DT: Prediction and functional analysis of native disorder in proteins from the three kingdoms of life. J Mol Biol 2004, 337: 635–645. 10.1016/j.jmb.2004.02.002

    Article  CAS  PubMed  Google Scholar 

  61. Tompa P, Dosztanyi Z, Simon I: Prevalent structural disorder in E. coli and S. cerevisiae proteomes. J Proteome Res 2006, 5(8):1996–2000. 10.1021/pr0600881

    Article  CAS  PubMed  Google Scholar 

  62. Xue B, Williams RW, Oldfield CJ, Dunker AK, Uversky VN: Archaic chaos: Intrinsically disordered proteins in Archaea. BMC Systems Biol 2010, 4(Suppl 1):S1. 10.1186/1752-0509-4-S1-S1

    Article  Google Scholar 

  63. Xie H, Vucetic S, Iakoucheva LM, Oldfield CJ, Dunker AK, Obradovic Z, Uversky VN: Functional anthology of intrinsic disorder. 3. Ligands, post-translational modifications, and diseases associated with intrinsically disordered proteins. J Proteome Res 2007, 6: 1917–1932. 10.1021/pr060394e

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  64. Xie H, Vucetic S, Iakoucheva LM, Oldfield CJ, Dunker AK, Uversky VN, Obradovic Z: Functional anthology of intrinsic disorder. 1. Biological processes and functions of proteins with long disordered regions. J Proteome Res 2007, 6: 1882–1898. 10.1021/pr060392u

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  65. Vucetic S, Xie H, Iakoucheva LM, Oldfield CJ, Dunker AK, Obradovic Z, Uversky VN: Functional anthology of intrinsic disorder. 2. Cellular components, domains, technical terms, developmental processes, and coding sequence diversities correlated with long disordered regions. J Proteome Res 2007, 6: 1899–1916. 10.1021/pr060393m

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  66. Vucetic S, Brown CJ, Dunker AK, Obradovic Z: Flavors of protein disorder. Proteins 2003, 52: 573–584. 10.1002/prot.10437

    Article  CAS  PubMed  Google Scholar 

  67. Williams RM, Obradovic Z, Mathura V, Braun W, Garner EC, Young J, Takayama S, Brown CJ, Dunker AK: The protein non-folding problem: amino acid determinants of intrinsic order and disorder. Proceedings of the Pac Symp Biocomput.:3–7 January 2001; Hawaii 2001, 89–100.

    Google Scholar 

  68. Uversky VN, Dunker AK: Understanding protein non-folding. Biochim. Biophys. Acta-Proteins and Proteomics 2010, 1804(6):1231–1264. 10.1016/j.bbapap.2010.01.017

    Article  CAS  Google Scholar 

  69. Radivojac P, Iakoucheva LM, Oldfield CJ, Obradovic Z, Uversky VN, Dunker AK: Intrinsic disorder and functional proteomics. Biophys J 2007, 92: 1439–1456. 10.1529/biophysj.106.094045

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  70. Vacic V, Uversky VN, Dunker AK, Lonardi S: Composition Profiler: a tool for discovery and visualization of amino acid composition differences. BMC Bioinformatics 2007, 8: 211. 10.1186/1471-2105-8-211

    Article  PubMed Central  PubMed  Google Scholar 

  71. Radivojac P, Obradovic Z, Smith DK, Zhu G, Vucetic S, Brown CJ, Lawson JD, Dunker AK: Protein flexibility and intrinsic disorder. Protein Sci 2004, 13: 71–80. 10.1110/ps.03128904

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  72. Zhang H, Zhang T, Chen K, Shen S, Ruan J, Kurgan L: On the relation between residue flexibility and local solvent accessibility in proteins. Proteins 2009, 76: 617–636. 10.1002/prot.22375

    Article  CAS  PubMed  Google Scholar 

  73. Lieutaud P, Canard B, Longhi S: MeDor: a metaserver for predicting protein disorder. BMC Genomics 2008, 9(Suppl 2):S25. 10.1186/1471-2164-9-S2-S25

    Article  PubMed Central  PubMed  Google Scholar 

  74. Bordoli L, Kiefer F, Schwede T: Assessment of disorder predictions in CASP7. Proteins 2007, 69(Suppl 8):129–136.

    Article  CAS  PubMed  Google Scholar 

  75. Noivirt-Brik O, Prilusky J, Sussman J: Assessment of disorder predictions in CASP8. Proteins 2009, 77(Suppl 9):210–216.

    Article  CAS  PubMed  Google Scholar 

  76. Sickmeier M, Hamilton JA, LeGall T, Vacic V, Cortese MS, Tantos A, Szabo B, Tompa P, Chen J, Uversky VN, et al.: DisProt: the database of disordered proteins. Nucleic Acids Res 2007, 35: D786–793. 10.1093/nar/gkl893

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  77. Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res 2000, 28: 235–242. 10.1093/nar/28.1.235

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  78. Wang G, Dunbrack RL Jr: PISCES: a protein sequence culling server. Bioinformatics 2003, 19: 1589–1591. 10.1093/bioinformatics/btg224

    Article  CAS  PubMed  Google Scholar 

  79. Sirota FL, Ooi HS, Gattermayer T, Schneider G, Eisenhaber F, Maurer-Stroh S: Parameterization of disorder predictors for large-scale applications requiring high specificity by using an extended benchmark dataset. BMC Genomics 2010, 11(Suppl 1):S15. 10.1186/1471-2164-11-S1-S15

    Article  PubMed Central  PubMed  Google Scholar 

  80. Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: 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  PubMed Central  CAS  PubMed  Google Scholar 

  81. Jones DT, Swindells MB: Getting the most from PSI-BLAST. Trends Biochem Sci 2002, 27: 161–164. 10.1016/S0968-0004(01)02039-4

    Article  CAS  PubMed  Google Scholar 

  82. McGuffin LJ, Bryson K, Jones DT: The PSIPRED protein structure prediction server. Bioinformatics 2000, 16: 404–405. 10.1093/bioinformatics/16.4.404

    Article  CAS  PubMed  Google Scholar 

  83. Faraggi E, Xue B, Zhou Y: Improving the prediction accuracy of residue solvent accessibility and real-value backbone torsion angles of proteins by fast guided-learning through a two-layer neural network. Proteins 2009, 74: 857–871. 10.1002/prot.22194

    Article  Google Scholar 

  84. Dor O, Zhou Y: Real-SPINE: an integrated system of neural networks for real-value prediction of protein structural properties. Proteins 2007, 68: 76–81. 10.1002/prot.21408

    Article  CAS  PubMed  Google Scholar 

  85. Plewczynski D, Slabinski L, Ginalski K, Rychlewski L: Prediction of signal peptides in protein sequences by neural networks. Acta Biochim Pol 2008, 55: 261–267.

    CAS  PubMed  Google Scholar 

  86. Shevade SK, Keerthi SS, Bhattacharyya C, Murthy KRK: Improvements to the SMO algorithm for SVM regression. IEEE Trans. Neural Networks 2000, 11(5):1188–1193. 10.1109/72.870050

    Article  CAS  PubMed  Google Scholar 

  87. Xu H, Yang L, Freitas MA: A robust linear regression based algorithm for automated evaluation of peptide identifications from shotgun proteomics by use of reversed-phase liquid chromatography retention time. BMC Bioinformatics 2008, 9: 347. 10.1186/1471-2105-9-347

    Article  PubMed Central  PubMed  Google Scholar 

  88. Gao J, Zhang T, Zhang H, Shen S, Ruan J, Kurgan L: Accurate prediction of protein folding rates from sequence and sequence-derived residue flexibility and solvent accessibility. Proteins 2010, 78(9):2114–2130.

    CAS  PubMed  Google Scholar 

  89. Jiang Y, Iglinski P, Kurgan L: Prediction of protein folding rates from primary sequences using hybrid sequence representation. J Comput Chem 2009, 30(5):772–83. 10.1002/jcc.21096

    Article  CAS  PubMed  Google Scholar 

  90. Wagner M, Adamczak R, Porollo A, Meller J: Linear regression models for solvent accessibility prediction in proteins. J Comput Biol 2005, 12(3):355–369. 10.1089/cmb.2005.12.355

    Article  CAS  PubMed  Google Scholar 

  91. Homaeian L, Kurgan L, Ruan J, Cios KJ, Chen K: Prediction of protein secondary structure content for the twilight zone sequences. Proteins 2007, 69(3):486–498. 10.1002/prot.21527

    Article  CAS  PubMed  Google Scholar 

  92. Sotriffer CA, Sanschagrin P, Matter H, Klebe G: SFCscore: scoring functions for affinity prediction of protein-ligand complexes. Proteins 2008, 73(2):395–419. 10.1002/prot.22058

    Article  CAS  PubMed  Google Scholar 

  93. Pan XY, Shen HB: Robust prediction of B-factor profile from sequence using two-stage SVR based on random forest feature selection. Protein Pept Lett 2009, 16(12):1447–1454. 10.2174/092986609789839250

    Article  CAS  PubMed  Google Scholar 

  94. Chang DT, Huang HY, Syu YT, Wu CP: Real value prediction of protein solvent accessibility using enhanced PSSM features. BMC Bioinformatics 2008, 9(Suppl 12):S12. 10.1186/1471-2105-9-S12-S12

    Article  PubMed Central  PubMed  Google Scholar 

  95. Song J, Tan H, Takemoto K, Akutsu T: HSEpred: predict half-sphere exposure from protein sequences. Bioinformatics 2008, 24(13):1489–1497. 10.1093/bioinformatics/btn222

    Article  CAS  PubMed  Google Scholar 

  96. Hall M, Frank E, Holmes G, Pfahringer B, Reutemann P, Witten IH: The WEKA Data Mining Software: an update. SIGKDD Explor 2009, 11: 10–18. 10.1145/1656274.1656278

    Article  Google Scholar 

  97. Uestuen B, Melssen WJ, Buydens LMC: Facilitating the application of Support Vector Regression by using a universal Pearson VII function based kernel. Chemometrics Intel. Lab. Sys 2006, 81: 29–40. 10.1016/j.chemolab.2005.09.003

    Article  CAS  Google Scholar 

  98. Hymowitz SG, O'Connell MP, Ultsch MH, Hurst A, Totpal K, Ashkenazi A, de Vos AM, Kelley RF: A unique zinc-binding site revealed by a high-resolution X-ray structure of homotrimeric Apo2L/TRAIL. Biochemistry 2000, 39(4):633–640. 10.1021/bi992242l

    Article  CAS  PubMed  Google Scholar 

  99. Whitby FG, Luecke H, Kuhn P, Somoza JR, Huete-Perez JA, Phillips JD, Hill CP, Fletterick RJ, Wang CC: Crystal structure of Tritrichomonas foetus inosine-5'-monophosphate dehydrogenase and the enzyme-product complex. Biochemistry 1997, 36(35):10666–10674. 10.1021/bi9708850

    Article  CAS  PubMed  Google Scholar 

  100. Brown CJ, Takayama S, Campen AM, Vise P, Marshall TW, Oldfield CJ, Williams CJ, Dunker AK: Evolutionary rate heterogeneity in proteins with long disordered regions. J Mol Evol 2002, 55(1):104–110. 10.1007/s00239-001-2309-6

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This work was sponsored in part by the Discovery grant from NSERC Canada to LK, the National Institutes of Health grant (R01 GM 085003) to YZ, and the Killam Memorial Scholarship to MJM. The funding agencies did not participate in the design, collection, analysis, and interpretation of the data.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Lukasz Kurgan.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

MJM contributed to the collection of the data, design of the prediction method, experimental validation, analysis and interpretation of the results, and implementation of the web server. TZ worked on the design and implementation of the prediction method and analysis and interpretation of the results. BX helped with the analysis and interpretation of the results. YZ and AKD helped with the conception of the project and analysis and interpretation of the results. VNU contributed to the analysis and interpretation of the results and writing of the manuscript. LK was responsible for the conception of the project, and contributed to the collection of the data, design and implementation of the prediction method, experimental validation, analysis and interpretation of the results, and writing of the manuscript. All authors have read, corrected and approved the manuscript.

Electronic supplementary material

12859_2010_4635_MOESM1_ESM.PDF

Additional file 1:Supplementary tables and figures. This file includes 4 supplementary tables and 3 supplementary figures. The tables summarize the input features and results obtained with alternative designs of the proposed predictor. The figures summarize correlation between the input features and the predictive target, the relation between the predictive quality and the input chain length, and the accuracy for the residue-level disorder predictions. (PDF 204 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://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

Mizianty, M.J., Zhang, T., Xue, B. et al. In-silico prediction of disorder content using hybrid sequence representation. BMC Bioinformatics 12, 245 (2011). https://doi.org/10.1186/1471-2105-12-245

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2105-12-245

Keywords