ReseaRch papeR

ReseaRch papeR

OncoImmunology 1:8, 1281–1289; November 2012; © 2012 Landes Bioscience

In silico prediction of tumor antigens derived from functional missense mutations of the cancer gene census
Jahan s. Khalili,1,†,* Russell W. hanson2,† and Zoltan szallasi2,3
1 Departments of Melanoma Medical Oncology and systems Biology; University of Texas M.D. anderson cancer center; houston, TX Usa; 2children’s hospital Informatics program; Division of health sciences and Technology; harvard–Massachusetts Institute of Technology; harvard Medical school; Boston, Ma Usa; 3Department of systems Biology; center for Biological sequence analysis; Technical University of Denmark; Lyngby, Denmark

†

These authors contributed equally to this work.

Keywords: cancer vaccines, computational biology, immunomics, immunotherapy, missense mutation, protein database, T cell therapy

antigen-specific immune responses against peptides derived from missense gene mutations have been identified in multiple cancers. The application of personalized peptide vaccines based on the tumor mutation repertoire of each cancer patient is a near-term clinical reality. These peptides can be identified for pre-validation by leveraging the results of massive gene sequencing efforts in cancer. In this study, we utilized NetMhc 3.2 to predict nanomolar peptide binding affinity to 57 human hLa-a and B alleles. all peptides were derived from 5,685 missense mutations in 312 genes annotated as functionally relevant in the cancer Genome project. Of the 26,672,189 potential 8–11 mer peptide-hLa pairs evaluated, 0.4% (127,800) display binding affinities < 50 nM, predicting high affinity interactions. These peptides can be segregated into two groups based on the binding affinity to hLa proteins relative to germline-encoded sequences: peptides for which both the mutant and wild-type forms are high affinity binders, and peptides for which only the mutant form is a high affinity binder. current evidence directs the attention to mutations that increase hLa binding affinity, as compared with cognate wild-type peptide sequences, as these potentially are more relevant for vaccine development from a clinical perspective. Our analysis generated a database including all predicted hLa binding peptides and the corresponding change in binding affinity as a result of point mutations. Our study constitutes a broad foundation for the development of personalized peptide vaccines that hone-in on functionally relevant targets in multiple cancers in individuals with diverse hLa haplotypes.

Introduction Missense mutations in cancer cells can generate unique T cell epitopes.1 Natural antigen-specific T-cell responses against missense mutant epitopes in acute myeloid leukemia (AML), melanoma, renal and lung cancer have been discovered using traditional discovery methods based on CD8 + tumor infiltrating lymphocytes (TILs).2–6 Such foreign epitopes have been shown is to facilitate immunosurveillance and cancer control in chemically-induced murine models of carcinogens.7 Furthermore, sequencing efforts have identified 50 mutated peptides that can serve as rejection antigens in transplanted B16 murine melanoma.8 The marked increase in tumor sequencing efforts has provided an opportunity to discover mutant epitopes that strongly bind to human HLA, compared with wild-type peptides, and may thus be effective mediators of T-cell responses. Thus, knowledge of a patient’s HLA type and missense mutation profile provides the opportunity to develop personalized peptide vaccines. Evidence supporting the efficacy of vaccination strategies based on mutant epitopes
*Correspondence to: Jahan S. Khalili; Email: jahkhalili@mdanderson.org Submitted: 06/18/12; Accepted: 07/17/12 http://dx.doi.org/10.4161/onci.21511 www.landesbioscience.com

has already been generated. For instance, peptides derived from codon 3 mutations in RAS family members can readily induce immunity in patients with pancreatic, lung, and colon cancer patients.1,9–14 However, RAS mutations do not generate optimal antigens in most cancer patients. Indeed, RAS is a relatively commonly mutated gene, yet is present only in a minority of patients. Models of peptide-HLA binding affinity can facilitate the identification of novel and “personal” targets for cancer vaccines. To the best of our knowledge, computational methods to discover mutant epitopes and the differential binding affinity to HLA were first applied by Segal et al. to a data set from breast and colorectal tumors.15 In this study, 1,152 peptides were interrogated for HLA-A*02:01 binding in silico. Similarly, Warren et al. subjected a general survey of mutations to multiple in silico HLA-binding algorithms in an effort to identify a polyvalent peptide vaccine optimized for prophylactic use.16 HLA allelic frequency in the United States, mutation frequency and tumor subtype frequency were given equal consideration to generate the proposed vaccine formulation.

OncoImmunology

1281

