Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Improved Method for Linear B-Cell Epitope Prediction Using Antigen’s Primary Sequence

  • Harinder Singh ,

    Contributed equally to this work with: Harinder Singh, Hifzur Rahman Ansari

    Affiliation Bioinformatics Center, Institute of Microbial Technology, Chandigarh, India

  • Hifzur Rahman Ansari ,

    Contributed equally to this work with: Harinder Singh, Hifzur Rahman Ansari

    Affiliation Bioinformatics Center, Institute of Microbial Technology, Chandigarh, India

  • Gajendra P. S. Raghava

    raghava@imtech.res.in

    Affiliation Bioinformatics Center, Institute of Microbial Technology, Chandigarh, India

Abstract

One of the major challenges in designing a peptide-based vaccine is the identification of antigenic regions in an antigen that can stimulate B-cell’s response, also called B-cell epitopes. In the past, several methods have been developed for the prediction of conformational and linear (or continuous) B-cell epitopes. However, the existing methods for predicting linear B-cell epitopes are far from perfection. In this study, an attempt has been made to develop an improved method for predicting linear B-cell epitopes. We have retrieved experimentally validated B-cell epitopes as well as non B-cell epitopes from Immune Epitope Database and derived two types of datasets called Lbtope_Variable and Lbtope_Fixed length datasets. The Lbtope_Variable dataset contains 14876 B-cell epitope and 23321 non-epitopes of variable length where as Lbtope_Fixed length dataset contains 12063 B-cell epitopes and 20589 non-epitopes of fixed length. We also evaluated the performance of models on above datasets after removing highly identical peptides from the datasets. In addition, we have derived third dataset Lbtope_Confirm having 1042 epitopes and 1795 non-epitopes where each epitope or non-epitope has been experimentally validated in at least two studies. A number of models have been developed to discriminate epitopes and non-epitopes using different machine-learning techniques like Support Vector Machine, and K-Nearest Neighbor. We achieved accuracy from ∼54% to 86% using diverse s features like binary profile, dipeptide composition, AAP (amino acid pair) profile. In this study, for the first time experimentally validated non B-cell epitopes have been used for developing method for predicting linear B-cell epitopes. In previous studies, random peptides have been used as non B-cell epitopes. In order to provide service to scientific community, a web server LBtope has been developed for predicting and designing B-cell epitopes (http://crdd.osdd.net/raghava/lbtope/).

Introduction

Identification of smallest regions in an antigen also called an antigenic region that can activate immune system is one of the major challenges in designing of a subunit or peptide-based vaccine. These antigenic regions, which stimulate B-cell response, are known as B-cell epitopes. Prediction of B-cell epitope is difficult but important for designing a peptide-based vaccine [1]. B-cell epitopes can be divided in two categories (i) continuous and (ii) discontinuous. The continuous or linear epitopes are made up of consecutive amino acids whereas the discontinuous or conformational epitopes constitute the spatially folded amino acids, which lie far away in the primary sequence. Linear B-cell epitope has vast application in the area of antibody production, immunodiagnostics; epitope-based vaccine design, selective de-immunization of therapeutic proteins and in autoimmunity [2], [3], [4]. Experimental methods for identification of B-cell epitopes are costly and time consuming.

In order to overcome limitations of experimental techniques, in the past several algorithms have been developed to predict linear B-cell epitopes [5], [6]. Due to variability in epitope length from 3–85 amino acids, the prediction of B-cell epitopes is much more complex than the prediction of T-cell epitopes. Recently, Kringelum et al, have analyzed conformational B-cell epitopes from antigen-antibody complexes and reported average length of conformation epitope as 15 residues [7]. Though the average length is 15, it does not mean that B-cell epitope core has 15 residues. To the best of author’s knowledge it is not known that what, is the minimum length of a conformational or continuous B-cell epitope. All methods developed so far have two major limitations, (i) they are based on a very limited number of epitopes (∼1000 epitopes and non-epitopes) and second (ii) these methods use random peptides as non B-cell epitope [8], [9], [10], [11], [12] (Table S1).

In this study, for the first time, we have exploited the availability of several thousands of experimentally verified epitopes and non-epitopes. We have derived five datasets from Immune Epitope Database (IEDB) called Lbtope_Fixed, Lbtope_Fixed_non_redundant, Lbtope_Variable, Lbtope_Variable_non_redundant and Lbtope_Confirm dataset. We developed various models on these datasets for discriminating B-cell epitopes from non-epitopes. A web server has been developed for predicting B-cell epitopes using best models developed on these datasets.

Materials and Methods

We have obtained experimentally validated 49694 B-cell epitopes and 50324 non B-cell epitopes from Immune Epitope Database (IEDB) in Jan 2012 [13]. These epitopes have 2 to 85 amino acids and belong to 3689 antigen sequences. We created five different datasets namely Lbtope_Fixed, Lbtope_Fixed_non_redundant, Lbtope_Variable, Lbtope_Variable_non_redundant, Lbtope_Variable and Lbtope_Confirm from this main dataset. The description of each dataset is as follows:

Lbtope_Fixed Dataset.

Most of the machine learning techniques commonly used for developing prediction or class discrimination need definite length patterns. Since B-cell epitopes have variable length, we used truncation and extension technique used in previous studies to generate definite length peptides (epitopes & non-epitopes) of 20 residues [8][11]. Following procedure has been adopted to generate fixed length epitopes; (i) all epitopes having less than five residues were removed, (ii) epitopes having more than 20 residues were trimmed from both ends to generate epitope of 20 residues from middle, (iii) epitopes containing less than 20 residues have been extended to 20 by adding an equal number of residues at both ends of the epitope, and (iv) finally, identical epitopes were removed. In order to extend an epitope, we mapped it on source antigen from where it has been derived and then we extended its length. In summary, Lbtope_Fixed dataset contains unique 19803 positive patterns or B-cell epitopes and 28329 negative patterns or non B-cell epitopes, where each pattern contains 20-residues. We also removed patterns common in both types of patterns. Our final Lbtope_fixed dataset contains 12063 B-cell epitopes and 20589 non-epitopes (Figure S3).

Lbtope_Fixed_non_redundant.

Using Lbtope_Fixed dataset, we have created an 80% non-redundant dataset using CD-HIT [14], [15]. The redundant dataset contains 7824 B-cell epitopes and 7853 non-epitopes.

Lbtope_Variable.

First, we removed all epitopes or non-epitopes having less than five residues or more than fifty residues. All epitopes common in B-cell epitopes and non B-cell epitopes were also removed. We found that majority of common epitopes are related to autoimmunity. Our final dataset Lbtope_Variable contains 14876 unique B-cell epitopes and 23321 unique non B-cell epitopes.

Lbtope_Variable_non_redundant.

We again created an 80% non-redundant Lbtope_Variable dataset using CD-HIT. We obtained 8011 B-cell epitopes and 10868 non-epitopes.

Lbtope_Confirm.

One of the challenges in creating dataset is its validity, though all epitopes, which we have extracted from IEDB, are experimentally tested. In order to improve the quality of epitopes/non-epitopes, we used only those epitopes/non-epitopes which reported in at least two studies. The final dataset Lbtope_Confirm contains 1042 unique B-cell epitopes and 1795 non B-cell epitopes.

Input Features

In this study, we generated and used various types of features of peptides that include binary profile or sparse matrix [16], physico-chemical profile (here only four properties tested) [17], dipeptide composition, Chen’s amino acid pair (AAP) propensities [18] and Composition-Transition-Distribution (CTD) profile [14]. All these features were already used in earlier methods (Text S1). Besides these features, we created modified AAP profile, termed as AAP* where instead of multiplication, each dipeptide value was assigned values to each dipeptide from a matrix, making a vector of 19 in place of 400 (Figure S1).

SVM and Weka Classifiers

In this study, we used SVM_light (http://svmlight.joachims.org/) package for implementing SVM technique. SVM has been used in several biological problems, including functional characterization of proteins [16], [19], [20]. Weka is a tool in which a number of algorithms like Baysian Network, SVMLib, Artificial Neural Network (ANN), Nearest Neighbor (IBk), Random Forest, etc. have been integrated with a user friendly, graphical front end. In weka, performance of all the methods can be compared on the single data set.

Cross Validation and Performance Measures

Although leave-one out or jackknife test is the best among cross-validation techniques, due to its time-consuming and heavy CPU requirements, n-fold cross validation is the optimum choice [21]. In this study, we have used five-fold cross validation on 90% data and remaining 10% data is used as independent dataset. We calculated sensitivity (Sen), Specificity (Spe), accuracy (Acc) and Matthew’s correlation coefficient (MCC) on an independent dataset. For detail description of these parameters see Ansari et al [20].

Results

Analysis of B-cell Epitopes

We analyzed B-cell epitopes to understand their charters tics. First, length wise distribution of B-cell epitopes was computed. As shown in Figure 1, most of the epitopes are in the range of 5–22 amino acid length. In order to understand the preference of residues in B-cell epitopes, we generated two-sample logo plot [22] using 20 mer epitope (upper panel) and non-epitope (lower panel). As shown in Figure 2, there is indeed elevated occurrence of surface accessible and flexible residue in the epitope region as compared to the non-epitope region. In addition, we observed high propensity of Proline and Glycine residue in the epitope region, which might be responsible for the creation of bends or flexibility in the epitope region.

thumbnail
Figure 1. Length-wise distribution of peptides (B-cell epitopes and non-epitopes), we divided peptides in different bins like peptides having a length less than five residues, having residues between 5 to 10 residues.

https://doi.org/10.1371/journal.pone.0062216.g001

thumbnail
Figure 2. Two-sample logo showing dominance of surface accessible residues in B-cell epitopes.

Yellow and black color residues indicate to surface accessible and non-accessible residues respectively.

https://doi.org/10.1371/journal.pone.0062216.g002

Performance of Binary Profile Based Models

We developed models for discriminating B-cell epitopes from non B-cell epitopes on Lbtope_Fixed dataset. SVM-based models have been developed using binary profile or sparse profile of patterns, which is represented by a vector length of Wx21 (W is window length, 20 in this study). Sparse matrix contains information for each position and each type of amino acids in the pattern. We achieved accuracy range from 37–67% with MCC of 0.03–0.22 and AUC of 0.65, which is better than random prediction (Table 1; Table S3).

thumbnail
Table 1. The performance of SVM models developed on Lbtope_Fixed dataset using various features.

https://doi.org/10.1371/journal.pone.0062216.t001

Performance of Models Based on Physico-chemical Properties

It is already known that physico-chemical properties of amino acids are responsible for structural and functional behavior of peptides and proteins. In our study, we have tested few topological properties (Table S2), which were shown to be a good index for B-cell epitope prediction such as relative connectivity, clustering coefficient, closeness and betweenness [23]. We developed SVM models using physico-chemical properties and achieved accuracy in the range of 43–64% with MCC of 0.06–0.13 and AUC of 0.58, which is poorer than models based on binary profile (Table 1; Table S4).

Performance of Composition Based Models

Besides understanding the positional effect of amino acids, we also computed and compared the overall composition of epitopes and non-epitopes (Figure S2). Similar to two-sample logo analysis, composition analysis revealed that it can be used to discriminate between epitopes and non-epitopes. Therefore, we have applied several distinctive types of models such as simple amino acid and dipeptide composition with different vector size. These models were trained and tested using SVM and IBk. While using SVM, simple amino acid composition performed best among binary and physico-chemical profiles with accuracy of 78%, MCC 0.53 and AUC 0.85. Dipeptide composition model performed better than Chen’s AAP with maximum accuracy of 81%, MCC 0.61 and AUC of 0.88, the highest among single feature models (Table S5, S6, S7, S8, S9). We tested different features of other algorithms implemented in Weka and found that IBk model performed best (results of other algorithms not shown). Dipeptide composition model performed better than AAP profile with accuracy of 81%, MCC 0.61 and AUC 0.86 (Table 2; Table S9).

thumbnail
Table 2. The performance of IBk models developed on Lbtope_Fixed dataset using various features.

https://doi.org/10.1371/journal.pone.0062216.t002

Since only composition-based method can be applied to variable data, we applied amino acid composition, CTD, AAP and dipeptide composition methods on Lbtope_Variable and Lbtope_Confirm datasets. It was observed that dipeptide-based method performed best among other methods with accuracy 75.89, 82.33 and MCC 0.51, 0.64 on Lbtope_Variable, Lbtope_Confirm dataset respectively (Table 3, 4, 5, 6; Table S10, S11, S12, S13, S14, S15, S16, S17). Performance of models was better in case of Lbtope_Confirm data as compare to Lbtope_Variable.

thumbnail
Table 3. The performance of SVM models developed on Lbtope_Variable dataset using various features.

https://doi.org/10.1371/journal.pone.0062216.t003

thumbnail
Table 4. The performance of IBk models developed on Lbtope_Variable dataset using various features.

https://doi.org/10.1371/journal.pone.0062216.t004

thumbnail
Table 5. The performance of SVM models developed on Lbtope_Confirm (epitope tested by at least two studies) dataset using various features.

https://doi.org/10.1371/journal.pone.0062216.t005

thumbnail
Table 6. The performance of IBk models developed on Lbtope_Confirm (epitope tested by at least two studies) dataset using various features.

https://doi.org/10.1371/journal.pone.0062216.t006

Performance on Non-redundant Peptide Dataset

Although we have considered unique epitopes, the redundancy could be expected among them similar to protein sequences. However, it is known that properties of peptide could change with a single amino acid variation. Nevertheless, to have an idea of redundancy and model performance, we created non-redundant dataset corresponding to both Lbtope_Fixed and Lbtope_Variable databases. We found that the number of peptides decreased as expected, but the performance remained significant.

We observed AUC of 0.61, 0.66 and 0.69 for simple amino acid, AAP and dipeptide composition respectively on Lbtope_Fixed_non_redundant dataset (Table S18, S19, S20, S21). We achieved better performance for Lbtope_Variable_non-redundant dataset. The AUC obtained was 0.60, 0.68 and 0.73 for simple amino acid, AAP and dipeptide composition respectively (Table S22, S23, S24, S25). In case of Lbtope_Fixed_non_redundant dataset, AUC sharply decreased from 0.88 to 0.69 and for Lbtope_Variable_non_redundant dataset AUC decreased from 0.83 to 0.73. Again, dipeptide based model performed better than other methods (Table 78; Table S18, S19, S20, S21, S22, S23, S24, S25).

thumbnail
Table 7. The performance of SVM models developed on Lbtope_Fixed_non_redundant dataset using various features.

https://doi.org/10.1371/journal.pone.0062216.t007

thumbnail
Table 8. The performance of SVM models developed on Lbtope_Variable_non_redundant dataset using various features.

https://doi.org/10.1371/journal.pone.0062216.t008

Benchmarking with Existing Methods

It is important to compare the newly developed algorithm with existing algorithms, which requires testing of all methods on same dataset. Unfortunately, our dataset is different than datasets used in previous studies. Thus one to one comparison is not feasible. In order to understand differences and similarities in our and existing models, we tested our models on datasets used in earlier methods. Similarly, we tested previously developed methods on our datasets. It was observed that earlier models failed on datasets used in this study, and our models failed on existing datasets (Table S26, S27, S28, S29, S30, S31, S32). Authors of ABCpred achieved sensitivity 57.14 and specificity 71.57 at the default threshold on their ABCpred dataset (Table 9). We assessed the performance of ABCpred at the default threshold on dataset Lbtope_Fixed used in this study and achieved sensitivity 54.55 and specificity 49.54 [11]. It was observed that sensitivity of ABCpred decreased slightly from 57.14 to 54.55 but specificity decreased drastically from 71.57 to 49.54 (Table 9). It suggested that ABCpred performance on B-cell epitopes decreased slightly but failed on non B-cell epitopes used in this study. We also evaluated models developed in this study on ABCpred and observed similar results. Models developed in this study failed on non B-cell epitopes (random peptides generated from proteins) used in ABCpred dataset (Table 9; Table S30, S31).

thumbnail
Table 9. The performance of our method LBtope on ABCpred dataset and performance of ABCpred on dataset Lbtope_Fixed (fixed length patterns used in this study).

https://doi.org/10.1371/journal.pone.0062216.t009

Similarly, we evaluated existing methods (BCPred and Chen’s method) on datasets used in this study. It was observed that these methods performed reasonably well on B-cell epitopes but failed on non B-cell epitopes (Table 10; Table S26, S27, S28, S29) [9], [10]. We also evaluated the performance of our dipeptide-based models on datasets used in previous studies (Table 10; Table S32) [14]. We observed similar trend, our models failed on negative patterns/example (random peptides used as non B-cell epitopes).

thumbnail
Table 10. The performance of BCPred, Chen on Lbtope_fixed and, Lbtope_variable dataset and performance of LBtope models on datasets used in previous studies.

https://doi.org/10.1371/journal.pone.0062216.t010

It can be suggested that existing models/methods perform reasonably fine on our B-cell epitopes, but failed on non B-cell epitopes. Similarly, our models failed on random peptides used in previous studies as non B-cell epitopes. This could be due to fact that our negative dataset comprised of experimentally verified non B-cell epitopes where as negative datasets of existing methods consist of random peptides generated from proteins.

To know the effect of using experimental proved non B-cell epitope instead of random peptides from Swiss-Prot in development of model. We created another dataset Lbtope-positive-fbcpred-negative, in which instead of experimental non B-cell epitopes, we used random peptides from FBCPred dataset as non B-cell epitopes. Next, we performed a five-fold cross validation on the above- dataset and obtained 85% sensitivity with 0.71 MCC. On the other hand, Lbtope_Confirm has achieved 81% sensitivity with 0.65 MCC, a bit poorer than Lbtope-positive-fbcpred-negative (Table 10; Table S33S34). Taken together all this results, it can be concluded that using experimental B-cell epitopes, the method can perform as good as using random peptides. In summary, performance of models depends upon the dataset used for training.

Implementation

In order to provide prediction service to scientific community, we have developed a user-friendly web server based on the model developed in this study. The server is developed using PHP 5.2.9, HTML and JavaScript as the front end and installed on a Red Hat Enterprise Linux 6 server environment. The server takes antigen primary amino acid sequence (s) in ‘FASTA’ format, generates 20 amino acids overlapping peptides for Lbtope_Fixed dataset model, 5–30 amino acids overlapping peptides for variable datasets model and predicts the linear epitopes. The non-redundant model is also implemented in case of very high specificity. The output is antigen sequence (s) mapped with B-cell epitopes with a probability scale of 20–80%. A higher score implies higher probability of peptide to be B-cell epitope. We have developed separate dedicated pages for antigen and peptide submission to avoid any complexity. In addition, we have developed a peptide mutation tool, which creates all possible single point mutations in given peptide and calculates the probability score based on the algorithm and also predicts the other properties. Using the mutation tool, user can design better epitopes or even choose fewer epitopic peptides for the de-immunization of therapeutic proteins. The web server is freely available at http://crdd.osdd.net/raghava/lbtope/.

Discussion

Epitope mapping is no doubt a very useful procedure, which has vast applications in the area of therapy and diagnostics. Experimental methods do exist, but they require time, resources and cannot handle the pace with which biological data is generated. Therefore, computer algorithms have been developed over the decades to predict the B- and T-cell epitopes from antigen sequence or structure if available. It’s observed that linear B-cell epitope prediction is more challenging than other epitope types like conformational B-cell or T-cell epitopes. This might be due to the reason that linear B-cell epitope posse’s variable length from 2–85 amino acids as compared to the almost fixed length core of the T-cell epitopes. This variability imposes several obstacles for algorithm developers. Besides variability, all the methods to date have been developed on very small data set with negatives examples obtained from randomly chosen UniProt peptides or same antigens, which are not experimentally validated. In the present study, for the first time we have used experimentally verified B-cell and non B-cell epitopes from IEDB database, which are much more in the number and rationally, created to previous methods. We created the 20 mer epitopes using corresponding ‘truncation-extension’ methodology and similar length, which were used in earlier methods. By using simple composition technique in combination with SVM and Weka implemented IBk; we came up with an algorithm, which is as good as existing tools. Performance of LBtope models decreased on non-redundant datasets, still performance remained as good as existing methods. It is also observed that model developed on Lbtope_Confirm dataset performed better than the models developed on Lbtope_Variable dataset. We have compared the performance of LBtope models on Lbtope and existing datasets. LBtope performs poor on negative dataset of existing methods, and they also performed poor on our negative dataset. It is because our negative dataset is experimentally verified B-cell epitope, whereas existing method, negative dataset were randomly generated from UniProt. We have implemented the algorithm in the form of a user-friendly web server: LBtope. The user can create the mutants of each peptide and test its epitopic or other desired probability using our server’s mutant tool. We hope that present model will aid the researchers in the field of linear B-cell epitope prediction.

Supporting Information

Figure S1.

Diagram showing calculation of Dipeptide composition, AAP and modified AAP (AAP*) from patterns.

https://doi.org/10.1371/journal.pone.0062216.s001

(TIF)

Figure S2.

Diagram showing % composition of B-cell epitopes and non-epitopes (LBtope data; 20 mers).

https://doi.org/10.1371/journal.pone.0062216.s002

(TIF)

Figure S3.

Flowchart showing preparation of LBtope dataset from IEDB database.

https://doi.org/10.1371/journal.pone.0062216.s003

(TIF)

Table S1.

Datasets used so far in the linear B-cell epitope prediction.

https://doi.org/10.1371/journal.pone.0062216.s004

(DOC)

Table S2.

Amino acid indices as obtained from Huang et al 2007.

https://doi.org/10.1371/journal.pone.0062216.s005

(DOC)

Table S3.

The performance of SVM models developed on Lbtope_Fixed dataset using binary profile. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s006

(DOC)

Table S4.

The performance of SVM models developed on Lbtope_Fixed dataset using Physico-chemical property (4 R indices). These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s007

(DOC)

Table S5.

The performance of SVM/IBK models developed on Lbtope_Fixed dataset using Amino acid composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s008

(DOC)

Table S6.

The performance of SVM/IBK models developed on Lbtope_Fixed dataset using Composition Transition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s009

(DOC)

Table S7.

The performance of SVM/IBK models developed on Lbtope_Fixed dataset using AAP profile. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s010

(DOC)

Table S8.

The performance of SVM/IBK models developed on Lbtope_Fixed dataset using AAA profile. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s011

(DOC)

Table S9.

The performance of SVM/IBK models developed on Lbtope_Fixed dataset using Dipeptide composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s012

(DOC)

Table S10.

The performance of SVM/IBK models developed on Lbtope_Variable dataset using Amino acid composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s013

(DOC)

Table S11.

The performance of SVM/IBK models developed on Lbtope_Variable dataset using Composition Transition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s014

(DOC)

Table S12.

The performance of SVM/IBK models developed on Lbtope_Variable dataset using AAP profile. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s015

(DOC)

Table S13.

The performance of SVM/IBK models developed on Lbtope_Variable dataset using Dipeptide composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s016

(DOC)

Table S14.

The performance of SVM/IBK models developed on Lbtope_Confirm (epitope tested by at least two studies) dataset using Amino acid composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s017

(DOC)

Table S15.

The performance of SVM/IBK models developed on Lbtope_Confirm (epitope tested by at least two studies) dataset using Composition Transition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s018

(DOC)

Table S16.

The performance of SVM/IBK models developed on Lbtope_Confirm (epitope tested by at least two studies) dataset using AAP profile. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s019

(DOC)

Table S17.

The performance of SVM/IBK models developed on Lbtope_Confirm (epitope tested by at least two studies) dataset using Dipeptide composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s020

(DOC)

Table S18.

The performance of SVM/IBK models developed on Lbtope_Fixed_non_redundant dataset using amino acid composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s021

(DOC)

Table S19.

The performance of SVM/IBK models developed on Lbtope_Fixed_non_redundant dataset using composition-transition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s022

(DOC)

Table S20.

The performance of SVM/IBK models developed on Lbtope_Fixed_non_redundant dataset using composition-transition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s023

(DOC)

Table S21.

The performance of SVM/IBK models developed on Lbtope_Fixed_non_redundant dataset using dipeptide composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s024

(DOC)

Table S22.

The performance of SVM/IBK model developed on Lbtope_Variable_non_redundant dataset using amino acid composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10%.

https://doi.org/10.1371/journal.pone.0062216.s025

(DOC)

Table S23.

The performance of SVM/IBK models developed on Lbtope_Variable_non_redundant dataset using composition-transition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s026

(DOC)

Table S24.

The performance of SVM/IBK models developed on Lbtope_Variable_non_redundant dataset using composition-transition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s027

(DOC)

Table S25.

The performance of SVM/IBK models developed on Lbtope_Variable_non_redundant dataset using dipeptide composition. These models were developed using 5-fold cross-validation on 90% data and tested on remaining 10% data.

https://doi.org/10.1371/journal.pone.0062216.s028

(DOC)

Table S26.

The performance of Chen’s AAP model on Lbtope_Fixed data.

https://doi.org/10.1371/journal.pone.0062216.s029

(DOC)

Table S27.

Performance of BCPred model (20 mer) on Lbtope_Fixed dataset.

https://doi.org/10.1371/journal.pone.0062216.s030

(DOC)

Table S28.

The performance of SVM models developed on Lbtope_Fixed dataset tested on Chen dataset.

https://doi.org/10.1371/journal.pone.0062216.s031

(DOC)

Table S29.

The performance of SVM models developed on Lbtope_Fixed dataset tested on BCPred dataset.

https://doi.org/10.1371/journal.pone.0062216.s032

(DOC)

Table S30.

Performance of ABCPred model (20 mer) on Lbtope_Fixed dataset.

https://doi.org/10.1371/journal.pone.0062216.s033

(DOC)

Table S31.

The performance of SVM models developed on Lbtope_Fixed and tested on ABCPred dataset.

https://doi.org/10.1371/journal.pone.0062216.s034

(DOC)

Table S32.

The performance of SVM models developed on Lbtope_Variable and tested on bcpred variable data.

https://doi.org/10.1371/journal.pone.0062216.s035

(DOC)

Table S33.

SVM results of five-fold cross validation on Lbtope_Confirm (epitope tested by at least two studies) dataset using dipeptide composition.

https://doi.org/10.1371/journal.pone.0062216.s036

(DOC)

Table S34.

SVM results of five-fold cross validation on Lbtope_positive_fbcpred_negative (B-cell epitope (positive) from Lbtope_Confirm dataset and random peptides (negative) from fbcpred) using dipeptide composition.

https://doi.org/10.1371/journal.pone.0062216.s037

(DOC)

Text S1.

Different B-cell epitope prediction techniques already used in earlier methods.

https://doi.org/10.1371/journal.pone.0062216.s038

(DOC)

Author Contributions

Conceived and designed the experiments: HS HRA GPSR. Performed the experiments: HS HRA. Analyzed the data: HS HRA GPSR. Contributed reagents/materials/analysis tools: HS HRA. Wrote the paper: HS HRA GPSR. Collected and performed the experiments: HS HRA.

References

  1. 1. Tomar N, De RK (2010) Immunoinformatics: an integrated scenario. Immunology 131: 153–168.
  2. 2. Dudek NL, Perlmutter P, Aguilar MI, Croft NP, Purcell AW (2010) Epitope discovery and their use in peptide based vaccines. Curr Pharm Des 16: 3149–3157.
  3. 3. Bryson CJ, Jones TD, Baker MP (2010) Prediction of immunogenicity of therapeutic proteins: validity of computational tools. BioDrugs 24: 1–8.
  4. 4. Steere AC, Drouin EE, Glickstein LJ (2011) Relationship between immunity to Borrelia burgdorferi outer-surface protein A (OspA) and Lyme arthritis. Clin Infect Dis 52 Suppl 3s259–265.
  5. 5. El-Manzalawy Y, Honavar V (2010) Recent advances in B-cell epitope prediction methods. Immunome Res 6 Suppl 2S2.
  6. 6. Xingdong Y, Xinglong Y (2009) An introduction to epitope prediction methods and software. Reviews in Medical Virology 19: 77–96.
  7. 7. Kringelum JV, Nielsen M, Padkjaer SB, Lund O (2013) Structural analysis of B-cell epitopes in antibody:protein complexes. Mol Immunol 53: 24–34.
  8. 8. Sweredoski M, Baldi P (2009) COBEpro: a novel system for predicting continuous B-cell epitopes. Protein Eng Des Sel 22: 113–120.
  9. 9. El-Manzalawy Y, Dobbs D, Honavar V (2008) Predicting linear B-cell epitopes using string kernels. J Mol Recognit 21: 243–255.
  10. 10. Chen J, Liu H, Yang J, Chou K (2007) Prediction of linear B-cell epitopes using amino acid pair antigenicity scale. Amino Acids 33: 423–428.
  11. 11. Saha S, Raghava G (2006) Prediction of continuous B-cell epitopes in an antigen using recurrent neural network. Proteins 65: 40–48.
  12. 12. Larsen JE, Lund O, Nielsen M (2006) Improved method for predicting linear B-cell epitopes. Immunome Res 2: 2.
  13. 13. Vita R, Zarebski L, Greenbaum JA, Emami H, Hoof I, et al. (2010) The immune epitope database 2.0. Nucleic Acids Res 38: D854–862.
  14. 14. EL-Manzalawy Y, Dobbs D, Honavar V (2008) Predicting flexible length linear B-cell epitopes. 7th International Conference on Computational Systems Bioinformatics: 121–131.
  15. 15. Li W, Godzik A (2006) Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22: 1658–1659.
  16. 16. Agarwal S, Mishra NK, Singh H, Raghava GP (2011) Identification of mannose interacting residues using local composition. PLoS One 6: e24039.
  17. 17. Ansari HR, Raghava GP (2010) Identification of conformational B-cell Epitopes in an antigen from its primary sequence. Immunome Res 6: 6.
  18. 18. Chen J, Liu H, Yang J, Chou KC (2007) Prediction of linear B-cell epitopes using amino acid pair antigenicity scale. Amino Acids 33: 423–428.
  19. 19. Yang ZR (2004) Biological applications of support vector machines. Brief Bioinform 5: 328–338.
  20. 20. Ansari HR, Raghava GP (2010) Identification of NAD interacting residues in proteins. BMC Bioinformatics 11: 160.
  21. 21. Chou KC, Zhang CT (1995) Prediction of protein structural classes. Crit Rev Biochem Mol Biol 30: 275–349.
  22. 22. Vacic V, Iakoucheva LM, Radivojac P (2006) Two Sample Logo: a graphical representation of the differences between two sets of sequence alignments. Bioinformatics 22: 1536–1537.
  23. 23. Huang J, Honda W, Kanehisa M (2007) Predicting B cell epitope residues with network topology based amino acid indices. Genome Inform 19: 40–49.