Next Article in Journal
Microbial Corrosion Resistance and Antibacterial Property of Electrodeposited Zn–Ni–Chitosan Coatings
Next Article in Special Issue
Novel pH Selective, Highly Lytic Peptides Based on a Chimeric Influenza Hemagglutinin Peptide/Cell Penetrating Peptide Motif
Previous Article in Journal
The Mutation of Conservative Asp268 Residue in the Peptidoglycan-Associated Domain of the OmpA Protein Affects Multiple Acinetobacter baumannii Virulence Characteristics
Previous Article in Special Issue
The Characteristics of PD-L1 Inhibitors, from Peptides to Small Molecules
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

ACPred: A Computational Tool for the Prediction and Analysis of Anticancer Peptides

by
Nalini Schaduangrat
1,
Chanin Nantasenamat
1,
Virapong Prachayasittikul
2 and
Watshara Shoombuatong
1,*
1
Center of Data Mining and Biomedical Informatics, Faculty of Medical Technology, Mahidol University, Bangkok 10700, Thailand
2
Department of Clinical Microbiology and Applied Technology, Faculty of Medical Technology, Mahidol University, Bangkok 10700, Thailand
*
Author to whom correspondence should be addressed.
Molecules 2019, 24(10), 1973; https://doi.org/10.3390/molecules24101973
Submission received: 1 April 2019 / Revised: 7 May 2019 / Accepted: 17 May 2019 / Published: 22 May 2019

Abstract

:
Anticancer peptides (ACPs) have emerged as a new class of therapeutic agent for cancer treatment due to their lower toxicity as well as greater efficacy, selectivity and specificity when compared to conventional small molecule drugs. However, the experimental identification of ACPs still remains a time-consuming and expensive endeavor. Therefore, it is desirable to develop and improve upon existing computational models for predicting and characterizing ACPs. In this study, we present a bioinformatics tool called the ACPred, which is an interpretable tool for the prediction and characterization of the anticancer activities of peptides. ACPred was developed by utilizing powerful machine learning models (support vector machine and random forest) and various classes of peptide features. It was observed by a jackknife cross-validation test that ACPred can achieve an overall accuracy of 95.61% in identifying ACPs. In addition, analysis revealed the following distinguishing characteristics that ACPs possess: (i) hydrophobic residue enhances the cationic properties of α-helical ACPs resulting in better cell penetration; (ii) the amphipathic nature of the α-helical structure plays a crucial role in its mechanism of cytotoxicity; and (iii) the formation of disulfide bridges on β-sheets is vital for structural maintenance which correlates with its ability to kill cancer cells. Finally, for the convenience of experimental scientists, the ACPred web server was established and made freely available online.

Graphical Abstract

1. Introduction

Cancer remains one of the leading causes of death worldwide, accounting for an estimated 18.1 million new cases and resulting in 9.6 million deaths in 2018 [1]. Of these, the most common cancer types causing mortality worldwide include lung, colorectal, stomach, liver and breast cancer, respectively [1,2]. The global burden of cancer has risen dramatically in the last few decades thereby making the prevention of cancer the most significant and challenging public health concern of the 21st century [3]. In spite of the many advanced clinical methods available for the treatment of cancer (e.g., chemotherapy, radiation therapy, hormonal therapy, etc.) their recurrence rate still remains high. In addition, anti-cancer drugs are known to exert their therapeutic effects by killing off normal cells and tissues, which exacerbates the immunodeficiency of patients. Thus, the discovery and development of novel anti-cancer drugs are crucial for decreasing premature death as well as for increasing the survival of many populations. In this regard, peptide-based therapeutics have gained much attention as a new drug class since they are considered to be relatively safe, highly selective, possess good tolerability, lower production cost, easy to modify and synthesize while also exhibiting promising pharmacological profiles [4,5,6,7]. In the last decade, more than 7000 naturally occurring peptides have been discovered and they are found to exhibit a wide range of bioactivities including tumor homing, antihypertensive, antiparasitic, antiviral, antiangiogenic, antibiofilm, antimicrobial and anticancer [8,9,10]. Currently, 60 peptide-based drugs have been FDA-approved [11] with preclinical and clinical studies in the pipeline for another 150 peptides [6].
Anticancer peptides (ACPs) are a subset of antimicrobial peptides (AMPs) that are have been shown to exhibit anticancer activities. These peptides are typically made up of short (5–30 amino acids) cationic residues that mostly adopts the ɑ-helical structure (e.g., BMAP-27, BMAP-28, Cecropin A and B, LL-37, Magainins etc.) or folds into a β-sheet (e.g., Defensins, Lactoferricin and Tachylepsin I) while extended linear forms (e.g., Tritrpticin and Indolicidin) have also been reported [12,13], as shown in Figure 1.
In comparison to normal cells, the cancer cell membrane exhibits different properties (i.e., larger surface area due to the higher number of microvili, high negative charge of the cell membrane, higher membrane fluidity, etc.) [14,15,16,17]. Thus, the inherent cationic property of ACPs is responsible for the electrostatic interactions between ACPs and the anionic cell membrane of cancer cells, and hence allow the selective killing of cancer cells [18]. Previously, many studies have reported that ACPs display several advantages over small molecule drugs. For instance, they exhibit a shorter half-life that consequently leads to decreased probability of resistance, lowered toxicity (not harming normal cells as much), higher specificity, good solubility as well as good tumor penetration capabilities thereby indicating their great potential in cancer therapy [19,20,21,22].
Accurate prediction of ACPs is of great importance for the exploration of their mechanism of action and for the development of therapeutic ACPs. As the experimental identification and development of novel ACPs is a very time-consuming and labor-intensive process, therefore bioinformatics tools are needed for effective analysis of the available big data on existing peptides so as to allow the identification of novel ACPs while also shedding light about their mechanism of action.
Recently, many sequence-based computational methods, for example, AntiCP [23], Hajisharifi et al.’s method [24], ACPP [25], iACP [26], Li and Wang’s method [27], iACP-GAEnsC [28], MLACP [29], SAP [30] and TargetACP [31], have been developed using a wide range of machine learning methods and peptide features as summarized in Table 1. In 2014, Hajisharifi et al. [24] established a benchmark dataset spanning 138 experimentally confirmed ACPs (positive dataset) and 227 non-AMPs (negative dataset). In their method, support vector machine (SVM) with pseudo amino acid composition were used for model building and an accuracy of 92.68% were attained. As noticed in Table 1, iACP [26] was the first computational method based on informative features. In this method, Chen et al. utilized a feature selection technique for obtaining an informative feature set subsequently followed by encoding them for constructing an ACP predictor known as iACP [26] via the use of SVM. This method yielded better prediction accuracy (95.06%) than the two previous methods [23,24]. Thusfar, state-of-the-art ACP predictors includes iACP-GAEnsC [28] and TargetACP [31] in which both afforded high prediction accuracies of 96.45% and 98.78%, respectively. The iACP-GAEnsC [28] method follows the concept of evolutionary intelligent genetic algorithm-based ensemble model for improving the true classification rate. iACP-GAEnsC was developed using various classification algorithms (e.g., random forest, k-nearest neighbor, support vector machine, generalized neural network and probabilistic neural network) and three types of peptide features (e.g., pseudo g-gap dipeptide composition, amphiphilic pseudo amino acid composition and reduced amino acid alphabet). More recently, Kabir et al. developed the TargetACP [31] for solving the problem of class imbalance presented in the previously described benchmark dataset [24] by integrating sequential and evolutionary-profile information as input features and using SVM as the classifier. In addition, the authors utilized the synthetic minority oversampling technique (SMOTE) to solve the imbalance phenomenon between minority (ACPs) and majority (non-ACPs) samples.
Although, iACP-GAEnsC and TargetACP yielded encouraging performances with reasonably high prediction accuracies, the overall utility of these two methods is limited in terms of interpretability and practical utility thereby warranting further improvements [32]. One of the major limitations of iACP-GAEnsC and TargetACP was that there was no web server provided for these two methods. Therefore, their utility is limited to researchers with informatics background who can develop in-house prediction models. In addition, the underlying mechanisms of the investigated bioactivity for these two methods affords limited interpretability for experimental scientists. Owing to the heterogeneous computational architecture and design of the study as well as differences in the benchmark dataset used, it is not easy to identify and assess which methods and features provided the greatest contribution to good prediction. Motivated by aforementioned issues presented in previous models [32,33,34,35,36], in order to establish an interpretable ACP predictor we propose a systematic effort for the prediction and analysis of anticancer activities of peptides called the ACPred. Firstly, ACPred was built using various types of peptide features and two powerful machine learning algorithms including random forest (RF) and support vector machine (SVM). Secondly, insights into the ACP mechanism of action were attained via the identified important interpretable features and subsequent in-depth analysis pertaining to the underlying biophysical and biochemical properties of anticancer activities of peptides [33,34,36,37,38]. Thirdly, if-then interpretable rules were extracted using the RF model (IR-ACP). Our prediction results on the benchmark dataset showed that the proposed ACPred achieved an accuracy and MCC of 95.61% and 0.91, respectively, which outperforms most of the existing sequence-based computational methods including Hajisharifi et al.’s method [24], iACP [26], Li and Wang’s method [27], iACP-GAEnsC [28], MLACP [29] and SAP [30]. As previously noted, although iACP-GAEnsC and TargetACP boasts good predictive performances but their interpretability and practical utility prompted the development of ACPred. Thus, ACPred is a user-friendly and publicly accessible web server that allow robust predictions to be made without the need to develop in-house prediction models.

2. Materials and Methods

The overall framework of ACPred is illustrated in Figure 2. It can be seen that four major steps are involved in the development of this method as follows: (i) data preparation, (ii) model construction (feature extraction and feature combination), (iii) feature importance analysis and (iv) web server construction. The first step prepares the benchmark dataset, the second step extracts and combines peptide features and the third step analyzes the important features governing ACP prediction via the RF model. Finally, the fourth step involves the construction of a publicly accessible web server using the best model.

2.1. Benchmark Dataset

A high-quality benchmark dataset can ensure the reliability of the proposed method. Therefore, the investigated dataset was retrieved from the work of Hajisharifi et al. [24], which contains 192 ACPs and 215 non-ACP. As to guarantee the quality of the benchmark dataset, the following steps were considered. Firstly, in order to avoid a dataset containing redundant peptides, the CD-HIT program was used to eliminate ACPs and non-ACPs with >90% identity. Secondly, ACPs and non-ACP containing special characters, such as X and U, were removed. After such screening procedures, a set of 138 ACPs and 205 non-ACPs was attained for prediction and analysis of ACPs.

2.2. Feature Representation

In order to build an effective prediction model, it is necessary to represent each peptide as numerical vectors that can encompass the perspective of its biophysical and biochemical properties. For peptide or protein sequences, the most classical and interpretable features are amino acid composition (AAC), dipeptide composition (DPC) and physicochemical property (PCP).
According to the classical definition, AAC and DPC are the proportions of each amino acid and dipeptide in a peptide sequence P that are expressed as fixed lengths of 20 and 400, respectively. Thus, in terms of AAC and DPC features, a peptide P can be expressed by vectors with 20D and 400D (dimension) spaces, respectively, as formulated by the following two equations:
P = [ a a 1 , a a 2 , ,   a a 20 ] T
P = [ d p 1 , d p 2 , ,   d p 400 ] T
where T is the transpose operator, while aa1, aa2, …, aa20 and dp1, dp2, …, dp400 are occurrence frequencies of the 20 and 400 native amino acids and dipeptides, respectively, in a peptide sequence P. PCP is one of the most intuitive features associated with biophysical and biochemical reactions. In fact, there are 544 PCPs of amino acids extracted from the amino acid index database (AAindex) [39], which is a collection of published literature as well as different biochemical and biophysical properties of amino acids. Each physicochemical property consists of a set of 20 numerical values for amino acids. After removing 13 PCPs with not applicable (NA) as their amino acid indices, a total of 531 PCPs were further used in this study.
As mentioned in the study of [40,41,42], the sequence order information of AAC and DPC would be lost. To deal with such a dilemma, the pseudo amino acid composition (PseAAC) and amphiphilic pseudo amino acid composition (Am-PseAAC) approaches were proposed. According to Chou’s PseAAC [41], the general form of PseAAC for a peptide P is formulated by: where the subscript Ω is an integer to reflect the feature’s dimension. The value of Ω and the component of Ψ u , where u = 1 , 2 , , Ω depends on the protein or peptide sequences. In this study, the parameters of PseAAC, i.e., the discrete correlation factor λ and weight of the sequence information ϖ , were set by default in the protr R package [43]. The dimension of PseAAC feature is 20 + λ × ϖ . Since the hydrophobic and hydrophilic properties of proteins play an important role in its folding and interaction of proteins, Am-PseAAC was introduced by Chou [40]. The dimension of Am-PseAAC feature is 20 + 2 λ whereby, the first 20 components are the 20 basic AAC ( p 1 , p 2 , ,   p 20 ) and the next 2 λ ones are a set of correlation factors that reveal the physicochemical properties such as hydrophobicity and hydrophilicity along a protein or peptide sequence as formulated by:
P = [ Ψ 1 , Ψ 2 , , Ψ u , ,   Ψ Ω ] T
where the subscript Ω is an integer to reflect the feature’s dimension. The value of Ω and the component of Ψ u , where u = 1 , 2 , , Ω depend on the protein or peptide sequences. In this study, the parameters of PseAAC, i.e., the discrete correlation factor λ and weight of the sequence information ϖ , were set by defaults using the protr R package [43]. The dimension of PseAAC feature is 20 + λ × ϖ . Since hydrophobic and hydrophilic properties play an important role in the folding and interaction of proteins, Am-PseAAC was introduced by Chou [40] for representing such features. The dimension of Am-PseAAC feature is 20 + 2 λ whereby the first 20 components are the 20 basic AAC ( p 1 , p 2 , ,   p 20 ) while the next 2 λ ones are a set of correlation factors that pertains to physicochemical properties namely hydrophobicity and hydrophilicity along a protein or peptide sequence as formulated by:
P = [ p 1 , p 2 , ,   p 20 , p 20 + λ , p 20 + λ + 1 , p 20 + 2 λ ] T
In this study, the five aforementioned features of peptide sequences were calculated using the protr package in the R software [43].