Figure 1. allelic distribution of tight binding mutated peptide-hLa pairs. allelic distribution of predicted tight binding (affinity score < 50 nM) mutated peptides from the cancer Gene census among Mhc class I alleles. (A) human hLa (n = 52), null data not shown for n = 5 hLa alleles (n = 5). (B) Murine h-2 alleles (n = 6). (C) Rhesus macaque Mamu alleles (n = 7). (D) chimpanzee patr alleles (n = 9).

Since the binding affinity of a peptide for HLA proteins is associated with the immunogenicity of the peptide, computational prediction is a convenient and practical first step toward the identification of optimal vaccine targets for cancer therapy.17 The ongoing development of peptide-MHC binding algorithms has benefitted from progressively more empirical binding data, allowing for increased prediction accuracy.18,19 By means of the artificial neural network-based prediction algorithm NetMHC3.2, developed by Lundegaard et al., highly accurate predictions of peptide binding to large numbers of HLA alleles are now possible.18 The prediction of peptide sequence-dependent antigen presentation, which entails multiple steps including proteosomal cleavage and TAP binding, is by far more complex (and hence less developed) than that of peptide-HLA interactions. Moreover, peptide processing is highly variable across multiple types of cancer cells and inflammatory states. The analysis performed in this study yields a potential functional class of immunotherapeutic targets that possess increased

HLA binding affinity upon mutation. In addition, we present a set of peptides predicted to bind HLA with similar affinity before and after a cancer-associated mutation. Thus, we offer a foundational database to support research of hypotheses related to the immunogenicity of peptides derived from missense cancerassociated mutations. Results Execution of MHC Class I binding peptide prediction with NetMHC 3.2. In order to identify peptides that may serve as tumor rejection antigens based on the predicted ability to bind human HLA, a database of missense mutation-derived peptides was assembled from the Catalogue of Somatic Mutations in Cancer (COSMIC) database. Mutations (n = 5,590) from 312 genes (Table S1) that are represented in the COSMIC database were used as a resource for the amino acid substitution affecting the wild-type sequence. Short peptide strings of 8, 9, 10 and 11

1282

OncoImmunology

Volume 1 Issue 8

Figure 2. Tight binding mutated peptides derived from facilitating mutations and human hLa alleles. (A) allelic distribution of predicted tightbinding peptides derived from facilitating mutations (mutated peptide affinity score < 50 nM aND wild-type peptide affinity score > 500 nM) from the cancer Gene census among hLa alleles (n = 52), null data not shown for n = 5 hLa alleles. (B) plot of tight-binding mutated peptides and peptides derived from faciliating mutations for each cognate hLa allele. Linear regression, y = 13.22(×) − 82.942, R2 as indicated.

amino acid lengths were generated from every mutation in the primary database. Each mutation resulted in 76 unique strings. In total, 26,672,189 peptides were generated for analysis. Each short peptide string was submitted to the NetMHC 3.2 HLAbinding algorithm. The results of this analysis were nanomolar (nM) binding affinity prediction scores for human HLA-A alleles and HLA-B alleles as indicated in the methods section.

From the 26,672,189 mutated peptides and 57 distinct human HLA alleles interrogated, 127,801 (0.3%) unique peptide-HLA pairs achieved a binding score of 50 nM of less (i.e., tight binding). This set of peptides was subjected to further analysis. Mutation-mediated alteration in HLA-binding affinity. For the purpose of comparison, a delta table was generated to display the change in predicted HLA affinity that occurred as a result of

www.landesbioscience.com

OncoImmunology

1283

