Development and validation of an epitope prediction tool for swine (PigMatrix) based on the pocket profile method

BackgroundT cell epitope prediction tools and associated vaccine design algorithms have accelerated the development of vaccines for humans. Predictive tools for swine and other food animals are not as well developed, primarily because the data required to develop the tools are lacking. Here, we overcome a lack of T cell epitope data to construct swine epitope predictors by systematically leveraging available human information. Applying the “pocket profile method”, we use sequence and structural similarities in the binding pockets of human and swine major histocompatibility complex proteins to infer Swine Leukocyte Antigen (SLA) peptide binding preferences.We developed epitope-prediction matrices (PigMatrices), for three SLA class I alleles (SLA-1*0401, 2*0401 and 3*0401) and one class II allele (SLA-DRB1*0201), based on the binding preferences of the best-matched Human Leukocyte Antigen (HLA) pocket for each SLA pocket. The contact residues involved in the binding pockets were defined for class I based on crystal structures of either SLA (SLA-specific contacts, Ssc) or HLA supertype alleles (HLA contacts, Hc); for class II, only Hc was possible. Different substitution matrices were evaluated (PAM and BLOSUM) for scoring pocket similarity and identifying the best human match. The accuracy of the PigMatrices was compared to available online swine epitope prediction tools such as PickPocket and NetMHCpan.ResultsPigMatrices that used Ssc to define the pocket sequences and PAM30 to score pocket similarity demonstrated the best predictive performance and were able to accurately separate binders from random peptides. For SLA-1*0401 and 2*0401, PigMatrix achieved area under the receiver operating characteristic curves (AUC) of 0.78 and 0.73, respectively, which were equivalent or better than PickPocket (0.76 and 0.54) and NetMHCpan version 2.4 (0.41 and 0.51) and version 2.8 (0.72 and 0.71). In addition, we developed the first predictive SLA class II matrix, obtaining an AUC of 0.73 for existing SLA-DRB1*0201 epitopes. Notably, PigMatrix achieved this level of predictive power without training on SLA binding data.ConclusionOverall, the pocket profile method combined with binding preferences from HLA binding data shows significant promise for developing T cell epitope prediction tools for pigs. When combined with existing vaccine design algorithms, PigMatrix will be useful for developing genome-derived vaccines for a range of pig pathogens for which no effective vaccines currently exist (e.g. porcine reproductive and respiratory syndrome, influenza and porcine epidemic diarrhea).

[1]  M. O. Dayhoff,et al.  22 A Model of Evolutionary Change in Proteins , 1978 .

[2]  M. O. Dayhoff A model of evolutionary change in protein , 1978 .

[3]  H. Takamatsu,et al.  Identification of novel foot-and-mouth disease virus specific T-cell epitopes in c/c and d/d haplotype miniature swine. , 2006, Virus research.

[4]  Leonard Moise,et al.  VennVax, a DNA-prime, peptide-boost multi-T-cell epitope poxvirus vaccine, induces protective immunity against vaccinia infection by T cell response alone. , 2011, Vaccine.

[5]  S. Altschul Amino acid substitution matrices from an information theoretic perspective , 1991, Journal of Molecular Biology.

[6]  K. Wiesmüller,et al.  Identification of T-cell epitopes in the structural and non-structural proteins of classical swine fever virus. , 2002, The Journal of general virology.

[7]  D. Andreu,et al.  Identification of T-Cell Epitopes in Nonstructural Proteins of Foot-and-Mouth Disease Virus , 2001, Journal of Virology.

[8]  Hau-San Wong,et al.  TEPITOPEpan: Extending TEPITOPE for Peptide Binding Prediction Covering over 700 HLA-DR Molecules , 2012, PloS one.

[9]  Jun Heon Lee,et al.  Sequence‐based characterization of swine leucocyte antigen alleles in commercially available porcine cell lines , 2009, International journal of immunogenetics.

[10]  Hans D. Mittelmann,et al.  MultiRTA: A simple yet reliable method for predicting peptide binding affinities for multiple class II MHC allotypes , 2010, BMC Bioinformatics.