2.3. Dataset Modelability

The prediction performances of QSAR models are directly influenced by the dataset and the types of feature used. Recently, the concept of dataset modelability has been proposed by Golbraikh et al. [32,35,44,45] for evaluating the feasibility of developing a robust and predictive model for a given dataset of interest. In this study, the MODI score for each type of peptide features and their combinations were calculated using an in-house developed R code. Prediction models are recommended for further analysis if the MODI score is greater than 0.65; otherwise, prediction models are not recommended. Further details pertaining to the calculation of the MODI score can be found elsewhere [44].

2.4. Classifier Selection

In this study, we employed two successful machine learning algorithms namely random forest (RF) and support vector machine (SVM). Previously, these approaches have been successfully used in the prediction of various functions and properties of peptides and proteins [33,36,38,46,47,48,49,50,51] as well as other biological or chemical entities [35,45,52,53,54,55]. Herein, the basic concepts and associated parameter optimizations for the two classifiers are briefly described hereafter.
RF models were constructed according to the described original RF algorithm [56]. These models were developed by growing many weak classification and regression tree (CART) classifiers where each classifier is generated using a random vector sampled independently from the input vector as to enhance the prediction performance of CART [56,57]. In the RF method, the out-of-bag (OOB) approach is used for evaluating the feature importance as follows: (1) two-thirds of the training data was utilized to construct the predictive classifier while the remaining was used for evaluating the performance of such classifier and (2) the feature importance of each feature can be evaluated by measuring the decrease in the prediction performance. It should be noted that the performance evaluation of the model can be either accuracy or Gini index. Herein, the RF classifier was established using the randomForest package in the R software [58]. To enhance the performance of the RF model, two parameters namely, ntree (i.e., the number of trees used for constructing the RF classifier) and mtry (i.e., the number of random candidate features) were determined using the caret package in the R software [59] via a 5-fold cross-validation (5-fold CV) approach. The search space of ntree and mtry were [100, 500] and [1, 10] using steps of 100 and 1, respectively.
SVM is a supervised learning model based on the principles of structure risk minimization and kernel method as proposed by Vapnik [60,61,62]. SVM model can deal with the problem of over-fitting arising with the use of a small training dataset by mapping the input samples to a higher dimension space and then searching for the maximum-margin hyperplane for constructing the classifier [63,64]. Previously, SVM models were used in many applications because of their predictive performance capabilities when applied to classification, prediction, and regression problems. In the optimization process, the regularization parameter C { 0.25 ,   0.50 ,   1 ,   2 ,   4 } was determined with a 5-fold CV approach using the caret package in the R software [59].

2.5. Identification of Important Features

The analysis and identification of feature importance can provide a better understanding on the underlying biophysical and biochemical properties governing anticancer activities of peptides. Herein, the efficient and effective built-in feature importance estimator of the RF method was used to reveal and characterize differences between ACPs and non-ACPs. As mentioned above, the RF method provides two measures for ranking feature importance, i.e., the mean decrease of the Gini index (MDGI) and the mean decrease of prediction accuracy. Since Calle and Urrea [65] demonstrated that the MDGI provided more robust results when compared with the mean decrease of prediction accuracy, we utilized the MDGI to rank the importance of interpretable features that included AAC, DPC and PCP. Until now, these three features have been used to characterize many peptides and proteins such as for predicting HIV-1 CRF01-AE co-receptor usage [46], predicting protein crystallization [51,66], predicting the oligomeric states of fluorescent proteins [38], predicting the bioactivity of host defense peptides [47], prediction of human leukocyte antigen gene [37,49], predicting antifreeze proteins [33], predicting the hemolytic activity of peptides [34] and predicting the antihypertensive activity of peptides [36].
The MDGI is an impurity measure that corresponds to the ability of each feature in discriminating the sample classes. The Gini index can be defined as follows:
1 c = 1 2 p 2 ( c | t )
where denotes the estimated class probability for node t in a tree classifier and c is the class label (i.e., either ACPs or non-ACPs). Features with the largest MDGI value is considered to be an important feature as it significantly contributes to the prediction performance.

2.6. Rule Extraction

This work presented an interpretable rule extraction of ACPs (IR-ACP) based on the RF method cooperating with AAC features for determining how the features work in combination within the learning process. A set of rules from an individual tree is derived from the root to the leaves. In this study, only 100 decision trees were used to extract if-then interpretable rules for explaining the prediction results by means of the RRF, inTrees and xtable packages in the R software [56,57,67]. More details of the rule extraction process can be found in previous related works [67,68].

2.7. Performance Evaluation

Three testing methods consisting of (i) sub-sampling test (2-, 5- or 10-fold cross-validation), (ii) jackknife test or also known as leave-one-out cross-validation (LOO-CV), and (iii) independent (or external) testing dataset are often used to evaluate the prediction performance. The sub-sampling test is one of the popular cross-validation methods for assess the predictive capability of a model. As elucidated in [41,69] and demonstrated by Equation (50) of [70], among the three cross-validation methods, the jackknife test is considered as one of the most rigorous and objective methods for cross-validation in statistics and can provide a unique result for a given benchmark dataset. Thus, in this study, the 5-fold cross-validation (5-fold CV) and jackknife test were used to evaluate the prediction performance of our models. The former set (5-fold CV) makes use of data from the training set where the data set is separated into five subsets. Practically, one subset from a total of five subsets is left out as the testing set while the remaining are used for training the model. This process is repeated iteratively until all data samples have had the chance to be left out as the testing set. During the jackknifing process, each sample in the dataset (N samples) is left out as an independent sample (one sample) and the remaining sequences (N-1 samples) are used for training the model. This process is repeated iteratively until all samples have had the chance to be left out as the independent sample.
In order to evaluate the prediction ability of the model, the following sets of four metrics are used:
A c = T P + T N ( T P + T N + F P + F N )
S n = T P ( T P + F N )
S p = T N ( T N + F P )
M C C = T P × T N F P × F N ( T P + F P ) ( T P + F N ) ( T N + F P ) ( T N + F N )
where Ac, Sn, Sp and MCC represent accuracy, sensitivity, specificity and Matthews coefficient correlation, respectively. TP, TN, FP, and FN represent the instances of true positive, true negative, false positive and false negative, respectively. Moreover, in order to evaluate the prediction performance of models using threshold-independent parameters, receiver operating characteristic (ROC) curves were plotted using the pROC package in the R software [71]. The area under the ROC curve (auAUC) was used to measure the prediction performance, where AUC values of 0.5 and 1 are indicative of perfect and random models, respectively.

2.8. Reproducible Research

To ensure the reproducibility of models proposed herein, all R codes and the benchmark dataset used in the construction of the predictive models, graphical figures and the ACPred web server are available on GitHub at https://github.com/Shoombuatong2527/acpred and https://github.com/chaninlab/acpred-webserver.

3. Results and Discussion

3.1. Prediction Performance

In this study, effort is directed toward the determination of feature types that are beneficial for the prediction of APCs. This was performed by comparing the performance of different feature types as constructed using SVM and RF by means of 5-fold CV and jackknife test. The five basic features consisting of AAC, DPC, PCP, PseAAC and Am-PseAAC were considered as well as their combinations: AAC + PseACC, AAC + Am-PseACC, PseACC + Am-PseACC and AAC + PseACC + Am-PseACC, which were also utilized to investigate the complementarity of the five aforementioned features.
Prior to the development of the predictive model, the assessment of modelability of the benchmark dataset was performed by computing the MODI index. This index helps modelers to estimate the feasibility of obtaining robust and reliable predictive models. For the binary classification problem, if the value of MODI index is greater than 0.65, the feature is considered to be reliable for classification modelling; otherwise, the feature is not recommended for classification modelling.
As shown in Table A1, all feature types and their combinations from the benchmark dataset met these criteria. Thus, it could be stated that our proposed features are reliable and efficient for constructing ACP predictor.
Table 2 and Table 3 shows the performance comparison amongst different feature types with RF and SVM on the benchmark dataset via 5-fold CV and jackknife test, respectively. In addition, Figure 3 shows the receiver operating characteristic (ROC) curve of the five considered feature sets obtained from RF (top and bottom left) and SVM (top and bottom right) as evaluated by 5-fold CV (top left and right) and jackknife test (bottom left and right). As seen in Table 2 and Table 3, SVM model with the Am-PseAAC feature afforded the highest accuracy of 95.03% as evaluated by the jackknife test. Meanwhile, RF model with PseACC feature and SVM model with AAC feature performed well with the second and third the highest accuracy of 93.28% and 92.98%, respectively. These results showed that ACC, PseACC and Am-PseAAC were effective for ACP prediction. Since using a combination of various features might yield better prediction performance, four combinations of the three effective features were also considered. The highest accuracy and MCC of 95.61% and 0.91, respectively, was achieved by using SVM model cooperating with the combination feature of AAC and Am-PseAAC, while using the combination of the three effective features, i.e., ACC, PseACC and Am-PseAAC, performed slightly worse than the combination feature of AAC and Am-PseAAC with an accuracy and MCC of 94.74% and 0.89, respectively. In addition, our prediction results were well consistent with previous studies [27,28,29] and related studies [33,34,36,72].
Furthermore, observations pertaining to the performance comparisons from Table 2 and Table 3 and Figure 3, it can be briefly summarized as follows. Table 2 and Table 3 shows that both RF and SVM models afforded improved prediction performances when using the Am-PseACC feature while the best prediction performance on the benchmark dataset via both 5-fold CV and jackknife test were achieved by the SVM model trained using the combination of AAC and Am-PseAAC features.
For the convenience of subsequent descriptions, we will refer to this method as ACPred as it represents the best model that will be used for further comparisons with other tools.

3.2. Comparison with Other Methods

To indicate the effectiveness of the proposed method, we compared the performance of our selected model (named ACPred) with other popular ACP predictors. Since Hajisharifi et al. [24] were the first group to establish the benchmark dataset, and iACP [26], iACP-GAEnsC [28] and TargetACP [31] were state-of-the-art ACP predictors that provided prediction results as assessed by the jackknife test, therefore these predictors were used for performance comparisons. Table 4 lists the performance comparisons of ACPred with the three ACP predictors. As noticed from Table 4, prediction results were obtained from two different experimental designs: (i) the prediction model was performed on the benchmark dataset as derived from the work of [24] as well as those of Hajisharifi et al.’s method, iACP, iACP-GAEnsC, TargetACP and ACPred, (ii) the prediction model was built using the balanced dataset including TargetACP and ACPred. The reported results from aforementioned methods as summarized in Table 4 was obtained directly from the work of iACP-GAEnsC [28].
As for performance comparisons on the benchmark dataset, it was found that ACPred was comparable with that of iACP-GAEnsC and TargetACP as indicated by four statistical parameters. Moreover, ACPred yielded a greater prediction performance than the method proposed by Hajisharifi et al. and iACP. To the best of our knowledge, iACP-GAEnsC is the only method that utilizes the synthetic minority oversampling technique (SMOTE) [73] technique for constructing a balanced dataset from the original benchmark dataset [24]. From the perspectives of machine learning, it is not fair to directly compare prediction results between TargetACP and ACPred because TargetACP was trained and tested on the balanced dataset consisting of 205 ACPs and 205 non-ACPs. In order to perform a fair comparison, the herein proposed ACPred method was applied on the balanced dataset (called ACPred-modified), which was generated by means of the SMOTE technique. Herein, the oversampling method was used to add synthetic samples for the minority class (i.e., ACPs) presented in the benchmark dataset by setting the parameters of the number of nearest neighbours (k), the number of extra cases (oversampling) to add to the minority class (perc.over) and the number of cases to reduce (undersampling) from the majority class (perc.under) were set to 9, 50 and 300, respectively. As seen in Table 4, ACPred-modified was found to outperform Hajisharifi et al.’s method, iACP and iACP-GAEnsC with an improvement of 2%–4% and 3%–7% on Ac and MCC, respectively. However, it was observed that ACPred-modified performed slightly worse than TargetACP by approximately 2%–3% where ACPred-modified and TargetACP provided 97.56% Ac/0.95 MCC and 98.78% Ac/0.97 MCC, respectively. Nevertheless, TargetACP was constructed with a sophisticated design thus, it is not easy to identify and assess which features offer the most contribution to the prediction improvement. On the other hand, ACPred was designed in a systematic manner for prediction and characterization of anticancer activities of peptides. In addition, a user-friendly web server ACPred was developed to facilitate high-throughput prediction of ACP. Therefore, the proposed ACPred model could become a practical tool for predicting and interpreting the anticancer activity of peptides, or at least as a complementary tool to existing methods in the field.
Thinking to the future, some possible improvements that could enhance the prediction performance for predicting anticancer activity of peptides are described hereafter. Firstly, it is worthy to explore the separation of peptides according to their sequence lengths (e.g., peptides having 4, 5 or 6 amino acids, peptides having 7–12 amino acids, and peptides having more than 12 amino acids) followed by the development of separate predictive models for each of the five sequence range [36,74]. Secondly, extract 5, 10 and 15 amino acids from the N or C terminus for the development of predictive QSAR models, which has successfully been demonstrated by the group of P.S. Raghava [23] Thirdly, explore the modelability of the dataset using the approach of Golbraikh et al. [44] such that robust models could be developed. Fourthly, the currently employed non-anticancer peptides [24] are not based on true experimentally-determined inactive peptides but are random peptides obtained from the UniProt database therefore, if such true inactive peptides could be determined they may aid in the development of robust models.
Aside from the aforementioned methodological improvements that have been witnessed from the literature in the development of anticancer peptide classifiers, however, the underlying origin of anticancer activity of investigated peptides as rationalized by peptide features is an area that deserves more attention. Current predictive models used for the prediction of anticancer activity of peptides are primarily based on sequence order-independent descriptors that does not consider the order of amino acids in a peptide sequence (e.g., AAC and DPC), which hampers the identification of motifs from peptides. Such limitation has been addressed by the group of KC Chou [41,75,76] in which they developed PseAAC and Am-PseAAC descriptors that take into account sequence order information based on physicochemical properties. However, the interpretability of these descriptors is limited and thus, existing approaches for rationalizing biological interpretation from these predictive models are resorted back to the sequence order-independent descriptors such as AAC and DPC. Analysis of the literature revealed that interpretable descriptors such as the z-scale descriptors [77,78] have been successful for predicting robust and interpretable models, however the flaw of this approach is that the peptides used in the development of the model should be of equal length so as to conduct sequence alignment. Methodological tweaks addressing such limitation have been made via the use of autocorrelation such that sequence alignment is not needed and would allow the models to be built on peptides of varied lengths. In spite of this, this method is not applicable in this study field in which the prepared dataset removes peptides having high similarity of greater than 90% therefore, rendering the removal of a large number of peptides that is normally needed when using z-scale descriptors. As such, there is ample room for the discovery of new and interpretable descriptors for the development of predictive models of peptides, which is an important area for further driving the field of anticancer peptide prediction.

