IL17eScan: A Tool for the Identification of Peptides Inducing IL-17 Response

IL-17 cytokines are pro-inflammatory cytokines and are crucial in host defense against various microbes. Induction of these cytokines by microbial antigens has been investigated in the case of ischemic brain injury, gingivitis, candidiasis, autoimmune myocarditis, etc. In this study, we have investigated the ability of amino acid sequence of antigens to induce IL-17 response using machine-learning approaches. A total of 338 IL-17-inducing and 984 IL-17 non-inducing peptides were retrieved from Immune Epitope Database. 80% of the data were randomly selected as training dataset and rest 20% as validation dataset. To predict the IL-17-inducing ability of peptides/protein antigens, different sequence-based machine-learning models were developed. The performance of support vector machine (SVM) and random forest (RF) was compared with different parameters to predict IL-17-inducing epitopes (IIEs). The dipeptide composition-based SVM-model displayed an accuracy of 82.4% with Matthews correlation coefficient = 0.62 at polynomial (t = 1) kernel on 10-fold cross-validation and outperformed RF. Amino acid residues Leu, Ser, Arg, Asn, and Phe and dipeptides LL, SL, LK, IL, LI, NL, LR, FK, SF, and LE are abundant in IIEs. The present tool helps in the identification of IIEs using machine-learning approaches. The induction of IL-17 plays an important role in several inflammatory diseases, and identification of such epitopes would be of great help to the immunologists. It is freely available at http://metagenomics.iiserb.ac.in/IL17eScan/ and http://metabiosys.iiserb.ac.in/IL17eScan/.

Background

Human body harbors complex microbial communities which may exist in planktonic forms or as higher order structures termed as biofilms (1). The interaction of the peripheral immune system with these microbes has an essential role in the pathophysiology of different diseases (2). One of the key components of the peripheral immune system is IL-17 family of cytokines, which play regulatory roles in host defense and during inflammatory diseases. They mediate pro-inflammatory responses via surface receptors on target cells and play several protective roles in host defense against pathogens at epithelial and mucosal barriers including skin, colon, and lung (3).

The induction of IL-17 by antigens present in gut commensal microbes and its relation with ischemic brain injury/stroke has been well established (2). The intestinal commensal microbes modulate the lymphocyte populations, which lead to various pathological conditions or dysbiosis. Similarly, in case of oral biofilms, the peptides Kgp467–477 of lysine-gingipain protein from Porphyromonas gingivalis induce IL-17 and further immunopathology in the case of periodontitis and gingivitis (4). On the other hand, the induction of IL-17 by peptide from agglutinin-like sequence protein in the case of oropharyngeal candidiasis makes it a suitable candidate for immunotherapeutics.

Similarly, there are reports of an increased level of gastric mucosal IL-17 level in response to Helicobacter pylori biofilm in mice (5, 6). The pneumococcal surface adhesin A231–268 (PsaA231–268), which is a highly conserved region in clinically relevant S. pneumonia strains, can induce an IL-17 response in mice upon infection (7). Furthermore, the Myelin basic protein 85–99 mimicking bacterial peptide can induce IL-17 in humanized transgenic mice (8). Likewise, myocarditogenic mimicry epitopes, such as BAC 25–40 peptide of Bacillus sp., induce IL-17 in autoimmune myocarditis in mouse model suggesting a role in its mediation (9). IL-17 secretion can also be triggered when CD4⁺ T-cells encounter viruses. For example, AA242–259 of rotaviral VP6 protein induces an IL-17 response in spleen cells from mice (10). Briefly, the induction of IL-17 in response to various antigens plays a pivotal role in initiation and/or development of several allergic inflammatory responses and autoimmune diseases such as multiple sclerosis (11), autoimmune encephalomyelitis (12), rheumatoid arthritis (13), systemic lupus erythematous (14), Behcet’s disease (15), and psoriasis (16). These evidences suggest that there is a peptide-sequence-specific induction of IL-17 through biofilms and planktonic microbial communities, which further leads to pro-inflammatory responses and pathogenesis. Further the role of selected residues in an epitope was demonstrated by a study carried out by mutating the key binding residues of epitopes and showed that the IL-17-producing CD8⁺ T cells were largely epitope specific (17). Similarly, five key residues essential for T cell activation were identified by replacing the residues with alanine amino acid in env_122–141 epitope of Friend murine leukemia virus (18).

Several studies have focused on the in silico prediction of different types of immune epitopes such as IL4-inducing peptides (19), IFN-gamma inducing major histocompatibility factor (MHC) binders (19), MHC binders (20), T cell epitopes (21, 22), B-cell epitopes (23, 24), and allergenicity (25, 26). However, there are no reports of any study in which the prediction of IL-17 induction by peptides was carried out. In this study, we have developed a classification method to predict the IL-17-inducing property of peptides using sequence-based features from experimentally validated IL-17-inducing and non-inducing epitopes.

Methods

Dataset

