Skip to main content
Download PDF

Incorporating support vector machine with sequential minimal optimization to identify anticancer peptides

BMC Bioinformatics volume 22, Article number: 286 (2021) Cite this article

Abstract

Background

Cancer is one of the major causes of death worldwide. To treat cancer, the use of anticancer peptides (ACPs) has attracted increased attention in recent years. ACPs are a unique group of small molecules that can target and kill cancer cells fast and directly. However, identifying ACPs by wet-lab experiments is time-consuming and labor-intensive. Therefore, it is significant to develop computational tools for ACPs prediction. Though some ACP prediction tools have been developed recently, their performances are not well enough and most of them do not offer a function to distinguish ACPs from antimicrobial peptides (AMPs). Considering the fact that a growing number of studies have shown that some AMPs exhibit anticancer function, this work tries to build a model for distinguishing AMPs from ACPs in addition to a model that predicts ACPs from whole peptides.

Results

This study chooses amino acid composition, N5C5, k-space, position-specific scoring matrix (PSSM) as features, and analyzes them by machine learning methods, including support vector machine (SVM) and sequential minimal optimization (SMO) to build a model (model 2) for distinguishing ACPs from whole peptides. Another model (model 1) that distinguishes ACPs from AMPs is also developed. Comparing to previous models, models developed in this research show better performance (accuracy: 85.5% for model 1 and 95.2% for model 2).

Conclusions

This work utilizes a new feature, PSSM, which contributes to better performance than other features. In addition to SVM, SMO is used in this research for optimizing SVM and the SMO-optimized models show better performance than non-optimized models. Last but not least, this work provides two different functions, including distinguishing ACPs from AMPs and distinguishing ACPs from all peptides. The second SMO-optimized model, which utilizes PSSM as a feature, performs better than all other existing tools.

Background

Cancer is a leading cause of death and the most important barrier to increasing life expectancy worldwide in this century [1]. This disease is caused by the growth and uncontrolled proliferation of abnormal cells. Conventional cancer treatments, including radiation therapy and chemotherapy, often have adverse effects on normal cells and thus not effective enough [2]. Moreover, some mechanisms also lead to drug resistance from the cancerous cells [3]. Therefore, a novel treatment that lacks adverse effects, targets specifically to cancer cells, and with a low possibility of drug resistance is in need urgently.

In recent years, a new group of small peptides, ACPs, has been discovered that can target and kill cancer cells specifically while not affecting healthy cells [4, 5]. The high selectivity and cancer-selective toxicity [6] of ACPs depend on multiple differences between cancer cells and normal cells, including membrane net charge and unique molecules on the membrane [6]. Due to their specificity and low toxicity, ACPs have attracted growing attention as a novel cancer treatment and have been considered to be promising [7]. For example, romidepsin (FK228), has been shown to have clinical effectiveness in patients with refractory cutaneous T-cell lymphoma [8]. To promote its application, it is of great significance to distinguish ACPs from all peptides. Nevertheless, finding anticancer peptides by experiments could be both time-consuming and labor-intensive [9]. To deal with this problem, computational identification before wet-lab experiments is necessary. Machine-learning-based methods could be of great help to classify and predict those special peptides. Moreover, some characteristics of cancer cells, such as the negative surface charge of their membrane, also shared by bacterial cells [10]. In fact, a hypothesis is proposed that ACPs share similar features with another group of small molecules that can specifically target and kill microbes, called AMPs [11]. Indeed, some AMPs are discovered to exhibit anticancer function according to recent studies [12]. Thus, distinguishing ACPs from AMPs may promote the discovery of ACPs more accurate, more convenient and faster.

To identify and predict ACPs, many computational tools for predicting ACPs have been designed, including Hajisharifi’s model [13], AntiCP [14], iACP [15], MLACP [16], mACPpred [17], ACPred [18], ACPred-Fuse [19] and ACPred-FL [20]. Hajisharifi et al. use physicochemical properties and PseAAC as characteristics of peptide sequences, and SVM as a machine learning method to identify ACPs. Their method is claimed to perform with an accuracy of 83.82% [13]. By analyzing the AAC of peptides and using SVM as a machine learning method, AntiCP offers two models that can distinguish ACPs from either AMPs or non-ACPs based on different datasets [14]. In MLACP, they analyze the AAC, dipeptide composition, atomic composition and physicochemical properties separately and hybridlike. Then they apply two machine learning methods: SVM and random forest to build models based on peptide characteristics. The performance of MLACP is claimed to be better than any other existing methods, with an accuracy of 87.5%. The deficiency of the MLACP study is that it does not offer a model that can distinguish ACPs from AMPs [16]. mACPpred, which achieves an accuracy of 88.5% in their independent test, uses SVM as the final classifier. ACPred also utilizes SVM and analyzes several different features, the accuracy of which is 97.56% according to their paper. ACPred-FL incorporated feature representation learning and feature selection with SVM. The prediction accuracy of this tool in their independent test is 85.7%. Similarly, ACPred-Fuse fuses a feature representation learning model that integrates 29 different features with random forest, and performs an 89% accuracy in their independent test.

