Table of Contents
ISRN Bioinformatics
Volume 2013, Article ID 671269, 11 pages
http://dx.doi.org/10.1155/2013/671269
Research Article

SUMOhunt: Combining Spatial Staging between Lysine and SUMO with Random Forests to Predict SUMOylation

National Institute of Biotechnology and Genetic Engineering, P.O. Box 577, Jhang Road, Faisalabad, Pakistan

Received 24 April 2013; Accepted 28 May 2013

Academic Editors: A. Lukas, A. Pulvirenti, J. T. L. Wang, and K. Yura

Copyright © 2013 Amna Ijaz. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Modification with SUMO protein has many key roles in eukaryotic systems which renders the identification of its target proteins and sites of considerable importance. Information regarding the SUMOylation of a protein may tell us about its subcellular localization, function, and spatial orientation. This modification occurs at particular and not all lysine residues in a given protein. In competition with biochemical means of modified-site recognition, computational methods are strong contenders in the prediction of SUMOylation-undergoing sites on proteins. In this research, physicochemical properties of amino acids retrieved from AAIndex, especially those involved in docking of modifier and target proteins and optimal presentation of target lysine, in combination with sequence information and random forest-based classifier presented in WEKA have been used to develop a prediction model, SUMOhunt, with statistics significantly better than all previous predictors. In this model 97.56% accuracy, 100% sensitivity, 94% specificity, and 0.95 MCC have been achieved which shows that proposed amino acid properties have a significant role in SUMO attachment. SUMOhunt will hence bring great reliability and efficiency in SUMOylation prediction.

1. Introduction

Posttranslational modifications on proteins offer spectacular diversity and functional variety to an organism’s otherwise constrained proteome. SUMOylation is one such PTM whose vast expanse of biological implications in organisms has brought it under attention; still till now many of its functional outcomes are not known. To name a few, SUMOylation is involved in transcriptional regulation [13], mRNA metabolism [4], apoptosis [5, 6], nuclear and subcellular transport [7, 8], protein trafficking [9], signal transduction [10], regulation of DNA damage and replication, cell-cycle progression, competition with other members of the ubiquitin family [2, 11, 12], prevention or promotion of deacetylation [13], chromosome segregation [14], structural integrity of chromatin and many proteins, and mitosis [15]. It has been reported to be involved in the perception of sound as well [16]. Also, it is known to participate in early developmental processes like cell differentiation, specification, division, and lineage commitment [17]. SUMOylation of a target protein can change its localization in a cell by altering its intermolecular and intramolecular interactions [18]. Hence, by determining whether a protein is SUMOylated or not, vital evidences can be gathered regarding its function and spatial association [19].

SUMO, a member of the ubiquitin family, is made up of 97 amino acids and is also called Smt3p, Pmt2p, PIC-1, GMP-1, Ubl1, and Sentrin [20]. It mostly modifies proteins present in the nucleus, cytoplasm, and sometimes the plasma membrane of eukaryotic cells [17]. SUMO proteins are highly conserved across eukaryotic proteomes. In mammals, including humans, there are four isoforms of SUMO called SUMO 1, SUMO 2, SUMO 3, and SUMO 4; in yeasts there is only one SUMO protein while plants produce at least eight SUMO isoforms [21]. SUMO 1, SUMO 2, and SUMO 3 are expressed throughout an organism’s body [2] with the latter two having greater sequence similarity as compared to SUMO 1 [22]. SUMO 4 dominates in lymph node, kidney, and spleen in mammals [23], having resemblance with SUMO 2 and 3 and a dominant occurrence in kidney [24].

SUMOylation can occur in either the cytoplasm or the nucleus depending on the locality of the target protein, though the modification may afterwards be responsible for regulation of production and change in its localization [2] mostly from other parts of cell to the nucleus. RanGAP1, which controls the transportation of ribonucleoproteins across the nuclear pore complex, was identified as the first SUMO target [25]. Before SUMOylation, this protein mainly resides in the cytoplasm and afterwards associates itself with the cytoplasmic fibers of the nuclear pore complex [25]. All essential components of SUMOylation are present at the nuclear pore complex which shows its involvement in nuclear import or even retention of any incoming proteins within the nucleus [26, 27]. However, TEL protein, modified at K99, dominates in the cytoplasm while TEL mutated and hence unable to be SUMOylated at this residue mostly resides in the nucleus, pointing towards possible role of SUMOylation in its export [28].

Since SUMOylation is a part of a very wide array of biological and cellular processes, even a minor dysfunction in the pathway can result in severe pathological conditions like cancer. Neurodegenerative diseases including Huntington’s disease [19], Alzheimer’s disease, Parkinson’s disease [29], and neuronal intranuclear inclusion disease [9] often have anomalies in the SUMOylation pathway involved in their onset. Also, type I diabetes [30] and familial amyotrophic sclerosis [13] have SUMOylation dysfunction as a frequent part of their beginning.