each missense mutation, compared with the wild-type sequence. This table is organized with a FASTA-formatted identifier for each mutation and HLA pair, reference to the position of the mutation in the peptide string, the peptide length, and whether the peptide referenced is the wild-type or mutant sequence. Following the FASTA-formatted identifier is the peptide sequence and the nM binding affinity score for the peptide and indicated HLA allele. This order is repeated for the mutated peptide sequence and each row ends with the calculated difference in nM affinity between the wildtype and mutated peptides for an individual HLA allele (Table S2). Some human HLA alleles present a broader repertoire of mutated peptides than others (Fig. 1A). The frequently studied human allele HLA-A0201 ranks 15th among 57 alleles in the breadth of its repertoire. The mean number of mutated peptides that human alleles may bind is 1,767 with a standard deviation of 2,377. The number of mutated peptides that any human allele
Table 1. Facilitating, mutated strong binding hLa-a 02:01 peptides FASTA wt > aLK-R401Q-wildtype-sequence1-hLa-a0201 > Bap1-h169Q-wildtype-sequence1-hLa-a0201 > BRaF-K475M-wildtype-sequence1-hLa-a0201 > cDK6-p199L-wildtype-sequence1-hLa-a0201 > cheK2-p536L-wildtypesequence-1-hLa-a0201 > eGFR-h773L-wildtype-sequence8-hLa-a0201 > FaNcF-p185L-wildtype-sequence1-hLa-a0201 > GNas-D141V-wildtype-sequence8-hLa-a0201 > ITK-G372V-wildtype-sequence-7hLa-a0201 > JaK1-e966V-wildtype-sequence8-hLa-a0201 > JaK2-K539L-wildtype-sequence8-hLa-a0201 > KRas-Q61L-wildtype-sequence10-hLa-a0201 > NOTch1-R1634L-wildtypesequence-1-hLa-a0201 > RB1-p515L-wildtype-sequence-1hLa-a0201 > Tp53-p47L-wildtype-sequence-8hLa-a0201 wt Peptide FRVaLeYI FhFVsYVpI GKWhGDVaV TpVDLWsV RpaVcaaV VMasVDNph RpaRFLssL sVMNVpDFD FVQeIGsG FLpsGsLKe hMNQMVFhK cLLDILDTaGQ KRaaeGWaa FpWILNVL aMDDLMLsp nM (wt) 15009 13143 15489 14601 20392 22464 22304 20809 19247 17955 18253 6354 21734 10736 10776

may bind with high affinity range from 0–11,005; the quartiles are bound by 0, 131, 768, 2,692 and 11,005. Analysis of mutated peptide binding to a limited number of MHC Class I alleles from mouse, rhesus macaque and chimpanzees demonstrates that these species may also be capable of presenting this set of antigens (Fig. 1B–D). Analysis of the number of tight-binding mutated peptides for which the corresponding wild-type peptide has a nonbinding score (greater than 500 nM) is reported in Figure 2. These peptide-MHC pairs numbered 7,982 (8% of the total tight-binding mutated peptides). The mutations leading to peptides that may exhibit greater binding to HLA than wild-type counterparts are referred to as facilitating mutations. A partial set of facilitating mutations for HLA-A02:01 is collected in Table 1 (n = 15). Neither the number of peptides binding each HLA molecule, nor facilitating mutations are equally distributed among the HLA alleles analyzed. However, the number of facilitating mutations

FASTA mt > aLK-R401Q-mutant-sequence-1hLa-a0201 > Bap1-h169Q-mutant-sequence-1hLa-a0201 > BRaF-K475M-mutant-sequence-1hLa-a0201 > cDK6-p199L-mutant-sequence-1hLa-a0201 > cheK2-p536L-mutant-sequence1-hLa-a0201 > eGFR-h773L-mutant-sequence-8hLa-a0201 > FaNcF-p185L-mutant-sequence-1hLa-a0201 > GNas-D141V-mutant-sequence-8hLa-a0201 > ITK-G372V-mutant-sequence-7hLa-a0201 > JaK1-e966V-mutant-sequence-8hLa-a0201 > JaK2-K539L-mutant-sequence-8hLa-a0201 > KRas-Q61L-mutant-sequence-10hLa-a0201 > NOTch1-R1634L-mutantsequence-1-hLa-a0201 > RB1-p515L-mutant-sequence-1hLa-a0201 > Tp53-p47L-mutant-sequence-8hLa-a0201

mt Peptide FQVaLeYI FQFVsYVpI GMWhGDVaV TLVDLWsV RLaVcaaV VMasVDNpL RLaRFLssL sVMNVpDFV FVQeIGsV FLpsGsLKV hMNQMVFhL cLLDILDTaGL KLaaeGWaa FLWILNVL aMDDLMLsL

nM (mt) 29 21 13 12 25 48 38 24 48 11 35 26 24 13 11

Delta nM 14980 13122 15476 14589 20367 22416 22266 20785 19199 17944 18218 6328 21710 10723 10765

Table of representative mutated peptides from the cancer Gene census predicted to be tight binders (affinity < 50 nM) to hLa-a02:01, for which cognate wild-type peptides are predicted to be non-binders (affinity > 500 nM). n = 15. 1284 OncoImmunology Volume 1 Issue 8

Figure 3. Facilitating mutation distribution for human hLa-a and hLa-B alleles. plot of tight binding mutated peptides (< 50 nM mutated peptide affinity score) from the cancer Gene census and corresponding wild-type peptides affinity score for each cognate hLa allele. Light green background indicates the threshold (500 nM) of predicted non-binding wild-type peptides.