[11]  I. Díaz,et al.  Predicted Peptides from Non-Structural Proteins of Porcine Reproductive and Respiratory Syndrome Virus Are Able to Induce IFN-γ and IL-10 , 2013, Viruses.

[12]  Hans-Georg Rammensee,et al.  MHC Ligands and Peptide Motifs , 1998, Molecular Biology Intelligence Unit.

[13]  J K Lunney,et al.  Nomenclature for factors of the swine leukocyte antigen class II system, 2005. , 2005, Tissue antigens.

[14]  M. Nielsen,et al.  Identification of peptides from foot-and-mouth disease virus structural proteins bound by class I swine leukocyte antigen (SLA) alleles, SLA-1*0401 and SLA-2*0401. , 2013, Animal genetics.

[15]  I. Díaz,et al.  In silico prediction and ex vivo evaluation of potential T-cell epitopes in glycoproteins 4 and 5 and nucleocapsid protein of genotype-I (European) of porcine reproductive and respiratory syndrome virus. , 2009, Vaccine.

[16]  O. Lund,et al.  NetMHCIIpan-3.0, a common pan-specific MHC class II prediction method including all three human MHC class II isotypes, HLA-DR, HLA-DP and HLA-DQ , 2013, Immunogenetics.

[17]  J K Lunney,et al.  Nomenclature for factors of the SLA class-I system, 2004. , 2005, Tissue antigens.

[18]  A. D. De Groot,et al.  Epitope-based vaccination against pneumonic tularemia. , 2009, Vaccine.

[19]  S. Buus,et al.  Analysis of Swine Leukocyte Antigen Peptide Binding Profiles and the Identification of T cell Epitopes by Tetramer Staining , 2012 .

[20]  Ronald E. Bontrop,et al.  Immunogenetics , 2005, Genes and Immunity.

[21]  F. Osorio,et al.  Location of T-cell epitopes in nonstructural proteins 9 and 10 of type-II porcine reproductive and respiratory syndrome virus. , 2012, Virus research.

[22]  Jean-Philippe Vert,et al.  Efficient peptide-MHC-I binding prediction for alleles with few known binders , 2008, Bioinform..

[23]  A. D. De Groot,et al.  Time for T? Immunoinformatics addresses vaccine design for neglected tropical and emerging infectious diseases , 2014, Expert review of vaccines.

[24]  Morten Nielsen,et al.  The PickPocket method for predicting binding specificities for receptors based on receptor pocket similarities: application to MHC-peptide binding , 2009, Bioinform..

[25]  Morten Nielsen,et al.  Porcine major histocompatibility complex (MHC) class I molecules and analysis of their peptide-binding specificities , 2011, Immunogenetics.

[26]  Julie A McMurry,et al.  In silico-accelerated identification of conserved and immunogenic variola/vaccinia T-cell epitopes. , 2009, Vaccine.

[27]  Qing Zhang,et al.  Immune epitope database analysis resource (IEDB-AR) , 2008, Nucleic Acids Res..

[28]  O. Lund,et al.  NetMHCpan, a Method for Quantitative Predictions of Peptide Binding to Any HLA-A and -B Locus Protein of Known Sequence , 2007, PloS one.

[29]  J. Neefjes,et al.  Towards a systems understanding of MHC class I and MHC class II antigen presentation , 2011, Nature Reviews Immunology.

[30]  Immunogenicity and T cell recognition in swine of foot-and-mouth disease virus polymerase 3D. , 2004, Virology.

[31]  A. D. De Groot,et al.  An interactive Web site providing major histocompatibility ligand predictions: application to HIV research. , 1997, AIDS research and human retroviruses.

[32]  A. D. De Groot,et al.  Peptide-pulsed dendritic cells induce the hepatitis C viral epitope-specific responses of naïve human T cells. , 2014, Vaccine.

[33]  Anne S De Groot,et al.  Reducing risk, improving outcomes: bioengineering less immunogenic protein therapeutics. , 2009, Clinical immunology.

[34]  Ora Schueler-Furman,et al.  Learning MHC I - peptide binding , 2006, ISMB.

[35]  Ron D. Appel,et al.  ExPASy: the proteomics server for in-depth protein knowledge and analysis , 2003, Nucleic Acids Res..