SUMO proteins attach covalently to lysine, which mostly lies in consensus motif ΨKXE [31], where Ψ is any hydrophobic amino acid, K is the target lysine residue, X is any of the twenty encoded amino acids, and E is glutamic acid, with the help of SUMO E2-conjugation enzyme, ubc9 [31]. In mammals, any one of the four isoforms of SUMO conjugates to target protein as an individual molecule or in the form of polymeric chains [32] at target lysine. The enzymes working the SUMOylation pathway are E1-activating enzyme, E2-conjugating enzyme, and E3 ligase. The SUMO pathway can be categorized into maturation, activation, conjugation, and finally ligation at target site (Figure 1). First of all, SUMO is processed by SUMO-specific proteases (SENPs) [33] to cleave its terminal into an exposed diglycine motif at the carboxyl end, converting it from immature to mature protein. The maturation process of all three major SUMO proteins is identical in mammals [34]. Secondly, activation of mature SUMO takes place through an ATP-dependent thioester bond formation between SUMO and UBA2/AOS1, a heterodimeric E1-activating enzyme [2]. Direct linkage occurs between SUMO and UBA2 subunit of the heterodimer. This enzyme is called SAE1/SAE2 or Sua1/hUba2 in humans [34]. The yeast homologue of UBA2 subunit is Uba2p, which together with Aos1 was the first activating enzyme to be discovered [35] while the human homologue is hUba2 [36]; both act in the same way. The third step in this pathway is catalyzed by E2-conjugating enzyme, Ubc9. Activated SUMO is transferred from E1-activating enzyme to a cysteine residue in Ubc9. In contrast to other enzymes involved in the process, Ubc9 is the only type of enzyme identified in its category [31, 37]. It is Ubc9 which identifies the consensus or nonconsensus sequence at target site for subsequent conjugation [38]. In the final step, SUMO is attached to target protein with isopeptide bond between the exposed diglycine on carboxyl terminal of SUMO and the -amino group of target lysine in protein with the help of E3-ligating enzyme. Three distinct characteristics of ligating enzymes in SUMOylation pathway have been sketched through research: (i) they should be able to directly or indirectly associate with the target protein, (ii) they should be able to bind with their preceding enzyme, ubc9, in the pathway, and (iii) they should be capable of transferring SUMO to target protein or another SUMO in case of poly-SUMOylation [34]. RanBP2 which is a nuclear pore protein [39], TOPORS [40, 41], PIAS proteins [42], Pc2 [43] which is polycomb group protein, and RNF4 [44] have all been identified to have E3 ligase activity in SUMOylation pathway.

671269.fig.001
Figure 1: The SUMOylation pathway: the SUMO pathway is categorized into maturation, activation, conjugation, and ligation. First of all, SENPs convert immature SUMO to mature SUMO by exposing diglycine motif at carboxyl terminal. This mature SUMO is activated by AOS1/UBA2 (E1—activating enzyme) with the help of ATP-dependent reaction. Activated SUMO conjugates with E2—conjugating enzyme ubc9—which finally transfers it to target lysine on substrate protein with the help of E3 ligase.

Nearly all research on SUMO attachment has pointed out the significance of ΨKXE motif, but on the other hand, SUMOylation has been reported to occur in regions outside of this consensus motif as well: in nonconsensus sites. Xu et al. [21] reported 26% SUMO occurrence in nonconsensus while Xue et al. [45] reported 23% such cases in bulk data used to develop their respective prediction servers. For example, there are four core histones, H2A, H2B, H3, and H4, that are frequently SUMOylated. However, none of the SUMOylated sites in these histones conform to the common consensus motif [46]. Several other consensus motifs have been proposed including NDSM and PDSM. NDSM proposes that negatively charged amino acids around the target lysine enhance SUMOylation [47] while PDSM is based on ΨKXEXXSP motif, both being only an extension of the original common motif [48].

An analysis of available PDB structures of several protein targets having 57 reported SUMOylation sites revealed 54 of the sites to be exposed on the surface while only 3 were buried within the proteins’ globular structure [21]. Also, research has led to the conclusion that SUMOylation is greatly enhanced when the target lysine is forced to adopt a favorable conformation [49]. From this behavior it can be inferred that the conjugation enzymes and ligases have sequence preference since they come into direct contact with target protein. SUMOylation pathway requires only three enzymes, specificity of subcellular localization and appropriate presentation of target residue on globular structure [2]. In case of other PTMs, a variety of enzymes with their target recognition and modification systems bring out varied site preference; typically they are not focused on any one type of residue. The major role in SUMOylation is played by only few enzymes discussed above; it suffers from the lack of efficient target recognition and modification systems, thus emphasizing the importance of motif and sequence information on the target protein as a device of recognition in the pathway.

Including sequence information as a principal contributor of computational prediction performance can provide rational computational tools, but focusing entirely on the consensus motif is not preferable as it can result in missing many true positives that lie in nonconsensus regions along with high false positive rate due to the many consensus sites that are actually not SUMOylated. In view of the occurrence of SUMO at both consensus and nonconsensus, sites it is proposed that other than the raw amino acid sequence around a SUMOylatable lysine residue there are factors of appropriate presentation and exposure of lysine and adjacent residues including steric hindrance, hydrophobicity, polarity, and entropy, playing a crucial role in determining whether a residue shall undergo SUMOylation or not. However, two residues downstream and one residue upstream of the target lysine play the most important role in SUMOylation [2, 50, 51].

We have employed different peptide lengths centered on lysine residues experimentally proved to undergo SUMOylation with sixteen amino acid properties (Table 1) from AAIndex to develop a prediction model named SUMOhunt. These amino acid properties were chosen on the basis of their contribution in increasing structural complementarity and association between incoming SUMO and target protein. Promising accuracy measures obtained on SUMOhunt developed by combining these properties with sequence information in the vicinity of target lysine and random-forest based algorithm presented in data mining software WEKA [52] have opened new paths for the development of an efficient prediction method.

tab1
Table 1: Physicochemical properties used to shape SUMOhunt. They were chosen because of their potential contribution in docking of SUMO with substrate protein at the target site.

2. Materials and Methods

2.1. Dataset

452 modified lysine instances (positive instances) were obtained from dbPTM [53], from training sets and sites supplemented for SUMOpre [21] and SUMOsp [45] and in the publication titled SUMO targets and Site Prediction: Combining Pattern Recognition and Phylogenetic Conservation by Xue et al. [54]. Primary sequences around these residues were retrieved from UniprotKB [55] in the form of 21-mer peptides: SUMOylated lysine residue had 10 residues upstream and 10 residues downstream of it. Standard single letter code was used for every amino acid residue.

Removal of redundant information gave 293 modified residues within 181 proteins (S1). The remaining 7346 lysine residues in these proteins, that is, lysine residues not reported to be modified, were assumed to be unmodified (negative instances) from a total lysine count of 7639 as calculated from MAPRes [56], and 21-mer peptides were generated for them as well. Such cases that were reported as both modified and unmodified were considered modified.