was highly associated with the total number of tight binding mutated peptides analyzed R = 0.97 (n = 57) (Fig. 2). There is no enrichment or depletion of facilitating mutations among peptidHLA pairs as represented in the Cancer Gene Census database, a collection of the observable functional missense mutations in human cancer (Fig. 3). The identified facilitating mutations are largely composed of peptides in which the mutated residue is located in a peripheral residue position (Fig. 4A and B). In the large majority of peptide-MHC pairs that remain in the total data set, mutations exert limited effect upon the HLAbinding scores. In this set of 79,733 neutral mutations, 79% of the tight-binding peptides retain < 50 nM binding scores with and without the point mutation numbers (Fig. 3). A partial set of neutral mutations for HLA-A02:01 is collected in Table 2 (n = 15). Discussion In this work, we have analyzed peptides that contain mutated amino acids from 5685 mutations in 312 genes for binding to HLA alleles in silico. All genes and mutations were selected from the Cancer Gene Census and annotated in the COSMIC database to have known functional roles in cancer.20 This curated data set of functional mutations was selected for analysis to facilitate vaccine efforts that may be able to exert immunological pressure on genes required for tumor cell survival. In addition to focusing our search on genes and mutations with important survival

properties, we have identified peptides that we predict will demonstrate enrich surface presentation as a result of enhanced HLAbinding affinity.21 Mutated peptides are a class of antigens with superior properties as immunological targets for cancer therapy. Self-protein derived peptides must overcome multiple and redundant mechanisms of T-cell suppression due to peripheral tolerance.22 Additionally, the frequency of autoreactive T cells is attenuated due to central tolerance. Mutated peptides are absent in the thymus during development and are not endogenously expressed in normal tissues. Moreover, altered protein folding of mutated proteins generally reduces their stability. Reduced stability in some cases is observed in the presence of stabilizing the heat-shock 90 KDa protein (HSP90) and in other cases by the increased proteosomal degradation accounted for by peptide and protein decay experiments.23 Loss-of-function of tumor suppressor genes can be facilitated by a rapid protein turnover due to destabilization. The increased source of peptides resulting from decreased stability can enrich these peptides as potential HLA binders as compared with their wild-type counterparts, as in the case of neutral effects on HLA-binding affinity.24 A simplistic mass action mechanism for increased peptide presentation may underlie the enhanced surface presentation of these antigens. However, if the HLA-binding affinity of a peptide is increased due to a mutation, this may translate into increased presentation

www.landesbioscience.com

OncoImmunology

1285

Figure 4. For figure legend, see page 1287.

1286

OncoImmunology

Volume 1 Issue 8

