- Open Access
Prediction of anti-inflammatory proteins/peptides: an insilico approach
Journal of Translational Medicine volume 15, Article number: 7 (2017)
The current therapy for inflammatory and autoimmune disorders involves the use of nonspecific anti-inflammatory drugs and other immunosuppressant, which are often accompanied with potential side effects. As an alternative therapy, anti-inflammatory peptides are recently being exploited as anti-inflammatory agents for treatment of various inflammatory diseases such as Alzheimer’s disease and rheumatoid arthritis. Thus, understanding the correlation between amino acid sequence and its potential anti-inflammatory property is of great importance for the discovery of novel and efficient anti-inflammatory peptide-based therapeutics.
In this study, we have developed a prediction tool for the classification of peptides as anti-inflammatory epitopes or non anti-inflammatory epitopes. The training was performed using experimentally validated epitopes obtained from Immune epitope database and analysis resource database. Different sequence-based features and their hybrids with motif information were employed for development of support vector machine-based machine learning models. Similarly, machine learning models were also constructed using random forest.
The composition and terminal residue conservation analysis of peptides revealed the dominance of leucine, serine, tyrosine and arginine residues in anti-inflammatory epitopes as compared to non anti-inflammatory epitopes. Similarly, the anti-inflammatory epitopes specific motifs were found to be rich in hydrophobic and polar residues. The hybrid of tripeptide composition-based support vector machine model and motif yielded the best performance on 10-fold cross validation with an accuracy of 78.1% and MCC of 0.58. The same displayed an accuracy of 72% and MCC of 0.45 on validation dataset, rejecting any possibility of over-fitting. The tripeptide composition-based random forest model displayed an accuracy of 0.8 and MCC of 0.59 on 10-fold cross validation, however, the accuracy (0.68) and MCC (0.31) was lower as compared to support vector machine model on validation dataset. Thus, the support vector machine model is implemented as the default model and an additional option of using the random forest model is provided.
The prediction models along with tools for epitope mapping and similarity search have been provided as a web server which is freely accessible at http://metagenomics.iiserb.ac.in/antiinflam/.
For maintenance of immune homeostasis and restricting the occurrence of aggravated inflammation and autoimmunity, the induction of immune tolerance is crucial . The effect of multiple mechanisms working in concert is necessary for the upkeep of the state of tolerance. Several endogenous peptides identified during inflammatory responses have emerged as anti-inflammatory agents that can be employed in new therapies for inflammatory and autoimmune conditions . The immunotherapeutic ability of these anti-inflammatory peptides in the inhibition of antigen-specific T(H)1-driven responses, and in the generation of regulatory T cells has numerous clinical applications .
Several synthetic and natural peptides are known to inhibit the signal transduction pathways for expression of inflammatory cytokines. For example, intranasal treatment of amyloid-beta (Abeta) peptide in mice with neuropathological symptoms of Alzheimer’s disease showed a decreased A-beta plaque burden along with anti-inflammatory cytokines IL-4, IL-10 and tumor growth factor-beta . Another neuropeptide: vasoactive intestinal peptide (VIP) is also known to reduce the inflammatory and autoimmune components of rheumatoid arthritis in the experimental model . Similarly, there are pieces of evidence showing production of pro-inflammatory cytokines and microglial activation during age-related cognitive deficits in the hippocampal region. A mimetic peptide FG Loop(FGL) derived from the neural cell adhesion molecule (NCAM) is shown to work as a novel anti-inflammatory agent in such age-related disorders . Furthermore, nonapeptides, Anti inflammin 1 and 2, can potentially reduce leukocyte trafficking in host defense and inflammation by inhibiting PMN adhesion to HCAEC by attenuating up-regulation of CD11/CD18 expression on leukocytes . Moreover, certain synthetic peptides are known to be involved in inhibition of early signal transduction pathways of inflammatory cytokines. One notable example is the role of intravesical RDP58 in interstitial cystitis (IC) inflammation-induced neurotrophic changes in the bladder. Mice models treated with RDP58 shows near complete abolition of LPS induced production of TNF-α in the bladder .
Taken together, the availability of experimental data allows for the computational correlation between amino acid sequences and their anti-inflammatory properties. To the best of the authors’ knowledge, there has been no study that investigates the sequence based signature for the anti-inflammatory behavior of target peptide sequences. Understanding the association between specific amino acid sequences and the anti-inflammatory properties conferred by them is likely to facilitate their application for therapeutic purposes. One such example is of altered peptide ligands (APLs) and their role in therapeutic modulation of T cell function . The aim of this study was to examine the amino acid sequences of experimentally validated anti-inflammatory epitopes (AIEs) and non anti-inflammatory epitopes (NAIEs) and develop a machine learning based classification method that incorporates the sequence based features to predict the anti-inflammatory nature of various peptides and proteins.
Preparation of dataset
An unambiguous and well-curated dataset is required for machine-learning based model building. Thus, the positive and negative assaying epitopes were retrieved from the immune epitope database . A peptide was considered as anti-inflammatory if it was shown to induce any one of the anti-inflammatory cytokines [IL-10, IL-13, IL-22, IL-1RA, TGFβ, IFN-α/β] . 863 epitopes were reported to show a positive assay, whereas, 1261 peptides tested negative for anti-inflammatory cytokines and were classified as NAIEs. The length of these epitopes ranged between 4 and 30 amino acids residues. Out of the total input dataset, 80% of both AIEs and NAIEs were assigned as training dataset and the remaining 20% constituted the validation dataset. Resultantly, a total of 690 AIEs (positive data) and 1009 NAIEs (negative data) constituted the training dataset, and the validation dataset contained 173 AIEs (positive data) and 253 NAIEs (negative data) (Fig. 1). Two datasets were constructed on which the various analysis were carried out in this manuscript, namely, (1) realistic dataset with all available positive and negative epitopes, and (2) balanced dataset including all positive epitopes (690 positives) and randomly chosen equal number of negative epitopes (690 negatives).
Input features for machine learning
Amino acid composition
Amino acid composition (AAC) is a widely used feature in composition-based peptide classification methods [12, 13]. It is the percentage of each amino acid in a given length of the peptide sequence. Each peptide/protein can be represented as a vector of 20 features (composition of 20 naturally occurring amino acids).
where, AAC is the amino acid composition of an amino acid (i) which is one of the 20 amino acids.
Several immune epitope prediction algorithms and numerous composition-based algorithms use dipeptide composition (DPC) for classification [14, 15]. DPC is the percentage composition of each of the 400 possible dipeptides formed by 20 amino acids. Unlike AAC, DPC carries the additional information of local arrangement in a peptide/protein since it is the pair of adjacently placed amino acids. DPC was calculated using the following equation.
where Dp(n) is a dipeptide n out of 400 dipeptides.
Tripeptide composition (TPC) is the percentage composition of each of the 8000 possible dipeptides formed by 20 amino acids and was calculated using the following formula. It can also be used for composition-based classification .
where TPC(n) is a tripeptide n out of 8000 tripeptides.
Identification of functional motifs in amino acid sequences is one of the widely used techniques in the functional annotation of peptides/proteins. Several previous immunoinformatics studies have examined the immunologically relevant motifs in peptides for various purposes [17–19]. To identify the motifs specific to AIEs by comparing it with NAIEs, MERCI software was used . The identification of motifs in positive data by comparing with negative data yields motifs exclusive for the positive data. Motifs specific for NAIEs were also identified by MERCI software in a similar manner. The maximum length of the motif was set to 9 amino acids because the core size of both MHC I and II equals nine amino acids [21, 22]. The gap length was set to the default length of 1.
Hybrid models of AAC, DPC, and TPC with MERCI-based motifs were also constructed in this study [23, 24]. The presence of anti-inflammatory and non anti-inflammatory specific motifs in every peptide in the training dataset was examined. If the epitope contained the anti-inflammatory motifs, a weightage of +1 was added to its AAC, DPC or TPC based SVM score. Similarly, if the epitope is positive for non anti-inflammatory motifs, the weight was added as −1 for the same.
Machine-learning-based prediction models
Caret package in R was used to evaluate the performance of various machine learning methods such as support vector machine (SVM), Random Forest (RF), k-nearest neighbors (kNN), linear discriminant analysis (LDA) and classification and regression tree (CART). Based on the results, SVM and RF were used as the machine learning algorithms for building predictive models through SVMlight package available at http://svmlight.joachims.org/ and randomForest package available at https://cran.r-project.org/web/packages/randomForest/index.html. Machine learning algorithms have also been used in previous studies for developing tools for immune epitope prediction with good performances [15, 25–28].
Performance evaluation of models
The comparison of different machine learning methods to select the best method is essential in predictive modelling. A 10-fold cross-validation technique was employed for different learning methods and evaluated on a separate unseen validation dataset. To carry out the 10-fold cross-validation, the data was divided into ten parts of which nine parts were used for training and the tenth segment was used for testing. The procedure is repeated ten times to test the learning methods on each part. The result of all ten predictions was taken together for evaluating the performance of various models built by different methods. For measuring the performance of each model, threshold dependent and threshold independent parameters were used. For all models, area under curve (AUC), a threshold independent parameter was calculated using PERF . Threshold-dependent parameters such as sensitivity (Sen), specificity (Spec), accuracy (Acc) and Matthews’s correlation coefficient (MCC) were also calculated for all models. The following equations were used for calculating the threshold dependent parameters.
In addition, a fivefold cross validation study was also carried out and performance was compared with the 10-fold cross validation study.
(TN = true negative, FP = false positive, TP = true positive, FN = false negative; Sens = sensitivity; Spec = specificity; Acc = accuracy)
To carry out the compositional analysis of AIEs and NAIEs, the frequencies of all amino acids (AAC), dipeptides (DPC), and tripeptides (TPC) were calculated in the peptides included in the two sets. The AAC analysis of AIEs revealed higher abundance of Leu, Ser, Tyr, Gln and Arg (Welch’s t test, p value ≤0.05) (Fig. 2). Similarly for NAIEs, Ala, Gly, Asp, Thr and Pro were observed in higher abundance (Welch’s t test, p value ≤0.05).The positional conservation of amino acids in AIEs and NAIEs was examined using a two sample logo (TSL) analysis (Fig. 3). Leu, Tyr, Ser, Arg and Glu were found as highly conserved at N-terminal of AIEs, whereas, Leu, Gln, Ser and Arg were found at C-terminal. In contrast, Ala, Gly, Val, Pro and Asp were highly conserved at N-terminal, whereas, Pro, Gly and Ala were found conserved at C-terminal of NAIEs.
In DPC, from the 400 possible dipeptides, 54 significant dipeptides (p value ≤0.01, Welch’s t-test) were found significantly different in AIEs as compared to NAIEs (Fig. 4). AIEs were found to be rich in Leu–Leu (LL), Ser–Leu (SL), Leu–Glu (LE), Leu–Lys (LK), Arg–Leu (RL) and Tyr–Leu (YL) dipeptides, whereas, NAIEs were rich in Ala–Ala (AA), Ala–Thr (AT), Gly–Ala (GA), Gly–Pro (GP), Pro–Gly (PG), Pro–Pro (PP), Leu–Pro (LP) and Thr–Gly (TG) dipeptides.
Similarly, the TPC analysis also revealed several tripeptides which were highly abundant in AIEs as compared to NAIEs. A total of 122 tripeptides were found significantly abundant (p value ≤0.01, Welch’s t-test) in AIEs (Additional file 1: Table S1). In particular, Ser–Leu–Ser (SLS), Leu–Lys–Leu (LKL), Leu–Arg–Leu (LRL) etc. are highly abundant in AIEs and their constituent dipeptides and amino acids are also abundant in AAC and DPC analysis.
HLA allele distribution
An analysis of the distribution of epitopes from IEDB database against different HLA alleles was also carried out (Fig. 5). Among the most contrasting results, most (50) of the peptides showed anti-inflammatory response against HLA-DRB1*04:01 in assays, whereas, only 14 peptides showed non anti-inflammatory response against the same HLA allele. Similarly, the assays against H2-IAd allele included 25 unique NAIEs compared to 10 unique AIEs tested against the same HLA allele.
Inflammatory responses arise due to the induction of T cells upon activation by antigen-presenting cells. In anti-inflammatory response, MHC molecules on antigen presenting cells present the bounded peptides and initiate the response. Peptide-MHC interaction occurs through specific amino acids in the form of motifs. Thus, the presence of these specific motifs can be utilized as a marker for better classification of peptides. MERCI software was used to predict specific motifs by comparing positive and negative datasets as explained in the “Methods” section. 13 different specific motifs for AIEs were observed as shown in Table 1. The “hydrophobic hydrophobic tiny F charged” motif was found in 22 unique AIEs and had the highest observed motif coverage among all motifs.
Machine learning-based classification
It is apparent from the earlier analysis that AIEs and NAIEs differ in several composition-based features. Therefore, these differences can be exploited for the classification of epitopes. Machine learning-based classification was implemented to exploit the sequence-specific features for the classification and prediction of AIEs and NAIEs.
Comparison of machine learning methods
Different machine learning algorithms such as SVM, RF, kNN, LDA and CART were implemented using caret package in R for comparing their performances on the same data. Performances of SVM and RF were found better as compared to other machine learning techniques, and hence, these two machine learning methods were optimized further using amino acid, dipeptide and tripeptide composition as input features (Additional file 1: Text S1, Additional file 2: Figure S1, Additional file 3: Figure S2, Additional file 4: Figure S3, Additional file 5: Figure S4, Additional file 6: Figure S5, Additional file 7: Figure S6, Additional file 8: Figure S7, Additional file 9: Figure S8, Additional file 10: Figure S9, Additional file 11: Figure S10 and Additional file 12: Figure S11).
Performance of SVM models
Amino acid composition-based models
Compositional analysis revealed that AIEs and NAIEs show compositional differences. These characteristics were utilized to classify peptides into anti-inflammatory or non anti-inflammatory using machine learning models. SVM model with gamma parameter (g) = 0.005, trade-off factor (c) = 1 and cost factor (j) = 1 produced the best classification in AAC-based models. It showed MCC of 0.37, accuracy of 68.16% and the area under curve (AUC) was found to be 0.74 (Fig. 6).
Dipeptide composition-based models
Dipeptide composition was also used for building the predictive machine learning model using 400 sized vectors representing all possible dipeptides. SVM model with gamma parameter (g) = 0.005, trade-off factor (c) = 1 and cost factor (j) = 3 was the best performing model with an accuracy of 70.98, MCC of 0.44 and an area under curve (AUC) of 0.79 (Fig. 6).
Tripeptide composition-based models
The local order in the form of DPC as feature yielded better predictive ability of SVM model as compared to the results obtained using AAC. However, to further evaluate the effect of local order in the form of tripeptide, TPC as input feature was used to build the SVM model. The best model performed with an accuracy of 75.04%, MCC of 0.51 and AUC of 0.81 on gamma parameter (g) = 0.001, trade-off factor (c) = 4 and cost factor (j) = 1 (Fig. 6).
Performance of RF models
RF models were constructed at optimized mtry and ntree for each input feature separately and performance of these models were evaluated using 10-fold cross validation and on the validation set (Additional file 1: Text S1). Performance of RF model was comparable to SVM models. For amino acid, dipeptide and selected tripeptide composition features SVM showed MCC of 0.37, 0.44 and 0.51, whereas, RF showed MCC of 0.39, 0.43 and 0.59 (Table 2 and Additional file 1: Table S2). The RF model constructed using selected tripeptide composition features displayed higher performance as compared to SVM, using 10-fold cross validation. However, on the validation set, the same RF model showed the MCC of 0.31 which is much lower as compared to the MCC (0.43) shown by SVM model (Table 3 and Additional file 1: Table S3). Hence, SVM models were implemented as the default model for the construction of hybrid models. Performance of RF models was also evaluated on the balanced dataset using 10-fold cross validation and on validation set (Additional file 1: Tables S4 and S5). Detailed description of machine learning comparison and RF models optimization and performance has been provided in Additional file 1: Text S1.
TPC-based SVM model yielded the best accuracy among the three composition-based models. However, to further improve the efficiency of models, hybrid models were constructed for AAC, DPC, and TPC based models by combining each one with the motif features. The hybrid of AAC-based model (AAC_HYB) with motif features performed with 71.1% accuracy, whereas, the hybrid of DPC-based model (DPC_HYB) with motif features showed an accuracy of 74.34%. The hybrid of TPC-based model with motif features (TPC_HYB) was the best performing among all hybrid models. It showed the highest accuracy of 78.1%, MCC of 0.58 and an area under curve (AUC) of 0.86 (Fig. 6).
Performance on validation dataset
The performance of all models was evaluated on the validation data (Table 3 and Fig. 6). This evaluation was carried out in addition to cross-validation to ensure that the obtained results were not due to over-optimization. TPC_HYB model which showed the highest accuracy in 10-fold cross validation also performed well on the validation dataset (72%). These results attest that the observed results were not due to over-optimization of data.
Furthermore, all the inputs such as amino acid, dipeptides, tripeptides etc., were processed for feature selection for composition and hybrid models by CfsSubsetEval algorithm of WEKA package. This feature selection is generally carried out to reduce the input data size without compromising on the performance. The final number of input features and performances of model with reduced input features are mentioned in Additional file 1: Table S6. It is apparent from Additional file 1: Table S7 that the performances of all models were reduced on validation dataset, with the reduction of feature vector size, however, in case of TPC_HYB model, the performance showed a marginal decline on the validation dataset, and thus, this model with reduced features was integrated in the webserver for efficient analysis.
In addition to 10-fold cross validation, models were also subjected to a fivefold cross validation analysis. The performance of models using fivefold cross validation also showed similar results. The TPC_HYB model displayed the maximum accuracy of 76% and also performed better than the other models for sensitivity and specificity values. It showed an MCC value of 0.56 and an AUC measure of 0.85 (Additional file 1: Table S8). Resultantly, the TPC_HYB model was incorporated in the webserver for epitope prediction.
Balanced dataset which contained 690 samples from both positive and negative dataset was also used to examine the performance of SVM models. The TPC_HYB model displayed the maximum accuracy of 77.8% and MCC value of 0.56. The sensitivity and specificity of this model was also higher than other models (Additional file 1: Table S9).
Web server and tools
A webserver was constructed for hosting a number of tools developed in this study for predicting the anti-inflammatory nature of peptide/protein antigens. These tools employ the validated best performing model developed in the study. The web server also incorporates several other analysis tools as described below.
This module is developed for predicting the anti-inflammatory nature of small length peptides. The user can submit a single peptide or a batch of peptides with lengths ranging from 4 to 30 amino acids. A threshold option is also provided which can be tuned to increase the stringency of positive prediction. Queries are passed through anti-inflammatory epitope prediction pipeline and the results are displayed as a table which shows the prediction of a peptide as anti-inflammatory or non anti-inflammatory based on SVM or RF scores.
Every peptide sequence in the results table is provided with option for its virtual screening and design which is an additional feature that screens the variants of the peptides for anti-inflammatory property. Variants are created by substituting each amino acid position by all the 19 other naturally occurring amino acids. This option enables the user to check the anti-inflammatory nature of peptides that are similar to the query peptides.
To determine the antigenic regions in a full-length protein which can induce anti-inflammatory responses, this module has been designed to allow the submission of query proteins with longer lengths. This tool divides the protein sequences into smaller peptides and runs through the same prediction pipeline. The user can also select the desired window length of peptides for prediction. Virtual screening option is also provided similar to the peptide prediction module.
In addition to the above tools, the user may wish to map the experimentally validated anti-inflammatory antigenic regions on a query protein. This module map experimentally validated AIEs from IEDB database on the protein query sequence. The results are displayed as an alignment with links to the related assay provided by IEDB database.
Unlike the epitope mapping tool, this module identifies the similar regions in the query peptide/protein and experimentally validated AIEs. Smith–Watermann search of query sequence in IEDB database is carried for displaying the top hits. Results are shown as alignments with links to related assay in IEDB database.
Anti-inflammatory drugs have widespread applications in the field of therapeutics because of their role in the reduction of swelling and inflammation. The novel peptides with anti-inflammatory properties may also be used as therapeutic agents against several diseases such as Alzheimer’s, rheumatoid arthritis as well as gut infections. Inflammatory response at the site of infection is caused due to the production of prostaglandin, thromboxane, etc. along with a release of histamine and kinin. Though it has been shown that nature of activating signal, timing, and sequence of cytokine action as well as experimental models affect their pro- or anti-inflammatory properties, IL4, IL-10, IL-13, IFN-alpha and transforming growth factor-beta are widely recognized as anti-inflammatory cytokines . Studies have shown that anti-inflammatory drugs function by down-regulating the expressions of 5-LOX, COX-2, MMPs and NF-κB . These drugs show low efficacy in curing the inflammation and leads to hypersensitivity reactions and deterioration of the immune system . Due to the limitation of small molecule drugs, peptide-based drugs appear as better candidates in curing inflammation . Thus, developing an in silico method for predicting the anti-inflammatory properties of newly synthesized or recently discovered amino acid sequence (peptide) is very valuable.
For the identification of anti-inflammatory targets, positive (anti-inflammatory) and negative (non anti-inflammatory) samples data was curated from immune epitope database which is a well-known reference database and contain the experimentally validated data. Using a good input quality data is a key for the application of machine learning techniques.
Compositional analysis revealed the presence of selected amino acids in abundance in AIEs, which is supported by other independent studies. For example, an amino acid stretch of 8-16 of SV-IV, which is rich in Ser and Gln residues, has been reported to exhibit anti-inflammatory properties . Motif based analysis helped to determine the patterns in binding sites of target peptides on MHC molecules. Annotation of MERCI motifs from unique AIEs revealed that ‘hydrophobic hydrophobic tiny F charged’ had the maximum coverage. From this study, it was also observed that certain HLA alleles such as HLA-DRB1*04:01 alleles were more frequent in anti-inflammatory positive assays. Interestingly, the HLA-DRB1*04:01 allele has been associated with an increased occurrence of rheumatoid arthritis in the Caucasian population .
Furthermore, three different types of compositional analysis including amino acid, dipeptide and tripeptide composition were carried out to obtain the information of signature sequences, of which TPC_HYB (tripeptide hybrid with motifs derived weightage) was selected and included as the final model for webserver as a default prediction model. TPC (tripeptide composition) based RF model was also included at the webserver as an additional prediction model. The availability of two machine learning models provides flexibility and options to the users to predict epitopes using two different methods.
The inherent ability of anti-inflammatory molecules to induce anti-inflammatory cytokines has the potential to play a major role in immunology and peptide therapeutic applications. Hence, an attempt has been made to identify sequence-based signatures of anti-inflammatory peptides and a machine learning tool has been developed for the prediction of anti-inflammatory epitopes. The observed accuracy values attest the performance of the tool and to our knowledge, it is the only available tool for the prediction of anti-inflammatory epitopes. Thus, the identification of anti-inflammatory target peptide sequences using this computational server provides valuable leads for experimental validations for the scientific community. The developed computational tools for investigation of anti-inflammatory sites have been made freely available for academic usage as a web server.
amino acid composition
area under curve
immune epitope database
human leukocyte antigen
Matthews correlation coefficient
major histocompatibility factor
non anti-inflammatory epitopes
support vector machine
two sample logo
linear discriminant analysis
classification and regression tree
Sun S-C, Chang J-H, Jin J. Regulation of nuclear factor-κB in autoimmunity. Trends Immunol. 2013;34(6):282–9.
Gonzalez-Rey E, Anderson P, Delgado M. Emerging roles of vasoactive intestinal peptide: a new approach for autoimmune therapy. Ann Rheumc Dis. 2007;66(suppl 3):iii70–6.
Delgado M, Ganea D. Anti-inflammatory neuropeptides: a new class of endogenous immunoregulatory agents. Brain Behav Immun. 2008;22(8):1146–51.
Weiner HL, Lemere CA, Maron R, Spooner ET, Grenfell TJ, Mori C, Issazadeh S, Hancock WW, Selkoe DJ. Nasal administration of amyloid-β peptide decreases cerebral amyloid burden in a mouse model of Alzheimer’s disease. Ann Neurol. 2000;48(4):567–79.
Delgado M, Abad C, Martinez C, Leceta J, Gomariz RP. Vasoactive intestinal peptide prevents experimental arthritis by downregulating both autoimmune and inflammatory components of the disease. Nat Med. 2001;7(5):563–8.
Downer EJ, Cowley TR, Cox F, Maher FO, Berezin V, Bock E, Lynch MA. A synthetic NCAM-derived mimetic peptide, FGL, exerts anti-inflammatory properties via IGF-1 and interferon-γ modulation. J Neurochem. 2009;109(5):1516–25.
Zouki C, Ouellet S, Filep JG. The anti-inflammatory peptides, antiflammins, regulate the expression of adhesion molecules on human leukocytes and prevent neutrophil adhesion to endothelial cells. FASEB J. 2000;14(3):572–80.
Gonzalez RR, Fong T, Belmar N, Saban M, Felsen D, Te A. Modulating bladder neuro-inflammation: RDP58, a novel anti-inflammatory peptide, decreases inflammation and nerve growth factor production in experimental cystitis. J Urol. 2005;173(2):630–4.
Bielekova B, Martin R. Antigen-specific immunomodulation via altered peptide ligands. J Mol Med. 2001;79(10):552–65.
Zhang Q, Wang P, Kim Y, Haste-Andersen P, Beaver J, Bourne PE, Bui H-H, Buus S, Frankild S, Greenbaum J. Immune epitope database analysis resource (IEDB-AR). Nucleic Acids Res. 2008;36(suppl 2):W513–8.
Marie C, Pitton C, Fitting C, Cavaillon J. Regulation by anti-inflammatory cytokines (IL-4, IL-10, IL-13, TGFβ) of interleukin-8 production by LPS-and/or TNFα-activated human polymorphonuclear cells. Mediat Inflamm. 1996;5(5):334–40.
Gupta A, Kapil R, Dhakan DB, Sharma VK. MP3: a software tool for the prediction of pathogenic proteins in genomic and metagenomic data. PLoS ONE. 2014;9(4):e93907.
Sharma AK, Gupta A, Kumar S, Dhakan DB, Sharma VK. Woods: a fast and accurate functional annotator and classifier of genomic and metagenomic sequences. Genomics. 2015;106(1):1–6.
Zhang J, Zhao X, Sun P, Gao B, Ma Z: Conformational B-cell epitopes prediction from sequences using cost-sensitive ensemble classifiers and spatial clustering. BioMed Res Int. 2014;2014.
Gupta S, Ansari HR, Gautam A, Raghava GP. Identification of B-cell epitopes in an antigen for inducing specific class of antibodies. Biol Direct. 2013;8(1):1.
Chaudhuri R, Ansari FA, Raghunandanan MV, Ramachandran S. FungalRV: adhesin prediction and immunoinformatics portal for human fungal pathogens. BMC Genom. 2011;12(1):1.
Sirskyj D, Diaz-Mitoma F, Golshani A, Kumar A, Azizi A. Innovative bioinformatic approaches for developing peptide-based vaccines against hypervariable viruses. Immunol Cell Biol. 2011;89(1):81–9.
Petrovsky N, Brusic V. Computational immunology: the coming of age. Immunol Cell Biol. 2002;80(3):248–54.
Gupta S, Madhu MK, Sharma AK, Sharma VK. ProInflam: a webserver for the prediction of proinflammatory antigenicity of peptides and proteins. J Transl Med. 2016;14(1):178.
Vens C, Rosso M-N, Danchin EG. Identifying discriminative classification-based motifs in biological sequences. Bioinformatics. 2011;27(9):1231–8.
Rammensee H-G, Friede T, Stevanović S. MHC ligands and peptide motifs: first listing. Immunogenetics. 1995;41(4):178–228.
Nielsen M, Lund O, Buus S, Lundegaard C. MHC class II epitope predictive algorithms. Immunology. 2010;130(3):319–28.
Dhanda SK, Gupta S, Vir P, Raghava G: Prediction of IL4 inducing peptides. Clin Dev Immunol. 2013;2013.
Dhanda SK, Vir P, Raghava GP. Designing of interferon-gamma inducing MHC class-II binders. Biol Direct. 2013;8(1):1.
Wang HW, Pai TW. Machine learning-based methods for prediction of linear B-cell epitopes. Immunoinformatics. 2014:217–36.
Liu Y. Active learning with support vector machine applied to gene expression data for cancer classification. J Chem Info Comput Sci. 2004;44(6):1936–41.
Mohammad TAS, Nagarajaram HA. Svm-based method for protein structural class prediction using secondary structural content and structural information of amino acids. J Bioinform Comput Biol. 2011;9(04):489–502.
Zavaljevski N, Stevens FJ, Reifman J. Support vector machines with selective kernel scaling for protein classification and identification of key amino acid positions. Bioinformatics. 2002;18(5):689–96.
Saito T, Rehmsmeier M. The precision-recall plot is more informative than the ROC plot when evaluating binary classifiers on imbalanced datasets. PLoS ONE. 2015;10(3):e0118432.
Cavaillon J-M. Pro-versus anti-inflammatory cytokines: myth or reality. Cell Mol Biol. 2001;47(4):695–702.
Vidal C, Gomez-Hernandez A, Sanchez-Galan E, Gonzalez A, Ortega L, Gomez-Gerique JA, Tunon J, Egido J. Licofelone, a balanced inhibitor of cyclooxygenase and 5-lipoxygenase, reduces inflammation in a rabbit model of atherosclerosis. J Pharmacol Exp Ther. 2007;320(1):108–16.
Janeway CA, Travers P, Walport M, Shlomchik MJ. The immune system in health and disease. 2001
Veljaca M. Anti-inflammatory peptides and proteins in inflammatory bowel disease. Curr Opinion Investig Drug. 2001;2(10):1387–94.
Ialenti A, Santagada V, Caliendo G, Severino B, Fiorino F, Maffia P, Ianaro A, Morelli F, Di Micco B, Cartenì M. Synthesis of novel anti-inflammatory peptides derived from the amino-acid sequence of the bioactive protein SV-IV. Eur J Biochem. 2001;268(12):3399–406.
Fries JF, Wolfe F, Apple R, Erlich H, Bugawan T, Holmes T, Bruce B. HLA–DRB1 genotype associations in 793 white patients from a rheumatoid arthritis inception cohort: frequency, severity, and treatment bias. Arthritis Rheum. 2002;46(9):2320–9.
SG developed SVM-based models. AKS developed RF models. SG and AKS developed web server. MKM, VS, SG, AKS & VKS analysed the data and drafted the manuscript. SG and VKS conceived the work, participated in the design of the study. All authors read and approved the final manuscript.
This work was supported by the intramural funding received from IISER Bhopal. We thank MHRD, Govt of India, funded Centre for Research on Environment and Sustainable Technologies (CREST) at IISER Bhopal for its support. However, the views expressed in this manuscript are that of the authors alone and no approval of the same, explicit or implicit, by MHRD should be assumed.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used for this study are provided at http://metagenomics.iiserb.ac.in/antiinflam/data.php.
This work was supported by the intramural funding received from IISER Bhopal.
Sudheer Gupta and Ashok K. Sharma contributed equally to this work
Additional file 2: Figure S1 Comparison of machine learning algorithms using amino acid composition as feature input at fivefold cross validation.
Additional file 3: Figure S2. Comparison of machine learning algorithms using amino acid composition as feature input at ten-fold cross validation.
Additional file 4: Figure S3 Optimization of mtry for random forest model using amino acid composition as feature input.
Additional file 5: Figure S4. Optimization of mtry and ntree for random forest model using amino acid composition as feature input.
Additional file 6: Figure S5. Performance of model constructed using amino acid composition as feature input at optimized parameters (mtry=2, ntree=2k) on validation set.
Additional file 7: Figure S6. Optimization of mtry for random forest model using dipeptide composition as feature input.
Additional file 8: Figure S7. Optimization of mtry and ntree for random forest model using dipeptide composition as feature input.
Additional file 9: Figure S8. Performance of model constructed using dipeptide composition as feature input at optimized parameters (mtry=11, ntree=2k) on validation set.
Additional file 10: Figure S9. Optimization of mtry for random forest model using tripeptide composition as feature input.
Additional file 11: Figure S10. Optimization of ntree for random forest model using tripeptide composition as feature input.
Additional file 12: Figure S11. Performance of model constructed using tripeptide composition as feature input at optimized parameters (mtry=24, ntree=1.5k) on validation set.
About this article
Cite this article
Gupta, S., Sharma, A.K., Shastri, V. et al. Prediction of anti-inflammatory proteins/peptides: an insilico approach. J Transl Med 15, 7 (2017). https://doi.org/10.1186/s12967-016-1103-6
- Support Vector Machine
- Random Forest
- Area Under Curve
- Support Vector Machine Model