It is better to use numerical data for the development of a predictor, so every amino acid was encoded furthermore with the coefficients given in numeric matrices of chosen amino acid properties (Table 1) in the AAIndex [57]. Each peptide could be represented by a 16*21 dimension feature vector; that is, it could have 336 possible feature dimensions vectors.

This data was made workable with WEKA [52] by converting all the information into CSV format followed by conversion to ARFF.

2.2. Training the Prediction Program

In WEKA [52], random forest algorithm was trained using the ARFF generated. With default settings, ten trees were trained to vote for the class of each given instance. The random forest was trained using a dataset that had equivalent amount of modified and unmodified instances; the unmodified instances were randomly selected from the larger bulk of 7346 sites.

The choice of 21-mer peptide length around the target lysine site was tested against the same sites being present in 11-mer and 7-mer peptides as well. Each peptide in the latter categories could be represented in 16*11 = 176 and 16*7 = 112 dimension feature vectors, respectively. Hence, three datasets made in exactly the same manner were obtained. The datasets were not “formally” divided into training or test sets for accuracy measurement; instead percentage split in WEKA [52] was employed for this purpose to optimize the size of train and test. The dictated percentage of the original data including the most ideal-to-train instances is extracted and used as training data while the rest is used to test (S2 and S3).

Accuracy measures were calculated using

2.3. Cross-Validation and Evaluation

To test our predictor’s power, we used the methods of Xu et al. [21], that is, self-consistency test, k-fold cross-validation, and jack knife (leave-one-out cross-validation). Specificity, sensitivity, accuracy, and correlation coefficient for these tests were computed (Figures 2 and 3).

671269.fig.002
Figure 2: Comparison of accuracy measures (-axis) at several -fold cross-validations, LOOCV and self-consistency test (-axis). In -fold cross-validation, after dividing dataset into sections, typical training procedure is conducted using sections while the remaining one is used as test. This is repeated times until every set has been used as test exactly once. For eight -fold cross-validations performed, the value of was kept from 3 to 10 having average AC at 81%, SN at 83%, SP at 79%, and MCC at 0.66. LOOCV is a type of -fold cross-validation in which is equal to the number of total instances. It has the similar average as -fold. Self-consistency is a type of test in which prediction of every instance is done using the rules of the training dataset itself. This is done for each and every instance in the training dataset. In this test, 100% result is achieved for all accuracy measures.
fig3
Figure 3: Area under the ROC (-axis: true positive rate (fraction of true positives out of total positives) and -axis: false positive rate (fraction of false positives out of total negatives)) of 10-fold cross-validation, self-consistency test, and LOOCV-ROC curve depicts the performance of a classifier by plotting true positive rate versus the false positive rate. The greater the area under the curve, higher is the performance of a classifier. For 10-fold cross-validation 89% area under the curve is obtained, for LOOCV it is 88%, while for self-consistency it is a full 100%. At the end, ROC for classification of test set has been done with an AUC of 1.

Self-consistency test: it is predicted whether a given instance is positive or negative using the rules of the training dataset itself. This is done for each and every instance in the training dataset.

-fold cross-validation: here the dataset is randomly divided into sections. Typical training procedure is conducted using sections while the remaining one is used as test. This is repeated times until every set has been used as test exactly once.

Jack knife cross-validation: this is also called leave-one-out cross-validation and is an extension of the -fold cross-validation, having equal to the exact number of instances in the dataset.

2.4. Algorithm

The random forest algorithm [58] in WEKA [52] is built on decision tree classification. The said number of decision trees is generated with each tree having paths and nodes. Every node then uses rules derived from patterns in the data to decide between two or more paths. A given instance is classified by the last rule. To develop and grow the decision trees, a random selection of inputs and features is done at each node.

Voting on the class for a given instance is then carried out by the trees. Significant increase in classification accuracies have been observed if assortment of classifier trees is used and allowed to vote for the most popular class. Often random vectors are generated which govern the growth of each tree in the assortment. To see how a random forest actually shapes its model we should know that inherently for the th tree a random vector is generated, independent of all the previous vectors produced, but it has the same distribution. Now a tree is grown using the training set provided and , which makes a classifier ; here is the input vector. The nature and dimension of depends on its use in the construction of the tree. After a large number of trees have been generated, they then vote for the most popular class.

It is difficult to interpret models developed through random forest [58]; but there are certain features which make it suitable for the prediction of PTMs: a mixture of discrete and continuous descriptors, binary, or multiclass data can be proficiently treated with random forest algorithm. This algorithm is successful even when there is a lot of disorder in the data [59].

3. Results

3.1. Frequency Analysis

The number of different types of amino acid residues prevalent around the modified sites was analyzed (Table 2) with all results substantiated through corresponding frequency plot (Figure 4) of the same dataset. Results confirm the prevalence of the ΨKXE motif; approximately only 24% modified sites in this research lack the consensus motif. At −1 position or in place of Ψ, it is not just any hydrophobic amino acid; the data based on 293 modified instances has a significantly higher occurrence of hydrophobic amino acids with aliphatic side chains: Val (73 occurrences), Leu (45), and Ile (92) as compared to those with aromatic side chains: Phe (9), Tyr (2), Trp (0), and His (1) and other hydrophobic residues. In place of X or at +1 position a dominance of polar amino acids including Glu (33), Gln (43), and Thr (32) is there. The bulky aromatic: Trp (3), Tyr (4), and His (4) and small-size amino acids: Asn (6) and Cys (3) are in a significantly lesser proportion. At +2 position, polar residues Glu and Asp are prevalent with 224 and 13 existences, respectively.