To ensure a clean and experimentally validated data, the epitope (peptide) sequences reported as IL-17 (IL-17 A or IL-17 F) inducing and non-inducing in different assays were downloaded from the Immune Epitope Database (IEDB) (27). The length of peptides in the epitope data was between 5 and 30 amino acids, and the longer peptides were not included in the study. A total of 338 IL-17-inducing unique epitopes (IIEs) were retrieved and labeled as positive data. The negative data comprised of 984 unique IL-17 non-inducing epitopes (INIEs) which do not elicit an IL-17 response. The peptides in the positive dataset which showed an exact match with the peptides present in the negative dataset were removed from the negative dataset (50 common peptides were removed from 1,034 peptides of negative data). Thus, the sequences of IIEs and INIEs were mutually exclusive with no overlapping peptides in the two groups. Of the total dataset, 80% of the sequences were randomly selected as the training dataset, and 20% data were kept as the validation dataset (Figure 1). The final training dataset contained 271 IIEs (positive data) and 786 INIEs (negative data), whereas the validation dataset consisted of 67 IIEs and 198 INIEs.

FIGURE 1

Figure 1. Flowchart showing steps involved in the development of prediction model and web server. The figure shows the steps involved in retrieval of data, training, validation, and construction of the IL17eScan tool.

To examine the positional amino acid conservation in terminal residues, five residues were cut from both the N′ and C′ terminals of the epitope sequences. The two sample logos (TSLs) were prepared with TSL software (http://www.twosamplelogo.org/) (28).

Input Features Model Development

Composition-Based Features

Amino Acid Composition (AAC)

Amino acid composition is the percentage of each amino acid in a peptide of given length. AAC has been widely used in binary classification problems in machine learning (29–31). Each peptide/protein can be represented by percentage composition of the 20 naturally occurring amino acids making a vector size of 20. AAC for each amino acid can be calculated as:

$\text{AAC}(\text{i})=\frac{\text{Total number of amino acid}(\text{i})}{\text{Total number of all possible amino acids}}\times 100,$

where AAC(i) is the AAC of the amino acid (i).

Dipeptide Composition (DPC)

Dipeptide composition is another widely used input feature for peptide/protein composition-based classification (23, 29, 31), which is calculated using the percentages of the 400 dipeptide combinations. Several immune epitope prediction algorithms have used the DPC-based classification (19, 23). Apart from the composition, DPC additionally provides information about the local arrangements of amino acids in a sequence. Percentage of all possible pairs of amino acids was calculated using the following equation:

$\text{DPC}(\text{i})=\frac{\text{Total number of dipeptides}(\text{i})}{\text{Total number of all possible dipeptides}}\times 100,$

where DPC(i) is the dipeptide frequency of dipeptide (i) and the dipeptide (i) is one out of 400 dipeptides.

Amino Acid Pair (AAP)

Amino acid pair can be defined as weighted DPC in which each pair carries a weight based on its propensity in the given dataset. The AAP-based feature has been used for the prediction of B-cell epitopes and IL4-inducing epitopes in the past by different authors (19, 23). The AAP feature was calculated as described in the earlier studies (19, 24, 32).

Machine Learning-Based Prediction Models

Support Vector Machine (SVM)

Support vector machine is a supervised machine-learning algorithm that can learn to classify positive and negative data by drawing an optimal hyperplane in high-dimensional feature space separating the two with the highest possible distance. This learning can be used for the classification of unlabeled data. It performs very well on biological data because of its ability to handle large feature spaces and avoid over-fitting, and thus, has been extensively implemented in several immune epitopes prediction tools (19, 33, 34), protein structure prediction (35) and genomic data (36). In this study, SVM^light package, available at http://svmlight.joachims.org/ was used for SVM-based predictive modeling. The linear, polynomial, and radial bias function (RBF) kernels were tested using various parameters.

Random Forest (RF)

Random forest is an ensemble-based classification and regression method in which a large number of independent decision trees are formed and are then combined to give the final decision. It was implemented in this study as it has a fast and robust algorithm. In this study, the randomForest package in R has been used for developing the classification model. Different mtry and ntrees were tested to build the models.

Performance Evaluation of Prediction Models

To evaluate and compare the machine-learning methods and prediction models, cross-validation technique was adopted. Cross-validation is a widely accepted method which involves division of the data into two segments. The first part is used to train the model and the other holdout or test data are used to test the model. A 10-fold cross-validation was carried out, where nine parts were used for training of the model, and the 10th one was used for testing the model. The process is iterated 10 times to test all the segments. Results obtained from all the 10 predictions are taken together for measuring the performance using threshold-dependent and threshold-independent parameters. The threshold-independent parameter, area under curve (AUC), was measured using PERF software. ACC, sensitivity (SEN), specificity (SPC), and Matthews correlation coefficient (MCC) were threshold-dependent parameters and were calculated as per the following equations:

$\text{ACC}=\frac{\text{TP}+\text{TN}}{\text{TP}+\text{FN}+\text{FP}+\text{TN}},$

$\text{SEN}=\frac{\text{TP}}{\text{TP}+\text{FN}},$

$\text{SPC}=\frac{\text{TN}}{\text{TN}+\text{FP}},$

$\text{MCC}=\frac{(\text{TP}\times \text{TN})-(\text{FP}\times \text{FN})}{\sqrt{(\text{TP}+\text{FP})(\text{TP}+\text{FN})(\text{TN}+\text{FP})(\text{TN}+\text{FN})}},$

where TP = True Positive, FP = False Positive, FN = False Negative, TN = True Negative.

Prediction of IL-17-Inducing Peptides in Microbes

To compare the distribution of IL-17-inducing epitopes (IIEs) in different microbes known to induce Th17 responses, or known to induce interleukins other than IL-17 and non-inducing saprophytic microbes (37, 38), the protein sequences of Segmented Filamentous Bacteria, Staphylococcus aureus, Candida albicans, Listeria monocytogenes, Mycobacterium tuberculosis, Acetobacter aceti and Propionibacterium acnes were retrieved from NCBI. Random synthetic peptides were generated in 10 different sets with 1,000 peptides (15-mers) in each set using in-house Perl scripts and were predicted for their IL-17-inducing property. The IIEs were predicted using the IL17eScan web server.

Results

Composition and Position-Based Conservation Analysis

The AAC analysis revealed Leu, Ser, Arg, and Asn as the most abundant amino acids in IIEs as compared with INIEs. Similarly, Ala, Asp, Gly, and Pro were found to be rich in INIEs (Figure 2; Data Sheet S1 in Supplementary Material). Furthermore, some dipeptides were found to be significantly abundant (Welch’s t-test, p < 0.01) in IIEs. The top 10 most abundant and significant dipeptides present in IIEs were LL, SL, LK, IL, LI, NL, LR, FK, SF, and LE, whereas top 10 most abundant and significant dipeptides present in INIEs were PG, GA, AA, GP, PA, PP, AG, GD, PE, and AP (Figure 3; Data Sheet S2 in Supplementary Material).

FIGURE 2

Figure 2. Amino acid composition analysis of IL-17-inducing epitopes (IIEs) and IL-17 non-inducing epitopes (INIEs). The proportions of amino acid in IIEs and INIEs along with the p-values (Welch’s t-test) are shown.

FIGURE 3

Figure 3. Distribution of dipeptide composition of IL-17-inducing epitopes (IIEs) and IL-17 non-inducing epitopes (INIEs). Significant dipeptides in IIEs and INIEs with p-value < 0.01 (Welch’s t test) are shown.

To explore the positional conservation of the amino acid residues, the first five residues from N′- and C′-terminal of epitopes were examined. The TSL analysis revealed the conservation and abundance of Leu residues at various positions (particularly at the N′-terminal), which was also observed as abundant in the compositional analysis of the positive dataset (Figure 4).

FIGURE 4

Figure 4. Two sample logo showing positional conservation of five residues at both the terminals (N′- and C′-) in IL-17-inducing epitopes and IL-17 non-inducing epitopes. At both N and C terminals, amino acid L was found abundant and conserved at different positions.

Human Leukocyte Antigen (HLA)-Allele Distribution Analysis

Antigenic epitopes are identified by HLA molecules in the host, and the presence of different HLA types is a key determinant of epitope’s action in IL-17 induction (39). The HLA-allele distribution analysis was carried out to examine the association of any specific allele with IIEs. The analysis revealed the association of HLA alleles such as HLADRB1*15:01, H2 s class II, HLAA*02:01, and HLADR with IL-17 induction. Similarly, HLA alleles, such as H2 Iab, H2 b class II, and H2 Iaq, showed association with INIEs (Figure 5; Data Sheet S3 in Supplementary Material). Some previous studies also suggested the association of HLADR alleles with the induction of IL-17, and thus, leading to autoimmune disease such as Rheumatoid arthritis (40).

FIGURE 5

Figure 5. Distribution of human leukocyte antigen (HLA) alleles among assays reporting IL-17-inducing epitopes (IIE) and IL-17 non-inducing epitopes (INIEs). The HLA-allele distribution analysis examines the association of any specific allele with IIEs. The analysis revealed IL-17 induction associated with HLA alleles such as HLADRB1*15:01, H2 s class II, HLAA*02:01, and HLADR. Similarly, HLA alleles, such as H2 Iab, H2 b class II, and H2 Iaq, showed association with INIEs.

Machine Learning-Based Classification

The compositional profiles of IIEs and INIEs were found to be different, and thus, could be exploited to classify the epitopes using machine learning-based algorithms. SVM- and RF-based models were developed and evaluated using 10-fold cross-validation. The performance of SVM- and RF-based models on different sequence-based features at various kernels and mtry, respectively are discussed (Tables 1 and 2; Figure 6). Since SVM emerged as the best classification method for IIE and INIE prediction, results of SVM-based models have been mentioned and discussed in the manuscript.

TABLE 1

Table 1. Performance of support vector machine-based models on different sequence-based features using various kernels.

TABLE 2

Table 2. Performance of random forest-based models on different sequence-based features using various mtry.

FIGURE 6

Figure 6. ROC plots of prediction models developed using SVM^light as machine-learning technique. The DPC_polynomial model (shown in purple) achieved highest area under curve (AUC = 0.91 shown in Table 1).

AAC-Based Models

Support vector machine-based classification using AAC showed the best performance with RBF kernel (t = 2), gamma parameter (g) = 0.005, trade-off factor (c) = 1 and a cost factor (j) of 5. This model performed with an accuracy (ACC) of 77.29% and MCC of 0.47 (Table 1). However, MCC at linear and polynomial kernel was found to be 0.35 and 0.43, respectively, which was lesser than the RBF kernel (Table 1; Figure 6).

Dipeptide-Based Models

Dipeptide composition was also used as input feature since it harbors more information because of the longer vector length (400). DPC-based models with polynomial kernel (t = 1) performed best with parameter d = 2. Unlike the AAC-based model which performed best at complex kernel (RBF), the DPC-based model could classify the IIPs from INIEs better with the simpler polynomial kernel. The ACC, MCC, and AUC of the model were found to be 82.4%, 0.62, and 0.91, respectively. Similarly, the models with linear and RBF kernel could only achieve MCC of 0.39 and 0.57, respectively (Tables 1 and 2). The best AUC value of 0.91 was obtained for DPC at polynomial kernel (t = 1) (Figure 6).

AAP-Based Models

To further improve the performance, weights were given to each dipeptide, and the AAP values were calculated from the DPC as discussed in the Methods section. The model constructed using RBF kernel (t = 2) showed the best performance with an ACC of 84.58 and MCC of 0.6. The optimized parameters included gamma parameter (g) = 0.01, trade-off factor (c) = 5 and a cost factor (j) = 1 for this model (Tables 1 and 2; Figure 6).

Performance on Validation Dataset

After the 10-fold cross-validation, the performance of different SVM- and RF-based models was evaluated on a validation dataset to ensure that there was no over-fitting and the achieved performance of the final model is not due to over-optimization. The performance on the validation dataset are summarized in Table 3 for SVM-based models and Table 4 for RF based models. As mentioned earlier, the best performing models for AAC-, DPC-, and AAP-based features achieved MCC of 0.47 (t = 2), 0.62 (t = 1), and 0.60 (t = 2), respectively, on training data. On the validation dataset, the same models displayed the MCC values of 0.5, 0.57, and 0.52 for AAC, DPC, and AAP, respectively.

TABLE 3

Table 3. Performance of different support vector machine-based models on validation dataset.

TABLE 4

Table 4. Performance of different random forest-based models on validation dataset.

IIEs in Biofilm-Forming Bacteria

To examine the epitopes which may modulate host immune system by inducing IL-17 in biofilm-forming microbes in various disease conditions (41), we extracted all the protein sequences of these microorganisms from SwissProt database and analyzed using the prediction pipeline. We identified several IIEs (15-mers) in different proteins belonging to different microorganisms. The top 10 proteins for every microbial species harboring the highest number of epitopes are mentioned in the Data Sheet S4 in Supplementary Material. Among the major predicted IL-17 inducers, “Probable sugar efflux transporter protein” is commonly found in Haemophilus influenzae as well as Klebsiella pneumonia. Similarly, “Na(+)/H(+) antiporter NhaB protein” from Proteus mirabilis and Pseudomonas aeruginosa were found to have a large number of IL-17-inducing peptides. DNA polymerase III subunit of P. mirabilis involved in urinary catheter cystitis was found to harbor IIE, which corroborates with a previous study on IL-17 induction by DNA polymerase of Human adenovirus 5 (42).

Prediction of IL-17-Inducing Peptides in Microbes

The IIEs were predicted in microbes known to induce IL-17 response, known to induce other interleukins and in saprophytes using IL17eScan web server. The IIEs were found enriched in the microbes known to induce Th17 responses (Data Sheet S4 in Supplementary Material). L. monocytogenes and M. tuberculosis, which promote Th1 responses showed a lower representation of IIEs in their proteins (37, 38). A similar lower representation of IIEs was also observed in the case of saprophytic microbes such as A. aceti and P. acnes (Data Sheet S4 in Supplementary Material). On increasing the threshold to 1, a notable reduction in the percentages of IL-17-inducing proteins was observed, where the percentage was highest (1%) in the case of IL-17-inducing bacteria and the lowest (0.1%) for the bacteria for which there are no reports of their role in IL-17 induction. To further validate the above predictions, random peptides were generated in 10 different sets with 1,000 peptides (15-mers) in each set and were predicted for their IL-17-inducing property at the threshold of 1. Interestingly, none of the synthetic peptides in any of the 10 datasets were predicted to be IL-17 inducing. These results attest the usability of IL17eScan to predict the IIEs in the real datasets.

Web Server and Tools

A web server “IL17eScan” is constructed to provide the tools for the prediction, virtual screening, and mapping of IIEs. These available modules for prediction incorporate the best performing algorithm (DPC-based model) as default, which runs the queries through a pipeline and classifies the query peptides into IIEs or INIEs. A peptide with a score higher than the threshold is predicted as IL-17 inducing. An increase in the threshold will increase the SPC, and the prediction will become more stringent. As a trade-off between SPC and SEN, an optimal threshold (0.5) is set as default on the web server. However, the user has the flexibility to increase or decrease this threshold and analyze the results as per the requirement. Also, the AAC-based model is provided in all the modules for handling large queries since AAC-based models are faster than DPC-based models due to smaller vector size (20).

PepPred

The module “PepPred” classifies one or more proteins/peptide sequence(s) of length ranging from 5 to 30 amino acids into IIEs or INIEs. The stringency of positive prediction can be set using a threshold value provided by the user. Also, the “virtual screening and designing” option has also been provided, which allows the user to select peptides based on their prediction score, modify the query peptides, and resubmit them for prediction. This option carries out substitution of each amino acid of the peptide with other amino acids. After the substitution, for the resubmitted peptides, the results are displayed in the same tabular format with prediction scores. It allows the users to predict the IL-17-inducing nature of the multiple variants of the query peptide, and thus, is useful in assessing the position-specific effects of each amino acid in modulating the IL-17-inducing activity of the peptide.

PepScan

In contrast to the “PepPred” module that deals with smaller peptides, the “PepScan” module predicts the antigenic regions in full-length proteins that can potentially induce an IL-17 response in a host. Users are allowed to provide a window length ranging from 5 to 30 peptides which determine the length of peptide sequences considered for prediction. Virtual screening and design option is also available for this module.

MetaGScan

To investigate IIEs in amino acid sequence data obtained from metagenomic studies, we have incorporated a separate module “MetaGScan.” This module requires raw translated reads (peptide orfs) from any metagenomic study and identifies the antigenic regions which may induce an IL-17 response. The peptide orf containing the positively predicted epitopes can be aligned for similarity search against the protein sequences present in SwissProt database using BLASTP. As an example, we have included metagenomic reads data from the gut of a diabetes type II patient (processed reads with annotation from https://www.ebi.ac.uk/metagenomics/projects/SRP008047/samples/SRS259434/runs/SRR341581/results/versions/1.0) in this module.

EpiScan

To examine the exact occurrence of IIEs on the protein of interest, the EpiScan tool is provided which allows the user to map experimentally validated IIEs from IEDB (27) on the query peptide or proteins. The results are also linked to the related assays available in IEDB.

SimSearch

Unlike EpiScan, which searches for exact matches, the “SimSearch” option maps the experimentally validated epitopes to their similar sequences in the query peptide/protein. This module implements Smith–Waterman search algorithm and displays the match along with the links to related assays in IEDB.

Discussion

Recent advances in metagenomic and high-throughput assay technologies have provided us with new insights into the diversity of human microbiome, and their interaction with host immune system in different inflammatory and autoimmune diseases. Among these interactions, induction of IL-17 is one of the most studied pro-inflammatory responses against pathogens (3, 43, 44). In this study, we have developed an in silico method to predict the IL-17-inducing ability of peptides/proteins based on the sequence-based features derived from a set of experimentally validated IIEs (positive set), and non-inducing epitopes (negative set) obtained from the IEDB. Although the IL-17 response can be defined as induction of any cytokine of IL-17 class, the epitope assay data in IEDB were limited only to IL-17 A and IL-17 F cytokines of IL-17 class. Thus, the present tool is aimed only at predicting the IIEs, which is one of the limitations of the tool. Further, the IIEs had lengths ranging from 5 to 30 amino acids except for a few longer epitopes, and thus, the length range of 5–30 amino acids was selected for training and prediction. The non-redundant dataset constructed from the IL-17-inducing and non-inducing peptides ensured no over-fitting or bias due to the presence of multiple instances of the same peptide. The IIEs belonged to 117 unique proteins from 54 different taxa, which further reduced the chances of any bias.

The compositional analysis and positional conservation of residues by TSL revealed that Leu is highly abundant in IIEs as compared with INIEs. The Leu-rich epitopes have also been shown to induce an IL-17 response in different autoimmune diseases such as NLRP3 (autoimmune encephalomyelitis) (45), FLRT2 (systemic lupus erythematosus) (46, 47), and LGI1 (limbic encephalitis) (48–50). A higher abundance of specific residues has been previously observed for epitopes inducing other interleukins and immune cells (17, 18, 21, 22, 51). These findings suggest that a few residues could be associated with IL-17 induction. However, determining the biological significance of these residues in IL-17 induction requires further studies and experimental validations.

The development of IL-17 prediction models was carried out after evaluating multiple machine-learning methods, and the best performing DPC-based SVM classification models with polynomial kernel was incorporated in the web server pipeline for the best results. The DPC-based model performed better than the AAC-based model perhaps due to the larger vector size. However, as a weighted DPC, AAP feature was not able to improve the performance. Given the large vector size and high performance, the models were also scrutinized for over-optimization by testing on a validation dataset. The validation of models on the validation dataset confirmed that the high performance of the models is not due to over-fitting.

Further, the performance of the tool on IL-17-inducing, non-inducing, and saprophytic microbes and on a random peptide set underscores its applicability on real biological datasets and reveals the differences in the percentage of such epitopes in IL-17-inducing and non-inducing organisms. The tool also provides a reliable and reproducible framework for epitope prediction in peptides or proteins from whole genomes and metagenomes. For any prediction-based method, setting an optimal threshold for the selection of hits is one of the limitations, where a lower threshold could result in a higher number of false positives, although it may improve the SEN and vice versa for a higher threshold. Thus, we have provided a default threshold to ensure optimal performance; however, the stringency of results should be adjusted by selecting an appropriate threshold by the user.

The availability of experimentally validated IIEs for all classes of IL-17 cytokines will help in further improving the applicability of the tool. The present tool will help in developing a better understanding of the IL-17-inducing property of the peptides and is anticipated to be widely used for the computational identification of IIEs from genomes and metagenomes.

Conclusion

The propensity of antigens to induce an IL-17 response is of significant importance in the initiation and development of several allergic inflammatory responses and autoimmune diseases. Therefore, the developed machine learning-based tool provides a useful resource for predicting the IL-17-inducing peptides by successfully utilizing the sequence-based signatures of experimentally validated IIEs. To the best of our knowledge, this is the only in silico based method available to predict the IIEs in genomic and metagenomic peptides/proteins, and the lead peptides may serve as potential candidates for immunotherapeutics. The IL17eScan is available freely as a web server for academic use.

Author Contributions

SG developed SVM-based models. PM developed RF-based models. SG, PM, and MM developed web server. SG, PM, MM, and VS analyzed the data and drafted manuscript. SG and VS conceived the work and participated in the design of the study. All the authors read and approved the final manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This work was supported by the intramural funding received from IISER Bhopal.

Supplementary Material

The Supplementary Material for this article can be found online at http://www.frontiersin.org/article/10.3389/fimmu.2017.01430/full#supplementary-material.

Abbreviations

AAC, amino acid composition; ACC, accuracy; AUC, area under curve; DPC, dipeptide composition; FN, false negative; FP, false positive; IEDB, Immune Epitope Database; IIE, IL-17-inducing epitope; INIE, IL-17 non-inducing epitope; HLA, human leukocyte antigen; MCC, Matthews correlation coefficient; MHC, major histocompatibility factor; RBF, radial bias function; RF, random forest; SEN, sensitivity; SPC, specificity; SVM, support vector machine; TN, true negative; TP, true positive; TSL, two sample logo.

References

1. Bjarnsholt T. The role of bacterial biofilms in chronic infections. APMIS Suppl (2013) 121:1–51. doi:10.1111/apm.12099

CrossRef Full Text | Google Scholar

2. Benakis C, Brea D, Caballero S, Faraco G, Moore J, Murphy M, et al. Commensal microbiota affects ischemic stroke outcome by regulating intestinal gammadelta T cells. Nat Med (2016) 22:516–23. doi:10.1038/nm.4068

CrossRef Full Text | Google Scholar

3. Jin W, Dong C. IL-17 cytokines in immunity and inflammation. Emerg Microbes Infect (2013) 2:e60. doi:10.1038/emi.2013.58

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Bittner-Eddy PD, Fischer LA, Costalonga M. Identification of gingipain-specific I-A(b) -restricted CD4+ T cells following mucosal colonization with Porphyromonas gingivalis in C57BL/6 mice. Mol Oral Microbiol (2013) 28:452–66. doi:10.1111/omi.12038

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Luzza F, Parrello T, Monteleone G, Sebkova L, Romano M, Zarrilli R, et al. Up-regulation of IL-17 is associated with bioactive IL-8 expression in Helicobacter pylori-infected human gastric mucosa. J Immunol (2000) 165:5332–7. doi:10.4049/jimmunol.165.9.5332

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Shiomi S, Toriie A, Imamura S, Konishi H, Mitsufuji S, Iwakura Y, et al. IL-17 is involved in Helicobacter pylori-induced gastric inflammatory responses in a mouse model. Helicobacter (2008) 13:518–24. doi:10.1111/j.1523-5378.2008.00629.x

CrossRef Full Text | Google Scholar

7. Singh R, Gupta P, Sharma PK, Ades EW, Hollingshead SK, Singh S, et al. Prediction and characterization of helper T-cell epitopes from pneumococcal surface adhesin A. Immunology (2014) 141:514–30. doi:10.1111/imm.12194

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Greene MT, Ercolini AM, DeGutes M, Miller SD. Differential induction of experimental autoimmune encephalomyelitis by myelin basic protein molecular mimics in mice humanized for HLA-DR2 and an MBP(85-99)-specific T cell receptor. J Autoimmun (2008) 31:399–407. doi:10.1016/j.jaut.2008.09.004

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Massilamany C, Gangaplara A, Steffen D, Reddy J. Identification of novel mimicry epitopes for cardiac myosin heavy chain-alpha that induce autoimmune myocarditis in A/J mice. Cell Immunol (2011) 271:438–49. doi:10.1016/j.cellimm.2011.08.013

CrossRef Full Text | Google Scholar

10. McNeal MM, Basu M, Bean JA, Clements JD, Choi AH, Ward RL. Identification of an immunodominant CD4+ T cell epitope in the VP6 protein of rotavirus following intranasal immunization of BALB/c mice. Virology (2007) 363:410–8. doi:10.1016/j.virol.2007.01.041

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Langrish CL, Chen Y, Blumenschein WM, Mattson J, Basham B, Sedgwick JD, et al. IL-23 drives a pathogenic T cell population that induces autoimmune inflammation. J Exp Med (2005) 201:233–40. doi:10.1084/jem.20041257

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Park H, Li Z, Yang XO, Chang SH, Nurieva R, Wang YH, et al. A distinct lineage of CD4 T cells regulates tissue inflammation by producing interleukin 17. Nat Immunol (2005) 6:1133–41. doi:10.1038/ni1261

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Gaffen SL. The role of interleukin-17 in the pathogenesis of rheumatoid arthritis. Curr Rheumatol Rep (2009) 11:365–70. doi:10.1007/s11926-009-0052-y

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Wong CK, Ho CY, Li EK, Lam CW. Elevation of proinflammatory cytokine (IL-18, IL-17, IL-12) and Th2 cytokine (IL-4) concentrations in patients with systemic lupus erythematosus. Lupus (2000) 9:589–93. doi:10.1191/096120300678828703

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Hamzaoui K, Hamzaoui A, Guemira F, Bessioud M, Hamza M, Ayed K. Cytokine profile in Behcet’s disease patients. Relationship with disease activity. Scand J Rheumatol (2002) 31:205–10. doi:10.1080/030097402320318387

CrossRef Full Text | Google Scholar

16. Hueber W, Patel DD, Dryja T, Wright AM, Koroleva I, Bruin G, et al. Effects of AIN457, a fully human antibody to interleukin-17A, on psoriasis, rheumatoid arthritis, and uveitis. Sci Transl Med (2010) 2:52ra72. doi:10.1126/scitranslmed.3001107

CrossRef Full Text | Google Scholar

17. Myoung J, Kang HS, Hou W, Meng L, Dal Canto MC, Kim BS. Epitope-specific CD8+ T cells play a differential pathogenic role in the development of a viral disease model for multiple sclerosis. J Virol (2012) 86:13717–28. doi:10.1128/JVI.01733-12

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Shimizu T, Uenishi H, Teramura Y, Iwashiro M, Kuribayashi K, Tamamura H, et al. Fine structure of a virus-encoded helper T-cell epitope expressed on FBL-3 tumor cells. J Virol (1994) 68:7704–8.

PubMed Abstract | Google Scholar

19. Dhanda SK, Vir P, Raghava GP. Designing of interferon-gamma inducing MHC class-II binders. Biol Direct (2013) 8:30. doi:10.1186/1745-6150-8-30

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Bhasin M, Raghava GP. A hybrid approach for predicting promiscuous MHC class I restricted T cell epitopes. J Biosci (2007) 32:31–42. doi:10.1007/s12038-007-0004-5

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bhasin M, Raghava GP. Analysis and prediction of affinity of TAP binding peptides using cascade SVM. Protein Sci (2004) 13:596–607. doi:10.1110/ps.03373104

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Bhasin M, Raghava G. Prediction of CTL epitopes using QM, SVM and ANN techniques. Vaccine (2004) 22:3195–204. doi:10.1016/j.vaccine.2004.02.005

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Gupta S, Ansari HR, Gautam A; Open Source Drug Discovery Consortium, Raghava GP. Identification of B-cell epitopes in an antigen for inducing specific class of antibodies. Biol Direct (2013) 8:27. doi:10.1186/1745-6150-8-27

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Singh H, Ansari HR, Raghava GP. Improved method for linear B-cell epitope prediction using antigen’s primary sequence. PLoS One (2013) 8:e62216. doi:10.1371/journal.pone.0062216

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Saha S, Raghava GP. AlgPred: prediction of allergenic proteins and mapping of IgE epitopes. Nucleic Acids Res (2006) 34:W202–9. doi:10.1093/nar/gkl343

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Dimitrov I, Flower DR, Doytchinova I. AllerTOP – a server for in silico prediction of allergens. BMC Bioinformatics (2013) 14(Suppl 6):S4. doi:10.1186/1471-2105-14-S6-S4

CrossRef Full Text | Google Scholar

27. Vita R, Overton JA, Greenbaum JA, Ponomarenko J, Clark JD, Cantrell JR, et al. The immune epitope database (IEDB) 3.0. Nucleic Acids Res (2015) 43:D405–12. doi:10.1093/nar/gku938

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Vacic V, Iakoucheva LM, Radivojac P. Two sample logo: a graphical representation of the differences between two sets of sequence alignments. Bioinformatics (2006) 22:1536–7. doi:10.1093/bioinformatics/btl151

PubMed Abstract | CrossRef Full Text | Google Scholar

29. 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:e93907. doi:10.1371/journal.pone.0093907

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Gupta S, Kapoor P, Chaudhary K, Gautam A, Kumar R, Raghava GP. Peptide toxicity prediction. Methods Mol Biol (2015) 1268:143–57. doi:10.1007/978-1-4939-2285-7_7

CrossRef Full Text | Google Scholar

31. 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–6. doi:10.1016/j.ygeno.2015.04.001

PubMed Abstract | CrossRef Full Text | Google Scholar

32. El-Manzalawy Y, Dobbs D, Honavar V. Predicting flexible length linear B-cell epitopes. Comput Syst Bioinformatics Conf (2008) 7:121–32. doi:10.1142/9781848162648_0011

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Saha S, Raghava GP. Prediction methods for B-cell epitopes. Methods Mol Biol (2007) 409:387–94. doi:10.1007/978-1-60327-118-9_29

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Desai DV, Kulkarni-Kale U. T-cell epitope prediction methods: an overview. Methods Mol Biol (2014) 1184:333–64. doi:10.1007/978-1-4939-1115-8_19

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Hu H-J, Pan Y, Harrison R, Tai PC. Improved protein secondary structure prediction using support vector machine with a new encoding scheme and an advanced tertiary classifier. IEEE Trans Nanobioscience (2004) 3:265–71. doi:10.1109/TNB.2004.837906

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Fletez-Brant C, Lee D, McCallion AS, Beer MA. kmer-SVM: a web server for identifying predictive regulatory sequence features in genomic data sets. Nucleic Acids Res (2013) 41:W544–56. doi:10.1093/nar/gkt519

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Yang Y, Torchinsky MB, Gobert M, Xiong H, Xu M, Linehan JL, et al. Focused specificity of intestinal TH17 cells towards commensal bacterial antigens. Nature (2014) 510:152–6. doi:10.1038/nature13279

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Sallusto F. Heterogeneity of human CD4(+) T cells against microbes. Annu Rev Immunol (2016) 34:317–34. doi:10.1146/annurev-immunol-032414-112056

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Mangalam AK, Taneja V, David CS. HLA class II molecules influence susceptibility versus protection in inflammatory diseases by determining the cytokine profile. J Immunol (2013) 190:513–8. doi:10.4049/jimmunol.1201891

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Choy E. Understanding the dynamics: pathways involved in the pathogenesis of rheumatoid arthritis. Rheumatology (Oxford) (2012) 51(Suppl 5):v3–11. doi:10.1093/rheumatology/kes113

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Aparna MS, Yadav S. Biofilms: microbes and disease. Braz J Infect Dis (2008) 12:526–30. doi:10.1590/S1413-86702008000600016

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Haveman LM, Bierings M, Legger E, Klein MR, de Jager W, Otten HG, et al. Novel pan-DR-binding T cell epitopes of adenovirus induce pro-inflammatory cytokines and chemokines in healthy donors. Int Immunol (2006) 18:1521–9. doi:10.1093/intimm/dxl085

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Yang XO, Chang SH, Park H, Nurieva R, Shah B, Acero L, et al. Regulation of inflammatory responses by IL-17F. J Exp Med (2008) 205:1063–75. doi:10.1084/jem.20071978

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Ishigame H, Kakuta S, Nagai T, Kadoki M, Nambu A, Komiyama Y, et al. Differential roles of interleukin-17A and -17F in host defense against mucoepithelial bacterial infection and allergic responses. Immunity (2009) 30:108–19. doi:10.1016/j.immuni.2008.11.009

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Gris D, Ye Z, Iocca HA, Wen H, Craven RR, Gris P, et al. NLRP3 plays a critical role in the development of experimental autoimmune encephalomyelitis by mediating Th1 and Th17 responses. J Immunol (2010) 185:974–81. doi:10.4049/jimmunol.0904145

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Nalbandian A, Crispin JC, Tsokos GC. Interleukin-17 and systemic lupus erythematosus: current concepts. Clin Exp Immunol (2009) 157:209–15. doi:10.1111/j.1365-2249.2009.03944.x

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Shirai T, Fujii H, Ono M, Nakamura K, Watanabe R, Tajima Y, et al. A novel autoantibody against fibronectin leucine-rich transmembrane protein 2 expressed on the endothelial cell surface identified by retroviral vector system in systemic lupus erythematosus. Arthritis Res Ther (2012) 14:R157. doi:10.1186/ar3897

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Komiyama Y, Nakae S, Matsuki T, Nambu A, Ishigame H, Kakuta S, et al. IL-17 plays an important role in the development of experimental autoimmune encephalomyelitis. J Immunol (2006) 177:566–73. doi:10.4049/jimmunol.177.1.566

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Lee JJ, Lee ST, Jung KH, Chu K, Lee SK. Anti-LGI1 limbic encephalitis presented with atypical manifestations. Exp Neurobiol (2013) 22:337–40. doi:10.5607/en.2013.22.4.337

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Zandi MS. Defining and treating leucine-rich glioma inactivated 1 antibody associated autoimmunity. Brain (2013) 136:2933–5. doi:10.1093/brain/awt256

CrossRef Full Text | Google Scholar

51. Nagpal G, Usmani SS, Dhanda SK, Kaur H, Singh S, Sharma M, et al. Computer-aided designing of immunosuppressive peptides based on IL-10 inducing potential. Sci Rep (2017) 7:42851. doi:10.1038/srep42851

PubMed Abstract | CrossRef Full Text | Google Scholar