Figure 4 (See opposite page). Facilitating mutations utilize peripheral anchor residues. plot of tight binding mutated 8, 9, 10 and 11 mer peptides (mutated peptide affinity score < 50 nM) and corresponding wild-type peptides affinity score for each cognate hLa-a (a) and hLa-B (B) allele. Light green background indicates the threshold (500 nM) of predicted non-binding wild-type peptides. coloring indicates the position of each mutation in the peptide string (starting from the c terminus): dark blue (1), light blue (2), orange (3), light orange (4), dark green (5), light green (6), red (7), pink (8), purple (9), light purple (10), brown (11). Table 2. Neutral, mutated strong binding hLa-a 02:01 peptides FASTA wt > aLK-R1275Q-wildtypesequence-3-hLa-a0201 > Bap1-h169Q-wildtypesequence-6-hLa-a0201 > BRaF-L618W-wildtypesequence-2-hLa-a0201 > cDK6-p199L-wildtypesequence-9-hLa-a0201 > cRLF2-F232c-wildtypesequence-4-hLa-a0201 > eGFR-D761G-wildtypesequence-2-hLa-a0201 > FBXW7-I563T-wildtypesequence-5-hLa-a0201 > GNa11-V223M-wildtypesequence-1-hLa-a0201 > IDh2-V294M-wildtypesequence-5-hLa-a0201 > JaK1-D660h-wildtypesequence-8-hLa-a0201 > JaK2-K607N-wildtypesequence-4-hLa-a0201 > KRas-Q22K-wildtypesequence-0-hLa-a0201 > NOTch1-s1598I-wildtypesequence-5-hLa-a0201 > ReT-D925h-wildtypesequence-3-hLa-a0201 > TcL1a-K23N-wildtypesequence-6-hLa-a0201 wt Peptide GMaRDIYRa RTMeaFhFV sILWMapeV VLLQssYaTpV KLsKFILI ILDeaYVMasV sLDTsIRV NVTsIMFLV LIDDMVaQV YLYGVcVRDV KLshKhLV QLIQNhFV FLReLsRV sLFDhIYTT RLWaWeKFV nM (wt) 43 9 20 12 25 11 23 48 18 23 44 42 40 13 11 FASTA mt > aLK-R1275Q-mutant-sequence3-hLa-a0201 > Bap1-h169Q-mutant-sequence6-hLa-a0201 > BRaF-L618W-mutant-sequence2-hLa-a0201 > cDK6-p199L-mutant-sequence9-hLa-a0201 > cRLF2-F232c-mutant-sequence4-hLa-a0201 > eGFR-D761G-mutant-sequence2-hLa-a0201 > FBXW7-I563T-mutant-sequence5-hLa-a0201 > GNa11-V223M-mutantsequence-1-hLa-a0201 > IDh2-V294M-mutant-sequence5-hLa-a0201 > JaK1-D660h-mutant-sequence8-hLa-a0201 > JaK2-K607N-mutant-sequence4-hLa-a0201 > KRas-Q22K-mutant-sequence0-hLa-a0201 > NOTch1-s1598I-mutantsequence-5-hLa-a0201 > ReT-D925h-mutant-sequence3-hLa-a0201 > TcL1a-K23N-mutant-sequence6-hLa-a0201 mt Peptide GMaQDIYRa RTMeaFQFV sIWWMapeV VLLQssYaTLV KLsKcILI ILGeaYVMasV sLDTsTRV NMTsIMFLV LIDDMMaQV YLYGVcVRhV KLshNhLV KLIQNhFV FLReLIRV sLFhhIYTT RLWaWeNFV nM (mt) 33 15 10 13 24 20 44 9 23 15 29 11 20 32 8 Delta nM 10 6 10 1 1 9 21 39 5 8 15 31 20 19 3

Table of representative mutated peptides from the cancer Gene census predicted to be tight binders (affinity < 50 nM) to hLa-a02:01, for which cognate wild-type peptides are also predicted to be tight binders (affinity < 50 nM). n = 15.

of the mutant peptide. During the loading of peptides on HLA in the ER and endosomes, competition can occur. Thus, facilitating mutations are subject to this advantage in loading and presentation.24 Once exposed on the cell surface, stable HLA/B2M/ peptide complexes are presented for a relatively long time and hence have a higher chance to activate T cells.21,24 Interestingly, in a study that determined the antigen specificity of T cells in an immunosurveillance model of murine sarcoma, the mutated peptide from spectrin β2 was shown to contain a facilitating mutation. Identified in silico as a MHC Class I (H-2-Db) binder after exome sequencing, the VAVVNQIAL peptide possesses a

binding score of 9 nM based on the methods used in this study. The cognate wildtype peptide VAVVNQIAR achieves only a score of 5304 nM.7 In part, the principle of HLA-specific peptide anchor residues motivates the search for these target peptides. However, the artificial neural network-trained binding approach of NetMHC 3.2 functions independently of the recognized dominant influence that some residues acquire in peripheral positions. Therefore it was not surprising to observe the high prevalence of peptides with mutations in the anchor position within the population of facilitating mutations. This result indicates that the machine-learning model embedded in NetMHC

www.landesbioscience.com

OncoImmunology

1287