tab2
Table 2: Frequency analysis (occurrence in numbers) of all amino acid residues around modified target lysine in experimentally proved positive dataset of 293 instances.
fig4
Figure 4: Frequency plot developed for all modified and unmodified instances with ten residues upstream and downstream of the target lysine; the plots are quite different from each other as expected. In the case of modified sites, there is a dominance of the conventional ΨKXE motif. Ψ is hydrophobic amino acid. Results show that aliphatic hydrophobic residues including valine, leucine, and isoleucine are dominant as compared to aromatic ones. K is the target lysine residue. X is any amino acid, which the results have shown it to be predominantly polar like glutamic acid, glutamine, and threonine. E is glutamic acid. In the dataset used here, about 76% sites follow the consensus motif.

The incidence of Trp and Cys is strikingly low at all twenty positions around the target lysine, the highest being only 11 for Cys at position +10. The other three aromatic amino acids Phe, Tyr, and His also do not prevail around the target lysine, supporting the imperative role of a catalytically favorable presentation of lysine in its SUMOylation. Aromatic molecules could potentially compromise this presentation leading towards the unavailability of target lysine to incoming SUMO.

3.2. Optimum Window Length for Prediction

In order to derive a good prediction model, the optimum window length of the peptides used for its training has to be determined. As shown in Figure 5, values for all four accuracy measures are highest for peptides of window size 7 as compared to others. For 93% split (discussed in Section 2.2) of original dataset to train (7% is test) at window size 21 lowest accuracy of 87.8% is encountered in comparison to 92.7% at 11 window size and 97.56% at 7 window size. MCC and sensitivity show significant ascensions: from 75% for 21-mer peptide to 95% for 7-mer peptide and from 85% for 21-mer to 100% for 7-mer, respectively. Specificity does not show any major improvement or decrease among the three peptide lengths. Interestingly, an evaluation conducted for the same accuracy measures against window size during development of SUMOpre [21] also resulted in the same window size being chosen as the best length. Hence, with all the information and comparisons, window size 7, which has 3 residues downstream and 3 residues upstream of the target lysine residue, is the most rational choice for our prediction model.

671269.fig.005
Figure 5: Comparison of accuracy measures at three window sizes: 21, 11, and 7 residue peptides; a good prediction model is dependent on optimum window size (-axis) around target lysine. Window size 7 having 3 residues upstream and downstream is the best with highest AC, MCC, SP, and SN (-axis) while window size 21 with 10 residues upstream and downstream has the lowest efficiency.
3.3. Prediction Accuracy and Stability of Model

Prediction accuracy was measured using varied sequence lengths, with and without physicochemical properties. The addition of information on physicochemical properties of amino acid residues produced powerful accuracies given in Table 3. From the several models generated, the best with the chosen window size was further substantiated by testing its stability through three procedures discussed in Section 2.3. In Figure 2, accuracy measures for several -fold cross-validations, jack-knife and self-consistency tests have been visualized. Results of the former two types of tests were all of nearly equal values with small deviations from mean while the latter varied significantly towards greater robustness; ROC for these tests can be observed in Figure 3. Average values of -fold tests and LOOCV for AC remained at 82%, SN at 86%, SP at 80%, and MCC at 0.66. The last type of test, self-consistency, gave a 100% result as shown in Figure 2.

tab3
Table 3: Accuracy measures on 21-mer, 11-mer, and 7-mer peptides containing target lysine with different combinations of data divided into training and test sets.

4. Discussion

The large range of biological processes and localizations populated by SUMO targets presents a great motivation to unravel information regarding SUMOylation and its targets by all possible means. Rigorous wet-lab experiments are frequently undertaken to isolate, identify, quantify, and report SUMOylation. SUMOylated proteins have been identified in yeast strains [51, 60] using mass spectrometry, chromatography proteolytic digestion, and so forth, on a trial and error basis to find target lysine residues [61]. On larger scale, proteomes having larger and more complex proteins with many SUMOylatable lysine residues are also analyzed through mutational analysis. However, these randomly executed experiments not only take up significant amount of time but also consume physical and chemical toils that are often futile as they result in discovering lysine residues that do not undergo SUMOylation. These approaches focus on identifying substrates rather than exact sites. Computational prediction of target sites has become mandatory before conducting experiments; this enables researchers to directly focus on residues which are potential candidates of SUMOylation. This dry-lab testing prior to corresponding wet-lab experimentation has gained much attention due to its cost effectiveness and power in proteomic data mining.

Till now, nine prediction models for SUMOylation have been proposed, out of which six have reported servers including SUMOplot [62] (web server), SUMOsp 1.0 and 2.0 [45] (downloadable), SUMOpre [21] (web server), Boshu Liu’s PSFS method [63] (web server), SSPFS [51] (available upon request), and seeSUMO [64] (web server).

SUMOplot [62] was the first step in development of computational server for the prediction of SUMOylation sites but had a bias of ΨKXE sites in data. SUMOsp, 2006 [45], was presented as SUMOsp 1.0 and SUMOsp 2.0 that were generated using GPS and MotifX (originally developed for phosphorylation) with sequence information only like SUMOplot [62]. Liu et al. developed a PSFS-based prediction model in 2007 [63] trained on hundreds of amino acid properties from which seven were selected as relevant with the help of sequential forward selection. SUMOpre, 2008 [21], was developed on the basis of sequence data only with probabilistic model of prediction. Its accuracy of 97.7% is impressive but at the cost of its sensitivity which is only 73.96%. SUMOsvm, 2008 [19], was developed using support vector machines trained on sequence information, solvent accessibility, secondary structure, and evolutionary profiles. FindSUMO, 2008 [65], was soon after developed by the PSSM system and with little progress in prediction efficiency. SSPFS, 2009 [51], was developed using mRMR and nearest neighbor algorithm trained on seven optimal amino acid properties selected from hundreds of amino acid properties. SUMOtr, 2010 [23], introduced the use of hydrophobicity, 3D structure, protein volume, and sequence to shape a model through tree classification algorithms. Recently, seeSUMO, 2011 [64], was introduced as a web server using random forest-based algorithm for training, but due to unavailability of its full publication, comparison and information about it have not been included here.