This research offers more functions and better performance. First of all, sequences of examined ACPs, non-ACPs, and AMPs without anti-cancer functions are collected. With these data, two different groups of datasets are constructed: (1) inspected ACPs as positive data and AMPs without anti-cancer function as negative data; (2) examined ACPs as positive data, simple non-ACPs as negative data. Then characteristics of those peptide sequences are analyzed, considering four features, amino acids composition (AAC), N5C5, k-space and PSSM, separately and also hybridized. It is the first time that PSSM is considered as a feature in ACP prediction studies. Based on the analysis of those features, several models are built based on two machine learning methods: SVM [21] and SMO [22]. Comparing the performance of those models, two best ones are chosen: SMO-1, which utilizes SMO to analyze AAC and k-space feature of a dataset (1), and SMO-2, which uses SMO as well and is based on analysis of AAC, N5C5 and k-space of the dataset (2). At last, the same testing dataset is applied to test the performance of SMO-1, SMO-2, AntiCP-1, AntiCP-2, mACPpred, ACPred, ACPred-Fuse and ACPred-FL. As for results, comparing to AntiCP-1, which is also designed to distinguish ACPs from AMPs, SMO-1 shows higher accuracy, specificity and Matthews Correlation Coefficient (MCC). Also, the performance of SMO-1 is of better balance. As for SMO-2, identifying ACPs from all kinds of peptides, performs better with consideration of accuracy, sensitivity, specificity and MCC, and shows relatively more balanced performance than AntiCP-2, mACPpred, ACPred, ACPred-Fuse and ACPred-FL do. In general, this research built two models with different functions: one is for predicting ACPs from AMPs, which share some similarities to ACPs, and another one is used to distinguish ACPs from all peptides. The second SMO-optimized model shows better performance than the unoptimized model and other existing tools.

Results

Characterization of the sequence-based features of ACPs

Comparing to AMPs but non-ACPs (peptides in negative dataset 1), K, L, A are much more frequent in ACPs, whereas N, Y, Q are dominant in negative dataset 1(with the lowest p-values). Similarly, comparing to non-ACPs in negative dataset 2, L, W, A are dominant in ACPs, whereas M, R, Q are dominant in non-ACPs (with the lowest p-values) (Fig. 1). Those significant differences in the frequency of each amino acid in different datasets contribute greatly to later classification.

Fig. 1
figure 1

AAC analysis of positive, negative 1, negative 2. The frequency percentage of each amino acid in positive, negative1 and negative 2 data groups is shown in blue, orange and gray respectively

Full size image

Also, some of these results agree with the physicochemical properties of amino acids. For example, K is the most predominant amino acid in the positive dataset while hydrophobic positively charged lysine-rich peptides which act as cationic peptides that can indeed interact with anionic membranes on cancer cells, disrupt the cell membrane integrity, penetrate into the membrane and thus serve as ACPs [23].

According to the result of the positive dataset, K, L are the two most dominant amino acids in N5C5 of ACPs. Taking position under consideration, K is dominant in the third position of C-terminal end, L is dominant in the first position of C-terminal end, and G is dominant in the first position of N-terminal end. G is also dominant in the first position of N-terminal end of the negative 1 group (AMPs but non-ACPs). Contrarily, M is the most frequent one in the first position of N-terminal end of non-ACPs in the negative dataset 2. Comparing positive dataset to negative 1 dataset, significant differences can be found: A, L, F, K are more dominant in the positive dataset while C is more dominant in the negative 1 dataset. On the contrary, comparing positive dataset to the negative 2 dataset, distributions of each amino acid in each position are more divergent, and less contrasts could be extracted (Fig. 2).

Fig. 2
figure 2

N5C5 analysis results. The three heatmaps on the left show frequency of each amino acid at each position in three different data groups, including positive, negative 1 and negative 2 data. The two heatmaps on the right show different values of each amino acid at each position comparing positive to negative 1 and negative 2 separately

Full size image

Finally, with X representing spacings between amino acids, the ten most diverse k-space pairs comparing the positive data to the negative 1 data are KXXXK, KXL, KXXK, LXK, LXXXXK, AK, KK, LXXXXXK, AXXXXK, KXXXXL. The ten most different k-space pairs comparing the positive data to the negative dataset 2 are KXXXK, LXK, KXL, LXXXXK, KXXK, LXXXXXK, LXXXL, KXA, AK, KXXXXA (Fig. 3). It should be noted that these results are roughly correspondent to the previous AAC and N5C5 results.