3.3. Biological Space

The identification of feature importance from AAC and DPC descriptors can provide a better understanding on the biochemical and biophysical properties of anticancer activities of peptides. Previously, AAC and DPC features have been analyzed as to further gain insights on how to characterize therapeutic peptides [29,34,36,47] and protein functions [33,51,79]. In this study, the value of MDGI was adopted to rank and estimate the importance of each AAC and DPC feature. Such information is derived from the analysis of the benchmark dataset that consists of 138 ACPs and 205 non-ACPs. Table 5 lists the percentage values of the 20 amino acids for both ACPs and non-ACPs along with amino acid composition difference between the two classes as well as their MDGI values. In addition, a heatmap showing the feature importance for DPC features is shown in Figure 4. From Table 5, it can be observed that the ten informative amino acids with the highest MDGI values are Lys, Cys, Gln, Arg, Ile, Leu, His, Pro, Glu and Gly (29.54, 19.71, 15.89, 11.42, 10.13, 7.84, 7.80, 7.37, 6.88 and 6.56, respectively). Meanwhile, Figure 4 shows that the 20 top-ranked dipeptides according to their MDGI value are KK, KI, CY, IK, AK, KV, IP, LK, TC, SC, CG, LF, FK, GL, FL, LR, PP, LG, IG and AQ. Interestingly, among the 5 top-ranked informative amino acids, 3 of these are polar (i.e., Lys, Gln and Arg), while 2 are non-polar (i.e., Cys and Ile) and hydrophobic residues. In addition, Lys and Arg are also positively-charged residues, which may support the anti-cancer cell penetrating properties of peptides (i.e., that is required for targeting tumors with specificity and low toxicity) [80]. Although the percentage compositions differentiating ACPs from non-ACPs as presented in Table 5 highlight that ACPs a have higher percentage of polar residues and lower percentage of basic residues as compared to non-ACPs, cell surface binding and internalization are critical for specific targeting of cells with anti-cancer activities. According to the membrane of cancer cells presenting anionic molecules such as phosphatidyl serine (PS), heparin sulfate and O-glycosylated mucins, the ACPs that contain positive charge are critical for endocytosis and selectively killing cancer cells. While hydrogen bonds that appear on ACPs were designed for improving the solubilization of hydrophobic molecules on ACPs. However, positively charged amino acids alone are not enough to completely neutralize cancer cells but rather aid in rapid internalization of cancer cells for selective killing with low toxicity to normal cells. Moreover, hydrogen bonded interactions to encapsulate anti-cancer drugs has been utilized for increasing the oral efficiency and drug-controlled release. Therefore, only hydrogen bond on ACPs are not adequate to resist cancer cells. Thus, the combination between crucial amino acids and physicochemical properties could provide a synergistic effect that enhances the efficiency of ACPs [81]. As such, results from the analysis of the top ranked informative AAC are discussed below.
To roughly analyze the characteristics of ACPs, a few studies have shed light on important amino acids and dipeptides commonly found in ACPs by using simple composition analysis approaches without the use of experimental methods. For instance, Tyagi et al. [23] reported results of residue preference at 10 N-terminus and 10 C-terminus by using the sequence logos. Their analysis revealed that Leu and Lys were typically found at the N-terminus while Cys, Leu and Lys were typically found at the C-terminus. Chen et al. [26] also confirmed that amino acids including Lys, Ile, Cys, Glu, and Gly were abundant in ACPs as compared to non-ACPs. Recently, Manavalan et al. [29] revealed that the three top-ranked informative dipeptides consisted of KK, AK and KL. As can be noticed in Table 5 and Figure 4, our analysis results derived from the MDGI values of AAC and DPC were quite consistent with the three aforementioned studies.
As mentioned above, amongst the five top-ranked informative amino acids, Lys had the highest MDGI value of 29.54. Such feature with the highest MDGI is considered to be the most important factor governing the anticancer activity. As a basic residue, Lys is highly conserved in the composition of therapeutic peptides as it enhances the formation of electrostatic interactions between the peptide and the plasma membrane owing to its cationicity [82]. The role of cationicity in ACPs have been previously investigated by the study of Gopal et al. [83]. The authors modified a synthetic AMP known as HPA3NT3 (FKRLKKLFKKIWNWK), which was derived from Helicobacter pylori by substituting Arg and Asn with Lys at positions 3 and 13, respectively. Gopal et al. reported that the Lys-modified peptide increased the peptide selectivity toward negatively-charged membrane surface, which could allow for endocytosis and specificity of the APC in killing cancer cells. In addition, Wang et al. [84] has shown that the cell-penetrating ability is critical for specifically targeting cationic cancer cells. The authors observed that a Lys-rich peptide, L-K6 (IKKILSKIKKLLK-NH2) could internalize into MCF-7 cells without significant cell surface disruption. Interestingly, the authors also found that the internalized L-K6 could induce MCF-7 cell death without significant cytoskeleton disruption and mitochondrial impairment. Moreover, these findings indicated that the negatively-charged phosphatidylserine (PS) reported to be abundantly exposed on cancer cell surfaces [85,86] might contribute to the preferential binding of L-K6 to MCF-7 cancer cells. Furthermore, defensins are a large subfamily of natural cationic human AMP showing anticancer properties. Their roles in the immunomodulation of the innate and adaptive immune system have extensively been studied in previous works [87,88]. Defensins are categorized as either α-defensins or β-defensins, depending on their sequence homology and Cys connectivity (Figure 1), they are found to be rich in cationic residues (e.g., Lys and Arg) that enhances their immunomodulatory effects [89]. Papo et al. [90] constructed a short 15 amino acid peptide consisting of d- or l- amino acids (d-K6L9; LKLLKKLLKKLLKLL-NH2) and observed that d-K6L9 specifically targets and lyses cancer cells as well as inhibiting the growth of primary and metastatic tumors in in vivo models. The ability of d-K6L9 in selectively targeting cancer cells can be attributed to its cationicity whereby the peptide can bind to highly enriched surface-exposed PS in cancer cells as compared to non-cancer cells [91]. Similarly, Arg, which is a positively-charged residue that differs from Lys by its side chain, was shown to be the fourth informative amino acid form our analysis thus, reiterating the importance of the cationic property of APCs.
The second most important amino acid was Cys that afforded an MDGI value of 19.71. Cys is a highly reactive amino acid in which it contains a thiol group that is oxidized to form a dimer thus creating a disulfide bridge between two spatially adjacent Cys residues. The significance of disulfide bridges are seen in their strong, hydrophobic nature which are extremely important for the overall structural fold of the peptide, while also increasing the peptide stability against proteolytic degradation [92]. Almost a decade ago, GO-201 ([R]9CQCRRKNYGQLDIFP) achieved in facilitating the inhibition of mucin-1 cytoplasmic domain (MUC1-CD) [93]. Thus far, many studies have established the function of the GO-201 peptide whereby its anticancer activity was demonstrated both in vitro and in vivo against human breast cancer, prostate cancer, chronic myelogenous leukemia and pancreatic cancer [93,94,95,96]. From these aforementioned studies, the researchers were able to determine key motifs responsible for the anticancer activity of GO-201 as CxC, thus reiterating the crucial role of Cys disulfide bridges [97]. Furthermore, Tyuryaeva et al. [98] explored the molecular mechanism of Cys-induced cytotoxicity and the resulting apoptosis of cancer cells by designing and testing peptide sequences containing the key CxC ([R]9KCGCFF; named DIL) or CxxC ([R]9FFCPHCYQ; named DOL) motifs. The results obtained from this study suggests that the disulfide oxidoreductase is responsible for the effective cytotoxtic action of Cys. Therefore, it can be inferred that peptides containing the CxC or CxxC motifs form disulfide bridges that pilot a partial loss of protein function and the rapid onset of apoptosis. It can therefore be anticipated that such peptide sequences could help to create a potentially valuable new class of ACPs. Similarly, Schroeder et al. [99] ascertained a novel peptide GKAKCCK having pronounced antimicrobial activity. The authors discovered that the efficiency of this peptide was critically dependent on the two continuous cysteines and that upon mutating these residues to Ala or Ser, its activity was completely abrogated. Hence, given the well-established connection between antimicrobial and anticancer peptides [13,100], it is likely that the GKAKCCK peptide is able to induce antitumorigenic effects.

3.4. Mechanistic Interpretation of Informative PCP

The physicochemical properties of amino acids play an essential role as effective features for identifying and characterizing protein or peptide functions from their primary sequences [33,34,36,51,79]. It is well known that PCPs [39], such as molecular volume, exposure or accessible surface, polarity (hydrophobicity/hydrophilicity), charge/pK, hydrogen-bonding potential and so forth are correlated with the structure and function of the amino acid sequence [101]. Herein, we have obtained the ten top-ranked informative PCPs corresponding to their MDGI values as shown in Table 6.
The most important PCP is seen to be ‘hydrophobicity’ with an MDGI value of 1.51, which corresponds to the AAindex ARGP820101. In addition, from the analysis, 4 of the 10 top-ranked informative PCPs are related to α-helix properties of peptides consisting of ARGP820102, BEGF750101, ISOY800106 and CHOP780201 with corresponding MDGI values of 1.40, 0.92, 0.64 and 0.61, respectively. Furthermore, 3 of 10 top-ranked informative PCPs including the AAindices of BEGF750102, BEGF750103 and CHOP780101 are related to β-sheet properties with MDGI values of 1.00, 0.87 and 0.83, respectively. Consequently, it could be inferred that the hydrophobicity of a peptide along with its structural orientation (i.e., α-helix or β-sheet) are significant in determining the anticancer properties of peptides. The important PCPs of ACPs are analyzed and discussed below.

3.4.1. Hydrophobic Residues on α-Helical Structure Enhances Cell Penetrating Properties of Anticancer Peptides

It is commonly known that the hydrophobicity property is critical for the 3D structure of a protein. Many descriptors have been developed for measuring the hydrophobicity of proteins, depending on the presence of hydrophobic residues at specific locations in the peptide sequence. These indices are usually based on the 3D crystal structure of proteins coupled with the physicochemical properties of their side chains [102,103]. For instance, Argos et al. [104] utilized the hydrophobic index as measured by Nozaki and Tanford [105] for the 20 canonical amino acids for determining the role of hydrophobicity on the structure of amino acids. As mentioned above, the property of ARGP820101 with an MDGI value of 1.51 pertained to the hydrophobicity property. This result suggests that the presence of hydrophobic amino acid residues on the α-helix might be important for governing the anticancer activity of peptides. Moreover, the property of ARGP820103 (the fourth most important PCPs) is also associated with hydrophobicity property. Many previous studies have described the need for symbiosis between hydrophobicity and cationicity of peptides in order to enhance its cell penetrating property and its ability to reach target cell as to exert its effects [106,107,108,109]. In addition, researchers have made efforts to assess the influence of hydrophobicity on anticancer effects. For instance, Huang et al. [108] systematically altered the hydrophobicity of a 26-residue α-helical peptide (peptide P; KWKSFLKTFKSAKKTVLHTALKAISS) by replacing Ala residues with the more hydrophobic Leu residues as to increase its hydrophobicity or by changing Leu residues to Ala residues as to decrease its hydrophobicity. On the basis of these results, the authors observed a correlation between higher hydrophobic variants whereby the Leu-substituted peptides (A12L and A20L) exhibited greater anticancer activity against cancer cells with higher IC50 values than peptide P. Thus, peptides with higher hydrophobicity are assumed to penetrate deeper into the hydrophobic core of the cell membrane, thereby causing pores or channels on the cancer cell membrane and exhibiting greater anticancer activity. From these results, it can be inferred that the presence of hydrophobic residues on the α-helical structure of ACPs can enhance its anticancer properties.