Physicochemical properties including hydrophobicity, buriability, isoelectric point, hydrophilicity, polarity, bulkiness, and molecular weight of residues control the spatial flexibility of target residue and hence can be very important in site attachment by developing complementarity between SUMO, enzymes in the pathway, and the target itself. Here, a computational system (S4) for the prediction of SUMOylation and investigation of its dependence on proposed properties (Table 1) has been developed using random forest-based classifier provided in WEKA [52]. Different programs developed for SUMOylation prediction are not directly comparable as they were developed using not only different datasets but also varied cross-validations and methods. Hence, the MCC, which is designed to assess predictive values for models from classes of different sizes, should be considered the primary measure for the purpose. It has a value between +1 and −1, with +1 being the highest level of prediction, 0 being the average level, and −1 showing inverse prediction. Table 4 shows the reported AC, SN, SP, MCC, and AUC with their proposed training and evaluation procedures.

tab4
Table 4: Comparison of SUMOhunt with previous models/servers.

In comparison to previous prediction methods that employed highly unbalanced datasets, equivalent amounts of modified and unmodified sites were used to train our model. Previous research (Table 4) has used uneven data with a very large part of the whole dataset being based on unmodified sites. Large number of unmodified instances yields a high specificity by making the correct prediction of nearly all such sites possible, at the cost of lowered sensitivity. Unbalanced dataset reduces performance and reliability so we have made a balanced set by sampling modified and unmodified sites in equal ratios. Moreover, approximately 24% of validated SUMOylation sites do not conform to the consensus motif and their representative peptides have been used in training of this model, and hence, the specificity of our prediction model is for both consensus and nonconsensus sequences, reducing the difficulty faced in prediction of the latter ones. These considerations ensure a uniformity of predictability when the predictor comes across either known or unknown sequences.

Prediction accuracy achieved in this research is significantly higher than all other prediction methods except that of SUMOpre [21] and SUMOsvm [19], to which it is fairly equal. However, our prediction model has exceeded all predictors with its high MCC, sensitivity, and AUC. Particularly, the sensitivity (correct classification of SUMOylated samples) and MCC (the measure of the overall performance of biased datasets) are higher than all others. On specificity (correct classification of non-SUMOylated samples) measures it retains a similar position even after the usage of an equivalent amount of modified and unmodified sites for training.

Based on the prediction performance, we believe that SUMOhunt can very well be used to implement a prediction server in the future that can assist as powerful and complementary tool for SUMOylation site identification and that the model will be available for the purpose. Correct analysis of this modification in all proteomes can greatly enhance our knowledge of the mechanism and working of many biological systems.

5. Conclusion

Generating reliable tools for identification of target sites of SUMOylation presents a great challenge. Computational methods of estimation can never replace experimental methods but can be of invaluable support to quicken and focus experimentation. In this research, based on experimental data, a prediction model has been developed that assures of robust computational method for highly accurate and specific SUMOylation-site prediction. Moreover, the physicochemical properties proposed to be playing crucial role in the appropriate presentation and hence rapid SUMOylation of target lysine have brought significant improvement in accuracy measures. This opens new paths to future work in analyzing the effect of these amino acid properties experimentally. Most importantly, based on the model presented here, it gives the possibility of building a server for prediction of SUMOylation sites in relation to the spatial properties of amino acid around them and sequence information.

Abbreviations

SUMO:Small ubiquitin-like modifier
NDSM: Negative-charge dependent SUMOylation motif
PDSM: Phosphorylation dependent SUMOylation motif
WEKA: Waikato environment of knowledge analysis
ROC: Receiver operator characteristic
LOOCV: Leave-one-out cross-validation
AC: Accuracy
SN: Sensitivity
SP: Specificity
MAPRes: Mining associationpatterns between preferred residues
MCC: Matthews correlation coefficient
AUC: Area under the curve
PSFS: Properties sequence forward selection
SSPFS: Sumoylation site prediction based on forward selection
PSSM: Position specific scoring matrix
mRMR:Minimum redundancy maximum relevance
dbPTM:Database of post translational modifications
CSV: Comma separated value format
ARFF: Attribute related file.

Conflict of Interests

The authors declared that there is no conflict of interests.

Acknowledgments

The author is grateful to Mr. Ishtiaq Ahmad Qureshy of GCUL for his invaluable support and guidance during the supervision of this project. The author also expresses a debt to all the authors who provided the lists of SUMOylated proteins in their works [21, 45, 53, 55].