Fig. 3
figure 3

10 most different k-space pairs comparing positive to negative 1, positive to negative 2. Blue bars on the left represent the percentage of the top 10 different k-space pairs comparing positive data to negative 1 data: KXXXK, KXL, KXXK, LXK, LXXXXK, AK, KK, LXXXXXK, AXXXXK, KXXXXL. Orange bars on the right show percentage of the ten most different k-space pairs in comparison with the positive data and negative 2 data: KXXXK, LXK, KXL, LXXXXK, KXXK, LXXXXXK, LXXXL, KXA, AK, KXXXXA

Full size image

Model performance

Characteristics of the peptide data are then utilized to build the models, using machine learning methods such as SVM and SMO. In general, AAC, N5C5, k-space range from 0 to 2 and PSSM are used to build the model separately. Then, AAC, N5C5, k-space = 0 are hybridized in pairs and all together to build some other models. It should be noted that in SVM models (Tables 1, 2), the weight of each training model is tried to be adjusted from 0.1 to 10 and the one which could obtain the best accuracy is chosen as the final weight. First, SVM is utilized to analyze both the positive dataset and the negative dataset 1. As mentioned above, to improve the performance of those models, the weight of each model is adjusted. Among all those models, the one hybridized AAC and k-space = 0 as the feature with a weight of 0.9 performs the best, whose testing accuracy is 79.5% (Table 1).

Table 1 The performance of models based on both positive dataset and negative 1 dataset using SVM as classifier
Full size table
Table 2 Performance of models based on both positive dataset and negative 2 dataset using SVM as classifier
Full size table

Similarly, SVM is used to analyze both positive dataset and negative dataset 2. After adjusting the weight of each model, the one with the highest testing accuracy which hybridizes AAC and N5C5 for analysis is chosen. With a weight of 0.4, the accuracy of that model reaches 95% (Table 2). Then SMO algorithm is used to compare positive dataset to negative dataset 1. Using AAC and k-space as the representative characteristics of peptides, the accuracy of that model ranks the top one among all the models in this group at 85.5% (Table 3). Finally, models are constructed using SMO as a machine learning method and negative dataset 2 as negative data. After the evaluation of performance, the model which hybridizes AAC, N5C5 and k-space shows the highest accuracy of 95.2% (Table 4).

Table 3 Performance of models based on both positive dataset and negative 1 dataset using SMO as classifier
Full size table
Table 4 Performance of models based on both positive dataset and negative 2 dataset using SMO as classifier
Full size table

Among all models, two models with the best performance are chosen as the final models of this research: using SMO method to analyze AAC and k-space = 0 feature of the positive dataset against the negative 1 dataset (named as SMO-1) and using SMO method to analyze AAC, N5C5 and k-space = 0 feature of the positive dataset against the negative dataset 2 (named as SMO-2).

Comparison with existing ACPs prediction tools in terms of performance

To show the significance and success of those two models, the testing dataset is also applied to test the existing models, including the two AntiCP models, mACPpred, ACPred, ACPred-Fuse and ACPred-FL. Testing data from positive and negative 1 dataset are applied on SMO-1, AntiCP-1, mACPpred, ACPred, ACPred-Fuse and ACPred-FL. The model constructed in this work, SMO-1, shows the highest accuracy and MCC (Table 5). Although AntiCP-1 performs with the highest sensitivity and ACPred-Fuse performs with the highest specificity in all models, SMO-1 performs a more balanced result. Similarly, testing data from the positive and negative 2 datasets are applied on SMO-2, AntiCP-2, mACPpred, ACPred, ACPred-Fuse and ACPred-FL, which are all models utilized to distinguish ACPs from all kinds of peptides. Considering accuracy, sensitivity, specificity and MCC, SMO-2 shows the best performance among all models (Table 5). In general, SMO-2 performs with the highest accuracy comparing to existing models (Fig. 4).

Table 5 Comparison of my models and some existing tools
Full size table
Fig. 4
figure 4

Accuracy comparison of SMO-1 (highlighted in orange), SMO-2 (highlighted in orange) and mACPpred, ACPred, ACPred-Fuse and ACPred-FL

Full size image

Discussion