3.4.2. Peptides Forming an Amphipathic α-Helix Contributes to Anticancer Activity

Interestingly, four of the ten top-ranked informative PCPs consisting of ARGP820102, BEGF750101, ISOY800106 and CHOP780201 described helical properties of peptides such as helical potential [104], inner helix [110], frequency of helix end [111] and frequency of α-helix [112], respectively, as shown in Table 6. This supports the notion that the helical structure plays a principal role in governing the anticancer property of peptides. In addition, the α-helical conformation as adopted by many ACPs in biological membranes, is now regarded as a key determinant of anticancer activity [113]. BMAP-27 (GRFKRFRKKFKKLFKKLSPVIPLLHL) and BMAP-28 (GGLRSLGRKILRAWKKY GPIIVPIIRI) are well-known peptides derived from bovine sources having a cationic NH2 terminal that forms an amphipathic α-helix. They exert their cytotoxic activity against human leukemia cells by enhancing membrane permeability, which allows for an influx of Ca++ ions, thereby leading to apoptosis via DNA fragmentation [114]. In addition, the structures of Cecropin A (KWKLFKKIEKVGQNIRDGIIKAGPAVAVVGQATQIAK) and Cecropin B (KWKVFKKIEKMGRNIRNGIVKAGPAIAVLGEAKAL), which were first derived from insects, primarily consisted of two α-helix [115,116] while its amphipathic N-terminal (i.e., that is capable of interacting with anionic membrane components) is responsible for mediating the cytotoxic activity against cancer cells. On the other hand, the C-terminal is hydrophobic and it is postulated to facilitate peptide entry into the membrane thereby enabling oligomerization to occur, which subsequently leads to leakage of cell components due to pore formation and eventual cell death [117]. In addition, Srisailam et al. [118] designed a custom lytic peptide by modifying the cecropin B (CB) to contain two identical hydrophobic segments on both the N and the C terminals, thus disrupting its normal α-helical structure. The resulting peptide CB-3, was unable to effectively lyse cell membranes on cancer cells as well as bacterial cells as compared to the control peptide, therefore highlighting the significance of the peptide structure on its function. Interestingly, Moore et al. [119] uncovered that CB exhibited anticancer activities when tested against colon adenocarcinoma cells in vivo and cytotoxicity against multidrug-resistant breast and ovarian cancer cell lines in vitro.
Similarly, Magainins isolated from the skin of Xenopus laevis is comprised of a α-helical secondary structure with separate cationic and hydrophobic faces. Magainin-2 (GIGKFLHSAKKFGKAFVGEIMNS) causes lysis of both hematopoietic and solid tumor cell lines as while observing selective cytotoxic activity against several human bladder cell lines [120]. In addition, it was perceived using florescence spectrometry that magainins could lyse tumor cells by forming ion-conducting α-helical channels in the cancer cell membrane [121]. Melittin (GIGAVLKVLTTGLPALISWIKRKRQQ) is another well-known ACP that is derived from the venom of European honeybee (Apis mellifera) [122]. The N-terminal region of melittin is largely hydrophobic whereas the C-terminal contains positively-charged amino acid residues. Melittin form channels in lipid bilayers and is lytic against both cancer cells and normal cells [123]. Self-association of amphipathic α-helical monomers of melittin is largely suggested to be responsible for perturbations of membrane integrity thereby leading to cellular lysis via the barrel stave mechanism [124]. Furthermore, the only human cathelicidin-derived cationic ACP to date, LL-37 (LLGDFFRKSKEKIGKEFKRIVQRIKDFLRNLVPRTES) is present throughout the body and assumes a α-helical conformation that is comprised of both cationic and hydrophobic faces oriented in a parallel manner with lipid membranes. As such these peptides are suggested to engage in a carpet-like mechanism for its cytotoxicity [125]. However, unlike other ACPs, LL-37 is toxic to eukaryotic cells at slightly higher concentrations and thus lacks the selectivity of an ACP. Nevertheless, Okumura et al. [126] showed that a C-terminal fragment of LL-37 (hCAP18) was able to induce mitochondrial depolarization and apoptosis in human oral squamous carcinoma cells but not in healthy human gingival fibroblasts and human keratinocyte cells. In summary, it is evident that the amphipathicity associated with certain α-helical structures is pivotal for the elicitation of peptide anticancer activities.

3.4.3. Peptides Forming β-Sheet Are Vital to Peptide Structure and Contribute to Anticancer Activity

As seen in Table 6, 3 out of the 20 top-ranked informative PCPs pertain to β-sheet structure (i.e., BEGF750102, BEGF750103 and CHOP78010) which correspond to properties such as conformational parameter of β-structure, conformational parameter of β-turn and normalized frequency of β-turn [110], respectively. Although, most of the well-known APCs are α-helical structures, some β-sheet ACPs have also been extensively studied. Many studies showed that defensins contain a group of closely related Cys-Arg rich cationic peptides. Out of these, human α- and β-defensins remain the most studied, comprising of six conserved Cys residues that form three intramolecular disulfide bridges between the N-terminal and the C-terminal regions of the peptide [127,128]. In the studies of [129,130], human neutrophil peptides (HNPs) -1 (ACYCRIPACIAGERRYGTCIYQGRLWAFCC), -2 (CYCRIPACIAGERRYGTCIYQGRLWAFCC) and -3 (DCYCRIPACIAGERRYGTCIY QGRLWAFCC) were shown to possess cytotoxic effects against several human and mouse tumor cells such as human B-lymphoma cells, human oral squamous carcinoma cells, and mouse teratocarcinoma cells. HNPs were also shown to kill cancer cells through membrane binding mediated by electrostatic interactions, followed by rapid collapse of the membrane potential and loss of membrane integrity [130,131]. In addition, Mader et al. [100,132] reported that bovine lactoferricin (LfcinB; FKCRRWQWRMKKLGAPSITCVRRAF) isolated from cow’s milk, represents a cationic ACP with amphipathic β-sheet configuration displaying anticancer activity against leukemia cells and various other carcinomas [100,132]. Moreover, LfcinB not only exerts its action by inducing apoptosis through direct disruption of the mitochondrial membrane, but is also capable of lysing the membrane itself, depending on the cancer cell type [133,134].
Furthermore, isolated from the hemocytes of horseshoe crab (Tachypleus tridentatus), the ACP tachyplesin I (KWCFRVCYRGICYRRCR) is arranged in two anti-parallel β-sheets that are held in place by two disulfide bonds [135]. Chen et al. [136] discovered that tachyplesin I killed cancer cells in a unique way by binding to hyaluronan on human prostate carcinoma cells, as well as to the C1q complement component in human serum, leading to the activation of the classical immune complement pathway. As a result, complement-mediated lysis of tachyplesin I-coated cancer cells was achieved. In a separate study, Adamia et al. [137] observed that the binding activity of tachyplesin I to C1q was dependent on its secondary structure which consists of two disulfide bonds formed by Cys residues at positions 3–16 and 7–12. Previously it has also been demonstrated that, the deletion of Cys residues from tachyplesin I altered its β-sheet structure to that of a linear form thus, disrupting its activity [138]. Therefore, from the aforementioned studies, it can be postulated that β-sheet is vital for maintaining the stability of peptide structures. As noticed in Table 5, our result of informative amino acids is consistent with the above-mentioned studies, where Cys is the second most important amino acid having a MDGI value of 19.71.
Nevertheless, further research is necessary to determine the key amino acids that could transform non-ACPs into ACPs and thus, allow the generation of a library of potentially novel ACPs. The design of effective ACPs from non-ACPs could possibly be based on the improvement of physicochemical properties, like hydrophobicity, positive charge, replacing the amino acid used or adding a functional motif with better activity towards the targeted cancer cells. In that regard, Tada et al. [139] revealed that the replacement of amino acid His to Arg in epidermal growth factor receptor (EGFR) binding peptide increased its anti-tumor activity via stronger binding affinity to the EGFR on cancer cells. Furthermore, the utilization of a motif containing the amino acids Trp, Met and Trp (WMW) was able to sensitize apoptosis and reduce the migratory effect on cancer cells which in turn could synergize the anti-cancer activity of peptides [140]. Therefore, it is of paramount importance to explore this area further.

3.5. Interpretable Rules Acquisition Obtained from the RF Model

Interpretable rules crucial for discriminating ACPs from non-ACPs were derived from the RF model built using AAC features on the benchmark dataset known as IR-ACP. The main advantages of these constructed rules are two folds: (i) to demonstrate which amino acid or which combination of amino acids are effective for ACP prediction, and (ii) to simply discriminate ACPs from non-ACPs without the need to go through the mathematical and computational details. Table 7 and Table A2 shows the twenty interpretable rules that were important for ACPs and eleven that were important for non-ACPs.
Table 7 presents eight interpretable rules that were important for identifying ACPs and each rule comprises a different number of criteria. If a query peptide meets all of the criteria in at least one of the eight rules, then it will be predicted as an ACP. As observed in Table 7, four rules afforded a prediction accuracy larger than 95%, namely rules #1, #2, #3 and #8, while the remaining rules provided acceptable results with accuracy ranging between 80.79% and 87.88%. These results indicated that the constructed rules are robust, both in terms of their accuracy and interpretability for predicting and characterizing ACPs.

3.6. ACPred Web Server

In order to make the prediction model presented herein a practical tool that can be widely used by the scientific community, we have constructed a web server called the ACPred using the best model as described in previous sections. The web interface has been established using the Shiny package under the R programming environment. The web server is freely accessible at http://codes.bio/acpred/. Screenshots of the ACPred web server is shown in Figure 5 in which panel A shows the web server prior to submission of input data and panel B shown the web server after the prediction has been made.
Briefly, a step-by-step guide on using the web server is given below:
  • Step 1. Enter the following URL into the web browser, http://codes.bio/acpred/.
  • Step 2. Either enter the query peptide sequence into the Input box or upload the sequence file by clicking on the “Choose file” button (i.e., found below the “Enter your input sequence(s) in FASTA format heading”).
  • Step 3. Click on the “Submit” button to initiate the prediction process.
  • Step 4. Prediction results are automatically displayed in the grey box found below the “Status/Output” heading. Typically, it takes a few seconds for the server to process the task. Users can also download the prediction results as a CSV file by pressing on the “Download CSV button”.
Additionally, users can also run a local copy of ACPred on their own computer by using the following one-line code in an R environment:
shiny::runGitHub(‘acpred-webserver’, ‘chaninlab’)
However, prior to running the aforementioned code, it is recommended that users first install the prerequisite R packages. This can be performed by using the following code:
install.packages(c(‘shiny’, ‘shinyjs’, ‘shinythemes’, ‘protr’, ‘seqinr’, ‘caret’, ‘markdown’))

4. Conclusions

There is increasing evidence indicating that ACPs could be a potential, valuable therapeutic agent for cancer treatment. With the avalanche of peptide sequences generated from many sources, there is significant demand for computational methods that can support the timely discrimination of ACPs from other peptides, as well as offer insights from analysis using only their primary sequences. In this study, we have developed a the ACPred model in a systematic manner by taking advantage of two powerful machine learning methods consisting of support vector machine (SVM) and random forest (RF) for the prediction and characterization of anticancer activities of peptides. The best model for further development of ACPred was achieved by using the SVM method coupled with the combination of amino acid composition and pseudo amino acid composition features. By comparing with the state-of-the-art method, the TargetACP, it was found that the performance of ACPred was comparable to that of TargetACP (95.61% vs. 96.22% accuracies) as assessed by the rigorous jackknife test. Moreover, the feature importance analysis and interpretable rules extraction were carried out by using the RF model to identify which feature was the most important in predicting ACPs while also providing a better understanding of the anticancer activities of investigated peptides. Results from model interpretation revealed that the hydrophobic residues on the α-helical structure, Cys residues on the β-sheet structure and the formation of amphipathic α-helix played a crucial role in the anticancer activity of the investigated peptides. Finally, a web server named ACPred was established and made freely available online at https://codes.bio/acpred/.

Author Contributions

W.S. conceived, designed, performed and analyzed the experiments. N.S. and C.N. analyzed the data. W.S., N.S. and V.P. drafted the manuscript. W.S. and C.N. contributed the code for constructing the web server. C.N. vetted the manuscript. All authors read and approved the manuscript.

Funding

This work is supported by the Office of Higher Education Commission and the Thailand Research Fund (No. MRG6180226); the New Researcher Grant (A31/2561) from Mahidol University; and the Center of Excellence on Medical Biotechnology (CEMB), the S&T Postgraduate Education and Research Development Office (PERDO) and the Office of Higher Education Commission (OHEC), Thailand.

Acknowledgments