3.2 for predicting peptide-HLA binding confirms experimentally informed theoretical expectations. This study provides an analysis of the Cancer Gene Census database revealing a subset of mutated peptides with superior HLA-binding affinity for the development of personalized peptide vaccines that target functionally relevant targets in diverse cancer and diverse HLA haplotypes. For this approach, substantial rationale exists for the continued use of immunostimulatory treatments such as anti-CTLA4 or anti-PD-1 monoclonal antibodies in combination with peptide vaccines.22 Still, in the clinical and experimental use of self peptides as tumor antigens, autoimmune responses are a source of significant toxicity when combined with anti-CTLA4 therapy.22,25 For this reason, as we and others previously stated, the discovery and characterization of tumor-specific epitopes is still a priority in immunotherapy. Exerting pressure on cancer cells with functional mutations that facilitate tumor survival may provide an additional level of protection, as the loss of such antigens may be detrimental to the tumorigenic phenotype. The generation of this evolutionary double bind is predicted to have superior effects in a mutating target cell population.26 This class of peptides displays conserved molecular homology and exerts lethal pressure on cells in relation to oncogenic properties, which underlies their theoretical value as clinical targets for T-cell therapy and as ligands for TCR-like antibodies.27–29 The results of this study are critical to the highthroughput experimental validation of mutated tumor epitopes in anticipation of the broad implementation of personalized mutation-specific cancer vaccines. Materials and Methods Mutated peptide library generation. Point substitutions observed in cancer cell lines and tumor specimens were collected from the Catalogue of Somatic Mutations in Cancer (COSMIC) database in February 2011 and February 2012.20 Mutant genes were taken from the Cancer Genes of the Cancer Gene Census. This expanding collection of validated cancer mutations is curated to maintain a collection of mutations extending functional properties to cancer cells.20 Gene mutations resulting in altered peptide sequences
References
1. Tuttle TM. Oncogene products represent potential targets of tumor vaccines. Cancer Immunol Immunother 1996; 43:135-41; PMID:9001566; http://dx.doi. org/10.1007/s002620050314. Gaudin C, Kremer F, Angevin E, Scott V, Triebel F. A hsp70-2 mutation recognized by CTL on a human renal cell carcinoma. J Immunol 1999; 162:1730-8; PMID:9973436. Echchakir H, Mami-Chouaib F, Vergnon I, Baurain JF, Karanikas V, Chouaib S, et al. A point mutation in the alpha-actinin-4 gene generates an antigenic peptide recognized by autologous cytolytic T lymphocytes on a human lung carcinoma. Cancer Res 2001; 61:4078-83; PMID:11358829. Corbière V, Chapiro J, Stroobant V, Ma W, Lurquin C, Lethé B, et al. Antigen spreading contributes to MAGE vaccination-induced regression of melanoma metastases. Cancer Res 2011; 71:1253-62; PMID:21216894; http://dx.doi.org/10.1158/0008-5472.CAN-10-2693. 5.

were used to generate a library of 8, 9, 10 and 11 amino acid long peptides for further in silico screening. Reference sequences for each gene were similarly acquired from the COSMIC database, and 8, 9, 10, and 11 amino acid wild-type peptides corresponding to each mutated peptide were also generated. In silico prediction of antigen presentation. Peptides were submitted to NetMHC 3.2 for prediction of binding affinity for all MHC Class I alleles available.18,30 Human HLA-A alleles: 01:01, 02:01, 02:02, 02:03, 02:06, 02:11, 02:12, 02:16, 02:19, 02:50, 03:01, 11:01, 23:01, 24:02, 24:03, 25:01, 26:01, 26:02, 26:03, 29:02, 30:01, 30:02, 31:01, 32:01, 33:01, 68:01, 68:02, 69:01, 80:01; HLA-B alleles; 07:02, 08:01, 08:02, 08:03, 15:01, 15:02, 15:03, 15:09, 15:17, 18:01, 27:03, 27:05, 35:01, 38:01, 39:01, 40:01, 40:02, 44:02, 44:03, 45:01, 46:01, 48:01, 51:01, 53:01, 54:01, 57:01, 58:01, and 73:01 were included in the analysis. Murine alleles: H-2-Db, H-2-Dd, H-2-Kb, H-2-Kd, H-2Kk, H-2Ld. Rhesus macaque alleles: Mamu-A01, Mamu-A02, Mamu-A11, Mamu-A2:201, Mamu-A26:01, Mamu-B01, Mamu-B03, Mamu-B08, Mamu-B17. Chimpanzees alleles: Patr-A01:01, Patr-A03:01, Patr-A04:01, Patr-A07:01, Patr-A09:01, Patr-B01:01, Patr-B24:01. Nanomolar affinity values for wild-type reference and mutation-derived peptides were ranked, with those less that 50 nM referred to “tight binders” and those less that 500 nM as “loose binders.” These thresholds were previously adopted by Istail et al.31 Data Analysis. Linear regression was analyzed with Graphpad Prism (GraphPad Software Inc.). Figures were generated with Graphpad Prism and Tableau Professional (Tableau Software).
Disclosure of Potential Conflicts of Interest

No potential conflicts of interest were disclosed.
Acknowledgments

Z.S. is supported by the Breast Cancer Research Foundation.
Supplemental Materials

