Predictor Of Naturally Disordered Regions

PONDR® Accuracy

Two types of PONDR® accuracies will be discussed. The first is the accuracy obtained during predictor training. The second is the accuracy observed when validating the predictor on out-of-sample data. If these two accuracies agree for all sets of out-of-sample data, the predictor is said to generalize well. As an example, the training of VL-XT is discussed here, since it is the most studied and characterized of the PONDR® algorithms.

There follows a brief discussion of improving accuracy as disorder prediction length increases. Finally a brief discussion of two sources of error is given, with a prominent example.

Predictor Training

During development of the algorithms, overtraining can be a problem. Overtraining of the NNP leads to memorization of the data by the NNP, rendering the algorithm useless. To avoid this, the following procedure was used. During training, only four-fifths of the data was used. The algorithm was periodically tested against the remaining 1/5 of the data. When the accuracy on the training set (the 4/5) matched that on the test set (the 1/5), the training was halted. This training was repeated 3 different times with different starting conditions. Then a different 1/5 of the data was held out as the test data and the process was started over. This was done a total of 5 times, yielding 15 predictors. (This method is called "5-cross validation".) The training accuracy of VL1 over these 15 cycles was 83 ± 3%. Finally, to generate one predictor to replace the 15 made during training, the data was merged and a new predictor was trained to the same level of accuracy (Romero et al., 2001).

In a similar manner, the training accurracies of the XN and XC predictors reached 75%, but varied by the position in the sequence (Li et al., 1999). The two single predictors, XN and XC, were developed by merging the data and training to the same levels of accuracy, just as above.

 

Accuracy as False Positive & False Negative

The generalized accuracy of the predictors may be assessed by looking at their false positive and false negative rates.
bullet

False positive means that the predictor indicates "disorder" when the residue or region is known to be ordered.

bullet

False negative means that the predictor indicates "order" when the residue or region is known to be disordered.

To measure these rates, the PONDR® algorithms were allowed to predict on two sets of out-of-sample data:
bullet

one nominally containing entirely ordered protein (for false positive rates)
bullet

Formed by deleting residues with missing backbone electron density from the set of proteins in PDB Select 25 as of August 3, 1999, giving O_PDB_S25. This set contains 1,111 proteins having 220,668 ordered residues. Since the starting set, PDB Select 25, contained one protein from each family in PDB, the O_PDB_S25 set spanned protein structure space as known at the time of construction.

bullet

one nominally containing entirely disordered protein (for false negative rates)
bullet

The disordered data was compiled from examples in the literature and contained 150 protein chains with a total of 17,417 disordered residues. The long disordered regions in this set were characterized by x-ray crystallography, by NMR, or by CD.

The results (on a per-residue basis) are summarized below.

PONDR® False
Negative
False
Positive
5-cross
Validation
VL-XT 40% 22% 75 - 83%
XL1 62% 19% 73±4%
CaN 39% 34% 83±5%

Of the three predictors, VL-XT is the most accurate. It has a 78% success rate on predicting order, in good agreement with the training accuracy of 75 - 83%. It has a 60% success rate on predicting disorder which suggests that it does not generalize well. However, note that both VL-XT and XL1 were trained to recognize long regions of disorder (> 39 residues, see below).

XL1 is as accurate as VL-XT on predicting order, but is worse at predicting disorder.

CaN does not generalize well at all, showing poor out of sample accuracy, but in some cases the contrast of its output with other predictors provides insight into binding regions of disordered sequences (Garner et al. 1999).

 

Accuracy over longer predicted disorder regions

It is important to realize that the above table shows accuracies on a per-residue basis. As prediction length (L) increases, the accuracy of the prediction increases dramatically, as the following table shows for VL-XT. (This particular algorithm is the most studied.) Note the sharp drop in the prediction errors as the length of the prediction increases. This is why we recommend focusing on long disorder predictions.

Minimum Number of Consecutive
Predictions of Disorder ( = L )
False Positive
per Residue
1 20.1%
10 15.4%
20 8.3%
30 3.9%
40 1.7%
50 0.5%
60 0.1%
False positive rates as a function of minimum disordered region length. This summary comes from predictions made on the August 1999 O_PDB_S25 set. These results differ slightly from the published data (Romero et al., 2001 and Dunker et al., 2001) because a different procedure was used for error calculation. The estimation used previously was "percentage of sequence windows, with length equal to L, that are falsely predicted to be disordered." In this table, the error estimate is "percentage of residues in predicted disordered regions of length greater than or equal to L that are falsely predicted to be disordered." Notice that the false positive rate decreases dramatically for increasing length of predicted disorder.

 

Sources of Error

There are several potential sources for prediction error. We will discuss two of them here. First, absolute assignment of residues in the training sets to be ordered or disordered is difficult, especially for the ones assigned to be disordered. Second, many protein regions in the PDB database classified as ordered are involved in protein-protein or protein-nucleic acid interactions or complexes (for example, see figure below). However, these same regions are disordered in the absence of their partners. By our definition, such regions should be classified as disordered, instead of as ordered in PDB. Both of these error sources lead to noisy data, which contributes to the observed error rates.

One example of this second source of error is the oligomerization domain of the tumor suppressor p53 (see the below figure). Even though this region is shown as ordered in PDB, our PONDR® algorithm predicts portions to be disordered. Specifically the C-terminal helix and a portion of the beta sheet immediately upstream of that helix are predicted to be disordered (the red portions in the below figure). Disorder in the monomers is one explanation for how this interweaving fold could occur when the monomers interact (undergoing a disorder-to-order transition).

Example of disorder/order prediction

Structure of tumor suppressor p53 oligomerization domain (PDB entry 3SAK). This representation includes only one pair from the tetramer present in the NMR structure. PONDR® predicted disorder and order is represented by red and blue ribbons, respectively. PONDR® predicts the carboxyl terminal helix and a portion of the beta sheet to be disordered.

A further discussion of sources of error (including mis-classification in both directions) can be found in Li et al., 2000.

 

Bibliography

Dunker AK, Lawson JD, Brown CJ, Williams RM, Romero P, Oh JS, Oldfield CJ, Campen AM, Ratliff CM, Hipps KW, Ausio J, Nissen MS, Reeves R, Kang C, Kissinger CR, Bailey RW, Griswold MD, Chiu W, Garner EC, and Obradovic Z. (2001) Intrinsically disordered protein, J. Mol. Graphics and Modeling, 19, 26-59.

Garner E, Guilliot S, Dunker AK, Romero P, and ObradovicZ. (1998) Predicting long disordered regions in protein from amino acid sequence, Biophysical Journal, 74, Part 2 of 2: A1 - A452.

Li X, Romero P, Rani M, Dunker AK, Obradovic Z. (1999) Predicting protein disorder for N-, C- and internal regions, Genome Informatics, 10, 30-40.

Li X, Obradovic Z, Brown CJ, Garner EC, Dunker AC. (2000) Comparing predictors of disordered protein, Genome Informatics, 11, 172-184.

Romero P, Obradovic Z, Li X, Garner E, Brown C, and Dunker AK. (2001) Sequence complexity of disordered protein, Proteins: Struct. Funct. Gen., 42, 38-48.