We thank the reviewers for their great comments.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Summary of MODI index as derived from various types of peptide features on the benchmark dataset.
Table A1. Summary of MODI index as derived from various types of peptide features on the benchmark dataset.
FeatureMODI
AAC0.868
DPC0.836
PCP0.722
PseAAC0.859
Am-PseAAC0.859
AAC + PseAAC0.859
AAC + Am-PseAAC0.859
PseAAC + Am-PseAAC0.856
AAC + PseAAC + Am-PseAAC0.856
Table A2. Twenty if-then rules for discriminating ACP from non-ACP using random forest and amino acid composition.
Table A2. Twenty if-then rules for discriminating ACP from non-ACP using random forest and amino acid composition.
% Covered SamplesRulePrediction Result
14.33R > 0.017 and C ≤ 0.1145 and H ≤ 0.0645 and L > 0.1605 and K ≤ 0.1035non-ACP
13.74G > 0.041 and I > 0.0615 and L ≤ 0.1715 and K > 0.0385 and M ≤ 0.027ACP
9.65E > 0.058 and G ≤ 0.068non-ACP
7.89R ≤ 0.0515 and Q ≤ 0.026 and K > 0.094ACP
7.89C ≤ 0.1145 and H ≤ 0.0645 and K ≤ 0.1055 and W > 0.028non-ACP
5.26C > 0.1145 and P ≤ 0.073ACP
4.97N > 0.019 and K ≤ 0.1705 and M > 0.037non-ACP
4.09Q > 0.098 and K ≤ 0.0575non-ACP
3.22C > 0.015 and C ≤ 0.043non-ACP
3.80L ≤ 0.093 and F > 0.0715 and S ≤ 0.152ACP
6.43R > 0.013 and R ≤ 0.1345 and C ≤ 0.1205 and H ≤ 0.0465 and F ≤ 0.1345non-ACP
2.63A ≤ 0.0145 and Q ≤ 0.026ACP
3.22I ≤ 0.0335 and P > 0.074non-ACP
2.34R ≤ 0.0665 and E ≤ 0.044 and H > 0.052ACP
1.75N ≤ 0.0565 and H ≤ 0.0415 and K≤ 0.0335non-ACP
1.46I > 0.108 and K > 0.055ACP
1.17E ≤ 0.0545 and G > 0.0365 and K > 0.108 and M ≤ 0.04ACP
2.05I ≤ 0.0665 and F ≤ 0.006non-ACP
2.05S ≤ 0.04non-ACP
2.05ElseACP

