iSUMOK-PseAAC: prediction of lysine sumoylation sites using statistical moments and Chou’s PseAAC

Sumoylation is the post-translational modification that is involved in the adaption of the cells and the functional properties of a large number of proteins. Sumoylation has key importance in subcellular concentration, transcriptional synchronization, chromatin remodeling, response to stress, and regulation of mitosis. Sumoylation is associated with developmental defects in many human diseases such as cancer, Huntington’s, Alzheimer’s, Parkinson’s, Spin cerebellar ataxia 1, and amyotrophic lateral sclerosis. The covalent bonding of Sumoylation is essential to inheriting part of the operative characteristics of some other proteins. For that reason, the prediction of the Sumoylation site has significance in the scientific community. A novel and efficient technique is proposed to predict the Sumoylation sites in proteins by incorporating Chou’s Pseudo Amino Acid Composition (PseAAC) with statistical moments-based features. The outcomes from the proposed system using 10 fold cross-validation testing are 94.51%, 94.24%, 94.79% and 0.8903% accuracy, sensitivity, specificity and MCC, respectively. The performance of the proposed system is so far the best in comparison to the other state-of-the-art methods. The codes for the current study are available on the GitHub repository using the link: https://github.com/csbioinfopk/iSumoK-PseAAC.


INTRODUCTION
Post-translational modifications are the chemical changes that occur during the structural and functional organization of a protein after the protein is synthesized by the translational process. These chemical reactions or changes that take place at certain amino acid residue after the translational process are known as post-translational modifications (PTMs). The post-translational modifications can be covalent or enzymatic. These modifications involve phosphorylation, glycosylation, ubiquitination, nitrosylation, methylation, acetylation, lipidation, and proteolysis have impacted all the details of cellular biology and pathogenesis. The post-translational modification also has involvement in protein for various functions with a minimal number of genes and can influence protein actions such as gene expressions. PTM modulates the cellular behavior and relates to any change in the amino acid chains of the protein after its alteration. They are moderated by enzymatic activity and can be reversed according to its change using enzymatic actions (Lu et al., 2010;Zhao et al., 2014;Beauclair et al., 2015).
Small ubiquitin-like modifier (SUMO) is like a group of tiny proteins that may be covalently similar to other proteins in the cells and alter their functions. Sumoylation is a post-translational modification that performs various cellular functions such as nuclear-cytosolic transport, transcriptional ordinance, protein reliability, and development throughout the cellular cycle. SUMO is a reversible PTM in which small ubiquitin proteins are covalently bound to a lysine residue in a process same as ubiquitylation. SUMO is constructed after the binding of the final four amino acids having C terminals and is responsible for an iso-peptide bond among the C terminal glycine residue of SUMO and lysine (Green, Dmochowski & Golshani, 2006;Geiss-Friedlander & Melchior, 2007).
SUMO modified proteins have peptide concord motifs Ψ, K, x, and D/E. Here 'Ψ' is the hydrophobic amino acid, 'k' is the lysine residue, 'x' is the amino acid and 'E' is the glutamic acid. The consensus motifs were examined with a percentage of 23% in the preliminary study. Furthermore, in recent studies (Rodriguez, Dargemont & Hay, 2001;Sampson, Wang & Matunis, 2001;Xue et al., 2006) ψ-K-x-E\D was the consensus motif and reveals that 26% of SUMO sites did not accompany the consensus motifs. In short from the current study on Sumoylation it was observed that 40% of sumoylation sites did not have the consensus motifs (Zhao et al., 2014), and the remaining used the consensus motifs. This lack of information has led to insufficient knowledge for Sumoylation.
Lysine is a post-translational modification residue. It is the most used amino acid of PTM having Sumoylation and also narrates for various sumoylation sites. Lysine sumoylation is the reversible most modulated PTM. It lies via the covalent link of a small ubiquitin-like modifier (SUMO) and is moderated by the generation of enzymes E1 link enzymes E2 and ligase E3. Lysine residue endures many PTM as there must be synchronization among them (Müller et al., 2001).
In past, biologists were expected to perform the traditional experiments with the utilization of costly equipment in order to identify the post-translation modification of a protein. In recent years, the improvement of bioinformatics permitted scientific community to solve complex protein problems by the combination of informatics, mathematics, and statistics. Due to the reversible behavior of sumoylation, it was observed that the different procedures or methods used for the deficiency of sumoylation site modification must be essential for the site prediction process. Green, Dmochowski & Golshani (2006) proposed parallel cascade identification method for the prediction of Sumoylation sites. The main purpose behind this method was to originate PCI based proteomics tools for the prediction of protein structure and function. Xue et al. (2006) introduced a novel computational method SUMOsp for the prediction of Sumoylation sites. SUMOsp was based on manually curated datasets using the integration of two methods, GPS and MotifX. SUMOsp used large datasets and is considered a vigorous tool for in vivo and in vitro sumoylation site prediction. Zhao et al. (2014) proposed a GPS-SUMO approach used for the identification of sumoylation sites and sumo relating motifs (SIM) in proteins and for examining the association among the sumoylation sites and SUMO relevant techniques. A webserver for GPS-SUMO is also available for the research community to further utilize the dataset and methods. Beauclair et al. (2015) proposed a scoring system based on position frequency matrix derived from the alignment of experimental sumoylation sites. The specified tool, JASSA revealed high gains in proportion of implementation. Therefore JASSA was considered as a costly tool for evaluating the ideal sites and providing useful information about Sumoylation in cellular biology. Chang et al. (2018) proposed SUMOgo using Random Forests, motif screening, and feature selections based on variant in combinations to originate a Sumoylation site prediction. They used sequence-based binary encoding, encoded chemical attributes, and encoded secondary structure details. They enhanced the prediction performance and obtained the MCC 0.51 in comparison to the other state-of-the-art methods. Ijaz (2013) also developed a prediction method that showed improvements in the prediction process of Sumoylation site and its significance. In recent past, López et al. (2020) implemented a computational prediction method using the Adaboost classifier. They utilized the sine and cosine of backbone torsion angles and the accessible surface area. To overcome the balancing issues in their training matrix, they applied NearMiss method with undersampling the majority class. C-iSUMO was effective in use of circular functions. C-iSUMO was effective in prediction of Sumoylation sites as compared to the other state of art systems from the past and achieved 74.6% accuracy with 0.494 MCC.
In this study, we propose a novel method iSUMOk-PseAAC, to predict the Sumoylation sites using relative position based features by integrating the Chou's Pseudo Amino Acid Composition (PseAAC) (Chou, 2001a). The results from the proposed system were compared with the other state of art prediction methods SUMOgo (Chang et al., 2018), GPS-SUMO (Zhao et al., 2014), SUMOsp2.0 (Xue et al., 2006), JASSA (Beauclair et al., 2015), and PCI-SUMO (Green, Dmochowski & Golshani, 2006). The proposed method was implemented using the Chou's five-step rule . This method has been utilized by various studies (Cheng, Xiao & Chou, 2017a, 2017b, 2017cCai et al., 2018;Chen et al., 2018aChen et al., , 2018b which follow these five steps: (i) Benchmark dataset construction, (ii) Formulation of Samples, (iii) Operational Algorithm, (iv) Tests using cross-validations, and (v) Implementation of a webserver. These steps are discussed in detail in further sections. The framework used in the proposed system used Chou's fivestep rule. This approach is shown in Fig. 1.

Benchmark dataset
In many studies, Chou's peptide formulation approach was used for the simplification of the dataset (Chou, 2001a). Universal resource of protein (UniProtKB) is a freely available as a central repository of protein sequences information. The benchmark dataset was constructed using the UniProtKB. The protein sequences with PTM processing annotations and features were used for the dataset construction of Sumoylation sites. A term named glycl lysine was used for Sumoylation domain to obtain the sequences of SUMO proteins. Only those sites were retrieved and included in datasets which were reviewed and annotated with experimental evidences. There were many redundant sequences or highly similar sequences that existed in the newly constructed dataset. For this purpose, those sequences were removed with the help of a tool CD-HIT (Fu et al., 2012). The similarity threshold value for removing redundant and ambiguous sequences was kept at 60%. Finally, the benchmark dataset was obtained after removing the similar sequences. The negative datasets were also obtained from the UniProtKB using the complement sequences of the searching criteria used for positive Sumoylation sites. The benchmark dataset consists of 4,987 positive and 5,000 negative samples. This dataset was divided into training and testing dataset using the 70:30 ratio. For training data, 3,487 positive and 3,500 negative random samples were selected. Furthermore, an independent testing dataset was constructed from the leftover samples after constructed the training dataset and 3,000 independent test samples were left in which 1,500 are positive and 1,500 are negative samples. In prediction models based on statistical analysis, construction of training and testing datasets is crucial. The benchmark datasets can be accordingly reduced to Eq. (1): (1) According to Eq. (1), here U is denoted as the union operator. K + represents the 4,987 positive samples and K − represents the 5,000 negative samples. The total sum K consists of 9,987 sample sequences which are provided in the Supplemental Information S1 for the convenience of the readers. The independent dataset is provided in the Supplemental Information S2. For the graphical representation of lysine amino acids, we analyzed the sequence samples by using the development of Web Logo (Chou, 2001a) (see Figs. 2 and 3). In this logo, the symbol 'X' was attached as the dummy code to keep the same length of all samples.

Sample formulation
As the biological sequences are increasing in biological databanks at exponential rate, the position of discrete models or vectors from the sequences of biological data has been an issue in many bioinformatics based research methods and its attributes for target analysis. There are many algorithms as discussed in the review study  for the formulation of sequence in vector form, but the machine learning algorithms are among the top-rated algorithms such as Random Forest (RF) algorithm (Lin et al., 2011;Jia et al., 2016), Covariance Discriminant (CD) (Chou & Elrod, 2002;Lin et al., 2012), etc. The pseudo amino acid composition (PseAAC) (Chou, 2001b) was proposed to overpower the sequence pattern-related information of protein which might be a loss if they are represented in the form of vector or discrete model. Chou's PseAAC has been utilized for many computational proteomics problems (Butt et al., 2016;Butt, Rasool & Khan, 2017, 2019Sabooh et al., 2018;Sankari & Manimegalai, 2018;Srivastava, Kumar & Kumar, 2018;Zhang & Liang, 2018;Zhao et al., 2018;Butt & Khan, 2020a, 2020b; Chou, 2001c). Because of its prevalence and significance in computational proteomics, there are three webservers for the facilitation of the research community called PseAAC-General (Du, Gu & Jiao, 2014), PseAAC-Builder (Du et al., 2012), and propy (Cao, Xu & Liang, 2013). The previous two were used for Chou's PseAAC by many researchers [68] but the initial one was used for functional domains, mods, gene ontology, and sequential evaluation of Chou's general PseAAC . After using PseAAC for many protein prediction problems, PseKNC was proposed for the computation of features vectors using DNA/RNA sequences (Liu et al., 2015a;Chen et al., 2016a;Feng et al., 2017). Nowadays, most competent and effective webserver named 'Pse in One' (Liu et al., 2015b) and 'Pse in One2.0' (Liu, Wu & Chou, 2017) was evolved for feature vector design from protein sequence, peptide sequence and DNA/RNA sequence. By applying the Chou's general PseAAC , protein site sequence 'S' can be represented as Eq. (2): Here N 1 in protein sample 'S' represents the first residue, N 2 as the second amino acid residue, and similarly, N P , as the last amino acid residue. P is the total length of the protein site sequence or the total number amino acids in a sequence. The following discrete model based on amino acid composition is used to represent protein S in Eq. (3): The extraction of useful features is very important from the relevant protein/peptide sequences which are explained in Eq. (3). Where F X (X = 1, 2, 3, … , 20) are the useful amino acid features and symbol T is the transpose of a site sequence of the protein S. Each sample site in dataset was represented as a peptide segment of length 41 with 20 amino acid residues upstream and 20 residues downstream of the amino acid residue 'K' lysine. The following Eq. (4) represents the modified sample sequence from Eq. (2): Here N 21 ¼ K which represents the target lysine residue and N P P ¼ 1; 2; 3; . . . ; 41; ð P 6 ¼ 21Þ can be any other amino acid or virtual X code. A virtual amino acid residue 'X' fills the positions where no sufficient residues exist to make sure that the length of the peptide segment remains unified. From now onwards, the numerical codes of amino acid used as per order of their name as alphabetical order according to their first letter, 1, 2, 3 …… 20 for 1 to all 20, one per amino acid and 21 for X.

Site vicinity vector
In the polypeptide chain, the chance of PTM is much higher. There are multiple reasons for that modification such as lysine ubiquitination sites, methylation sites, sumoylation sites, and numerous other lysine PTM sites. Observing the potential site along with neighboring residues where PTM is also an important aspect (Lo & Don, 1989). A small sequence that contains the potential PTM residue site from the main sequence is called Site Vicinity Vector (SVV).
Let Ωi be the variable that shows the possible PTM site containing neighboring residues in the main sequence are given below in Eq. (5) The SVV can be defined as a subsequence of the primary sequence as given below as Eq. (6): The symbol i holds the minimum value of 20. In subsequence of protein where every residue symbolizes a specific amino acid out of known 20 amino acids. In SVV to assigned the unique numerical code to each residue position from the range of 1-20 of amino acid and one value 21 for the virtual amino acid 'X'. The value of i and j is the size of the window, which is chosen after extensively probing different values from a given range (10-100). The values that outcome best performance in results are selected and used accordingly.

Statistical moments
In most problems of pattern recognition, the statistical moment based features have been widely used for the qualitative measures of benchmarks datasets. Due to sensitive information of the sequence of amino acid residue in the benchmark dataset, a statistical approach was applied to manage the order of sequence samples of protein. The collection of numerous types of information was derived from a variety of statistical moments from the dataset of protein samples, to use the evaluation of data size and some other orientation and eccentricity. Various kinds of moments were used by the mathematicians and statisticians (Khan, Ahmad & Anwar, 2012;Khan, Ahmed & Khan, 2014) and these moments are depended on the distribution functions and polynomials.
Moments are a collection of statistical parameters used to understand the characteristics of a function and to capture its significance features; therefore Hahn, Raw and Central moments were calculated for the iSUMOk-PseAAC prediction model. Hahn moments are scale and position variant (Khan, 2014a) and calculated using the Hahn polynomials. Raw moments are also scaled and position variant calculated using the probability distribution of the dataset. Raw moments are very helpful for the calculation of asymmetry of samples of proteins, variance, and mean of benchmark dataset. Furthermore, central moments were also calculated which are scale variant and vicinity invariant (Butt et al., 2016;Butt, Rasool & Khan, 2017) and based on a centroid. It also has the property for the computation of variance, asymmetry, and mean of benchmark dataset.
The proposed method used a two-dimensional matrix Ṕhaving (n × n) dimensions that contain all the samples of proteins of the dataset. In this study, each method characterizes a benchmark dataset that represents quantified values and by passes the scale variant moments (Khan et al., 2014b).
Á Á Á R n!n 2 6 6 6 6 6 6 6 6 4 3 7 7 7 7 7 7 7 7 5 The benefit of the transformation of 2D matrix P (see Eq. (7)) is the fast calculation by using Hahn moments calculation. The transformation matrix was done into P using a hold function discussed in recent study (Akmal, Rasool & Khan, 2017) and all moments are computed up to third-degree by using the elements of P. Two-dimension matrix used only orthogonal input, therefore we computed Hahn moments. To restore the benchmark dataset used the inverse function of orthogonal Hahn moment. By computing the Hahn moments to 'N' order, the following Eq. (8) was used: According to Eq. (8), the pochammer symbol and the Gamma operator are discussed in the review study (Akmal, Rasool & Khan, 2017). The following Eq. (9) is used for the computation of normalized orthogonal Hahn moments.
c qr g h u;x q i; P ð Þ g h u;x r j; P ð Þ; n ¼ 0; 1; . . . ; P À 1 Vital information has been stored in protein samples by calculating the central moments which are related to mean, variance, and asymmetry. The following Eq. (10) is given for the calculation of the central moment.
To store the essential protein samples of the dataset by calculating raw moments using the property of probability distribution as followed in the given Eq. (11): The raw moments are calculated by the third degree which is r + S and M00, M01, M09, M10, M19, M20, M29, M30, and M03.

Frequency vector (FV)
In a benchmark dataset, one of the important parameters is frequency distribution, which gives valuable information about the distribution of the dataset for each sample of protein. Frequency is store in the form of a vector, this frequency is calculated for every amino acid, and this is called the frequency vector. The frequency vector tells us about the distribution and composition of the protein sample sequence. The computation of the frequency vector is defined in Eq. (12) as follows: where r i explaining the frequency of each single amino acid residue in alphabetical order.

Position relative incidence matrix (PRIM)
The relative position of amino acid is directly related to the physical characteristics of the protein. The position relative incidence matrix (PRIM) is the relative position of amino acid in the polypeptide chains. The input query of protein formation is the basic step for feature extraction and the size of protein sequences builds PRIM and frequency matrix (FM). The matrix is very helpful for the computation of moments by which the feature vectors are formed. A matrix (21×21) is formed called Z PRIM and it shows the protein sequences with relative associated information of residue through Eq. (13) as follow: In the Z PRIM matrix, element Z ij is the sum value of i th residue by calculating the relative position with the first happening of j th residue. PRIM gives up to 441 coefficients which is a huge number, therefore computed statistical moments using PRIM as the input for the minimization of coefficient and it attains thirty coefficients.

Reverse position relative incidence matrix (R-PRIM)
The PRIM matrix has information regarding the relative position of amino acids along with the polypeptide chains. Sample of protein sequences with reverse input is calculated by using Z PRIM for finding out the unclear hidden features from the potential sequences of protein containing the similar position uncertainty of protein. Below Z R-PRIM is calculated as Eq. (14): For decreasing the dimensionality of Z R-PRIM , by applying the same method to perform the coefficient reduction and statistical moments which gives back the 441 coefficients similar to Z PRIM , and obtained thirty coefficients.

Accumulative absolute position incidence vector (AAPIV)
The frequency matrix gives the information about the unclear features of samples of protein in the polypeptide chains but the accumulative absolute position incidence vector (AAPIV) gives the information of relative positional information regarding amino acid residues. Having the 21 amino acid that each single amino acid used for the sum of original values by the residues that may exist in the primary structure as Eq. (15): The computation of the arbitrary 'u i ' element of AAPIV as given below as Eq. (16):

Reverse accumulative absolute position incidence vector (R-AAPIV)
RAAPIV is formed by reversing the primary sequence of the protein and computed with the similar technique of AAPIV to know the obscure properties regarding relative positional data. Calculated the R-AAPIV by utilizing the reverse proteins sequences samples as Eq. (17):

Classification algorithm
The human brain comprises billions of neurons and they all have individual features and perform various tasks, get knowledge, take actions, and exchange information. When the neurons are operated, the Brain gains the information and performs the tasks according to conditions without having any observation or experience. The artificial neural network mechanism is constructed from the brain like the system to learn from the pattern, and past knowledge for various issues. It comprises particular neurons that perform various functions. They get the details from neurons, perform actions, and then utilize the design or model from various examples and experiences. The ANN acts on 2 models, first one is training where ANN is trained on specific data and the targeted attribute is given from the dataset to learn the pattern. The second one is testing in which we test our predicted model by giving the unseen data to its input and find the best result using trained examples and accessible information (Jiang et al., 2016) as shown in  Fig. 4), in short IFM was used for getting input and OM was used for the output and compute errors during learning the model through the backpropagation technique (Jiang et al., 2016).
The scikit-learn library was used to implement the neural network. The "max_iter" parameter was increased from the default parameter value of 200 to 437. The "max_iter" parameter value was optimized to 437 using hyper parameter tuning methods and optimal value for the parameter was searched using the successive halving technique in scikit-learn. The searching space for the parameters "hidden_layer_sizes", "learning_rate" and "momentum" was (10-50), (0.0001-0.1) and (0-1) respectively. These parameters were optimized to 23, 0.001 and 0.7 for the parameters "hidden_layer_sizes", "learning_rate" and "momentum", respectively, after successful halving technique. One of the key findings observed during the experimentation process was that "max_iter" with more than 500 iterations minimally contributed to the accuracy of the classifier substantially.