Some problems of conventional anticancer treatments, such as drug resistance and toxicity to other normal cells, make it necessary and urgent to discover other novel anticancer treatments [2, 3]. Among those promising treatments, anticancer peptides have attracted broad attention and interest. Due to the special structure of ACPs and its specific interaction with cancer cells, this special group of molecules can target and kill cancer cells without destroying other normal cells [4]. Before the wet-lab experiment, a computational predictive tool will definitely be helpful for the identification of ACPs. Moreover, in consideration of the similarity between ACPs and AMPs [8], it is regarded as a more efficient way of ACPs identification by searching from AMPs, because there are more examined sequence data of AMPs which could be obtained.

Nevertheless, most of the existing tools only provide a function that identifies ACPs from all kinds of peptides [13,14,15,16]. Therefore, this work creates a tool with more functions and better performance of prediction. To achieve this goal, several innovative efforts or improvements have been made. In this study, we create tools with ACPs comparing to previous studies. In total, 1492 positive and 7068 negative (4433 for negative 1 and 2635 for negative 2) data are gathered from seven different sources. Then, balanced datasets with 463 sequences in each training dataset and 100 sequences, which are independent of training data, in each testing dataset are constructed. Another improvement in this research is that new features are chosen for characterization, including N5C5, k-space and PSSM. The hybridization of some of those features greatly enhances the performance.

In the model construction step, a better machine-learning algorithm, SMO [22], is chosen and applied for classification, and increases the accuracy by 7.55% and 0.2% comparing to the SVM models. This performance suggests that SMO is a better choice than SVM in this case, and shows the success of SMO in text classification, proteomics projects, and analysis of high-dimensional data. Models built in this research is further compared with previous ACP prediction tools using an independent testing dataset. SMO-1 performs better than other tools considering accuracy and MCC value, and shows more balanced results. As for SMO-2, it performs better than all other tools in general.

Even though most of the accessible data of examined ACPs are collected in this study, the amount is still not adequate. As a result, this research may have some limitation, and could be improved in the future with more sequence data.

Conclusions

This research presents a new scheme for the identification of ACPs, including utilizing a new important feature, PSSM, and a new helpful algorithm, SMO, for optimizing SVM for classification. Also, this work offers two functions: (1) distinguishing ACPs from AMPs and (2) distinguishing ACPs from all kinds of peptides. With the help of SMO, optimized models perform better than ordinary models and other existing tools.

Methods

The process of this research is extracted and shown as a flowchart in Fig. 5. Details of the process will be explained in the following sections.

Fig. 5
figure 5

Flowchart of this work. Four major steps are involved: data preparation, feature investigation, model construction and evaluation, and comparison with existing tools. In the first step, data preparation, two datasets are constructed, which are then separated into training and testing data after CD-HIT. Then four features are investigated: amino acid composition, N5C5, k-space and position-specific scoring matrix. The third step involves model learning, cross-validation, parameter optimization and evaluation. Finally, two models proposed are compared with other existing tools

Full size image

Dataset preparation

In this research, three datasets are constructed: the positive dataset, negative dataset 1 and negative dataset 2. A positive dataset refers to anticancer peptides that are examined by experiments. They are collected from LEE dataset (total: 422) [16], Tyagi dataset (total: 450) [14], APD (total: 225) [24] and CancerPPD (total: 422) [25]. Negative dataset 1 is a collection of anti-microbial peptides without anti-cancer function. They are adapted from dbAMP dataset (total: 4057) [26] and Tyagi dataset (total: 1372). Peptides in negative dataset 2 are peptides without anti-microbial nor anti-cancer functions, which are collected from UniProt (total:281,665). Since anticancer peptides have been proved to be effective small molecules (< 50 amino acids) [27], peptides longer than 50 amino acids are removed out of datasets. Also, peptide contains artificial amino acids are removed. After this filtration step, 1492 peptide sequences in the positive dataset, 4433 peptide sequences in negative dataset 1 and 2635 peptide sequences in negative dataset 2 are obtained. To reduce identical or similar peptides sequence, CD-HIT program [28] is utilized in this research.

100% sequence-identity cut-off is applied on positive and negative 1 datasets using a Python program. Then the processed positive dataset is compared with processed negative dataset 1 using CD-HIT-2D [28]. It identifies and removes sequences in negative datasets that are similar to ones in positive dataset above a threshold of 40% (Table 6). To balance datasets, some of the peptide sequences in negative 1 dataset are removed randomly. Ultimately, both positive and negative 1 datasets have 563 peptide sequences. Each dataset is then divided randomly into two subsets, the one that contained 463 peptides is utilized as a training dataset and the other one which contained 100 peptides is used as a testing dataset. Considering the fact that normal peptides are much more abundant than ACPs in nature, the negative dataset 2 is constructed with 963 peptides collected randomly from the original 2635 peptides, and then be randomly divided into positive and negative dataset 2, with 463 and 500 in each (Table 7).