References

  1. WHO Cancer. Available online: https://www.who.int/news-room/fact-sheets/detail/cancer (accessed on 30 January 2019).
  2. Bray, F.; Ferlay, J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; 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, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Ferlay, J.; Colombet, M.; Soerjomataram, I.; Mathers, C.; Parkin, D.M.; Pineros, M.; Znaor, A.; Bray, F. Estimating the global cancer incidence and mortality in 2018: GLOBOCAN sources and methods. Int. J. Cancer 2018, 144, 1941–1953. [Google Scholar] [CrossRef] [Green Version]
  4. Craik, D.J.; Fairlie, D.P.; Liras, S.; Price, D. The future of peptide-based drugs. Chem. Biol. Drug Des. 2013, 81, 136–147. [Google Scholar] [CrossRef] [PubMed]
  5. Vlieghe, P.; Lisowski, V.; Martinez, J.; Khrestchatisky, M. Synthetic therapeutic peptides: Science and market. Drug Discov. Today 2010, 15, 40–56. [Google Scholar] [CrossRef]
  6. Lau, J.L.; Dunn, M.K. Therapeutic peptides: Historical perspectives, current development trends, and future directions. Bioorg. Med. Chem. 2018, 26, 2700–2707. [Google Scholar] [CrossRef] [PubMed]
  7. Fosgerau, K.; Hoffmann, T. Peptide therapeutics: Current status and future directions. Drug Discov. Today 2015, 20, 122–128. [Google Scholar] [CrossRef]
  8. O’Brien-Simpson, N.M.; Hoffmann, R.; Chia, C.S.B.; Wade, J.D. Editorial: Antimicrobial and Anticancer Peptides. Front. Chem. 2018, 6, 13. [Google Scholar] [CrossRef] [Green Version]
  9. Jin, G.; Weinberg, A. Human antimicrobial peptides and cancer. Semin. Cell Dev. Biol. 2018, 88, 156–162. [Google Scholar] [CrossRef] [PubMed]
  10. Karpinski, T.M.; Adamczak, A. Anticancer Activity of Bacterial Proteins and Peptides. Pharmaceutics 2018, 10, 54. [Google Scholar] [CrossRef]
  11. Usmani, S.S.; Bedi, G.; Samuel, J.S.; Singh, S.; Kalra, S.; Kumar, P.; Ahuja, A.A.; Sharma, M.; Gautam, A.; Raghava, G.P.S. THPdb: Database of FDA-approved peptide and protein therapeutics. PLoS ONE 2017, 12, e0181748. [Google Scholar] [CrossRef]
  12. Deslouches, B.; Di, Y.P. Antimicrobial peptides with selective antitumor mechanisms: Prospect for anticancer applications. Oncotarget 2017, 8, 46635–46651. [Google Scholar] [CrossRef]
  13. Gaspar, D.; Veiga, A.S.; Castanho, M.A. From antimicrobial to anticancer peptides. A review. Front. Microbiol. 2013, 4, 294. [Google Scholar] [CrossRef] [Green Version]
  14. Ran, S.; Downes, A.; Thorpe, P.E. Increased exposure of anionic phospholipids on the surface of tumor blood vessels. Cancer Res. 2002, 62, 6132–6140. [Google Scholar]
  15. Yoon, W.H.; Park, H.D.; Lim, K.; Hwang, B.D. Effect of O-glycosylated mucin on invasion and metastasis of HM7 human colon cancer cells. Biochem. Biophys. Res. Commun. 1996, 222, 694–699. [Google Scholar] [CrossRef]
  16. Sok, M.; Sentjurc, M.; Schara, M. Membrane fluidity characteristics of human lung cancer. Cancer Lett. 1999, 139, 215–220. [Google Scholar] [CrossRef]
  17. Chaudhary, J.; Munshi, M. Scanning electron microscopic analysis of breast aspirates. Cytopathology 1995, 6, 162–167. [Google Scholar] [CrossRef]
  18. Dobrzynska, I.; Szachowicz-Petelska, B.; Sulkowski, S.; Figaszewski, Z. Changes in electric charge and phospholipids composition in human colorectal cancer cells. Mol. Cell. Biochem. 2005, 276, 113–119. [Google Scholar] [CrossRef]
  19. Domalaon, R.; Findlay, B.; Ogunsina, M.; Arthur, G.; Schweizer, F. Ultrashort cationic lipopeptides and lipopeptoids: Evaluation and mechanistic insights against epithelial cancer cells. Peptides 2016, 84, 58–67. [Google Scholar] [CrossRef]
  20. Gaspar, D.; Freire, J.M.; Pacheco, T.R.; Barata, J.T.; Castanho, M.A. Apoptotic human neutrophil peptide-1 anti-tumor activity revealed by cellular biomechanics. Biochim. Biophys. Acta 2015, 1853, 308–316. [Google Scholar] [CrossRef] [Green Version]
  21. Figueiredo, C.R.; Matsuo, A.L.; Massaoka, M.H.; Polonelli, L.; Travassos, L.R. Anti-tumor activities of peptides corresponding to conserved complementary determining regions from different immunoglobulins. Peptides 2014, 59, 14–19. [Google Scholar] [CrossRef]
  22. 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, 766–781. [Google Scholar] [CrossRef] [PubMed]
  23. Tyagi, A.; Kapoor, P.; Kumar, R.; Chaudhary, K.; Gautam, A.; Raghava, G. In silico models for designing and discovering novel anticancer peptides. Sci. Rep. 2013, 3, 2984. [Google Scholar] [CrossRef]
  24. Hajisharifi, Z.; Piryaiee, M.; Beigi, M.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. [Google Scholar] [CrossRef] [PubMed]
  25. Vijayakumar, S.; Lakshmi, P. ACPP: A web server for prediction and design of anti-cancer peptides. Int. J. Pept. Res. Ther. 2015, 21, 99–106. [Google Scholar] [CrossRef]
  26. Chen, W.; Ding, H.; Feng, P.; Lin, H.; Chou, K.-C. iACP: A sequence-based tool for identifying anticancer peptides. Oncotarget 2016, 7, 16895. [Google Scholar] [CrossRef] [PubMed]
  27. Li, F.-M.; Wang, X.-Q. Identifying anticancer peptides by using improved hybrid compositions. Sci. Rep. 2016, 6, 33910. [Google Scholar] [CrossRef]
  28. Akbar, S.; Hayat, M.; Iqbal, M.; Jan, M.A. iACP-GAEnsC: Evolutionary genetic algorithm based ensemble classification of anticancer peptides by utilizing hybrid feature space. Artif. Intell. Med. 2017, 79, 62–70. [Google Scholar] [CrossRef]
  29. Manavalan, B.; Basith, S.; Shin, T.H.; Choi, S.; Kim, M.O.; Lee, G. MLACP: Machine-learning-based prediction of anticancer peptides. Oncotarget 2017, 8, 77121. [Google Scholar] [CrossRef]
  30. Xu, L.; Liang, G.; Wang, L.; Liao, C. A Novel Hybrid Sequence-Based Model for Identifying Anticancer Peptides. Genes 2018, 9, 158. [Google Scholar] [CrossRef]
  31. Kabir, M.; Arif, M.; Ahmad, S.; Ali, Z.; Swati, Z.N.K.; Yu, D.-J. Intelligent computational method for discrimination of anticancer peptides by incorporating sequential and evolutionary profiles information. Chemom. Intell. Lab. Syst. 2018, 182, 158–165. [Google Scholar] [CrossRef]
  32. Shoombuatong, W.; Schaduangrat, N.; Nantasenamat, C. Unraveling the bioactivity of anticancer peptides as deduced from machine learning. EXCLI J. 2018, 17, 734. [Google Scholar]
  33. Pratiwi, R.; Malik, A.A.; Schaduangrat, N.; Prachayasittikul, V.; Wikberg, J.E.; Nantasenamat, C.; Shoombuatong, W. CryoProtect: A Web Server for Classifying Antifreeze Proteins from Nonantifreeze Proteins. J. Chem. 2017, 2017. [Google Scholar] [CrossRef]
  34. Win, T.S.; Malik, A.A.; Prachayasittikul, V.; S Wikberg, J.E.; Nantasenamat, C.; Shoombuatong, W. HemoPred: A web server for predicting the hemolytic activity of peptides. Future Med. Chem. 2017, 9, 275–291. [Google Scholar] [CrossRef]
  35. Shoombuatong, W.; Schaduangrat, N.; Nantasenamat, C. Towards understanding aromatase inhibitory activity via QSAR modeling. EXCLI J. 2018, 17, 688. [Google Scholar]
  36. Win, T.S.; Schaduangrat, N.; Prachayasittikul, V.; Nantasenamat, C.; Shoombuatong, W. PAAP: A web server for predicting antihypertensive activity of peptides. Future Med. Chem. 2018, 10, 1749–1767. [Google Scholar] [CrossRef]
  37. Shoombuatong, W.; Mekha, P.; Chaijaruwanich, J. Sequence based human leukocyte antigen gene prediction using informative physicochemical properties. Int. J. Data Min. Bioinform. 2015, 13, 211–224. [Google Scholar] [CrossRef]
  38. Simeon, S.; Shoombuatong, W.; Anuwongcharoen, N.; Preeyanon, L.; Prachayasittikul, V.; Wikberg, J.E.; Nantasenamat, C. osFP: A web server for predicting the oligomeric states of fluorescent proteins. J. Cheminform. 2016, 8, 72. [Google Scholar] [CrossRef]
  39. Kawashima, S.; Kanehisa, M. AAindex: Amino acid index database. Nucleic Acids Res. 2000, 28, 374. [Google Scholar] [CrossRef]
  40. Chou, K.-C. Using amphiphilic pseudo amino acid composition to predict enzyme subfamily classes. Bioinformatics 2004, 21, 10–19. [Google Scholar] [CrossRef]
  41. Chou, K.-C. Some remarks on protein attribute prediction and pseudo amino acid composition. J. Theor. Biol. 2011, 273, 236–247. [Google Scholar] [CrossRef]
  42. Shen, H.-B.; Chou, K.-C. PseAAC: A flexible web server for generating various kinds of protein pseudo amino acid composition. Anal. Biochem. 2008, 373, 386–388. [Google Scholar] [CrossRef]
  43. Xiao, N.; Cao, D.-S.; Zhu, M.-F.; Xu, Q.-S. protr/ProtrWeb: R package and web server for generating various numerical representation schemes of protein sequences. Bioinformatics 2015, 31, 1857–1859. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Golbraikh, A.; Muratov, E.; Fourches, D.; Tropsha, A. Data set modelability by QSAR. J. Chem. Inf. Model. 2014, 54, 1–4. [Google Scholar] [CrossRef]
  45. Shoombuatong, W.; Prathipati, P.; Prachayasittikul, V.; Schaduangrat, N.; Malik, A.A.; Pratiwi, R.; Wanwimolruk, S.; Wikberg, J.E.; Gleeson, M.P.; Spjuth, O. Towards Predicting the Cytochrome P450 Modulation: From QSAR to Proteochemometric Modeling. Curr. Drug Metab. 2017, 18, 540–555. [Google Scholar] [CrossRef]
  46. Shoombuatong, W.; Hongjaisee, S.; Barin, F.; Chaijaruwanich, J.; Samleerat, T. HIV-1 CRF01_AE coreceptor usage prediction using kernel methods based logistic model trees. Comput. Biol. Med. 2012, 42, 885–889. [Google Scholar] [CrossRef]
  47. Simeon, S.; Li, H.; Win, T.S.; Malik, A.A.; Kandhro, A.H.; Piacham, T.; Shoombuatong, W.; Nuchnoi, P.; Wikberg, J.E.; Gleeson, M.P. PepBio: Predicting the bioactivity of host defense peptides. RSC Adv. 2017, 7, 35119–35134. [Google Scholar] [CrossRef]
  48. Shoombuatong, W.; Prachayasittikul, V.; Prachayasittikul, V.; Nantasenamat, C. Prediction of aromatase inhibitory activity using the efficient linear method (ELM). EXCLI J. 2015, 14, 452. [Google Scholar]
  49. Shoombuatong, W.; Mekha, P.; Waiyamai, K.; Cheevadhanarak, S.; Chaijaruwanicha, J. Prediction of human leukocyte antigen gene using k-nearest neighbour classifier based on spectrum kernel. ScienceAsia 2013, 39, 42–49. [Google Scholar] [CrossRef]
  50. Shoombuatong, W.; Traisathit, P.; Prasitwattanaseree, S.; Tayapiwatana, C.; Cutler, R.; Chaijaruwanich, J. Prediction of the disulphide bonding state of cysteines in proteins using Conditional Random Fields. Int. J. Data Min. Bioinform. 2011, 5, 449–464. [Google Scholar] [CrossRef]
  51. Charoenkwan, P.; Shoombuatong, W.; Lee, H.-C.; Chaijaruwanich, J.; Huang, H.-L.; Ho, S.-Y. SCMCRYS: Predicting protein crystallization using an ensemble scoring card method with estimating propensity scores of P-collocated amino acid pairs. PLoS ONE 2013, 8, e72368. [Google Scholar] [CrossRef] [PubMed]
  52. Anuwongcharoen, N.; Shoombuatong, W.; Tantimongcolwat, T.; Prachayasittikul, V.; Nantasenamat, C. Exploring the chemical space of influenza neuraminidase inhibitors. PeerJ 2016, 4, e1958. [Google Scholar] [CrossRef]
  53. Shoombuatong, W.; Prachayasittikul, V.; Anuwongcharoen, N.; Songtawee, N.; Monnor, T.; Prachayasittikul, S.; Prachayasittikul, V.; Nantasenamat, C. Navigating the chemical space of dipeptidyl peptidase-4 inhibitors. Drug Des. Dev. Ther. 2015, 9, 4515. [Google Scholar]
  54. Suvannang, N.; Preeyanon, L.; Malik, A.A.; Schaduangrat, N.; Shoombuatong, W.; Worachartcheewan, A.; Tantimongcolwat, T.; Nantasenamat, C. Probing the origin of estrogen receptor alpha inhibition via large-scale QSAR study. RSC Adv. 2018, 8, 11344–11356. [Google Scholar] [CrossRef] [Green Version]
  55. Simeon, S.; Anuwongcharoen, N.; Shoombuatong, W.; Malik, A.A.; Prachayasittikul, V.; Wikberg, J.E.; Nantasenamat, C. Probing the origins of human acetylcholinesterase inhibition via QSAR modeling and molecular docking. PeerJ 2016, 4, e2322. [Google Scholar] [CrossRef]
  56. Breiman, L. Random forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  57. Breiman, L. Classification and Regression Trees; Routledge: London, UK, 2017. [Google Scholar]
  58. Breiman, L. randomForest: Breiman and Cutler’s Random Forests for Classification and Regression. R Package Version. 2006. Available online: http://stat-www. berkeley. edu/users/breiman/RandomForests (accessed on 30 January 2019).
  59. Kuhn, M. Building predictive models in R using the caret package. J. Stat. Softw. 2008, 28, 1–26. [Google Scholar] [CrossRef]
  60. Cortes, C.; Vapnik, V. Support-vector networks. Mach. Learn. 1995, 20, 273–297. [Google Scholar] [CrossRef]
  61. Drucker, H.; Burges, C.J.; Kaufman, L.; Smola, A.J.; Vapnik, V. Support Vector Regression Machines. In Advances in Neural Information Processing Systems; MIT Press: Cambridge, MA, USA, 1997; pp. 155–161. [Google Scholar]
  62. Joachims, T. Svmlight: Support Vector Machine. SVM-Light Support. Vector Machine. University of Dortmund, 1999, 19. Available online: http://svmlight. joachims. org/ (accessed on 30 January 2019).
  63. Lameski, P.; Zdravevski, E.; Mingov, R.; Kulakov, A. SVM parameter tuning with grid search and its impact on reduction of model over-fitting. In Rough Sets, Fuzzy Sets, Data Mining, and Granular Computing; Springer: Berlin, Germany, 2015; pp. 464–474. [Google Scholar]
  64. Tung, C.-W.; Ziehm, M.; Kämper, A.; Kohlbacher, O.; Ho, S.-Y. POPISK: T-cell reactivity prediction using support vector machines and string kernels. BMC Bioinform. 2011, 12, 446. [Google Scholar] [CrossRef]
  65. Calle, M.L.; Urrea, V. Letter to the editor: Stability of random forest importance measures. Brief. Bioinform. 2010, 12, 86–89. [Google Scholar] [CrossRef]
  66. Shoombuatong, W.; Huang, H.-L.; Chaijaruwanich, J.; Charoenkwan, P.; Lee, H.-C.; Ho, S.-Y. Predicting Protein Crystallization Using a Simple Scoring Card Method. In Proceedings of the IEEE Symposium on Computational Intelligence in Bioinformatics and Computational Biology (CIBCB), Singapore, 16–19 April 2013; pp. 23–30. [Google Scholar]
  67. Hasan, M.M.; Yang, S.; Zhou, Y.; Mollah, M.N.H. SuccinSite: A computational tool for the prediction of protein succinylation sites by exploiting the amino acid patterns and properties. Mol. Biosyst. 2016, 12, 786–795. [Google Scholar] [CrossRef]
  68. Worachartcheewan, A.; Shoombuatong, W.; Pidetcha, P.; Nopnithipat, W.; Prachayasittikul, V.; Nantasenamat, C. Predicting metabolic syndrome using the random forest method. Sci. World J. 2015, 2015, 58150. [Google Scholar] [CrossRef] [PubMed]
  69. Chou, K.-C.; Shen, H.-B. Cell-PLoc: A package of Web servers for predicting subcellular localization of proteins in various organisms. Nat. Protoc. 2008, 3, 153. [Google Scholar] [CrossRef] [PubMed]
  70. Chou, K.-C.; Shen, H.-B. Recent progress in protein subcellular location prediction. Anal. Biochem. 2007, 370, 1. [Google Scholar] [CrossRef]
  71. Robin, X.; Turck, N.; Hainard, A.; Tiberti, N.; Lisacek, F.; Sanchez, J.-C.; Müller, M. pROC: An open-source package for R and S+ to analyze and compare ROC curves. BMC Bioinform. 2011, 12, 77. [Google Scholar] [CrossRef]
  72. Liaw, C.; Tung, C.-W.; Ho, S.-Y. Prediction and analysis of antibody amyloidogenesis from sequences. PLoS ONE 2013, 8, e53235. [Google Scholar] [CrossRef]
  73. Chawla, N.V.; Bowyer, K.W.; Hall, L.O.; Kegelmeyer, W.P. SMOTE: Synthetic minority over-sampling technique. J. Artif. Intell. Res. 2002, 16, 321–357. [Google Scholar] [CrossRef]
  74. Kumar, R.; Chaudhary, K.; Chauhan, J.S.; Nagpal, G.; Kumar, R.; Sharma, M.; Raghava, G.P. An in silico platform for predicting, screening and designing of antihypertensive peptides. Sci. Rep. 2015, 5, 12512. [Google Scholar] [CrossRef] [PubMed]
  75. Shen, H.-B.; Chou, K.-C. Gneg-mPLoc: A top-down strategy to enhance the quality of predicting subcellular localization of Gram-negative bacterial proteins. J. Theor. Biol. 2010, 264, 326–333. [Google Scholar] [CrossRef]
  76. Chou, K.-C.; Shen, H.-B. Plant-mPLoc: A top-down strategy to augment the power for predicting plant protein subcellular localization. PLoS ONE 2010, 5, e11335. [Google Scholar] [CrossRef]
  77. Hellberg, S.; Sjoestroem, M.; Skagerberg, B.; Wold, S. Peptide quantitative structure-activity relationships, a multivariate approach. J. Med. Chem. 1987, 30, 1126–1135. [Google Scholar] [CrossRef]
  78. Sandberg, M.; Eriksson, L.; Jonsson, J.; Sjöström, M.; Wold, S. New chemical descriptors relevant for the design of biologically active peptides. A multivariate characterization of 87 amino acids. J. Med. Chem. 1998, 41, 2481–2491. [Google Scholar] [CrossRef]
  79. Huang, H.-L. Propensity scores for prediction and characterization of bioluminescent proteins from sequences. PLoS ONE 2014, 9, e97158. [Google Scholar] [CrossRef] [PubMed]
  80. Borrelli, A.; Tornesello, A.L.; Tornesello, M.L.; Buonaguro, F.M. Cell Penetrating Peptides as Molecular Carriers for Anti-Cancer Agents. Molecules 2018, 23, 295. [Google Scholar] [CrossRef]
  81. Benabdelouahab, Y.; Muñoz-Moreno, L.; Frik, M.; de la Cueva-Alique, I.; El Amrani, M.A.; Contel, M.; Bajo, A.M.; Cuenca, T.; Royo, E. Hydrogen Bonding and Anticancer Properties of Water-Soluble Chiral p-Cymene RuII Compounds with Amino-Oxime Ligands. Eur. J. Inorg. Chem. 2015, 2015, 2295–2307. [Google Scholar] [CrossRef]
  82. Carlberg, C. Current understanding of the function of the nuclear vitamin D receptor in response to its natural and synthetic ligands. Recent Results Cancer Res. 2003, 164, 29–42. [Google Scholar]
  83. Gopal, R.; Park, S.C.; Ha, K.J.; Cho, S.J.; Kim, S.W.; Song, P.I.; Nah, J.W.; Park, Y.; Hahm, K.S. Effect of Leucine and Lysine substitution on the antimicrobial activity and evaluation of the mechanism of the HPA3NT3 analog peptide. J. Pept. Sci. 2009, 15, 589–594. [Google Scholar] [CrossRef]
  84. Wang, C.; Dong, S.; Zhang, L.; Zhao, Y.; Huang, L.; Gong, X.; Wang, H.; Shang, D. Cell surface binding, uptaking and anticancer activity of L-K6, a lysine/leucine-rich peptide, on human breast cancer MCF-7 cells. Sci. Rep. 2017, 7, 8293. [Google Scholar] [CrossRef] [PubMed]
  85. Riedl, S.; Rinner, B.; Asslaber, M.; Schaider, H.; Walzer, S.; Novak, A.; Lohner, K.; Zweytick, D. In search of a novel target—Phosphatidylserine exposed by non-apoptotic tumor cells and metastases of malignancies with poor treatment efficacy. Biochim. Biophys. Acta 2011, 1808, 2638–2645. [Google Scholar] [CrossRef] [PubMed]
  86. Wang, C.; Chen, Y.W.; Zhang, L.; Gong, X.G.; Zhou, Y.; Shang, D.J. Melanoma cell surface-expressed phosphatidylserine as a therapeutic target for cationic anticancer peptide, temporin-1CEa. J. Drug Target. 2016, 24, 548–556. [Google Scholar] [CrossRef] [PubMed]
  87. Yang, D.; Biragyn, A.; Hoover, D.M.; Lubkowski, J.; Oppenheim, J.J. Multiple roles of antimicrobial defensins, cathelicidins, and eosinophil-derived neurotoxin in host defense. Annu. Rev. Immunol. 2004, 22, 181–215. [Google Scholar] [CrossRef] [PubMed]
  88. Yang, D.; Oppenheim, J.J. Antimicrobial proteins act as “alarmins” in joint immune defense. Arthritis Rheum. 2004, 50, 3401–3403. [Google Scholar] [CrossRef]
  89. Zou, G.; de Leeuw, E.; Li, C.; Pazgier, M.; Li, C.; Zeng, P.; Lu, W.Y.; Lubkowski, J.; Lu, W. Toward understanding the cationicity of defensins. Arg and Lys versus their noncoded analogs. J. Biol. Chem. 2007, 282, 19653–19665. [Google Scholar] [CrossRef]
  90. Papo, N.; Seger, D.; Makovitzki, A.; Kalchenko, V.; Eshhar, Z.; Degani, H.; Shai, Y. Inhibition of tumor growth and elimination of multiple metastases in human prostate and breast xenografts by systemic inoculation of a host defense-like lytic peptide. Cancer Res. 2006, 66, 5371–5378. [Google Scholar] [CrossRef]
  91. Zwaal, R.F.; Comfurius, P.; Bevers, E.M. Surface exposure of phosphatidylserine in pathological cells. Cell. Mol. Life Sci. 2005, 62, 971–988. [Google Scholar] [CrossRef]
  92. Tanabe, H.; Ayabe, T.; Maemoto, A.; Ishikawa, C.; Inaba, Y.; Sato, R.; Moriichi, K.; Okamoto, K.; Watari, J.; Kono, T.; et al. Denatured human alpha-defensin attenuates the bactericidal activity and the stability against enzymatic digestion. Biochem. Biophys. Res. Commun. 2007, 358, 349–355. [Google Scholar] [CrossRef]
  93. Raina, D.; Ahmad, R.; Joshi, M.D.; Yin, L.; Wu, Z.; Kawano, T.; Vasir, B.; Avigan, D.; Kharbanda, S.; Kufe, D. Direct targeting of the mucin 1 oncoprotein blocks survival and tumorigenicity of human breast carcinoma cells. Cancer Res. 2009, 69, 5133–5141. [Google Scholar] [CrossRef]
  94. Joshi, M.D.; Ahmad, R.; Yin, L.; Raina, D.; Rajabi, H.; Bubley, G.; Kharbanda, S.; Kufe, D. MUC1 oncoprotein is a druggable target in human prostate cancer cells. Mol. Cancer Ther. 2009, 8, 3056–3065. [Google Scholar] [CrossRef] [Green Version]
  95. Yin, L.; Ahmad, R.; Kosugi, M.; Kawano, T.; Avigan, D.; Stone, R.; Kharbanda, S.; Kufe, D. Terminal differentiation of chronic myelogenous leukemia cells is induced by targeting of the MUC1-C oncoprotein. Cancer Biol. Ther. 2010, 10, 483–491. [Google Scholar] [CrossRef] [Green Version]
  96. Banerjee, S.; Mujumdar, N.; Dudeja, V.; Mackenzie, T.; Krosch, T.K.; Sangwan, V.; Vickers, S.M.; Saluja, A.K. MUC1c regulates cell survival in pancreatic cancer by preventing lysosomal permeabilization. PLoS ONE 2012, 7, e43020. [Google Scholar] [CrossRef]
  97. Raina, D.; Ahmad, R.; Rajabi, H.; Panchamoorthy, G.; Kharbanda, S.; Kufe, D. Targeting cysteine-mediated dimerization of the MUC1-C oncoprotein in human cancer cells. Int. J. Oncol. 2012, 40, 1643–1649. [Google Scholar]
  98. Tyuryaeva, I.I.; Lyublinskaya, O.G.; Podkorytov, I.S.; Skrynnikov, N.R. Origin of anti-tumor activity of the cysteine-containing GO peptides and further optimization of their cytotoxic properties. Sci. Rep. 2017, 7, 40217. [Google Scholar] [CrossRef] [Green Version]
  99. Schroeder, B.O.; Wu, Z.; Nuding, S.; Groscurth, S.; Marcinowski, M.; Beisner, J.; Buchner, J.; Schaller, M.; Stange, E.F.; Wehkamp, J. Reduction of disulphide bonds unmasks potent antimicrobial activity of human beta-defensin 1. Nature 2011, 469, 419–423. [Google Scholar] [CrossRef]
  100. Hoskin, D.W.; Ramamoorthy, A. Studies on anticancer activities of antimicrobial peptides. Biochim. Biophys. Acta 2008, 1778, 357–375. [Google Scholar] [CrossRef] [Green Version]
  101. Tsai, C.S. Biomacromolecules: Introduction to Structure, Function and Informatics; John Wiley & Sons: Hoboken, NJ, USA, 2007. [Google Scholar]
  102. Roseman, M.A. Hydrophilicity of polar amino acid side-chains is markedly reduced by flanking peptide bonds. J. Mol. Biol. 1988, 200, 513–522. [Google Scholar] [CrossRef]
  103. Kyte, J.; Doolittle, R.F. A simple method for displaying the hydropathic character of a protein. J. Mol Biol 1982, 157, 105–132. [Google Scholar] [CrossRef] [Green Version]
  104. Argos, P.; Rao, J.K.; Hargrave, P.A. Structural prediction of membrane-bound proteins. Eur. J. Biochem. 1982, 128, 565–575. [Google Scholar] [CrossRef]
  105. Nozaki, Y.; Tanford, C. The solubility of amino acids and two glycine peptides in aqueous ethanol and dioxane solutions. Establishment of a hydrophobicity scale. J. Biol. Chem. 1971, 246, 2211–2217. [Google Scholar]
  106. Mader, J.S.; Hoskin, D.W. Cationic antimicrobial peptides as novel cytotoxic agents for cancer treatment. Expert Opin. Investig. Drugs 2006, 15, 933–946. [Google Scholar] [CrossRef]
  107. Chen, C.; Yang, C.; Chen, Y.; Wang, F.; Mu, Q.; Zhang, J.; Li, Z.; Pan, F.; Xu, H.; Lu, J.R. Surface Physical Activity and Hydrophobicity of Designed Helical Peptide Amphiphiles Control Their Bioactivity and Cell Selectivity. ACS Appl. Mater. Interfaces 2016, 8, 26501–26510. [Google Scholar] [CrossRef]
  108. Huang, Y.B.; Wang, X.F.; Wang, H.Y.; Liu, Y.; Chen, Y. Studies on mechanism of action of anticancer peptides by modulation of hydrophobicity within a defined structural framework. Mol. Cancer Ther. 2011, 10, 416–426. [Google Scholar] [CrossRef]
  109. Park, J.H.; Kim, H.A.; Cho, S.H.; Lee, M. Characterization of hydrophobic anti-cancer drug-loaded amphiphilic peptides as a gene carrier. J. Cell. Biochem. 2012, 113, 1645–1653. [Google Scholar] [CrossRef]
  110. Chou, P.Y.; Fasman, G.D. Conformational parameters for amino acids in helical, beta-sheet, and random coil regions calculated from proteins. Biochemistry 1974, 13, 211–222. [Google Scholar] [CrossRef]
  111. Isogai, Y.; Nemethy, G.; Rackovsky, S.; Leach, S.J.; Scheraga, H.A. Characterization of multiple bends in proteins. Biopolymers 1980, 19, 1183–1210. [Google Scholar] [CrossRef] [PubMed]
  112. Chou, P.Y.; Fasman, G.D. Empirical predictions of protein conformation. Annu. Rev. Biochem. 1978, 47, 251–276. [Google Scholar] [CrossRef] [PubMed]
  113. Vermeer, L.S.; Lan, Y.; Abbate, V.; Ruh, E.; Bui, T.T.; Wilkinson, L.J.; Kanno, T.; Jumagulova, E.; Kozlowska, J.; Patel, J.; et al. Conformational flexibility determines selectivity and antibacterial, antiplasmodial, and anticancer potency of cationic alpha-helical peptides. J. Biol. Chem. 2012, 287, 34120–34133. [Google Scholar] [CrossRef] [PubMed]
  114. Risso, A.; Zanetti, M.; Gennaro, R. Cytotoxicity and apoptosis mediated by two peptides of innate immunity. Cell. Immunol. 1998, 189, 107–115. [Google Scholar] [CrossRef]
  115. Holak, T.A.; Engstrom, A.; Kraulis, P.J.; Lindeberg, G.; Bennich, H.; Jones, T.A.; Gronenborn, A.M.; Clore, G.M. The solution conformation of the antibacterial peptide cecropin A: A nuclear magnetic resonance and dynamical simulated annealing study. Biochemistry 1988, 27, 7620–7629. [Google Scholar] [CrossRef]
  116. Hung, S.C.; Wang, W.; Chan, S.I.; Chen, H.M. Membrane lysis by the antibacterial peptides cecropins B1 and B3: A spin-label electron spin resonance study on phospholipid bilayers. Biophys. J. 1999, 77, 3120–3133. [Google Scholar] [CrossRef]
  117. Merrifield, R.B.; Merrifield, E.L.; Juvvadi, P.; Andreu, D.; Boman, H.G. Design and synthesis of antimicrobial peptides. Ciba Found. Symp. 1994, 186, 5–20, discussion 20-6. [Google Scholar]
  118. Srisailam, S.; Kumar, T.K.; Arunkumar, A.I.; Leung, K.W.; Yu, C.; Chen, H.M. Crumpled structure of the custom hydrophobic lytic peptide cecropin B3. Eur. J. Biochem. 2001, 268, 4278–4284. [Google Scholar] [CrossRef]
  119. Moore, A.J.; Devine, D.A.; Bibby, M.C. Preliminary experimental anticancer activity of cecropins. Pept. Res. 1994, 7, 265–269. [Google Scholar] [PubMed]
  120. Lehmann, J.; Retz, M.; Sidhu, S.S.; Suttmann, H.; Sell, M.; Paulsen, F.; Harder, J.; Unteregger, G.; Stockle, M. Antitumor activity of the antimicrobial peptide magainin II against bladder cancer cell lines. Eur. Urol. 2006, 50, 141–147. [Google Scholar] [CrossRef]
  121. Cruciani, R.A.; Barker, J.L.; Zasloff, M.; Chen, H.C.; Colamonici, O. Antibiotic magainins exert cytolytic activity against transformed cell lines through channel formation. Proc. Natl. Acad. Sci. USA 1991, 88, 3792–3796. [Google Scholar] [CrossRef] [PubMed]
  122. Gauldie, J.; Hanson, J.M.; Shipolini, R.A.; Vernon, C.A. The structures of some peptides from bee venom. Eur. J. Biochem. 1978, 83, 405–410. [Google Scholar] [CrossRef] [PubMed]
  123. Tosteson, M.T.; Alvarez, O.; Hubbell, W.; Bieganski, R.M.; Attenbach, C.; Caporales, L.H.; Levy, J.J.; Nutt, R.F.; Rosenblatt, M.; Tosteson, D.C. Primary structure of peptides and ion channels. Role of amino acid side chains in voltage gating of melittin channels. Biophys. J. 1990, 58, 1367–1375. [Google Scholar] [CrossRef] [Green Version]
  124. Sui, S.F.; Wu, H.; Guo, Y.; Chen, K.S. Conformational changes of melittin upon insertion into phospholipid monolayer and vesicle. J. Biochem. 1994, 116, 482–487. [Google Scholar] [CrossRef]
  125. Henzler Wildman, K.A.; Lee, D.K.; Ramamoorthy, A. Mechanism of lipid bilayer disruption by the human antimicrobial peptide, LL-37. Biochemistry 2003, 42, 6545–6558. [Google Scholar] [CrossRef]
  126. Okumura, K.; Itoh, A.; Isogai, E.; Hirose, K.; Hosokawa, Y.; Abiko, Y.; Shibata, T.; Hirata, M.; Isogai, H. C-terminal domain of human CAP18 antimicrobial peptide induces apoptosis in oral squamous cell carcinoma SAS-H1 cells. Cancer Lett. 2004, 212, 185–194. [Google Scholar] [CrossRef]
  127. Dhople, V.; Krukemeyer, A.; Ramamoorthy, A. The human beta-defensin-3, an antibacterial peptide with multiple biological functions. Biochim. Biophys. Acta 2006, 1758, 1499–1512. [Google Scholar] [CrossRef] [Green Version]
  128. Selsted, M.E.; Harwig, S.S. Determination of the disulfide array in the human defensin HNP-2. A covalently cyclized peptide. J. Biol. Chem. 1989, 264, 4003–4007. [Google Scholar]
  129. Lichtenstein, A.; Ganz, T.; Selsted, M.E.; Lehrer, R.I. In vitro tumor cell cytolysis mediated by peptide defensins of human and rabbit granulocytes. Blood 1986, 68, 1407–1410. [Google Scholar] [Green Version]
  130. Lichtenstein, A.K.; Ganz, T.; Nguyen, T.M.; Selsted, M.E.; Lehrer, R.I. Mechanism of target cytolysis by peptide defensins. Target cell metabolic activities, possibly involving endocytosis, are crucial for expression of cytotoxicity. J. Immunol. 1988, 140, 2686–2694. [Google Scholar]
  131. Lichtenstein, A. Mechanism of mammalian cell lysis mediated by peptide defensins. Evidence for an initial alteration of the plasma membrane. J. Clin. Investig. 1991, 88, 93–100. [Google Scholar] [CrossRef]
  132. Mader, J.S.; Salsman, J.; Conrad, D.M.; Hoskin, D.W. Bovine lactoferricin selectively induces apoptosis in human leukemia and carcinoma cell lines. Mol. Cancer Ther. 2005, 4, 612–624. [Google Scholar] [CrossRef]
  133. Furlong, S.J.; Ridgway, N.D.; Hoskin, D.W. Modulation of ceramide metabolism in T-leukemia cell lines potentiates apoptosis induced by the cationic antimicrobial peptide bovine lactoferricin. Int. J. Oncol. 2008, 32, 537–544. [Google Scholar] [CrossRef] [Green Version]
  134. Eliassen, L.T.; Berge, G.; Leknessund, A.; Wikman, M.; Lindin, I.; Lokke, C.; Ponthan, F.; Johnsen, J.I.; Sveinbjornsson, B.; Kogner, P.; et al. The antimicrobial peptide, lactoferricin B, is cytotoxic to neuroblastoma cells in vitro and inhibits xenograft growth in vivo. Int. J. Cancer 2006, 119, 493–500. [Google Scholar] [CrossRef]
  135. Nakamura, T.; Furunaka, H.; Miyata, T.; Tokunaga, F.; Muta, T.; Iwanaga, S.; Niwa, M.; Takao, T.; Shimonishi, Y. Tachyplesin, a class of antimicrobial peptide from the hemocytes of the horseshoe crab (Tachypleus tridentatus). Isolation and chemical structure. J. Biol. Chem. 1988, 263, 16709–16713. [Google Scholar]
  136. Chen, J.; Xu, X.M.; Underhill, C.B.; Yang, S.; Wang, L.; Chen, Y.; Hong, S.; Creswell, K.; Zhang, L. Tachyplesin activates the classic complement pathway to kill tumor cells. Cancer Res. 2005, 65, 4614–4622. [Google Scholar] [CrossRef]
  137. Adamia, S.; Maxwell, C.A.; Pilarski, L.M. Hyaluronan and hyaluronan synthases: Potential therapeutic targets in cancer. Curr. Drug Targets Cardiovasc. Haematol. Disord. 2005, 5, 3–14. [Google Scholar] [CrossRef]
  138. Rao, A.G. Conformation and antimicrobial activity of linear derivatives of tachyplesin lacking disulfide bonds. Arch. Biochem. Biophys. 1999, 361, 127–134. [Google Scholar] [CrossRef]
  139. Tada, N.; Horibe, T.; Haramoto, M.; Ohara, K.; Kohno, M.; Kawakami, K. A single replacement of histidine to arginine in EGFR-lytic hybrid peptide demonstrates the improved anticancer activity. Biochem. Biophys. Res. Commun. 2011, 407, 383–388. [Google Scholar] [CrossRef] [PubMed]
  140. Barras, D.; Chevalier, N.; Zoete, V.; Dempsey, R.; Lapouge, K.; Olayioye, M.A.; Michielin, O.; Widmann, C. A WXW motif is required for the anticancer activity of the TAT-RasGAP317–326 peptide. J. Biol. Chem. 2014, 289, 23701–23711. [Google Scholar] [CrossRef] [PubMed]