RESULTS
Evaluation of the new machine learning prediction model is the essential step for the success rate and completeness of the model . For the best examination of the model, we must focus on two parts of the research study of the prediction model, the type of metrics we can use to represent the model prediction standard and the test methods that are best for scoring the metrics.

Evaluation metrics
For the correctness and the efficiency of the proposed model, we must consider four main points' metrics: (1) Find out the accuracy Acc of the proposed method, (2) Find out sensitivity Sn of the prediction model, (3) Find out specificity Sp of a predictor, and (4) Find out MCC for the strength of a predictor. The traditional metrics are mostly used in math to validate the accuracy of the prediction model, but it may be a very complex job for the biologist to perceive and use these metrics. For this reason, Mathew's correlation coefficient (MCC) is the best matrices used for the prediction model reliability and strength. Most of the symbols used in protein signal peptide by Chou (Chou), a collection of four built-in equations were worked (Chou, 2001b(Chou, , 2015Xu et al., 2013;Arif, Hayat & Jan, 2018) as follow in Eq. (18).
As explained in above Eq. (18) that when the value of И þ À ¼ 0 which simply means that the predictor predicts all the SUMOk sites correctly and no one reaming to predict the sites, therefore we have the sensitivity Sn = 1. On the other side if И þ À ¼ И þ which means the predictor predicts all the SUMOk sites incorrectly, therefore we have the sensitivity Sn = 0. Similarly, when И À þ ¼ 0 which means the predictor predicts all the non-SUMOk sites correctly and no one reaming to predict, therefore we get specificity Sp = 1. Further, if we have И À þ ¼ И À which means that the predictor predicts all the non-SUMOk sites incorrectly, therefore we get the specificity Sp = 0. Furthermore, when И þ À ¼ И À þ which means that the predictor predicts correctly for all positive dataset as well as for all negative dataset and no remaining sites for prediction, that gives the overall accuracy Acc = 0 and MCC = 1. On the other hand, if И þ À ¼ И þ and И À þ ¼ И À which means that the predictor predict all the values for the positive dataset and negative dataset incorrectly and we get MCC = −1 and Acc = 0; whereas if И þ À ¼ И þ 2 and И À þ ¼ И À 2 it gives the MCC = 0 and Acc = 0.5, which means nothing, and a random guess is better. Therefore Eq. (18) provides the detailed meaning of stability, comprehensive accuracy, specificity, and sensitivity for a better understanding, as discussed by many investigators (Chen et al., 2016b;Xiao et al., 2016).
In recent publications (Lin et al., 2014;Xu et al., 2014;Chen et al., 2016c;Zhang et al., 2016;Ehsan et al., 2018), these set of metrics were utilized for research in state-of-the-art methods. According to Eq. (18), SUMOk sites or non-SUMOk sites are applicable only for the binary classification data. Multi-label class problem is not applicable for this kind of solution, because this problem is different in biomedicine  and biological (Xiao, Wu & Chou, 2011;Chou, Wu & Xiao, 2012;Lin et al., 2013), therefore a different set of metrics required for this problem as explained in this study .

Self-consistency tests
Self-consistency testing was implemented for iSUMOk-PseAAC sites, which mean we have to train and test on the benchmark dataset . It is mostly used when the results are already known and saved in dataset; here we used true positive results for validation of our model. Results of our validation are shown in Table 1, which show all the predicted and actual classification. This tells us about the overall performance of our model.

Independent tests
Independent dataset testing plays a very important role in the evaluation of machine learning models because this testing dataset is divided into two parts one is the training dataset and the other is the testing dataset. The dividing ratio of the dataset can be different. In this research 70:30 ratio is used for training and testing and measure the highest accuracy which is given below in Table 2.

Receiver operating characteristics (ROC)
The performance measure is the important step of every machine learning model, therefore we draw Area under the curve receiver operating characteristic (AUC-ROC) for classification problem. It represents the model's accuracy in the form of true positive rate and false-positive rate. If the AUC is higher or near to 1 which means the model performed well and separate the class correctly. If the AUC is worst or near to 0 which means the performance of the model is poor and fails to separate the class. If AUC is 0.5, it means the performance of the model is neutral. This curve is plotted with TPR against the  Fig. 8.

10-Fold cross-validation test
Normally a dataset for prediction model is needed which should be labeled, in some cases which are not easy to obtain, so this testing of the dataset is taken from the original dataset, which was also used for training, this set is very useful to test the performance of the model. Especially in cross-validation, a dataset is divided into K parts where K can be any number and then one part of that data is kept for testing while other parts are used for training purposes. Next, a new part of data is chosen as the testing sample while others are used for training, this is repeated until all parts are tested. Finally, the average of all the results is calculated. In our case, the value of K = 10 and the result was average of all the accuracy for each fold. Our results are illustrated in Table 3. Figures 9-11 show the 10 Fold cross validation results ROCs.

DISCUSSIONS
There have been many predictors to predict the Sumoylation site, but their results not good as compared to our iSUMOk-PseAAC predictor. Recently Chang et al. (2018) proposed the SUMOgo predictor, which shows the highest accuracy of other Sumoylation site prediction tools through independent testing. To detect the overall accuracy of a prediction   (Zhao et al., 2014), and PCI-SUMO (Green, Dmochowski & Golshani, 2006). The overall accuracy in terms of ACC, Sp, Sn, and MCC is much higher in the proposed Sumoylation tool because iSumok-PseAAC performed well under statistical moments based features extractions as shown in Table 4 of their accuracy metrics. The PseAAC methods were utilized by the position relative features and statistical calculation for prediction of sumoylation sites. First of all 20 residues of the amino acid are surrounding left and 20 rights for targeted residue, after that computed statistical moments and for reduction of dimensionality constructed the frequency vector, site vicinity vector, PRIM, RPRIM, AAPIV, and RAAPIV. Finally, the composition of sequence and position relative feature is given to predictor for prediction. The Fig. 12 represents the ROCs analysis for state-of-the-art methods.

WEB-SERVER DEVELOPMENT
The final step of Chou's five-step rule  is the Web Server, which plays an important role in various computational analyses and findings. In recent research studies Cheng, Xiao & Chou, 2018c;Xiao et al., 2018), the effect of computational biology on medical science as effectively enhanced the availability, and also this server provides user-friendly environment (Chou & Shen, 2009) for the easiness of user as well as biologists, which drives the medical science (Chou, Figure 12 Comparative analysis ROCs of state-of-the-art methods. Full-size  DOI: 10.7717/peerj.11581/ fig-12 2015) into an exceptional revolution . The software code for the current study is provided at GitHub: https://github.com/csbioinfopk/iSumoK-PseAAC. The webserver to the current study will be provided for the research community in near future.

CONCLUSIONS
The presence of SUMO modification has occurred for over 10 years. Sumoylation plays an essential role in the regulation of various cellular functions. Sumoylation progress towards various biological processes, diseases, and medication; therefore, it is considered one of the significant aspects of cellular functions along with transcriptional ordinance, protein reliability, and the development throughout the cell cycle. Sumoylation relations with many different types of diseases like cancer, diabetes, inherit heart flaws, and most important neurodegenerative diseases which are directly linked to Sumoylation synchronization and modulation. Consequently, the classification of potential Sumoylation sites is beneficial, for this purpose we propose a SUMOk prediction tool, which is more accurate and efficient for easy experimental results followed by the five-step rule. To test the accuracy of iSUMOk-PseAAC, 10-fold cross-validation was implemented with the help of metrics. We achieved the result of 10 fold cross-validation with 94.51% accuracy, 94.24% sensitivity, 94.79% specificity, and 0.8903% MCC. For that reason, iSUMOk-PseAAC predictors are very helpful for predicting Sumoylation sites in an accurate and precise manner, although the results of the proposed model would be better by grouping the growing number of SUMO sites in a sequence of the protein.

ADDITIONAL INFORMATION AND DECLARATIONS Funding
The authors received no funding for this work.