Table 6 CD-HIT results of datasets
Full size table
Table 7 Number of peptides in each dataset
Full size table

Features investigation

To utilize machine learning methods analyzing peptide sequences, features of sequences have to be extracted. In this research, 4 features are considered: amino acids composition (AAC), N5C5, k-space and PSSM.

AAC

The AAC is the proportion of each amino acid in a given peptide sequence. It summarizes the peptide information in a vector of 20 dimensions. The AAC method has been successfully and widely applied in sequence-based classifications [29,30,31,32].

N5C5

Five amino acids from both the N-terminal and C-terminal end of a given peptide are cut off and then connected as a novel sequence. Then the proportion of each amino acid in those new N5C5 sequences is calculated. Furthermore, to better analyze N5C5 sequences and visualize analysis results, heatmaps that show frequencies of each amino acid in each position are generated.

K-space

The K-space method extracts pairs of amino acids that have k (k = 0, 1, …) spacing from a given peptide sequence. In total, (N-k-1) pairs are selected from a peptide sequence which consists of N amino acids. After gathering all amino-acid-pairs, the frequency of each kind of pair is counted. To explore k-space diversity between the positive dataset and those two negative datasets, the difference value of k-space frequency in the positive dataset and that in the negative datasets is then calculated. At last, those difference values of amino-acid-pairs are sorted, and ten pairs with the highest difference values are listed.

PSSM

PSSM is generated from a group of sequences previously aligned according to structural or sequence similarity. A PSSM for a given protein is an N 20 matrix P = {Pij: i = 1… N and j = 1 … 20}, where N is the length of the protein sequence. It assigns a score Pij for the j-th amino acid in the i-th position of the query sequence. A large value indicates a highly conserved position while a small value indicates a weakly conserved position [33].

Model construction by machine learning techniques

In this study, a supervised learning technique should be applied on text data for classification. Therefore, SVM [21] is utilized in cooperation with SMO [22]. For model construction, WEKA software (version 3.8.4) [22], and packages including LIBSVM (version 3.24) [23] and SMO package (using default parameters) within WEKA are utilized.

SVM is a data-driven supervised algorithm that constructs separating hyperplanes in high-dimensional space and selects the maximum-margin one for classification [34]. Based on its solid theoretical foundations, SVM has been successfully applied in various recognition and classification studies, including text classification [35], which is utilized in this research. SVM has also been successfully and widely used for high-dimensional biological data, including examination of gene expression profiles [36], mass spectra and genomics projects [37]. Comparing to other classifiers, such as artificial neural networks, SVM shows higher accuracy, particularly when the numbers of features are large [37]. Furthermore, to improve the performance of the SVM model, a program is designed to determine the optimum value of the weight vector for each model in this research. As for adjusting gamma and cost value, a program in the LIBSVM package [38] is applied to each model.

However, SVM does have some problems, including complexity and slow training speed for large-scale data. To solve these problems, another algorithm, SMO, is also applied for classification and shows both faster speed and better performance. SMO is a new algorithm for training SVMs, which breaks large quadratic programming (QP) optimization problem, a significant obstacle in the original SVM algorithm, into a series of smallest possible QP problem. By solving those smaller QP problems analytically, a time-consuming numerical QP optimization as an inner loop could be circumvented, and thus the computational time is shortened. The SVM maximization problem is as:

$$\mathop {\max }\limits_\lambda \sum\limits_{j = 1}^m {{\lambda _j}} - \frac{1}{2}\sum\limits_{j = 1}^m {\sum\limits_{k = 1}^n {{\lambda _j}} } {\lambda _k}{y_j}{y_k}{x_j}{x_k},\quad 0 \leq {\lambda _j} \leq C,\;\quad {\forall _j},\;\sum\limits_{j = 1}^m {{y_j}} {\lambda _j} = 0$$

where λ is the Lagrange multiplier, x is the input data and y represent the class label. In SMO, two Lagrange multipliers \({\lambda_1},{\lambda_2}\) are optimized while all the other multipliers are kept constant using this equation [39]:

$${\lambda_1}{y_1} + {\lambda_2}{y_2} = - \mathop \sum \limits_{j = 3}^m {\lambda_j}{y_j} = c.$$

Moreover, since SMO only utilizes linear amount of memory, it can handle very large training sets [22], which is perfectly aligned with the need in the biological data analysis. To compute a linear SVM, only one weight vector needs to be stored. The stored weight vector can be easily updated to reflect new Lagrange multiplier values by:

$${\vec w^{{\text{new}}}} = \vec w + {y_1}\left( {\alpha_1^{new} - {\alpha_1}} \right)\overrightarrow {x_1} + {y_2}\left( {\alpha_2^{new,clipped} - {\alpha_2}} \right)\overrightarrow {x_2} l.\,[23]$$