Sample Availability: Not applicable.
Figure 1. Overview of the structural diversity of three classes of anticancer peptides. Each structure is labeled by its common name followed by the Protein Data Bank (PDB ID) in parenthesis on the subsequent line. In cases where the PDB ID was not available, the SWISS-MODEL server (available at: https://swissmodel.expasy.org/) was used to construct the structure.
Figure 1. Overview of the structural diversity of three classes of anticancer peptides. Each structure is labeled by its common name followed by the Protein Data Bank (PDB ID) in parenthesis on the subsequent line. In cases where the PDB ID was not available, the SWISS-MODEL server (available at: https://swissmodel.expasy.org/) was used to construct the structure.
Molecules 24 01973 g001
Figure 2. Schematic framework of ACPred.
Figure 2. Schematic framework of ACPred.
Molecules 24 01973 g002
Figure 3. ROC curve of RF (top and bottom left) and SVM (top and bottom right) models as assessed by 5-fold cross-validation (top left and right) and jackknife test or leave-one-out cross-validation (bottom left and right).
Figure 3. ROC curve of RF (top and bottom left) and SVM (top and bottom right) models as assessed by 5-fold cross-validation (top left and right) and jackknife test or leave-one-out cross-validation (bottom left and right).
Molecules 24 01973 g003
Figure 4. Heat map of the mean decrease of Gini index of dipeptide compositions.
Figure 4. Heat map of the mean decrease of Gini index of dipeptide compositions.
Molecules 24 01973 g004
Figure 5. Screenshots of the ACPred web server before (A) and after (B) submission of sequence data for prediction.
Figure 5. Screenshots of the ACPred web server before (A) and after (B) submission of sequence data for prediction.
Molecules 24 01973 g005
Table 1. Summary of existing methods for predicting anticancer peptides.
Table 1. Summary of existing methods for predicting anticancer peptides.
Method (Year)Classifier aSequence Features bInterpretableWeb Server
AntiCP (2013)SVMAAC, DPC, BPNo
Hajisharifi et al. (2014)SVMPseACC, LAKNo
ACPP (2015)SVMPRMNoc
iACP (2016)SVMg-gap DPCYes
Li and Wang (2016)SVMAAC, RACC, acACSNo
iACP-GAEnsC (2017)Ensemble methodPse-g-gap DPC, Am-PseAAC, RACCNo
MLACP (2017)RFAAC, ATC, DPC, PCPYes
SAP (2018)SVMg-gap DPCNo
TargetACP (2018)SVMCPSR, SAAC, PsePSSMNo
ACPred (this study)SVMAAC, DPC, PCP, PseAAC, Am-PseAACYes
a RF: random forest, SVM: support vector machine. b AAC: amino acid composition, ATC: atomic composition, acACS: auto covariance of the average chemical shift, Am-PseAAC: amphiphilic pseudo amino acid composition, BP: binary profile, CPDR: composite protein sequence representation, DPC: dipeptide composition, g-gap DPC: G-Gap dipeptide composition, LAK: local alignment kernel, PCP: Physicochemical properties, PseACC: Pseudo amino acid composition, Pse-g-gap DPC: Pseudo G-Gap dipeptide composition, PRM: protein relatedness measure, RACC: reduce amino acid composition, SAAC: split amino acid composition. c The web server is not accessible.
Table 2. Performance comparison of SVM and RF with various types of sequence features over five-fold cross-validation.
Table 2. Performance comparison of SVM and RF with various types of sequence features over five-fold cross-validation.
FeatureClassifierAc (%)Sn (%)Sp (%)MCCauROC
AACSVM92.6983.9498.540.8500.977
RF91.2392.8090.320.8170.958
DPCSVM83.92100.0078.850.6870.942
RF87.1491.8984.850.7330.944
PCPSVM84.8063.5099.020.6980.938
RF83.6390.1080.910.6610.872
PseAACSVM92.9884.6798.540.8560.990
RF92.1193.6591.200.8350.959
Am-PseAACSVM95.0387.59100.000.8990.995
RF92.40100.0088.750.8480.974
AAC + PseAACSVM93.5785.4099.020.8690.991
RF92.9896.6990.950.8550.964
AAC + Am-PseAACSVM95.3289.0599.510.9040.994
RF93.2897.5090.990.8620.969
PseAAC + Am-PseAACSVM94.4486.8699.510.8870.994
RF92.6998.2889.820.8510.967
AAC + PseAAC + Am-PseAACSVM94.1586.1399.510.8810.993
RF92.9898.2990.220.8570.972
Table 3. Performance comparison of SVM and RF with various types of sequence features over jackknife test.
Table 3. Performance comparison of SVM and RF with various types of sequence features over jackknife test.
FeatureClassifierAc (%)Sn (%)Sp (%)MCCauROC
AACSVM92.9883.9499.020.8570.978
RF91.2392.8090.320.8170.959
DPCSVM85.0998.8680.320.7060.941
RF86.8489.6685.400.7250.947
PCPSVM84.8063.5099.020.6980.937
RF83.6388.5781.440.6590.868
PseAACSVM92.9884.6798.540.8560.990
RF93.2897.5090.990.8620.960
Am-PseAACSVM95.0387.59100.000.8990.995
RF92.4099.1289.080.8470.969
AAC + PseAACSVM93.5785.4099.020.8690.990
RF93.2898.3190.630.8630.962
AAC + Am-PseAACSVM95.6189.7899.510.9100.994
RF93.5798.3291.030.8690.967
PseAAC + Am-PseAACSVM93.8685.4099.510.8750.992
RF93.5799.1590.670.8700.959
AAC + PseAAC + Am-PseAACSVM94.7487.5999.510.8930.994
RF92.9899.1389.870.8580.973
Table 4. Performance comparison of the proposed ACPred model with existing methods.
Table 4. Performance comparison of the proposed ACPred model with existing methods.
Method aAc (%)Sn (%)Sp (%)MCC
Hajisharifi et al. b92.6889.7085.180.78
iACP b95.0689.8698.540.90
iACP-GAEnsC b96.4595.3697.540.91
TargetACP b96.2294.2097.570.92
TargetACP c98.7899.0298.540.97
ACPred b95.6189.7899.510.91
ACPred-modified c97.5696.0899.020.95
a Results were reported from the work of TargetACP. b Results were performed on the benchmark dataset consisting of 138 ACPs and 205 non-ACPs. c Results were performed on the balanced dataset consisting of 205 ACPs and 205 non-ACPs by using the SMOTE technique on the benchmark dataset.
Table 5. Amino acid compositions (%) of anticancer and non-anticancer peptides along with their difference as well as MDGI values. The rank of each amino acid amongst their 20 amino acids are shown in parenthesis for AAC difference and MDGI.
Table 5. Amino acid compositions (%) of anticancer and non-anticancer peptides along with their difference as well as MDGI values. The rank of each amino acid amongst their 20 amino acids are shown in parenthesis for AAC difference and MDGI.
Amino AcidACP (%)Non-ACP (%)Differencep-ValueMDGI
A-Ala7.62311.005−3.383 (7)<0.056.41 (11)
C-Cys3.9068.015−4.109 (5)<0.0519.71 (2)
D-Asp2.4173.418−1.002 (15)<0.053.85 (15)
E-Glu1.7073.523−1.816 (12)<0.056.88 (9)
F-Phe6.8232.414.413 (2)<0.056.12 (12)
G-Gly1.9754.123−2.148 (11)<0.056.56 (10)
H-His1.5365.798−4.262 (4)<0.057.80 (7)
I-Ile10.0726.983.092 (9)<0.0510.13 (5)
K-Lys2.5421.6510.892 (17)0.08629.54 (1)
L-Leu8.0993.7394.36 (3)<0.057.84 (6)
M-Met9.83113.888−4.057 (6)<0.054.34 (14)
N-Asn11.4973.9647.533 (1)<0.054.98 (13)
P-Pro0.9052.224−1.319 (13)<0.057.37 (8)
Q-Gln5.3852.7112.674 (10)<0.0515.89 (3)
R-Arg4.2117.495−3.283 (8)<0.0511.42 (4)
S-Ser6.5375.8320.705 (19)0.2453.50 (16)
T-Thr3.7814.704−0.923 (16)0.0983.13 (17)
V-Val2.2581.5600.698 (20)0.0832.96 (18)
W-Trp2.2441.4230.821 (18)<0.052.51 (20)
Y-Tyr6.655.5391.111 (14)0.0912.90 (19)
Table 6. Ten top-ranked physiocochemical properties from the AAindex having the highest MDGI values.
Table 6. Ten top-ranked physiocochemical properties from the AAindex having the highest MDGI values.
RankAAindexCategorized PropertyDescriptionMDGI
1ARGP820101HydrophobicityHydrophobicity index (Argos et al., 1982)1.51
2ARGP820102HydrophobicitySignal sequence helical potential (Argos et al., 1982)1.40
3BHAR880101HydrophobicityAverage flexibility indices (Bhaskaran-Ponnuswamy, 1988)1.08
4ARGP820103HydrophobicityMembrane-buried preference parameters (Argos et al., 1982)1.04
5BEGF750102Beta propensityConformational parameter of beta-structure (Beghin-Dirkx, 1975)1.00
6BEGF750101Alpha and turn propensitiesConformational parameter of inner helix (Beghin-Dirkx, 1975)0.92
7BIGC670101Physicochemical propertiesResidue volume (Bigelow, 1967)0.91
8BEGF750103Alpha and turn propensitiesConformational parameter of beta-turn (Beghin-Dirkx, 1975)0.87
9BIOV880102HydrophobicityInformation value for accessibility; average fraction 23% (Biou et al., 1988)0.85
10ISOY800107HydrophobicityNormalized relative frequency of double bend (Isogai et al., 1980)0.83
Table 7. Eight if-then rules for the prediction of anticancer peptides using random forest and amino acid composition.
Table 7. Eight if-then rules for the prediction of anticancer peptides using random forest and amino acid composition.
No.RuleCovered SamplesMisclassified SampleAc (%)
1G > 0.041 and I > 0.0615 and L ≤ 0.1715 and K > 0.0385 and M ≤ 0.02748197.92
2R ≤ 0.0515 and Q ≤ 0.026 and K > 0.09448197.92
3C > 0.1145 and P ≤ 0.073350100.00
4L ≤ 0.093 and F > 0.0715 and S ≤ 0.15233487.88
5A ≤ 0.0145 and Q ≤ 0.02626484.62
6R ≤ 0.0665 and E ≤ 0.044 and H > 0.05221480.95
7I > 0.108 and K > 0.05542685.71
8E ≤ 0.0545 and G > 0.0365 and K > 0.108 and M ≤ 0.0449197.96

Share and Cite

MDPI and ACS Style

Schaduangrat, N.; Nantasenamat, C.; Prachayasittikul, V.; Shoombuatong, W. ACPred: A Computational Tool for the Prediction and Analysis of Anticancer Peptides. Molecules 2019, 24, 1973. https://doi.org/10.3390/molecules24101973

AMA Style

Schaduangrat N, Nantasenamat C, Prachayasittikul V, Shoombuatong W. ACPred: A Computational Tool for the Prediction and Analysis of Anticancer Peptides. Molecules. 2019; 24(10):1973. https://doi.org/10.3390/molecules24101973

Chicago/Turabian Style

Schaduangrat, Nalini, Chanin Nantasenamat, Virapong Prachayasittikul, and Watshara Shoombuatong. 2019. "ACPred: A Computational Tool for the Prediction and Analysis of Anticancer Peptides" Molecules 24, no. 10: 1973. https://doi.org/10.3390/molecules24101973

APA Style

Schaduangrat, N., Nantasenamat, C., Prachayasittikul, V., & Shoombuatong, W. (2019). ACPred: A Computational Tool for the Prediction and Analysis of Anticancer Peptides. Molecules, 24(10), 1973. https://doi.org/10.3390/molecules24101973

Article Metrics

Back to TopTop