Supplemental materials may be found here: w w w.landesbioscience.com /journa ls /oncoimmunolog y/ article/21511
9. Toubaji A, Achtar M, Provenzano M, Herrin VE, Behrens R, Hamilton M, et al. Pilot study of mutant ras peptide-based vaccine as an adjuvant treatment in pancreatic and colorectal cancers. Cancer Immunol Immunother 2008; 57:1413-20; PMID:18297281; http://dx.doi.org/10.1007/s00262-008-0477-6. 10. Meyer RG, Korn S, Micke P, Becker K, Huber C, Wölfel T, et al. An open-label, prospective phase I/ II study evaluating the immunogenicity and safety of a ras peptide vaccine plus GM-CSF in patients with non-small cell lung cancer. Lung Cancer 2007; 58:8894; PMID:17599645; http://dx.doi.org/10.1016/j. lungcan.2007.05.003. 11. Fenton RG, Keller CJ, Hanna N, Taub DD. Induction of T-cell immunity against Ras oncoproteins by soluble protein or Ras-expressing Escherichia coli. J Natl Cancer Inst 1995; 87:1853-61; PMID:7494229; http://dx.doi.org/10.1093/jnci/87.24.1853.

2.

6.

3.

7.

4.

8.

Robbins PF, El-Gamil M, Li YF, Kawakami Y, Loftus D, Appella E, et al. A mutated beta-catenin gene encodes a melanoma-specific antigen recognized by tumor infiltrating lymphocytes. J Exp Med 1996; 183:1185-92; PMID:8642260; http://dx.doi. org/10.1084/jem.183.3.1185. Greiner J, Ono Y, Hofmann S, Schmitt A, Mehring E, Götz M, et al. Mutated regions of nucleophosmin 1 (NPM1) elicit both CD4+ and CD8+ T cell responses in patients with acute myeloid leukemia. Blood 2012; Epub ahead of print; PMID:22592607; http://dx.doi. org/10.1182/blood-2011-11-394395. Matsushita H, Vesely MD, Koboldt DC, Rickert CG, Uppaluri R, Magrini VJ, et al. Cancer exome analysis reveals a T-cell-dependent mechanism of cancer immunoediting. Nature 2012; 482:400-4; PMID:22318521; http://dx.doi.org/10.1038/nature10755. Castle JC, Kreiter S, Diekmann J, Löwer M, van de Roemer N, de Graaf J, et al. Exploiting the mutanome for tumor vaccination. Cancer Res 2012; 72:1081-91; PMID:22237626; http://dx.doi.org/10.1158/00085472.CAN-11-3722.

1288

OncoImmunology

Volume 1 Issue 8

12. Gjertsen MK, Bakka A, Breivik J, Saeterdal I, GeddeDahl T 3rd, Stokke KT, et al. Ex vivo ras peptide vaccination in patients with advanced pancreatic cancer: results of a phase I/II study. Int J Cancer 1996; 65:450-3; PMID:8621226; http://dx.doi.org/10.1002/ (SICI)1097-0215(19960208)65:4<450::AIDIJC10>3.0.CO;2-E. 13. Gjertsen MK, Bakka A, Breivik J, Saeterdal I, Solheim BG, Søreide O, et al. Vaccination with mutant ras peptides and induction of T-cell responsiveness in pancreatic carcinoma patients carrying the corresponding RAS mutation. Lancet 1995; 346:1399-400; PMID:7475823; http://dx.doi.org/10.1016/S01406736(95)92408-6. 14. Gjertsen MK, Saeterdal I, Thorsby E, Gaudernack G. Characterisation of immune responses in pancreatic carcinoma patients after mutant p21 ras peptide vaccination. Br J Cancer 1996; 74:1828-33; PMID:8956801; http://dx.doi.org/10.1038/bjc.1996.638. 15. Segal NH, Parsons DW, Peggs KS, Velculescu V, Kinzler KW, Vogelstein B, et al. Epitope landscape in breast and colorectal cancer. Cancer Res 2008; 68:889-92; PMID:18245491; http://dx.doi.org/10.1158/00085472.CAN-07-3095. 16. Warren RL, Holt RA. A census of predicted mutational epitopes suitable for immunologic cancer control. Hum Immunol 2010; 71:245-54; PMID:20035814; http:// dx.doi.org/10.1016/j.humimm.2009.12.007. 17. Keogh E, Fikes J, Southwood S, Celis E, Chesnut R, Sette A. Identification of new epitopes from four different tumor-associated antigens: recognition of naturally processed epitopes correlates with HLAA*0201-binding affinity. J Immunol 2001; 167:78796; PMID:11441084. 18. Lundegaard C, Lund O, Nielsen M. Prediction of epitopes using neural network based methods. J Immunol Methods 2011; 374:26-34; PMID:21047511; http:// dx.doi.org/10.1016/j.jim.2010.10.011.