[36]  G. Gao,et al.  Crystal Structure of Swine Major Histocompatibility Complex Class I SLA-1*0401 and Identification of 2009 Pandemic Swine-Origin Influenza A H1N1 Virus Cytotoxic T Lymphocyte Epitope Peptides , 2011, Journal of Virology.

[37]  M. O. Dayhoff,et al.  Atlas of protein sequence and structure , 1965 .

[38]  M F del Guercio,et al.  Several common HLA-DR types share largely overlapping peptide binding repertoires. , 1998, Journal of immunology.

[39]  G. Allan,et al.  T lymphocyte epitope mapping of porcine circovirus type 2. , 2007, Viral immunology.

[40]  L. E. Pedersen,et al.  Identification of swine influenza virus epitopes and analysis of multiple specificities expressed by cytotoxic T cell subsets , 2014, Virology Journal.

[41]  G. Chelvanayagam A roadmap for HLA-A, HLA-B, and HLA-C peptide binding specificities , 1996, Immunogenetics.

[42]  Bjoern Peters,et al.  HLA class I supertypes: a revised and updated classification , 2008, BMC Immunology.

[43]  Julie A McMurry,et al.  Epitope-driven TB vaccine development: a streamlined approach using immuno-informatics, ELISpot assays, and HLA transgenic mice. , 2007, Current molecular medicine.

[44]  Immunoinformatics prediction of linear epitopes from Taenia solium TSOL18 , 2011, Bioinformation.

[45]  Don C. Wiley,et al.  Crystal structure of the human class II MHC protein HLA-DR1 complexed with an influenza virus peptide , 1994, Nature.

[46]  J. Lunney,et al.  Molecular characterization of swine leucocyte antigen class II genes in outbred pig populations. , 2009, Animal genetics.

[47]  Judy Lieberman,et al.  Mapping cross-clade HIV-1 vaccine epitopes using a bioinformatics approach. , 2003, Vaccine.

[48]  O. Lund,et al.  NetMHCpan, a method for MHC class I binding prediction beyond humans , 2008, Immunogenetics.

[49]  U. Şahin,et al.  Generation of tissue-specific and promiscuous HLA ligand databases using DNA microarrays and virtual HLA class II matrices , 1999, Nature Biotechnology.

[50]  L. Stern,et al.  HLA-DR: molecular insights and vaccine design. , 2009, Current pharmaceutical design.

[51]  G Chelvanayagam,et al.  A roadmap for HLA-DR peptide binding specificities. , 1997, Human immunology.

[52]  V. Nene,et al.  T cell epitope identification for bovine vaccines: an epitope mapping method for BoLA A-11. , 2003, International journal for parasitology.

[53]  Leonard Moise,et al.  HelicoVax: epitope-based therapeutic Helicobacter pylori vaccination in a mouse model. , 2011, Vaccine.

[54]  B. Ober,et al.  Vaccine-Induced, Pseudorabies Virus-Specific, Extrathymic CD4+CD8+ Memory T-Helper Cells in Swine , 1998, Journal of Virology.

[55]  Yanjun Zhou,et al.  Identification of immunodominant T-cell epitopes in membrane protein of highly pathogenic porcine reproductive and respiratory syndrome virus. , 2011, Virus research.

[56]  Adeeba Kamarulzaman,et al.  AIDS Res Hum Retroviruses , 2006 .

[57]  T. Goldberg,et al.  Identification of immunodominant T-cell epitopes present in glycoprotein 5 of the North American genotype of porcine reproductive and respiratory syndrome virus. , 2008, Vaccine.

[58]  D. Madden The three-dimensional structure of peptide-MHC complexes. , 1995, Annual review of immunology.

[59]  Shanfeng Zhu,et al.  MHC2SKpan: a novel kernel based approach for pan-specific MHC class II peptide binding prediction , 2013, BMC Genomics.

[60]  Shen Jean Lim,et al.  Immunoinformatics: how in silico methods are re-shaping the investigation of peptide immune specificity. , 2012, Frontiers in bioscience.

[61]  J. Lunney,et al.  Molecular characterization of swine leucocyte antigen class I genes in outbred pig populations. , 2009, Animal genetics.