This algorithm has shown success in some biological applications, such as metabolism studies [40, 41], genomics [42] and molecular studies [43, 44].

Parameter optimization

In the SVM model training process, the probability distribution of positive prediction and negative prediction are listed as P+ and P−. Then a weight, ranging from 0.1 to 10 with 0.1 as the interval, is multiplied to both P+ and P−. Classification is redone according to:

$$Result = \left\{ {\begin{array}{*{20}{c}} {positive,\,if\,weight \times {P^+ } \geq weight \times {P^- }} \\ {negative,\,\;if\,weight \times {P^+ } < weight \times {P^- }} \end{array}} \right.$$

In this way, the performance of the training model is changed. The weight that contributes to the highest accuracy of the training model is chosen as the final parameter and then be applied on the testing dataset, which leads to the final testing performance.

Performance evaluation

To evaluate the performance of machine learning models, four indexes are calculated: accuracy, specificity (SP), sensitivity (SN) and Matthews correlation coefficient (MCC). Details of these metrics are shown as the following equations:

$${\text{SP}} = \frac{TN}{{TN + FP}}$$
$${\text{SN}} = \frac{TP}{{TP + FN}}$$
$${\text{ACC}} = \frac{TP + TN}{{TP + TN + FP + FN}}$$
$${\text{MCC}}\frac{TP \times TN - FP \times FN}{{\sqrt {\left( {TP + FP} \right)\left( {TP + FN} \right)\left( {TN + FP} \right)\left( {TN + FN} \right)} }}$$

where TP-true positive-represents the number of correctly predicted positive labels, TN-true negative-refers to the number of corrected predict negative labels, FP-false positive-represents the number of negative labels that are wrongly predicted as positive, and FN-false negative-refers to the number of positive labels that are wrongly predicted as negative by the classifier. In addition to those evaluation metrics, the receiver operating characteristic (ROC) curve (Fig. 6) is also generated in the step of weight adjustment to visualize the relationship of true positive rate and false positive rate, and used for comparison of performance.

Fig. 6
figure 6

ROC curve for training models using SVM

Full size image

Cross-validation and independent testing sets

To test potential overfitting and evaluate the model, ten-fold cross-validation is applied in the model training step. Also, to evaluate the model built in this research and compare its performance with that of other existing tools, independent testing datasets are constructed in the dataset preparation step.

Availability of data and materials

The dataset generated and analyzed during the current study are available at: https://github.com/georginawan/proteomics_project.

Abbreviations

AAC:

Amino acid composition

ACP:

Anticancer peptides

AMP:

Antimicrobial peptides

PSSM:

Position specific scoring matrix

SMO:

Sequential minimal optimization

SVM:

Support vector machine

References

  1. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  Google Scholar 

  2. Palumbo MO, Kavan P, Miller WH Jr, Panasci L, Assouline S, Johnson N, Cohen V, Patenaude F, Pollak M, Jagoe RT, et al. Systemic cancer therapy: achievements and challenges that lie ahead. Front Pharmacol. 2013;4:57.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  3. Gatti L, Zunino F. Overview of tumor cell chemoresistance mechanisms. Methods Mol Med. 2005;111:127–48.

    CAS  PubMed  Google Scholar 

  4. Gaspar D, Veiga AS, Castanho MA. From antimicrobial to anticancer peptides. A review. Front Microbiol. 2013;4:294.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Chung CR, Jhong JH, Wang Z, Chen S, Wan Y, Horng JT, Lee TY. Characterization and identification of natural antimicrobial peptides on different organisms. Int J Mol Sci. 2020;21(3):1–26.

    Article  CAS  Google Scholar 

  6. Schweizer F. Cationic amphiphilic peptides with cancer-selective toxicity. Eur J Pharmacol. 2009;625(1–3):190–4.

    Article  CAS  Google Scholar 

  7. Riedl S, Zweytick D, Lohner K. Membrane-active host defense peptides–challenges and perspectives for the development of novel anticancer drugs. Chem Phys Lipids. 2011;164(8):766–81.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Kitagaki J, Shi GB, Miyauchi S, Murakami S, Yang YL. Cyclic depsipeptides as potential cancer therapeutics. Anticancer Drugs. 2015;26(3):259–71.

    Article  CAS  PubMed  Google Scholar 

  9. Thundimadathil J. Cancer treatment using peptides: current therapies and future prospects. J Amino Acids. 2012;2012:967347.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  10. Hoskin DW, Ramamoorthy A. Studies on anticancer activities of antimicrobial peptides. Biochim Biophys Acta. 2008;1778(2):357–75.

    Article  CAS  PubMed  Google Scholar 

  11. van Zoggel H, Carpentier G, Dos Santos C, Hamma-Kourbali Y, Courty J, Amiche M, Delbe J. Antitumor and angiostatic activities of the antimicrobial peptide dermaseptin B2. PLoS ONE. 2012;7(9):e44351.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  12. Chung CR, Kuo TR, Wu LC, Lee TY, Horng JT. Characterization and identification of antimicrobial peptides with different functional activities. Brief Bioinform. 2019;21:1098–114.

    Article  CAS  Google Scholar 

  13. Hajisharifi Z, Piryaiee M, Mohammad Beigi M, Behbahani M, Mohabatkar H. Predicting anticancer peptides with Chou’s pseudo amino acid composition and investigating their mutagenicity via Ames test. J Theor Biol. 2014;341:34–40.

    Article  CAS  PubMed  Google Scholar 

  14. Tyagi A, Kapoor P, Kumar R, Chaudhary K, Gautam A, Raghava GP. In silico models for designing and discovering novel anticancer peptides. Sci Rep. 2013;3:2984.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Chen W, Ding H, Feng P, Lin H, Chou KC. iACP: a sequence-based tool for identifying anticancer peptides. Oncotarget. 2016;7(13):16895–909.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Manavalan B, Basith S, Shin TH, Choi S, Kim MO, Lee G. MLACP: machine-learning-based prediction of anticancer peptides. Oncotarget. 2017;8(44):77121–36.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Boopathi V, Subramaniyam S, Malik A, Lee G, Manavalan B, Yang DC. mACPpred: a support vector machine-based meta-predictor for identification of anticancer peptides. Int J Mol Sci. 2019;20(8):1–14.

    Article  CAS  Google Scholar 

  18. Schaduangrat N, Nantasenamat C, Prachayasittikul V, Shoombuatong W. ACPred: a computational tool for the prediction and analysis of anticancer peptides. Molecules. 2019;24(10):1–28.

    Article  CAS  Google Scholar 

  19. Rao B, Zhou C, Zhang G, Su R, Wei L. ACPred-Fuse: fusing multi-view information improves the prediction of anticancer peptides. Brief Bioinform. 2020;21(5):1846–55.

    Article  PubMed  Google Scholar 

  20. Wei LY, Zhou C, Chen HR, Song JN, Su R. ACPred-FL: a sequence-based predictor using effective feature representation to improve the prediction of anti-cancer peptides. Bioinformatics. 2018;34(23):4007–16.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Osuna EE. Support vector machines: training and applications. Ph.D. Massachusetts Institute of Technology, Sloan School of Management; 1998.

  22. Platt JC. Fast training of support vector machines using sequential minimal optimization. Adv Kernel Methods: Support Vector Learn. 1999; 185–208. Isbn: 0-262-19416-3.

  23. Dai Y, Cai X, Shi W, Bi X, Su X, Pan M, Li H, Lin H, Huang W, Qian H. Pro-apoptotic cationic host defense peptides rich in lysine or arginine to reverse drug resistance by disrupting tumor cell membrane. Amino Acids. 2017;49(9):1601–10.

    Article  CAS  PubMed  Google Scholar 

  24. Wang G, Li X, Wang Z. APD3: the antimicrobial peptide database as a tool for research and education. Nucl Acids Res. 2016;44(D1):D1087-1093.

    Article  CAS  PubMed  Google Scholar 

  25. Tyagi A, Tuknait A, Anand P, Gupta S, Sharma M, Mathur D, Joshi A, Singh S, Gautam A, Raghava GP. CancerPPD: a database of anticancer peptides and proteins. Nucl Acids Res. 2015;43(Database issue):D837–43.

    Article  CAS  PubMed  Google Scholar 

  26. Jhong JH, Chi YH, Li WC, Lin TH, Huang KY, Lee TY. dbAMP: an integrated resource for exploring antimicrobial peptides with functional activities and physicochemical properties on transcriptome and proteome data. Nucl Acids Res. 2019;47(D1):D285–97.

    Article  CAS  PubMed  Google Scholar 

  27. Li FM, Wang XQ. Identifying anticancer peptides by using improved hybrid compositions. Sci Rep. 2016;6:33910.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Huang Y, Niu B, Gao Y, Fu L, Li W. CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics. 2010;26(5):680–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Usmani SS, Bhalla S, Raghava GPS. Prediction of antitubercular peptides from sequence information using ensemble classifier and hybrid features. Front Pharmacol. 2018;9:954.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  30. Kao HJ, Nguyen VN, Huang KY, Chang WC, Lee TY. SuccSite: incorporating amino acid composition and informative k-spaced amino acid pairs to identify protein succinylation sites. Genomics Proteomics Bioinform. 2020;18(2):208–19.

    Article  Google Scholar 

  31. Huang CH, Su MG, Kao HJ, Jhong JH, Weng SL, Lee TY. UbiSite: incorporating two-layered machine learning method with substrate motifs to predict ubiquitin-conjugation site on lysines. BMC Syst Biol. 2016;10(Suppl 1):6.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Chen SA, Lee TY, Ou YY. Incorporating significant amino acid pairs to identify O-linked glycosylation sites on transmembrane proteins and non-transmembrane proteins. BMC Bioinform. 2010;11:536.

    Article  CAS  Google Scholar 

  33. Gribskov M, McLachlan AD, Eisenberg D. Profile analysis: detection of distantly related proteins. Proc Natl Acad Sci USA. 1987;84(13):4355–8.

    Article  CAS  PubMed  Google Scholar 

  34. Noble WS. What is a support vector machine? Nat Biotechnol. 2006;24(12):1565–7.

    Article  CAS  PubMed  Google Scholar 

  35. Tong S, Koller D. Support vector machine active learning with applications to text classification. J Mach Learn Res. 2002;2(1):45–66.

    Google Scholar 

  36. Golub TR, Slonim DK, Tamayo P, Huard C, Gaasenbeek M, Mesirov JP, Coller H, Loh ML, Downing JR, Caligiuri MA, et al. Molecular classification of cancer: class discovery and class prediction by gene expression monitoring. Science. 1999;286(5439):531–7.

    Article  CAS  PubMed  Google Scholar 

  37. Byvatov E, Schneider G. Support vector machine applications in bioinformatics. Appl Bioinform. 2003;2(2):67–77.

    Google Scholar 

  38. Chih-Chung C, Chih-Jen L. LIBSVM : a library for support vector machines. ACM Trans Intell Syst Technol. 2011;2(3):27:1–27.

    Google Scholar 

  39. Naveed H, Khan G, Khan AU, Siddiqi A, Khan MUG. Human activity recognition using mixture of heterogeneous features and sequential minimal optimization. Int J Mach Learn Cybern. 2019;10(9):2329–40.

    Article  Google Scholar 

  40. Chen L, Chu C, Feng K. Predicting the types of metabolic pathway of compounds using molecular fragments and sequential minimal optimization. Comb Chem High Throughput Screen. 2016;19(2):136–43.

    Article  CAS  PubMed  Google Scholar 

  41. Lee TY, Hsu JB, Chang WC, Wang TY, Hsu PC, Huang HD. A comprehensive resource for integrating and displaying protein post-translational modifications. BMC Res Notes. 2009;2:111.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Huang T, Shu Y, Cai YD. Genetic differences among ethnic groups. BMC Genomics. 2015;16:1–10.

    Article  CAS  Google Scholar 

  43. Periwal V, Rajappan JK, Jaleel AU, et al. Predictive models for anti-tubercular molecules using machine learning on high-throughput biological screening datasets. BMC Res Notes. 2011;4:1–10.

    Article  CAS  Google Scholar 

  44. Huang KY, Lee TY, Kao HJ, Ma CT, Lee CC, Lin TH, Chang WC, Huang HD. dbPTM in 2019: exploring disease association and cross-talk of post-translational modifications. Nucl Acids Res. 2019;47(D1):D298–308.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

I would like to thank the support from Undergraduate Research Award, CUHKSZ.

Funding

This work was supported by the Warshel Institute for Computational Biology, School of Life and Health Sciences, The Chinese University of Hong Kong, Shenzhen, China and National Natural Science Foundation of China (Grant No. 32070659). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations

  1. School of Life and Health Sciences, The Chinese University of Hong Kong, Shenzhen, Shenzhen, 518172, Guangdong, People’s Republic of China

    Yu Wan & Tzong-Yi Lee

  2. Warshel Institute for Computational Biology, The Chinese University of Hong Kong, Shenzhen, Shenzhen, 518172, Guangdong, People’s Republic of China

    Zhuo Wang & Tzong-Yi Lee

Authors
  1. Yu Wan
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Zhuo Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Tzong-Yi Lee
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

TYL contributed to the study concept and design. YW carried out the data analysis and drafted the manuscript. ZW and YW helped conceive the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Tzong-Yi Lee.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wan, Y., Wang, Z. & Lee, TY. Incorporating support vector machine with sequential minimal optimization to identify anticancer peptides. BMC Bioinformatics 22, 286 (2021). https://doi.org/10.1186/s12859-021-03965-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12859-021-03965-4

Keywords

BMC Bioinformatics

ISSN: 1471-2105

Contact us