19. Zhang GL, Ansari HR, Bradley P, Cawley GC, Hertz T, Hu X, et al. Machine learning competition in immunology - Prediction of HLA class I binding peptides. J Immunol Methods 2011; 374:1-4; PMID:21986107; http://dx.doi.org/10.1016/j.jim.2011.09.010. 20. Forbes SA, Bindal N, Bamford S, Cole C, Kok CY, Beare D, et al. COSMIC: mining complete cancer genomes in the Catalogue of Somatic Mutations in Cancer. Nucleic Acids Res 2011; 39(Database issue):D945-50; PMID:20952405; http://dx.doi. org/10.1093/nar/gkq929. 21. van der Burg SH, Visseren MJ, Brandt RM, Kast WM, Melief CJ. Immunogenicity of peptides bound to MHC class I molecules depends on the MHC-peptide complex stability. J Immunol 1996; 156:3308-14; PMID:8617954. 22. Sanderson K, Scotland R, Lee P, Liu D, Groshen S, Snively J, et al. Autoimmunity in a phase I trial of a fully human anti-cytotoxic T-lymphocyte antigen-4 monoclonal antibody with multiple melanoma peptides and Montanide ISA 51 for patients with resected stages III and IV melanoma. J Clin Oncol 2005; 23:741-50; PMID:15613700; http://dx.doi. org/10.1200/JCO.2005.01.128. 23. Grbovic OM, Basso AD, Sawai A, Ye Q, Friedlander P, Solit D, et al. V600E B-Raf requires the Hsp90 chaperone for stability and is degraded in response to Hsp90 inhibitors. Proc Natl Acad Sci U S A 2006; 103:5762; PMID:16371460; http://dx.doi.org/10.1073/ pnas.0609973103. 24. Busch DH, Pamer EG. MHC class I/peptide stability: implications for immunodominance, in vitro proliferation, and diversity of responding CTL. J Immunol 1998; 160:4441-8; PMID:9574549. 25. Amos SM, Duong CP, Westwood JA, Ritchie DS, Junghans RP, Darcy PK, et al. Autoimmunity associated with immunotherapy of cancer. Blood 2011; 118:499-509; PMID:21531979; http://dx.doi. org/10.1182/blood-2011-01-325266.

26. Gatenby RA, Brown J, Vincent T. Lessons from applied ecology: cancer control using an evolutionary double bind. Cancer Res 2009; 69:7499-502; PMID:19752088; http://dx.doi.org/10.1158/00085472.CAN-09-1354. 27. Klechevsky E, Gallegos M, Denkberg G, Palucka K, Banchereau J, Cohen C, et al. Antitumor activity of immunotoxins with T-cell receptor-like specificity against human melanoma xenografts. Cancer Res 2008; 68:6360-7; PMID:18676861; http://dx.doi. org/10.1158/0008-5472.CAN-08-0928. 28. Michaeli Y, Denkberg G, Sinik K, Lantzy L, ChihSheng C, Beauverd C, et al. Expression hierarchy of T cell epitopes from melanoma differentiation antigens: unexpected high level presentation of tyrosinase-HLA-A2 Complexes revealed by peptide-specific, MHC-restricted, TCR-like antibodies. J Immunol 2009; 182:6328-41; PMID:19414786; http://dx.doi. org/10.4049/jimmunol.0801898. 29. Denkberg G, Lev A, Eisenbach L, Benhar I, Reiter Y. Selective targeting of melanoma and APCs using a recombinant antibody with TCR-like specificity directed toward a melanoma differentiation antigen. J Immunol 2003; 171:2197-207; PMID:12928363. 30. Lundegaard C, Hoof I, Lund O, Nielsen M. State of the art and challenges in sequence based T-cell epitope prediction. Immunome Res 2010; 6(Suppl 2):S3; PMID:21067545; http://dx.doi.org/10.1186/17457580-6-S2-S3. 31. Istrail S, Florea L, Halldórsson BV, Kohlbacher O, Schwartz RS, Yap VB, et al. Comparative immunopeptidomics of humans and their pathogens. Proc Natl Acad Sci U S A 2004; 101:13268-72; PMID:15326311; http://dx.doi.org/10.1073/pnas.0404740101.

www.landesbioscience.com

OncoImmunology

1289