References

  1. D. W. H. Girdwood, M. H. Tatham, and R. T. Hay, “SUMO and transcriptional regulation,” Seminars in Cell and Developmental Biology, vol. 15, no. 2, pp. 201–210, 2004. View at Publisher · View at Google Scholar · View at Scopus
  2. R. T. Hay, “SUMO: a history of modification,” Molecular Cell, vol. 18, no. 1, pp. 1–12, 2005. View at Publisher · View at Google Scholar · View at Scopus
  3. P. J. Hamard, M. Boyer-Guittaut, B. Camuzeaux et al., “Sumoylation delays the ATF7 transcription factor subcellular localization and inhibits its transcriptional activity,” Nucleic Acids Research, vol. 35, no. 4, pp. 1134–1144, 2007. View at Publisher · View at Google Scholar · View at Scopus
  4. T. Li, E. Evdokimov, R. F. Shen et al., “Sumoylation of heterogeneous nuclear ribonucleoproteins, zinc finger proteins, and nuclear pore complex proteins: a proteomic analysis,” Proceedings of the National Academy of Sciences of the United States of America, vol. 101, no. 23, pp. 8551–8556, 2004. View at Publisher · View at Google Scholar · View at Scopus
  5. M. S. Y. Huen and J. Chen, “The DNA damage response pathways: at the crossroad of protein modifications,” Cell Research, vol. 18, no. 1, pp. 8–16, 2008. View at Publisher · View at Google Scholar · View at Scopus
  6. T. Li, R. Santockyte, R. F. Shen et al., “Expression of SUMO-2/3 induced senescence through p53- and pRB-mediated pathways,” Journal of Biological Chemistry, vol. 281, no. 47, pp. 36221–36227, 2006. View at Publisher · View at Google Scholar · View at Scopus
  7. C. Fu, K. Ahmed, H. Ding et al., “Stabilization of PML nuclear localization by conjugation and oligomerization of SUMO-3,” Oncogene, vol. 24, no. 35, pp. 5401–5413, 2005. View at Publisher · View at Google Scholar · View at Scopus
  8. F. Melchior, M. Schergaut, and A. Pichler, “SUMO: ligases, isopeptidases and nuclear pores,” Trends in Biochemical Sciences, vol. 28, no. 11, pp. 612–618, 2003. View at Publisher · View at Google Scholar · View at Scopus
  9. V. Dorval and P. E. Fraser, “Small ubiquitin-like modifier (SUMO) modification of natively unfolded proteins tau and α-synuclein,” Journal of Biological Chemistry, vol. 281, no. 15, pp. 9919–9924, 2006. View at Publisher · View at Google Scholar · View at Scopus
  10. M. Liang, F. Melchior, X. H. Feng, and X. Lin, “Regulation of Smad4 sumoylation and transforming growth factor-β signaling by protein inhibitor of activated STAT1,” Journal of Biological Chemistry, vol. 279, no. 22, pp. 22857–22865, 2004. View at Publisher · View at Google Scholar · View at Scopus
  11. M. B. Kroetz, “SUMO: a ubiquitin-like protein modifier,” The Yale Journal of Biology and Medicine, vol. 78, no. 4, pp. 197–201, 2005. View at Google Scholar · View at Scopus
  12. J. S. Seeler and A. Dejean, “Nuclear and unclear functions of sumo,” Nature Reviews Molecular Cell Biology, vol. 4, no. 9, pp. 690–699, 2003. View at Publisher · View at Google Scholar · View at Scopus
  13. J. Zhao, “Sumoylation regulates diverse biological processes,” Cellular and Molecular Life Sciences, vol. 64, no. 23, pp. 3017–3033, 2007. View at Publisher · View at Google Scholar · View at Scopus
  14. F. Z. Watts, “The role of SUMO in chromosome segregation,” Chromosoma, vol. 116, no. 1, pp. 15–20, 2007. View at Publisher · View at Google Scholar · View at Scopus
  15. M. Dasso, “Emerging roles of the SUMO pathway in mitosis,” Cell Division, vol. 3, article 5, 2008. View at Publisher · View at Google Scholar · View at Scopus
  16. F. Zhou, Y. Xue, H. Lu, G. Chen, and X. Yao, “A genome-wide analysis of sumoylation-related biological processes and functions in human nucleus,” FEBS Letters, vol. 579, no. 16, pp. 3369–3375, 2005. View at Publisher · View at Google Scholar · View at Scopus
  17. H. Lomelí and M. Vázquez, “Emerging roles of the SUMO pathway in development,” Cellular and Molecular Life Sciences, vol. 68, no. 24, pp. 4045–4064, 2011. View at Publisher · View at Google Scholar · View at Scopus
  18. R. Giess-Friedlander and F. Melchoir, “Concepts in sumoylation: a decade on,” Nature Reviews Molecular Cell Biology, vol. 8, no. 12, pp. 947–956, 2007. View at Google Scholar
  19. D. C. Bauer, F. A. Buske, and M. Bod’en, “Predicting SUMOylation sites,” in Pattern Recognition in Bioinformatics, Lecture Notes in Computer Science Series, pp. 28–40, Springer, Berlin, Germany, 2008. View at Google Scholar
  20. J. S. Seeler and A. Dejean, “Sumo: of branched proteins and nuclear bodies,” Oncogene, vol. 20, no. 49, pp. 7243–7249, 2001. View at Publisher · View at Google Scholar · View at Scopus
  21. J. Xu, Y. He, B. Qiang, J. Yuan, X. Peng, and X. M. Pan, “A novel method for high accuracy sumoylation site prediction from protein sequences,” BMC Bioinformatics, vol. 9, article 8, 2008. View at Publisher · View at Google Scholar · View at Scopus
  22. D. C. Schwartz and M. Hochstrasser, “A superfamily of protein tags: ubiquitin, SUMO and related modifiers,” Trends in Biochemical Sciences, vol. 28, no. 6, pp. 321–328, 2003. View at Publisher · View at Google Scholar · View at Scopus
  23. A. S. Yavuz and U. Sezerman, “SUMOtr: SUMOylation site prediction based on 3D structure and hydrophobicity,” in Proceedings of the 5th International Symposium on Health Informatics and Bioinformatics (HIBIT '10), pp. 93–97, IEEE, Antalya, Turkey, April 2010. View at Publisher · View at Google Scholar · View at Scopus
  24. K. M. Bohren, V. Nadkarni, J. H. Song, K. H. Gabbay, and D. Owerbach, “A M55V polymorphism in a novel SUMO gene (SUMO-4) differentially activates heat shock transcription factors and is associated with susceptibility to type I diabetes mellitus,” Journal of Biological Chemistry, vol. 279, no. 26, pp. 27233–27238, 2004. View at Publisher · View at Google Scholar · View at Scopus
  25. M. J. Matunis, E. Coutavas, and G. Blobel, “A novel ubiquitin-like modification modulates the partitioning of the Ran-GTPase-activating protein RanGAP1 between the cytosol and the nuclear pore complex,” Journal of Cell Biology, vol. 135, no. 6, pp. 1457–1470, 1996. View at Publisher · View at Google Scholar · View at Scopus
  26. H. Zhang, H. Saitoh, and M. J. Matunis, “Enzymes of the SUMO modification pathway localize to filaments of the nuclear pore complex,” Molecular and Cellular Biology, vol. 22, no. 18, pp. 6498–6508, 2002. View at Publisher · View at Google Scholar · View at Scopus
  27. A. Pichler and F. Melchior, “Ubiquitin-related modifier SUMO1 and nucleocytoplasmic transport,” Traffic, vol. 3, no. 6, pp. 381–387, 2002. View at Publisher · View at Google Scholar · View at Scopus
  28. L. D. Wood, B. J. Irvin, G. Nucifora, K. S. Luce, and S. W. Hiebert, “Small ubiquitin-like modifier conjugation regulates nuclear export of TEL, a putative tumor suppressor,” Proceedings of the National Academy of Sciences of the United States of America, vol. 100, no. 6, pp. 3257–3262, 2003. View at Publisher · View at Google Scholar · View at Scopus
  29. Y. Shinbo, T. Niki, T. Taira et al., “Proper SUMO-1 conjugation is essential to DJ-1 to exert its full activities,” Cell Death and Differentiation, vol. 13, no. 1, pp. 96–108, 2006. View at Publisher · View at Google Scholar · View at Scopus
  30. M. Li, D. Guo, C. M. Isales et al., “SUMO wrestling with type 1 diabetes,” Journal of Molecular Medicine, vol. 83, no. 7, pp. 504–513, 2005. View at Publisher · View at Google Scholar · View at Scopus
  31. J. M. P. Desterro, J. Thomson, and R. T. Hay, “Ubch9 conjugates SUMO but not ubiquitin,” FEBS Letters, vol. 417, no. 3, pp. 297–300, 1997. View at Publisher · View at Google Scholar · View at Scopus
  32. A. C. O. Vertegaal, “Small ubiquitin-related modifiers in chains,” Biochemical Society Transactions, vol. 35, no. 6, pp. 1422–1423, 2007. View at Publisher · View at Google Scholar · View at Scopus
  33. E. van Damme, K. Laukens, T. H. Dang, and X. van Ostade, “A manually curated network of the pml nuclear body interactome reveals an important role for PML-NBs in SUMOylation dynamics,” International Journal of Biological Sciences, vol. 6, no. 1, pp. 51–67, 2010. View at Google Scholar · View at Scopus
  34. R. S. Hilgarth, L. A. Murphy, H. S. Skaggs, D. C. Wilkerson, H. Xing, and K. D. Sarge, “Regulation and function of SUMO modification,” Journal of Biological Chemistry, vol. 279, no. 52, pp. 53899–53902, 2004. View at Publisher · View at Google Scholar · View at Scopus
  35. E. S. Johnson, I. Schwienhorst, R. J. Dohmen, and G. Blobel, “The ubiquitin-like protein Smt3p is activated for conjugation to other proteins by an Aos1p/Uba2p heterodimer,” EMBO Journal, vol. 16, no. 18, pp. 5509–5519, 1997. View at Publisher · View at Google Scholar · View at Scopus
  36. T. Okuma, R. Honda, G. Ichikawa, N. Tsumagari, and H. Yasuda, “In vitro SUMO-1 modification requires two enzymatic steps, E1 and E2,” Biochemical and Biophysical Research Communications, vol. 254, no. 3, pp. 693–698, 1999. View at Publisher · View at Google Scholar · View at Scopus
  37. L. Gong, T. Kamitani, K. Fujise, L. S. Caskey, and E. T. H. Yeh, “Preferential interaction of sentrin with a ubiquitin-conjugating enzyme, Ubc9,” Journal of Biological Chemistry, vol. 272, no. 45, pp. 28198–28201, 1997. View at Publisher · View at Google Scholar · View at Scopus
  38. M. S. Rodriguez, C. Dargemont, and R. T. Hay, “SUMO-1 conjugation in vivo requires both a consensus modification motif and nuclear targeting,” Journal of Biological Chemistry, vol. 276, no. 16, pp. 12654–12659, 2001. View at Publisher · View at Google Scholar · View at Scopus
  39. A. Pichler, A. Gast, J. S. Seeler, A. Dejean, and F. Melchior, “The nucleoporin RanBP2 has SUMO1 E3 ligase activity,” Cell, vol. 108, no. 1, pp. 109–120, 2002. View at Publisher · View at Google Scholar · View at Scopus
  40. S. Weger, E. Hammer, and R. Heilbronn, “Topors acts as a SUMO-1 E3 ligase for p53 in vitro and in vivo,” FEBS Letters, vol. 579, no. 22, pp. 5007–5012, 2005. View at Publisher · View at Google Scholar · View at Scopus
  41. P. Pungaliya, D. Kulkarni, H. J. Park et al., “TOPORS functions as a SUMO-1 E3 ligase for chromatin-modifying proteins,” Journal of Proteome Research, vol. 6, no. 10, pp. 3918–3923, 2007. View at Publisher · View at Google Scholar · View at Scopus
  42. H. Yamamoto, M. Ihara, Y. Matsuura, and A. Kikuchi, “Sumoylation is involved in β-catenin-dependent activation of Tcf-4,” EMBO Journal, vol. 22, no. 9, pp. 2047–2059, 2003. View at Publisher · View at Google Scholar · View at Scopus
  43. M. H. Kagey, T. A. Melhuish, and D. Wotton, “The polycomb protein Pc2 is a SUMO E3,” Cell, vol. 113, no. 1, pp. 127–137, 2003. View at Publisher · View at Google Scholar · View at Scopus
  44. M. H. Tatham, M. C. Geoffroy, L. Shen et al., “RNF4 is a poly-SUMO-specific E3 ubiquitin ligase required for arsenic-induced PML degradation,” Nature Cell Biology, vol. 10, no. 5, pp. 538–546, 2008. View at Publisher · View at Google Scholar · View at Scopus
  45. Y. Xue, F. Zhou, C. Fu, Y. Xu, and X. Yao, “SUMOsp: a web server for sumoylation site prediction,” Nucleic Acids Research, vol. 34, pp. W254–W257, 2006. View at Publisher · View at Google Scholar · View at Scopus
  46. D. Nathan, K. Ingvarsdottir, D. E. Sterner et al., “Histone sumoylation is a negative regulator in Saccharomyces cerevisiae and shows dynamic interplay with positive-acting histone modifications,” Genes and Development, vol. 20, no. 8, pp. 966–976, 2006. View at Publisher · View at Google Scholar · View at Scopus
  47. S. H. Yang, A. Galanis, J. Witty, and A. D. Sharrocks, “An extended consensus motif enhances the specificity of substrate modification by SUMO,” EMBO Journal, vol. 25, no. 21, pp. 5083–5093, 2006. View at Publisher · View at Google Scholar · View at Scopus
  48. X. J. Yang and S. Grégoire, “A recurrent phospho-sumoyl switch in transcriptional repression and beyond,” Molecular Cell, vol. 23, no. 4, pp. 779–786, 2006. View at Publisher · View at Google Scholar · View at Scopus
  49. M. J. Matunis and C. M. Pickart, “Beginning at the end with SUMO,” Nature Structural and Molecular Biology, vol. 12, no. 7, pp. 565–566, 2005. View at Publisher · View at Google Scholar · View at Scopus
  50. M. Mann and O. N. Jensen, “Proteomic analysis of post-translational modifications,” Nature Biotechnology, vol. 21, no. 3, pp. 255–261, 2003. View at Publisher · View at Google Scholar · View at Scopus
  51. L. Lu, X. H. Shi, S. J. Li et al., “Protein sumoylation sites prediction based on two-stage feature selection,” Molecular Diversity, vol. 14, no. 1, pp. 81–86, 2010. View at Publisher · View at Google Scholar · View at Scopus
  52. E. Frank, M. Hall, L. Trigg, G. Holmes, and I. H. Witten, “Data mining in bioinformatics using Weka,” Bioinformatics, vol. 20, no. 15, pp. 2479–2481, 2004. View at Publisher · View at Google Scholar · View at Scopus
  53. T. Y. Lee, H. D. Huang, J. H. Hung, H. Y. Huang, Y. S. Yang, and T. H. Wang, “dbPTM: an information repository of protein post-translational modification,” Nucleic Acids Research, vol. 34, pp. D622–627, 2006. View at Google Scholar · View at Scopus
  54. Y. Xue, F. Zhou, H. Lu, G. Chen, and X. Yao, “SUMO substrates and site prediction: combining pattern recognition and phylogenetic conservation,” http://arxiv.org/ftp/q-bio/papers/0409/0409011.pdf.
  55. E. Boutet, D. Lieberherr, M. Tognolli, M. Schneider, and A. Bairoch, “UniProtKB/Swiss-Prot: the manually annotated section of the UniProt KnowledgeBase,” Methods in Molecular Biology, vol. 406, pp. 89–112, 2007. View at Google Scholar · View at Scopus
  56. I. Ahmad, D. C. Hoessli, W. M. Qazi et al., “MAPRes: an efficient method to analyze protein sequence around post-translational modification sites,” Journal of Cellular Biochemistry, vol. 104, no. 4, pp. 1220–1231, 2008. View at Publisher · View at Google Scholar · View at Scopus
  57. S. Kawashima, H. Ogata, and M. Kanehisa, “AAindex: amino acid index database,” Nucleic Acids Research, vol. 27, no. 1, pp. 368–369, 1999. View at Publisher · View at Google Scholar · View at Scopus
  58. L. Breiman, “Random forests,” Statistics Department, University of California, 2001, http://www.stat.berkeley.edu/~breiman/RandomForests/cc_home.htm.
  59. S. E. Hamby and J. D. Hirst, “Prediction of glycosylation sites using random forests,” BMC Bioinformatics, vol. 9, article 500, 2008. View at Publisher · View at Google Scholar · View at Scopus
  60. J. A. Wohlschlegel, E. S. Johnson, S. I. Reed, and J. R. Yates III, “Global analysis of protein sumoylation in Saccharomyces cerevisiae,” Journal of Biological Chemistry, vol. 279, no. 44, pp. 45662–45668, 2004. View at Publisher · View at Google Scholar · View at Scopus
  61. C. W. Tung and S. Y. Ho, “Computational identification of ubiquitylation sites from protein sequences,” BMC Bioinformatics, vol. 9, article 310, 2008. View at Publisher · View at Google Scholar · View at Scopus
  62. SUMOplot: Abgent. SUMOplotAnalysis Program, http://www.abgent.com/tools/sumoplot.
  63. B. Liu, S. Li, Y. Wang, L. Lu, Y. Li, and Y. Cai, “Predicting the protein SUMO modification sites based on Properties Sequential Forward Selection (PSFS),” Biochemical and Biophysical Research Communications, vol. 358, no. 1, pp. 136–139, 2007. View at Publisher · View at Google Scholar · View at Scopus
  64. S. Teng, H. Luo, and L. Wang, “Random forest-based prediction of protein sumoylation sites from sequence features,” in Proceedings of the ACM International Conference on Bioinformatics and Computational Biology (BCB '10), pp. 120–126, New York, NY, USA, August 2010. View at Publisher · View at Google Scholar · View at Scopus
  65. C. Friedline, X. Zhang, Z. Zehner, and Z. Zhao, “FindSUMO: a PSSM-based method for sumoylation site prediction,” in Advanced Intelligent Computing Theories and Applications with Aspects of Artificial Intelligence, Lecture Notes in Computer Science Series, pp. 1004–1011, Springer, Berlin, Germany, 2008. View at Google Scholar