Biochemical and functional characterization of mutant KRAS epitopes validates this oncoprotein for immunological targeting

Bear, Adham S.; Blanchard, Tatiana; Cesare, Joseph; Ford, Michael J.; Richman, Lee P.; Xu, Chong; Baroja, Miren L.; McCuaig, Sarah; Costeas, Christina; Gabunia, Khatuna; Scholler, John; Posey, Avery D.; O’Hara, Mark H.; Smole, Anze; Powell, Daniel J.; Garcia, Benjamin A.; Vonderheide, Robert H.; Linette, Gerald P.; Carreno, Beatriz M.

doi:10.1038/s41467-021-24562-2

Download PDF

Article
Open access
Published: 16 July 2021

Biochemical and functional characterization of mutant KRAS epitopes validates this oncoprotein for immunological targeting

Nature Communications volume 12, Article number: 4365 (2021) Cite this article

21k Accesses
48 Citations
58 Altmetric
Metrics details

Subjects

Abstract

Activating RAS missense mutations are among the most prevalent genomic alterations observed in human cancers and drive oncogenesis in the three most lethal tumor types. Emerging evidence suggests mutant KRAS (mKRAS) may be targeted immunologically, but mKRAS epitopes remain poorly defined. Here we employ a multi-omics approach to characterize HLA class I-restricted mKRAS epitopes. We provide proteomic evidence of mKRAS epitope processing and presentation by high prevalence HLA class I alleles. Select epitopes are immunogenic enabling mKRAS-specific TCRαβ isolation. TCR transfer to primary CD8⁺ T cells confers cytotoxicity against mKRAS tumor cell lines independent of histologic origin, and the kinetics of lytic activity correlates with mKRAS peptide-HLA class I complex abundance. Adoptive transfer of mKRAS-TCR engineered CD8⁺ T cells leads to tumor eradication in a xenograft model of metastatic lung cancer. This study validates mKRAS peptides as bona fide epitopes facilitating the development of immune therapies targeting this oncoprotein.

FOXO1 is a master regulator of memory programming in CAR T cells

Article Open access 10 April 2024

Alexander E. Doan, Katherine P. Mueller, … Evan W. Weber

Personalized neoantigen vaccine and pembrolizumab in advanced hepatocellular carcinoma: a phase 1/2 trial

Article Open access 07 April 2024

Mark Yarchoan, Edward J. Gane, … Niranjan Y. Sardesai

The present and future of bispecific antibodies for cancer therapy

Article 06 March 2024

Christian Klein, Ulrich Brinkmann, … Roland E. Kontermann

Introduction

The RAS family of GTPases (KRAS, NRAS, HRAS) is mutated in approximately 20% of all human malignancies^1,2. The vast majority of RAS genomic alterations are the result of missense mutations at codon positions G12, G13, or Q61 that involve the RAS protein GTP-binding domain, resulting in the activation of downstream effector substrates ERK and PI3-K leading to dysregulated cell growth and survival³. Among cancers in which KRAS mutations predominate, including pancreatic ductal adenocarcinomas (PDA), lung adenocarcinomas (LAC), and colorectal carcinomas (CRC), over 75% of amino acid substitutions occur at the G12 codon position⁴. The high frequency of G12 amino acid substitutions makes this codon position an ideal drug target; however, no small molecule inhibitors of G12 variants have been successfully developed aside from KRAS G12C that has only recently demonstrated clinical promise^5,6.

Somatic gene mutations within cancer cells may be translated into peptides that are processed and presented on the surface of tumor cells^7,8. These mutated peptides can serve as foreign epitopes, or neoantigens, that may be recognized by αβ T cells of the host immune system. Neoantigen-specific T cell responses have been well documented in patients for whom immune checkpoint blockade therapy has been successful, and they are believed to be the key mediators of anti-tumor activity⁹. Neoantigen-specific T cells can be isolated from the peripheral blood or tumor tissue of antigen-exposed cancer patients¹⁰, induced and expanded from the peripheral blood of cancer patients following neoantigen vaccination¹¹, and generated in vitro utilizing the naive T-cell receptor (TCR) repertoire of healthy donors¹². These observations have garnered interest in the development of neoantigen-targeted cancer vaccines and adoptive T cell therapies.

Neoantigen-based treatment strategies are inherently highly personalized as somatic tumor mutations are rarely shared between patients^13,14. The high prevalence and conserved mutational profile of KRAS affords a unique opportunity to develop a neoantigen-targeted therapy with broad generalizability. Mutant KRAS (mKRAS) has been previously studied as a potential target of cancer vaccines, and clinical studies have demonstrated the successful generation of CD4⁺ and CD8⁺ αβ T cell responses with reactivity against allogeneic or autologous mKRAS tumor cell lines^15,16,17,18. More recently, mKRAS-specific T cells have been isolated and characterized from the peripheral blood of patients with mKRAS epithelial cancers^19,20, and they have been induced in vaccinated HLA-transgenic mice²¹. The therapeutic potential of targeting mKRAS as a cancer neoantigen was highlighted in a case report demonstrating clinical benefit in a patient with KRAS G12D metastatic CRC following the adoptive transfer of KRAS G12D-specific T cells restricted to HLA-C*08:02²². However, mKRAS as an immunological target remains poorly characterized, and there is a dearth of evidence regarding tumor cell processing and presentation of mKRAS-derived epitopes to guide the development of targeted immunotherapies.

In this work, we aim to validate mKRAS as an immunological target. Utilizing computational epitope prediction, biochemical assays, and proteomic analysis, we predict and identify high-affinity and/or high-stability mKRAS G12 peptides to HLA-A*03:01, HLA-A*11:01, and HLA-B*07:02, and we confirm these epitopes as constituents of their corresponding HLA class I ligandome. The induction of mKRAS-reactive T cells from healthy donors by select peptides both confirms peptide immunogenicity and enables the isolation of mKRAS-specific TCRs directed against the following three epitopes: G12V/HLA-A*03:01, G12V/HLA-A*11:01, and G12R/HLA-B*07:02. Expression of these mKRAS-specific TCRαβ pairs in a NFAT-inducible reporter T cell line or TCRαβ^null primary CD8⁺ T cells serves as sensitive probes to characterize the presentation of mKRAS G12 epitopes by tumor cell lines. TCR-redirected T cells display cytotoxic function against mKRAS tumor cell lines of various histologies without reactivity to wild-type (WT) KRAS peptides, thereby authenticating mKRAS G12V and G12R peptide ligands as bona fide neoantigens. Finally, the adoptive transfer of TCR redirected CD8⁺ human T cells specific for G12V/HLA-A*03:01 or G12V/HLA-A*11:01 leads to tumor eradication and prolongs survival in a mouse xenograft model of metastatic lung cancer. These findings provide further validation and strengthen the nomination of mKRAS as an immunological target.

Results

Identification and characterization of mKRAS neoepitopes

Although predictive peptide-HLA (p-HLA) binding affinity and stability are important criteria for neoantigen nomination^23,24,25, only a small percentage of the human proteome is presented by HLA molecules²⁶. Therefore, only a fraction of candidate neoantigens are constituents of the tumor HLA ligandome²⁷, and predicting high-affinity/stability p-HLA complexes is not indicative of intracellular processing and presentation²⁸. To improve the identification of bona fide mKRAS neoepitopes, we employed a comprehensive multi-omic approach consisting of (1) computational epitope prediction, (2) biochemical assessment of p-HLA binding / stability, (3) proteomic validation of antigen processing and presentation, and (4) characterization of immunogenicity (Supplementary Fig. 1).

We utilized the neoantigen prediction tool antigen.garnish²⁴ to perform a comprehensive identification of candidate HLA class I allele-restricted mKRAS G12 epitopes with recurrent amino acid substitutions (G12C, G12D, G12R, G12V) embedded within the 24 amino acid sequence derived from the WT KRAS protein. Epitopes derived from these 24-mer amino acid sequences were evaluated for predicted binding to HLA class I alleles to nominate HLA-A, HLA-B, and HLA-C nonamer (9-mer) and/or decamer (10-mer) mKRAS peptides (Supplementary Data 1, candidates <500 nM affinity). Predicted epitopes restricted to HLA-A*02:01, HLA-A*03:01, HLA-A*11:01, and HLA-B*07:02 were selected for further study given the high prevalence of these HLA class I alleles within the USA population (Fig. 1a). Peptides encoding these candidate epitopes are designated by their KRAS protein amino acid positions followed by the G12 amino acid substitution (Fig. 1b). The predicted HLA class I affinities and p-HLA complex stabilities of candidate epitopes are presented in Supplementary Tables 1 and 2.

**Fig. 1: Identification, characterization, and validation of candidate mKRAS epitopes.**

We experimentally assessed p-HLA binding affinity and complex stability of candidate peptides as these two properties improve the performance of T cell epitope identification²⁹. Furthermore, p-HLA complex stability positively correlates with immunogenicity¹². To evaluate p-HLA binding affinity, we employed competitive peptide-binding assays using a FITC-labeled reference peptide and purified soluble HLA class I molecules to quantitate the binding capacity of non-labeled mKRAS peptide candidates. Peptide affinities are reported as concentrations of unlabeled peptide that result in 50% binding inhibition of FITC-labeled reference peptide (Supplementary Fig. 2 and Supplementary Table 1)³⁰. To evaluate p-HLA complex stability, dissociation assays were performed where the stability of p-HLA complexes is expressed as the half-life of beta-2-microglobulin dissociation from the HLA heavy chain (Supplementary Table 2)²⁹. Both assays incorporated validated HLA class I-restricted T cell epitopes as positive controls (Supplementary Table 3). Affinity binding values determined in competitive assays are reported as log₁₀[IC₅₀(nM)] and IC₅₀(nM) and are defined according to a previously reported scale detailed in Supplementary Tables 1 and 3. For HLA-A*02:01 epitopes, only peptide 5-14V exhibited strong class I binding affinity (log₁₀[IC₅₀] = 3.469 nM) (Fig. 1c). However, all p-HLA-A*02:01 class I complexes exhibited short half-lives (t_1/2 < 1 h) compared to control peptides (t_1/2 1.3–10.1 h, Fig. 1d). For HLA-A*03:01 epitopes, peptides 7-16R, 8-16V, and 7-16V exhibited strong binding affinity (log₁₀[IC₅₀] ≤ 3.589 nM) (Fig. 1c) and formed p-HLA complexes with equivalent stability to positive controls (t_1/2 > 1 h). All 9-mer (8-16C-V) and 10-mer (7-16C-V) peptides bound with strong affinity to HLA-A*11:01 (log₁₀[IC₅₀] ≤ 3.121 nM) (Fig. 1c), and 10-mer peptides formed more stable p-HLA-A*11:01 complexes compared to 9-mer peptides (t_1/2 6–12 h vs 0.9–9.2 h) (Fig. 1d). For HLA-B*07:02 epitopes, peptides 10-19C, 10-19D, and 10-19R bound with high affinity (log₁₀[IC₅₀] ≤ 3.460 nM) (Fig. 1c), but only 10-19C and 10-19R exhibited p-HLA complex stability equivalent to positive control peptides (t_1/2 > 1 h) (Fig. 1d). A comparison of predicted and experimentally determined p-HLA affinities showed a strong correlation for 8-16C-V and 7-16C-V peptides restricted by HLA-A*03:01 and HLA-A*11:01 (R² > 0.9). However, no correlation was observed for either 5-14C-V peptides restricted by HLA-A*02:01 (R² = 0.1033) or 10-19C-V peptides restricted by HLA-B*07:02 (R² = 0.6945) (Supplementary Fig. 3 and Supplementary Table 1). A similar analysis of predicted and experimentally determined p-HLA stabilities (Supplementary Table 2) showed no correlation among p-HLA complexes studied. Experimental stabilities for all mKRAS/HLA-A*11:01 complexes yielded higher binding affinities than were predicted. In summary, we experimentally identified 8 mKRAS G12 epitopes restricted to HLA-A*03:01 (n = 3), HLA-A*11:01 (n = 3), or HLA-B*07:02 (n = 2) that display high affinity and formed high-stability p-HLA complexes equivalent to positive control peptides.

Validation of mKRAS antigen processing and presentation

To validate epitope processing and presentation, we performed targeted mass spectrometry (MS) on peptides isolated from single HLA class I allele/tandem minigene construct (TMC)-expressing cell lines. HLA class I negative cell lines were engineered to express a GFP-tagged HLA class I/β2-microglobulin single-chain dimer (HLA-SCD) construct and a mCherry-tagged TMC (Supplementary Fig. 4a). Use of monoallelic cells lines eliminates the ambiguity that arises from the coexpression of multiple HLA class I alleles³¹. The TMC consisted of nine minigenes encoding WT and mKRAS long 25-mer peptides (G12C, G12D, G12R, G12V) as well as validated viral epitopes restricted by HLA-A*02:01 (NLV), HLA-A*03:01 (ILR), HLA-A*11:01 (IVT), and HLA-B*07:02 (TPR) as positive controls (Supplementary Fig. 4b and Supplementary Table 3). TMC-expressing monoallelic cell lines were phenotyped for reporter and HLA class I expression (Supplementary Fig. 4c–e) then subjected to HLA class I immunoprecipitation and peptide elution. Peptide sequence identities were determined by quadrupole—orbitrap tandem MS using data-dependent acquisition (DDA) and targeted parallel reaction monitoring (PRM)³². Synthetic mKRAS peptides served as “beacons” to zero-in on TMC-encoded peptides that were naturally loaded onto HLA class I molecules. MS/MS fragmentation pattern comparison of eluted and synthetic peptides ensured accurate assignment of AA sequence identity of mKRAS epitopes (Fig. 1e and Supplementary Figs 5–7). In this manner, peptides 7-16D, 8-16V, 7-16V, and 8-16R were eluted in the context of HLA-A*03:01 (Supplementary Fig. 5), and peptides 8-16D, 7-16D, 8-16V, 7-16V, and 8-16R were eluted in the context HLA-A*11:01 (Supplementary Fig. 6). Among these peptides, only 8-16V and 7-16V form p-HLA complexes with HLA-A*03:01 and HLA-A*11:01 with comparable stability to positive controls. Lastly, peptides 10-19D and 10-19R were also identified in the context of HLA-B*07:02 (Supplementary Fig. 7). To further support the identity assignments between eluted and synthetic peptides, we determined precursor ion similarity via isotopic dot product (idotp, value >0.9) and dot product (dotp, value >0.5) measurements³³. No mKRAS G12 peptides presented by HLA-A*02:01 were observed despite the detection of the validated NLV viral epitope (Supplementary Fig. 8). In addition, no cysteine-containing epitopes were detected. Free cysteines are prone to post-translational modifications and are underrepresented in MS-based HLA-peptide datasets^34,35. Future analysis accounting for post-translation modifications will be performed to address the processing of mKRAS cysteine-containing peptides. Importantly, neither 8-16C nor 7-16C in the context of HLA-A*03:01 was observed by Wang et al.³⁶ suggesting these peptides are not part of this HLA class I peptidome. In sum, we identified 11 processed and presented mKRAS epitopes restricted to HLA-A*03:01 (n = 4), HLA-A*11:01 (n = 5), or HLA-B*07:02 (n = 2) class I alleles (Fig. 1e).

mKRAS peptide immunogenicity and mKRAS-specific TCRs isolation

To evaluate the immunogenicity of selected mKRAS peptides, CD8⁺ T cells isolated from healthy donors were cocultured with autologous monocyte-derived mature dendritic cells (mDC) pulsed with candidate epitopes. Following two stimulations, cultures were screened for CD8⁺ T cell responses by IFN-γ ELISPOT assay (Fig. 2a). No responses were observed to peptides 5-14C-V among HLA-A*02:01⁺ donors. Responses were observed to peptide 7-16V in HLA-A*03:01⁺ (n = 1/4) and HLA-A*11:01⁺ (n = 4/4) donors, and responses to 7-16C (n = 4/4) and 7-16D (n = 1/4) were also noted in HLA-A*11:01⁺ donors. Responses to peptide 10-19R (n = 2/5) in HLA-B*07:02⁺ donors were also observed.

**Fig. 2: Assessment of mKRAS neoantigen immunogenicity and identification of mKRAS-specific TCR sequences.**

Given the high sensitivity of TCRs for p-HLA complexes, we sought to isolate TCRs directed against 7-16V and 10-19R epitopes in order to characterize their expression by tumor cell lines and validate these epitopes as neoantigens. To isolate TCRα/β pairs, cell cultures containing CD8⁺ T cell responses to 7-16V/HLA-A*03:01, 7-16V/HLA-A*11:01, and 10-19R/HLA-B*07:02 were expanded using peptide-pulsed artificial antigen presenting cells (APC). Antigen specificity was confirmed by p-HLA multimer staining (Fig. 2b–d), and multimer⁺ cells were flow cytometrically sorted to >99% purity. TCRα/β sequences were determined by next-generation DNA and RNA sequencing³⁷. A total of five TCRα/β pairs were identified, three of which were selected for further studies: TCRA3V (7-16V/HLA-A*03:01), TCRA11V (7-16V/HLA-A*11:01), and TCRB7R (10-19R/HLA-B*07:02) (Table 1).

Table 1 Summary table of mKRAS-TCRs.

Full size table

We initially characterized the expression and function of these mKRAS-TCRs using a cell reporter system developed in our laboratory designated J^ASP90 reporter cells. The J^ASP90 reporter cells are derived from CD8⁺/TCRαβ null Jurkat E6.1 cells engineered to express NFAT-inducible eGFP as a readout for TCR signaling. mKRAS-TCRα/β chain pairs were cloned into lentiviral vectors and expressed in J^ASP90 reporter cells. Custom p-HLA multimer staining was performed to validate TCRA3V, TCRA11V, and TCRB7R expression as well as confirm antigen specificity and HLA restriction (Fig. 3a).

**Fig. 3: Validation and functional characterization of mKRAS-specific TCRs.**

To further assess antigen specificity and measure TCR avidity, TCR-engineered J^ASP90 reporter cells were cocultured with HLA class I matched monoallelic K562 cells pulsed with WT or mKRAS synthetic peptides. TCRA3V, TCRA11V, and TCRB7R-expressing J^ASP90 reporter cells demonstrated specific reactivity to mKRAS cognate epitopes without reactivity to WT KRAS peptide-pulsed or unpulsed monoallelic K562 cells (Fig. 3b). Using titrated peptide concentrations, the functional avidities of each TCR were determined based on the mean concentration required to achieve 50% NFAT activation (Fig. 3b). The avidities of TCRA3V, TCRA11V, and TCRB7R were determined to be 1.64, 0.15, and 3.38 nM, respectively. Furthermore, TCRA3V and TCRB7R cells exhibited no cross-reactivity to other mKRAS G12 epitopes, while TCRA11V cells demonstrated cross-reactivity to peptide 7-16C, albeit with decreased avidity as compared to the cognate peptide. These TCRs constituted high-affinity probes for the characterization of naturally processed and presented mKRAS epitopes by human tumor cells.

mKRAS-specific TCRs confer lytic activity against human tumor cells

To probe for p-HLA complexes on tumor cell lines, the anti-tumor activity of mKRAS TCRs was evaluated using engineered TCRαβ^null primary CD8⁺ T cells. To this end, healthy donor primary CD8⁺ T cells were gene edited using CRISPR/Cas9 and single guide RNAs (sgRNAs) to knock out endogenous TCRα and β expression (Supplementary Fig. 9a–c) to eliminate TCRαβ mispairing and enhance receptor cell surface expression^38,39. Engineered TCR expression was assessed by cell surface CD3 and TCRαβ staining (Supplementary Fig. 9e). Secondary expansion of mKRAS-TCR T cells using peptide-pulsed artificial APCs led to enrichment of antigen-specific CD8⁺ T cells (Supplementary Fig. 9f), yielding populations with high (60–90%) expression of TCRA3V, TCRA11V, or TCRB7R as assessed by p-HLA multimer staining (Fig. 4a). In 4 h ⁵¹Cr-release assays, TCRA3V, TCRA11V, and TCRB7R-expressing CD8⁺ T cells kill HLA class I matched K562 cells pulsed with exogenous cognate mKRAS peptide but not WT KRAS peptide (Fig. 4b). In addition, recognition of processed and presented mKRAS antigen by TCRA3V, TCRA11V, and TCRB7R-expressing CD8⁺ T cells was demonstrated using TMC-expressing K562 cells as targets (Fig. 4b). Lack of recognition of processed and presented WT KRAS was demonstrated using a WT KRAS-expressing cell line (Fig. 4b, c). Altogether, peptide-pulsing experiments validate antigen specificity while TMC experiments, provide support for proteomic results (Fig. 1e), and validate presentation of cell surface mKRAS p-HLA complexes.

**Fig. 4: mKRAS TCRαβ gene transfer confers HLA-restricted lytic activity against human tumor cells.**

In order to evaluate the recognition of p-HLA complexes on the surface of cancer cells by TCRs, we characterized the expression of HLA class I alleles in a panel of mKRAS G12 cancer cell lines of various tissue origins. Cell lines with reported KRAS WT (TPM 23), G12V (TPM 56-308), and G12R (TPM 16-130) expression, as well as cell lines expressing target HLA class I alleles but alternative KRAS mutations were selected for further study (Supplementary Table 4). Aside from G12V/HLA-A*03:01, no cell lines expressing endogenous G12V/HLA-A*11:01 or G12R/HLA-B*07:02 were identified. When applicable, KRAS G12V⁺ cell lines were modified to express either HLA-A*03:01 or HLA-A*11:01, and KRAS G12R⁺ lines were modified to express HLA-B*07:02. As shown in Fig. 4c, d, TCRA3V and TCRA11V-expressing CD8⁺ T cells exhibited specific cytotoxicity, irrespective of tissue of origin, against all HLA-A*03:01 or HLA-A*11:01-matched KRAS G12V⁺ cancer cell lines, respectively. Of note, neither TCR conferred cytotoxic activity against HLA matched BxPC-3 tumor cells, a KRAS WT pancreatic adenocarcinoma cell line, validating the specificity of each TCR for the cognate G12V peptide. These observations were replicated and expanded upon using the J^ASP90 reporter cell system. Upon co-culture with mKRAS/HLA-matched cell lines, TCRA3V-, TCRA11V-, and TCRB7R-J^ASP90 cells demonstrated NFAT-induced GFP expression indicative of TCR sensing of p-HLA complexes presented on the surface of tumor cells (Supplementary Fig. 10a–c). No NFAT-induced GFP expression was observed in J^ASP90 reporter cells cultured with HLA matched KRAS WT-expressing BxPC-3 cells. The graded pattern of NFAT activation induced in J^ASP90 reporter cells may suggest varied surface expression of p-HLA complexes and/or differing patterns of coreceptor expression by cancer cells. Altogether, these results strongly support the hypothesis that mKRAS epitopes are processed and presented by tumor cells in the context of HLA-A*03:01, HLA-A*11:01 and HLA-B*07:02 at sufficient levels to activate recognition by mKRAS-TCR-engineered T cells.

Quantitation and recognition of mKRAS peptides presented by human tumor cells

To gain further insights into the processing and presentation of 8-16V and 7-16V peptides on the surface of cancer cell lines, we performed targeted MS and absolute peptide quantitation to enumerate p-HLA complexes on the surface of CORL23 lung tumor cells engineered to express HLA-A*03:01 (CORL23-A3) or HLA-A*11:01 (CORL23-A11)⁴⁰. CORL23 was chosen for further studies given its ability to recapitulate a model of metastatic lung cancer (see below)⁴¹. In the context of HLA-A*03:01, we detected 8 complexes of 8-16V and 27 complexes of 7-16V per cell as compared to 43 complexes of 8-16V and 78 complexes of 7-16V per cell in the context of HLA-A*11:01 (Fig. 5a and Supplementary Fig. 11). This result confirms processing and presentation data obtained using TMC-engineered HLA class I monoallelic cell lines (Supplementary Figs 5 and 6) and suggests that KRAS G12V⁺ tumors present both 9-mer (8-16V) and 10-mer (7-16V) epitopes in the context of HLA-A*03:01 and HLA-A*11:01. Although these cells were genetically engineered to express the restricting HLA class I allele, HLA-A*03:01 and HLA-A*11:01 cell surface levels were noted to be within the range of expression comparable to other non-engineered HLA-A*03:01⁺ and HLA-A*11:01⁺ mKRAS tumor cells lines (Supplementary Fig. 12a, b). Importantly, enumerated p-HLA complexes are within the ranges reported for other mKRAS/HLA complexes expressed by non-engineered cancer cell lines³⁶.

**Fig. 5: Quantitation and recognition of KRAS G12V peptide ligands presented by human tumor cell lines.**

The cytotoxic activity of TCRA3V- and TCRA11V-engineered CD8⁺ T cells against CORL23 tumor cell lines was then characterized via live cell imaging and cellular impedance. These complimentary measurements of cell death allow for visualization of the loss of GFP-labeled tumor cells and the accumulation of Annexin V-CF594 dye (Supplementary Fig. 13 and Supplementary Movie 1), as well as cell shrinkage and fragmentation measurements over a 6-day period. TCRA3V cells promoted rapid apoptotic cell death of CORL23-A3 tumor cells at a 10:1 E:T ratio compared to control T cells, with slower killing kinetics at lower E:T ratios of 3:1 and 1:1 (Fig. 5b). TCRA11V cells promoted rapid cell death of CORL23-A11 cells at all E:T ratios (Fig. 5c). The kinetics to reach 50% cytolysis (KT50) was significantly faster for TCRA11V cells compared to TCRA3V cells at all E:T ratios (Fig. 5d), a finding consistent with the higher avidity of TCRA11V (Fig. 3b) and the higher number of 7-16V complexes detected in the context of HLA-A*11:01 expressed on the surface of CORL23-A11 cells.

We next evaluated the sensitivity of mKRAS TCRs to recognize p-HLA complexes in a non-engineered cell line. We identified the LAC cell line NCI-H441 as being both KRAS G12V⁺ and HLA-A*03:01⁺ (Supplementary Table 4), highlighting it as a potential target of TCRA3V. We confirmed this cell line retained HLA class I expression and expressed low levels of HLA-A*03:01 (Supplementary Fig. 12a). In both 4 h ⁵¹Cr-release and cellular impedance assays, TCRA3V cells exhibited specific lysis against parental (non-engineered) NCI-H441 cells, while TCRA11V cells, which recognizes the identical 7-16V epitope but restricted by HLA-A*11:01, exhibited no cytotoxic activity (Fig. 5e, f). Together, we conclude that TCRA3V and TCRA11V can “sense” low numbers of p-HLA complexes per cell and represent valuable reagents to immunologically target mKRAS.

In vivo anti-tumor activity of mKRAS-specific TCR+ CD8+ T cells

CORL23-A3 and CORL23-A11 tumor cells engineered to express click beetle red (CBR) luciferase were infused into NOD/scid/γcnull (NSG) mice via tail vein injection, and the establishment of pulmonary tumors was confirmed by bioluminescence imaging. Four days post tumor inoculation, mice were treated with TCRA3V- or TCRA11V-engineered CD8⁺ T cells. As controls, tumor-bearing mice received either no treatment (Mock) or TCRαβ^null (TCR KO) CD8⁺ T cells. Mice treated with TCRA3V- or TCRA11V- engineered CD8⁺ T experienced complete eradication of CORL23-A3 (Fig. 6a, c) or CORL23-A11 (Fig. 6b, d) tumors, respectively, relative to control groups (p < 0.01). The treatment effect was long-lasting with no evidence of tumor outgrowth by day 39 in mice treated with mKRAS-TCR therapy, while control mice experienced persistent tumor growth. Consequently, mice treated with mKRAS-TCR therapy experienced prolonged survival relative to control groups (p < 0.0002, Fig. 6e, f). Altogether, these results demonstrate that the transfer of T cells engineered with high avidity mKRAS TCRs promotes tumor rejection and prolongs survival.

**Fig. 6: Adoptive transfer of mKRAS-TCR T cells leads to in vivo eradication of KRAS G12V⁺ tumor cells in a xenograft model of metastatic lung cancer.**

Discussion

Epitope discovery is an essential step in designing immunotherapies such as cancer vaccines, bi-specific engagers, and TCR therapies. The ubiquity of mKRAS along with its role as a recurrent clonal driver makes this oncoprotein an ideal target^42,43,44. To identify mKRAS epitopes, we undertook a comprehensive multi-omics approach (e.g., bioinformatic, biochemical, proteomic and immunological) to address the intricate processes involved in the generation of immunogenic peptides. By subjecting each candidate peptide to a series of assays that interrogate critical properties of antigen processing and presentation, we have identified a curated, albeit partial set of mKRAS peptides that represent candidate epitopes. Using this approach, we identified three epitopes (7-16V/HLA-A*03:01, 7-16V/HLA-A*11:01, and 10-19R/HLA-B*07:02) that fulfilled criteria of strong measured p-HLA binding affinity, high complex stability, and proteomically validated antigen processing and presentation. Furthermore, we determined that these epitopes were immunogenic allowing for the identification of mKRAS TCR sequences that were used as sensitive probes to validate epitope presentation by human cancer cell lines. Direct quantification of mKRAS p-HLA complexes on CORL23 cells highlights the sensitivity of mKRAS TCRs to recognize low abundance p-HLA complexes and eliminate mKRAS tumors cells in a xenograft mouse model of metastatic cancer.

Most mKRAS epitopes described herein have not been previously reported as constituents of the HLA peptidome. These epitopes include the following: HLA-A*03:01 – 8-16V, 8-16R, and 7-16V, HLA-A*11:01 – 8-16D, 8-16R, 8-16V, 7-16D, and 7-16V, HLA-B*07:02 – 10-19D and 10-19R. These epitopes exhibit high affinity for the restricting HLA class I allele and relatively high complex stability as each forms p-HLA complexes with t_1/2 > 0.5 h. In agreement with our data, the processing and presentation of the 7-16D peptide in the context of HLA-A*03:01 has also recently been reported by Wang et al.³⁶ despite exhibiting low measured affinity and complex stability. Our study and those of others³⁶ using cell lines expressing mKRAS have failed to identify 5-14 mKRAS peptides as constituents of the HLA-A*02:01 peptidome. Furthermore, biochemical assessment of 5-14C-V peptide-binding and stability properties for HLA-A*02:01 suggest that these ligands fail to form stable p-HLA complexes, and only the 5-14V peptide demonstrated high HLA class I affinity in competitive binding assays. By contrast, Mishto et al.⁴⁵ reported mKRAS linear (5-14V) and spliced (5-6/8-14V, lacking the V at position 7 of the canonical peptide) peptides as products of in vitro proteasome digests and have demonstrated these peptides to have high affinity for HLA-A*02:01. Given the high prevalence of HLA-A*02:01 in the population, further studies are needed to evaluate the discordant results between processing of the 5-14V mKRAS peptide by cell lines and in vitro proteasome digest assays.

To investigate the immunogenicity of mKRAS epitopes, immune responses were “outsourced” from healthy donors as naive TCR repertoires may recognize tumor neoantigens¹². In addition to 7-16V and 10-19R epitopes described, we also identified CD8⁺ T cell responses to 7-16C and 7-16D in HLA-A*11:01⁺ donors. Of these epitopes, only 7-16D- and 7-16V-specific TCRs isolated from immunized transgenic HLA-A*11:01 mice have been previously reported²¹. We did not observe immune responses directed against 5-14C-V in HLA-A*02:01⁺ donors despite previous studies reporting HLA-A2-restricted 5-14D and 5-14V CD8⁺ T cell responses^17,46. Given that these studies were reported over 15 years ago, it is likely high-resolution HLA typing was not performed, and it is possible responding patients expressed less prevalent HLA-A2 sub-alleles (e.g., HLA-A*02:03 or HLA-A*02:06) that may exhibit higher binding affinity to 5-14V and 5-14D peptides relative to HLA-A*02:01. Clearly, this issue related to mKRAS epitopes restricted by HLA-A*02:01 requires further investigation.

The availability of TCRs directed against 7-16V/HLA-A*03:01 (TCRA3V), 7-16V/HLA-A*11:01 (TCRA11V), and 10-19R/HLA-B*07:02 (TCRB7R) permitted their use as sensitive probes demonstrating the presence of these p-HLA complexes on the cell surface of various cancer cell lines from multiple histologies. Our observation that NCI-H441 can be recognized by TCRA3V-engineered CD8⁺ human T cells is noteworthy given this KRAS G12V mutant cell line naturally expresses the HLA-A*03:01 restricting allele. Further insights were gained by the quantification of p-HLA complexes on CORL23 cell lines demonstrating expression of both 9-mer and 10-mer peptides in the context of HLA-A*03:01 and HLA-A*11:01. As TCRA3V and TCRA11V recognize the same 7-16V peptide, we hypothesized differences in TCR avidity as well as p-HLA abundance by CORL23 lines may contribute to the distinct kinetics of in vitro killing. Indeed, higher expression of 7-16V/A*11:01 complexes (78 complexes/cell) relative to 7-16V/A*03:01 (27 complexes/cell) correlates with the more rapid elimination of CORL23-A11 relative to CORL23-A3 tumor cells. Importantly, the adoptive T cell transfer of either TCRA3V or TCRA11V engineered CD8⁺ human cells efficiently eliminated CORL23 tumor cells in vivo and proved curative for a majority of treated mice.

In summary, epitope identification and validation strategies reliant on a contemporary toolbox of bioinformatic, biochemical, genomic, and immunological assays have recently been employed by several groups in the development of personalized melanoma vaccines targeting passenger mutations^11,47,48. Here, we have employed a similar strategy to characterize epitopes presented by high prevalence HLA class I alleles targeting a recurrent clonal driver—mKRAS. Our findings provide further evidence highlighting this oncogenic protein as a clinical target of immune-based therapies.

Methods

Primary cells

Peripheral blood mononuclear cells and purified CD8⁺ T cells were provided by the Human Immunology Core after cell isolation from apheresis products of HLA class I- and class II-typed healthy donors enrolled on the Institutional Review Board-approved research protocol 705906 at the University of Pennsylvania after providing informed consent.

Cell lines

The following cell lines were cultured in RPMI media supplemented with 10% fetal bovine serum (FBS), 2 mM l-glutamine (CORNING 25-005-CI), and 1X penicillin/streptomycin (CORNING 30-002-CI): 721.221 (Fred Hutchinson Cancer Research Center International Histocompatibility Working Group IHW00001), K562 (American Type Culture Collection (ATCC CCL-243)), Jurkat E6-1 (ATCC TIB-152), BxPC-3 (ATCC CRL-1687), CORL23 (Sigma-Aldrich 92031919), NCI-H441 (ATCC HTB-174), SW620 (ATCC CCL-227), YAPC (German Collection of Microorganisms and Cell Cultures (DSMZ) ACC 382), KP-2 (Japanese Collection of Research Bioresources Cell Bank (JCRB) JCRB0181), QGP-1 (JCRB JCRB0183), PSN-1 (ATCC CRL-3211), NCI-H2030 (ATCC CRL-5914), NCI-H358 (ATCC CRL-5807), HuCCT1 (CellBank Australia JCRB0425), and RERF-LC-Ad1 (JCRB JCRB1020). The following cell lines were cultured in DMEM media supplemented with 10% FBS, 2 mM l-glutamine, and 1X penicillin/streptomycin: SK-CO-1 (ATCC HTB-39) and HuP-T3 (DSMZ ACC 259). The following cell lines were cultured in EMEM media supplemented with 10% FBS, l-glutamine, and penicillin/streptomycin: CAL-62 (DSMZ ACC 448) and PANC1 (Laboratory of Dr Gregory Beatty at the University of Pennsylvania).

Antibodies and flow cytometry

The following primary antibodies along with vendor, catalog number, and working dilution used for flow cytometry in this manuscript are as follows: anti-human HLA-A, B, C (W6/32)—APC (BioLegend, 311410, 1:50), anti-human HLA-A*02:01—PE (BD, 558570, 1:50), anti-human HLA-A*03:01—Biotin (One Lambda, BIH0269, 1:20), anti-human HLA-A*11:01—Biotin (One Lambda, BIH0084, 1:20), anti-human HLA-B*07:02—PE (BioLegend, 372403, 1:50), anti-human CD3—FITC (Invitrogen, MHCD0301, 1:100), anti-human CD8—APC (ThermoFisher Scientific, MHCD0805, 1:100), anti-human CD8 (SK1)—APC (BioLegend, 344722, 1:100), anti-human TCRαβ (IP26)—APC (BioLegend, 306717, 1:20), and anti-human TCRαβ (IP26)—Biotin (BioLegend, 306703, 1:20). The following custom p-HLA multimers along with vendor and working dilution information are as follows: VVVGAVGVGK/HLA-A*03:01—PE (Immudex, Custom, 1:100), VVVGAVGVGK / HLA-A*11:01—PE (Immudex, Custom, 1:100), GARGVGKSAL/HLA-B*07:02—PE (Immudex, Custom, 1:100). Cells were stained at the specified antibody dilutions and washed in FACS buffer. Data were acquired with a BD LSRFortessa flow cytometer using BD FACSDiva software (version 8.0.2) and analyzed using FlowJo software (version 10.0.7). ELISPOT assays utilized anti-human IFN-γ 1-D1K (Mabtech, 3420-3-250, 1:100), and anti-human IFN-γ 7-B6-1—Biotin (Mabtech, 3420-6-1000, 1:1000). For rapid expansion of sorted T cells, anti-human CD3 (OKT3) (LEAF, 317326, 30 ng/ml) was used. For HLA class I immunoprecipitation, anti-human HLA-A, B, C (W6/32)—AC (Santa Cruz Biotech, sc-32235 AC, 10 μg per 1 mg protein lysate) was used.

Computational prediction of mKRAS G12 neoepitopes

We utilized the tumor neoantigen prediction and multi-parameter quality analysis tool antigen.garnish²⁴ to identify HLA class I binding peptide ligands of mKRAS G12. The input file consisted of the amino acids 1–24 of WT KRAS as well as the most prevalent mKRAS somatic G12 variants (G12C, G12D, G12R, and G12V) screened against HLA class I alleles with a prevalence >1% of the USA population. Nonamer and decamer peptide ligands with ensemble predicted binding affinities <500 nM to target HLA class I molecules were designated candidate mKRAS binding ligands. HLA prevalence among Caucasians, African Americans, Hispanics, and Asians represented in the USA population were gathered from the Allele Frequency Net Database⁴⁹. WT and mKRAS G12 somatic variant frequencies were collated from The International Cancer Genome Consortium Data Portal utilizing pancreatic adenocarcinoma (Project: PAAD-US, n = 177; PACA-CA, n = 268; PACA-AU, n = 391), colorectal cancer (Project: COAD-US, n = 402; READ-US, n = 321; COCA-CN, n = 321), and LAC (Project: LUAD-US, n = 516) patient cohorts.

Peptides and biochemical peptide-MHC binding assays

Custom peptides were synthesized by New England Peptide (Gardner, MA) to >95% purity. Lyophilized peptides were dissolved in 10% DMSO (Mylan Cryoserv 67457-178-50) in sterile water pH 7.4 and passed through a 0.2 uM Centrex filter (10467013). Peptide binding to HLA class I alleles was evaluated using a fluorescence polarization-based competition assay (Pure Protein, L.L.C.) to determine the half maximal inhibitory concentration (IC₅₀) of test peptides as a measure of immunogenic potential. The affinity scale for this assay is as follows: High = log₁₀[IC₅₀] < 3.7 nM, Intermediate = log₁₀[IC₅₀] = 3.7–4.7 nM, Low = log₁₀[IC₅₀] = 4.7–5.5 nM, and Very Low = log₁₀[IC₅₀] = 5.5–6.0 nM³⁰. p-HLA complex stability measurements were determined using a scintillation proximity assay with association/dissociation curves fitted using GraphPad Prism version 7. In a random sample of 384 peptides, the scintillation proximity assay could discriminate among peptides with affinities <500 nM, stable complexes with half-life between 0.1 and 46 h⁵⁰.

HLA class I genotyping and phenotyping of cell lines

Tumor cell line HLA class I typing data were obtained using the TRON Cell Line Portal⁵¹ when available. HLA class I typing for CORL23 and CAL-62 cell lines was performed by standard molecular typing (Hospital of the University of Pennsylvania, Philadelphia, PA). HLA haplotypes of cell lines are presented in Supplementary Table 4. Tumor cell lines were transduced with lentiviral particles expressing HLA class I/β2-microglobulin single-chain dimer constructs (HLA-SCD, eGFP⁺) to generate mKRAS G12/HLA class I-matched cell lines. HLA class I cell surface expression was assessed by flow cytometry using APC-conjugated anti-human HLA-A, B, C (clone W6/32). HLA-allele-specific antibodies conjugated to PE, APC, or biotin were used to assess target HLA expression.

Lentiviral vector constructs and production

All lentiviral constructs were generated using the third generation lentiviral transfer vector pTRPE-eGFP-T2A-mCherry (generously provided by Dr Michael C. Milone, University of Pennsylvania, Philadelphia, PA). The HLA-A*02:01-SCD was PCR amplified from a HLA-A*02:01-SCD plasmid⁵² and used to replace the mCherry moiety on pTRPE-eGFP-T2A-mCherry vector backbone. All other HLA class I-SCD constructs were generated by exchanging HLA-A*02:01 with HLA-A*03:01, HLA-A*11:01 or HLA-B*07:02. To generate the TMC, eGFP within the lentiviral transfer vector backbone was replaced with ubiquitin to generate pTRPE-Ubiquitin-T2A-mCherry. A synthetic DNA vector construct was synthesized (TWIST Bioscience, San Francisco, CA) encoding the first 25 amino acids of WT and mKRAS (G12C, G12D, G12R, G12V) expressed in tandem without linker sequences (Supplementary Fig. 4). The TMC also included viral epitopes embedded within natural 20 amino acid viral protein sequences restricted to HLA-A*02:01 (NLVPMVATV), HLA-A*03:01 (ILRGSVAHK), HLA-A*11:01 (IVTDFSVIK), and HLA-B*07:02 (TPRVTGGGAM). The TMC was PCR amplified and inserted between ubiquitin and T2A sequences. Synthetic TCR DNA vector constructs were synthesized (TWIST Bioscience, San Francisco, CA) to include the TCRα chain followed by the TCRβ chain separated by a T2A sequence. TRAC and TRBC regions of synthetic TCRs were codon altered to make them resistant to Cas9 protein riboprobes targeting endogenous TRAC and TRBC1/TRBC2. TCR constructs were PCR amplified and cloned into a lentiviral vector backbone. High titer lentiviral vector production was performed^53,54.

Monoallelic KRAS–tandem minigene construct (TMC) cell lines

HLA class I monoallelic cell lines were generated by transduction of K562 or 721.221 cell lines using lentiviral particles encoding HLA-SCD/eGFP⁵². Cells were subsequently transduced with lentiviral particle encoding TMC/mCherry then sorted to >99% purity by flow cytometry based on coexpression of GFP, mCherry, and HLA class expression as measured by APC-conjugated anti-human HLA-A, B, C (clone W6/32) staining.

Proteomic validation of mKRAS antigen processing and presentation

Sample preparation

Monoallelic KRAS-TMC K562 and 721.221 cell lines were expanded to 1–2 × 10⁸ total cells and HLA class I immunopurification (IP) was performed³¹ using MHC class I (W6/32) antibody non-covalently linked to agarose beads (Santa Cruz Biotechnology, Dallas, TX). Peptides were eluted from HLA class I molecules using 0.1% trifluoroacetic acid (TFA). Immunoprecipitation eluent was passed through a 10,000 Da Amicon molecular weight cut off filter (Merck Millipore) via centrifugation at 10,000 × g for 10 min. Filtered eluent was dried and resuspended in 100 μl of 0.1% TFA. Stage tip C18 columns (Harvard Apparatus) were equilibrated with 200 μl of acetonitrile and 200 μl of 0.1% TFA, and samples were loaded onto columns for desalting. Samples were then washed with three rounds of 200 μl of 0.1% TFA then eluted in 70% acetonitrile in 0.1% formic acid (FA) and dried.

Data-dependent acquisition (DDA)

DDA analysis of monoallelic KRAS-TMC K562 and 721.221 cell lines was performed by the Quantitative Proteomics Resource Core at the University of Pennsylvania. DDA experiments were carried out using half of each enriched sample by nano LC-MS/MS using an Orbitrap Fusion (Thermo Scientific) coupled to an Easy-nLC system (Thermo Scientific). Samples were resuspended in 10 μl of 0.1% TFA and 5 μl was injected for each analysis. Samples were separated with an in-house packed column with ReproSil-Pur C18 AQ 3 μm resin with dimensions 75 μm × 20 cm (Dr Maisch GmbH, Ammerbuch, Germany) at a flow rate of 400 nl/min. Using 0.1% FA as buffer A and 80% ACN 0.1% FA as buffer B, peptides were eluted with a gradient of 4% buffer B to 35% buffer B in 50 min, then to 65% buffer B in 5 min, and finally a 5 min wash at 95% buffer B. Peptides were ionized in a Nanospray Flex Ion Source (Thermo Scientific) at 2.3 kV. An MS1 scan was acquired using a resolution of 120,000 FWHM, AGC target of 4e5, and maximum inject time of 50 ms. Precursors were isolated with 1.6 m/z isolation window and fragmented with an HCD collision energy of 30 eV. An MS2 scan was acquired using a resolution of 30,000 FWHM, AGC target of 5e4, and maximum ion inject time of 54 ms in 3 s cycle time scan mode. DDA data were processed with PEAKS X+ version 10.5 with 10 ppm parent mass tolerance, 0.02 Da fragment mass error tolerance, and no enzyme specificity. Data were searched against the SwissProt Human proteome appended with the custom KRAS peptide sequences.

Parallel reaction monitoring (PRM)

PRM analysis of monoallelic KRAS-TMC K562 and 721.221 cell lines was performed by the Quantitative Proteomics Resource Core at the University of Pennsylvania. Targeted Analysis by PRM was performed on the remaining half of each sample to increase the sensitivity of mKRAS epitope detection. Samples were analyzed on a Q-Exactive HF-X (Thermo Scientific) coupled to an Ultimate 3000 nano UHPLC system (Thermo Scientific) using a PRM strategy³². Samples were resuspended in 10 μl of 0.1% TFA and 5 μl was injected for each analysis. Samples were separated with an in-house packed column with ReproSil-Pur C18 AQ 3 μm resin with dimensions 75 μm × 20 cm (Dr Maisch GmbH, Ammerbuch, Germany) at a flow rate of 400 nl/min. Using 0.1% FA as buffer A and 80% ACN 0.1% FA as buffer B, peptides were eluted with a gradient of 4% buffer B to 30% buffer B in 42 min, then to 65% buffer B in 6 min, followed by a 7 min wash at 95% buffer B, and a 5 min equilibration at 4% buffer B. Peptides were ionized in a Nanospray Flex Ion Source (Thermo Scientific) at 2.3 kV. An MS1 scan was acquired using a resolution of 120,000 FWHM, AGC target of 1e5, and maximum inject time of 50 ms. An inclusion list of 33 unique ions was used to sequentially isolate and fragment each target peptide with an isolation window of 1.6 m/z and an HCD collision energy of 28 eV. An MS2 scan was acquired using a resolution of 15,000 FWHM, AGC target of 1e5, maximum ion inject time of 100 ms, and loop count of 10. Results were analyzed with Skyline⁵⁵ (version 20.1) and reference spectra of synthetic peptides (New England Peptide, Gardner, MA) were used to validate fragmentation patterns. Spectra were manually annotated with IPSA⁵⁶. Synthetic peptides were analyzed on a Q-Exactive HF (Thermo Scientific) coupled to an Ultimate 3000 nano UHPLC system (Thermo Scientific) with DDA. Each synthetic peptide was injected at 1 pmol/μl and analyzed with similar liquid chromatography conditions. An MS1 scan was acquired at a resolution of 120,000 FWHM, AGC target of 3e6, maximum inject time of 32 ms. The top 20 intense ions were isolated and fragmented with a dynamic exclusion of 45 s. For each peptide, an MS2 scan was acquired using a resolution of 15,000 FWHM, AGC target of 2e5, maximum ion inject time of 32 ms, and isolation window of 1.4 m/z. Ions were filtered for charges 2–5. Raw files were searched using Proteome Discoverer 2.2 (Thermo Scientific) against a database of targets with non-trypic digestion, precursor mass tolerance of 10 ppm, and fragment mass tolerance of 0.02 Da. Search results were filtered with the target decoy approach⁴⁰. The validation of the precursor ion identity was based on isotopic dot product (idotp) and dot product (dotp) values generated by Skyline³³. Both values provide a measurement of similarity between eluted and expected or reference peptide. The idotp value provides a measure to assess the precursor isotope distribution and its correlation between observed (eluted) and expected (synthetic) peptide while the dotp value provides a measure of similarity between observed and library spectrum (Proteome Discoverer Database) peak areas. Values range from a best of 1.0 to a worst of 0.0.

Proteomic quantitation of mKRAS peptides presented by human tumor cells

Sample preparation

CORL23-A3 and CORL23-A11 cells were expanded to 1–2 × 10⁸ total cells. HLA class I IP was performed by Cayman Chemical (Ann Arbor, MI). Cells were resuspended in MHC lysis buffer consisting of DPBS (0.25 % Sodium deoxycholate, 200 μM iodoacetamide, 1% N-Octyl-β-D- thioglucoside, 1 mM EDTA) containing protease inhibitor (1 ml of lysis buffer per 1 × 10⁸ cells). Cell pellets were resuspended by pipetting and vortexing six times for 3–4 s each at 5-min intervals. Lysate was centrifuged at 800 × g for 5 min and the supernatant was transferred to a 50-ml tube. Additional lysis buffer was added to the supernatants (to double the volume) prior to centrifugation at 20,000 × g for 60 min. The supernatants were transferred to a new 15 ml conical tube. Then, 200 μl from each tube was saved for analysis (Pre-IP). The cell lysate was then mixed with 100 μl of W6/32 resin and incubated overnight at 4 °C with gentle rotation. The IA resin (W6/32 antibody DMP crosslinked to protein A resin) was washed with lysis buffer before use. The mixture was centrifuged at 800 × g for 5 min and the supernatant was transferred to a new tube (500 μl of the supernatant (Post-IP) was retained for analysis). The IA resin pellets were washed in 2.5 ml of lysis buffer, centrifuged at 800 × g for 5 min at 4 °C, and the supernatant was aspirated. Two additional washes were carried out using 2.5 ml each of wash buffers 2 (20 mM Tris-HCL, 400 mM Nacl, pH 8.0) and 3 (20 mM Tris-HCL, 150 mM Nacl, pH 8.0). The IA resin was transferred to a low protein binding Eppendorf tube for a final wash. 0.75 ml of wash buffer 4 (20 mM Tris-HCL, pH 8.0) was added, mixed, centrifuged at 800 × g for 5 min and the supernatant was aspirated. Peptides were eluted from MHC class I molecules using 1 ml of MHC-I elution buffer (0.1 M Acetic acid, 0.1% TFA), incubated for 5 min at 37 °C and centrifuged at 800 × g for 5 min at 4 °C. The eluate was collected into fresh low protein binding tubes, and was centrifuged at 20,000 × g for 2 min to clear the sample of any remaining resin. The eluate was collected into a new low protein binding tube, flash frozen in liquid nitrogen, and stored at −80 °C. A quality assurance ELISA assay was performed to quantify the enrichment of HLA class I protein in Pre- and Post-IP samples.

Data-dependent acquisition (DDA)

DDA analysis of CORL23 cell lines was performed by MS Bioworks (Ann Arbor, MI). DDA experiments were carried out using half of each enriched sample by nano LC-MS/MS using a Waters M-Class system interfaced to a ThermoFisher Fusion Lumos mass spectrometer. Peptides were loaded on a trapping column and eluted over a 75-μm analytical column at 350 nl/min; both columns were packed with Luna C18 resin (Phenomenex). A 2 h reverse phase gradient was employed. The mass spectrometer was operated in a combined data-dependent EThcD/CID mode, with MS and MS/MS performed in the Orbitrap at 60,000 FWHM resolution and 15,000 FWHM resolution, respectively. The instrument was run with a 3-s cycle for MS and MS/MS. DDA data were processed with PEAKS X+ version 10.5 with 10 ppm parent mass tolerance, 0.02 Da fragment mass error tolerance, and no enzyme specificity. Data were searched against the SwissProt Human proteome appended with the custom KRAS peptide sequences. Peptide 8-16V (VVGAVGVGK) and 7-16V (VVVGAVGVGK) identifications were confirmed with analysis of synthetic peptide standards (New England Peptide).

Quantification of mKRAS epitopes by parallel reaction monitoring (PRM)

The quantification of mKRAS epitopes expressed by HLA-engineered CORL23 cell lines was performed by targeted MS using PRM (MS Bioworks, Ann Arbor, MI). Peptides were enriched as described above. Synthetic stable labeled peptides VVGAVGVGK^ and VVVGAVGVGK^ were purchased from New England Peptide, where ^ is Lysine (¹³C₆¹⁵N₂). In all, 200 fmols of each stable labeled peptide was added to the enriched samples for analysis. Each enriched sample was analyzed in duplicate (50% of the sample per injection). PRM was performed with a Waters M-Class HPLC system interfaced to a ThermoFisher Fusion Lumos mass spectrometer. Peptides were loaded on a trapping column and eluted over a 75-μm analytical column at 350 nl/min; both columns were packed with Luna C18 resin (Phenomenex). The mass spectrometer was operated in PRM mode with the Orbitrap operating at 17,500 FWHM resolution. Collision-induced dissociation data were collected for the (M+2H)²⁺ charge state ions of the target peptides VVGAVGVGK, VVGAVGVGK^, VVVGAVGVGK, and VVVGAVGVGK^. Extracted ion chromatograms for each of the target peptides were generated manually using the XCalibur QualBrowser software version 4.1.31.9 (ThermoFisher) with a 20 ppm mass tolerance for product ions. For the calculation of 8-16V and 7-16V epitopes expressed by HLA-A*03:01 and HLA-A*11:01 complexes, eluted and internal standard peptide peak area data were used to calculate the number of moles of peptide present in the sample. This number was doubled (half the IP was analyzed in a single injection) and converted to molecules by multiplying by Avogadro’s number. The result was divided by the number of input cells (164M and 163M for the CORL23-A3 and CORL23-A11 cell lines, respectively) to give number of peptide molecules/cell.

Generation of mKRAS-reactive T cell lines

The generation of monocyte-derived mDC and antigen-specific T cells was performed⁵⁷. Briefly, DC were generated from monocyte cultures in GM-CSF (100 ng/ml, Berlex) and IL-4 (20 ng/ml, Miltenyi Biotec) and matured using CD40L-expressing K562 cells, IFN-γ (100 u/ml, CellGenix), poly I:C (Invivogen, Inc), and R848 (Invivogen, Inc.) for 16 h to generate mDC. CD8⁺ T cells were isolated from peripheral blood mononuclear cells, using a CD8 negative selection kit (Miltenyi Biotech, Auburn, CA). Purified CD8⁺ T cells (5 × 10⁶ cells/ml) were cultured at a 20:1 ratio with irradiated (2500 Rads) autologous mDC pulsed with peptide (40 μg per 1 × 10⁶ DC/ml) in 24-well trays in Optimizer CST media (Gibco) supplemented with 5% pooled human sera. Human IL-7 (10 ng/ml), IL-15 (5 ng/ml), and IL-12 (10 ng/ml) were added on day 0. Fresh media supplemented with IL-7 (10 ng/ml) and IL-15 (5 ng/ml) was added on day 7. Fourteen days after primary mDC stimulation, T cell cultures were harvested and re-stimulated with irradiated (2500 Rads) peptide-pulsed mDC. Cell culture media was supplemented with 50 U/ml IL-2 (Chiron, Emeryville, CA) starting day 2 then every 48 h following secondary stimulation. On days 10–14 of secondary mDC stimulation, antigen-specific T cell responses were identified by IFN-γ ELISPOT assay and flow cytometry by p-HLA multimer staining.

IFN-γ ELISPOT assay

CD8⁺ T cell reactivity to peptide antigen was assessed by interferon-γ (IFN-γ) ELISPOT assay⁵⁸. The spot number was determined in an independent blinded fashion (ZellNet Consulting, New York, NY) using the high-resolution automated KS ELISPOT reader (Zeiss, Thornwood, NY) and KS ELISPOT 4.9.16 software with reading parameters established per International Harmonization Guidelines⁵⁹. A positive response was recorded if the number of spots in the peptide-exposed wells was two times or more higher than the number of spots in the unstimulated wells and if there was a minimum of 20 (after subtraction of background spots) peptide-specific spots per 5 × 10⁵ CD8⁺ cells.

p-HLA multimer assay

mKRAS-specific CD8⁺ T cell frequencies were determined by staining with PE-conjugated p-HLA dextramers (Immudex), followed by addition of APC-conjugated anti-human CD8 antibody (Invitrogen). Cells were washed and resuspended in FACS buffer. One million events in the CD8⁺ gate were collected using a hierarchical gating strategy. Data were acquired using BD FACSDiva software (version 8.0.2) and analyzed using FlowJo software (version 10.7.2).

Isolation and validation of KRAS-specific TCRs

mKRAS-specific CD8⁺ T cell cultures underwent antigen-specific expansion using irradiated (10,000 Rads) HLA-SCD/4-1BBL expressing K562 artificial APC at a 1:1 ratio. Cultures were supplemented with IL-2 (500 U/ml) 48 h after initiation then every 48 h thereafter until culture termination. T cells were expanded for 12–14 days, sorted to 97–99% purity by CD8 and p-HLA multimer coexpression.

TCRα/TCRβ sequencing

DNA was extracted using the Gentra Puregene cell kit following the manufacturer’s directions (Qiagen, Valencia, CA, Cat. No. 158388). Bulk DNA TCR Vβ and Vα libraries were prepared for sequencing on the Illumina MiSeq platform using a cocktail of 23 Vβ families from framework region 2 (FR2) forward primers, and 13 Jβ region reverse primers, modified from the BIOMED2 primer series⁶⁰ and a cocktail of 39 Vα from FR3 primers and 50 Jα region reverse primers, respectively. Libraries were generated using the QIAGEN Multiplex PCR and Illumina Nextera XT index kits⁶¹. RNA was extracted using RNeasy Plus Mini Kit following the manufacturer’s directions (Qiagen, Cat. No. 74134). Bulk RNA TCR Vβ and Vα libraries were prepared using SMARTer^® Human TCR α/β Profiling Kit following the manufacturer’s directions (Takara Bio USA, Inc., Mountain View, CA, Cat. No. 635014). Single cells were sorted into a 96-well plate and single cell sequencing libraries were prepared using SMARTer^® Human scTCR a/b Profiling Kit following the manufacturer’s directions (Takara Bio USA, Inc., Cat. No. 634431). Libraries were sequenced (2 × 300 bp paired end reads, MiSeq Reagent Kit v3-600 cycle, Illumina, San Diego, CA, USA; Cat. No. 102-3003) on an Illumina MiSeq instrument in the Human Immunology Core Facility at the University of Pennsylvania.

Data analysis

Raw sequences were quality filtered⁶² and clone assemblies were processed with MiXCR (v. 3.0.7)⁶³ and VDJtools (v1.2.1)⁶⁴.

Jurkat reporter system to assess TCR antigen specificity and avidity

Using Cas9 protein and riboprobes as guides against TRAC (5’-TGTGCTAGACATGAGGTCTA-3’) and TRBC1/TRBC2 (5’-GGAGAATGACGAGTGGACCC-3’)³⁸, a Jurkat E6.1 TCRα/β-negative cell line was generated. To generate J^ASP90 reporter cells, TCR-negative Jurkat E6.1 cells were then transduced with lentiviral particles encoding human CD8αβ heterodimer (generously provided by Dr Jim L. Riley, University of Pennsylvania, Philadelphia, PA) and the Uni-Vect reporter system that features constitutive mCherry and NFAT-inducible eGFP expression as a means to assess TCR signaling (generously provided by Dr Daniel J. Powell Jr., University of Pennsylvania, Philadelphia, PA). J^ASP90 reporter cells were flow cytometrically sorted to purity based on TCR^– (anti-human TCRα/β antibody clone IP26, BioLegend, San Diego, CA), CD8⁺ (anti-CD8 antibody, ThermoFisher Scientific, Waltham, MA), mCherry⁺ expression, and low basal NFAT (eGFP⁺) signal. The TCRα and TCRβ chains of TCRA3V, TCRA11V, and TCRB7R were expressed via lentivirus in J^ASP90 using the pTRPE vector to generate the TCRA3V, TCRA11V, and TCRB7R J^ASP90 reporter cell lines. These cell lines were sorted to yield a unimodal p-HLA multimer⁺ cell population and expanded for use in functional assays. Sorted TCRA3V, TCRA11V, and TCRB7R J^ASP90 reporter cell lines were mixed 1:1 with HLA-SCD-expressing K562 cells pulsed with titrated peptide concentrations (10 μM–1 pM). After 16–20 h, cells were analyzed by flow cytometry to determine the percentage of eGFP⁺ cells in each sample. Cells activated with PMA (50 ng/ml) and Ionomycin (750 ng/ml) were included as positive controls (%GFP_max), and J^ASP90 reporter cells cultured in media alone were used as negative controls (%GFP_min). Data were fitted to a dose-response curve by linear non-regression analysis using GraphPad Prism version 7.0.

TCRs as probes for recognition of mKRAS epitopes

Tumor cell lines (2.5 × 10⁵ cells/well) were cocultured with TCR transduced J^ASP90 cells at a 1:1 ratio in 48-well flat bottom plates in a final volume of 500 μl of media (RPMI, 10% FBS, l-Glutamine, Penicillin/Streptomycin). Each tumor cell line was tested in its parental (non-HLA matched) and HLA-engineered (HLA-matched) form using lentiviral particles encoding HLA-SCD constructs of interest and flow cytometric sorting to >99% purity based on the coexpression of the target HLA class I allele and eGFP. G12R⁺ tumor cell lines were pre-cultured with IFN-γ (500 U/ml) for 48 h prior to co-culture with TCRB7R transduced J^ASP90 cells. After 16–20 h of incubation at 37 °C, cells were analyzed by flow cytometry to determine the percentage of live, CD8⁺mCherry⁺eGFP⁺ cells in each sample. Cells activated with PMA and Ionomycin were included as positive controls.

Expression of mKRAS-specific TCRs in primary CD8⁺ T cells

Primary T cell expansion

On day 0, human primary CD8⁺ T cells were activated and expanded using anti-CD3/CD28 antibody–conjugated paramagnetic microbeads (Life Technologies). On day 1 the cells were transduced with lentiviral vector particles containing TCRA3V, TCRA11V, or TCRB7R genes at an MOI of 5. Cell cultures and maintained in media supplemented with IL-7 (5 ng/ml; R&D Systems) and IL-15 (5 ng/ml; R&D Systems). On day 5, T cell cultures were harvested and paramagnetic microbeads were magnetically removed prior to electroporation.

Cas9 RNP assembly and electroporation

Successful knock out of the endogenous TCR was carried out using single guide RNAs targeting TRAC (5’-TGTGCTAGACATGAGGTCTA-3’) and TRBC1/TRBC2 (5’-GGAGAATGACGAGTGGACCC-3’)³⁸ obtained from Integrated DNA Technologies. For Cas9 and sgRNA delivery, the ribonucleoprotein (RNP) complex was first pre-formed by incubating TrueCut Cas9 Protein V2 (10 ug, ThermoFisher Scientific) with sgRNA (5 ug) for 10 min at RT. Cells were spun down at 300 × g for 10 min and resuspended at a concentration of 5 × 10⁶ cells/100 μl in the P3 buffer. Both TCR α/β RNP complexes and 100 μl of resuspended cells were combined and electroporated (P3 Primary Cell 4-D Kit, Lonza:PBP3-02250). TCRαβ^null T cells were then cultured in media supplemented with IL-7 and IL-15 for a total of 14 days.

Secondary T cell expansion

T cell cultures were harvested and re-stimulated with peptide-pulsed (500 ng/ml) irradiated (10,000 Rads) HLA-SCD/4-1BBL expressing K562 artificial APCs at a 1:1 ratio. Cultures were supplemented with IL-7 and IL-15 every 48 h until culture termination. Engineered TCR expression was then assessed by CD3, CD8, and TCRαβ expression as well as p-HLA multimer binding then cryopreserved for use in subsequent assays.

Cytotoxicity assays

⁵¹Cr-release assay

K562 cells and KRAS-mutated tumor cell lines were labeled with 25 μCi ⁵¹Cr in the presence or absence of peptide (10 μM) for 1 h at 37 °C, washed and tested as targets in a standard 4 h ⁵¹Cr-release assay. Effector cells consisted of primary gene-edited (TCRα-/TCRβ-) mKRAS-specific TCR CD8⁺ T cells, designated TCRA3V, TCRA11V, TCRB7R. Transgenic TCR expression was assessed by p-HLA multimer assay. Assays were performed, in triplicate, at various effector: target ratios. Data were collected using a MicroBeta2 LumiJET Microplate Counter (PerkinElmer). Data are represented as percent-specific lysis reported as mean ± standard deviation (SD).

Real-time apoptotic cell death analysis

Real-time apoptotic cell death analysis (live cell imaging with cellular impedance) was performed to assess extended cytotoxic activity using the xCELLigence Real Time Cell Analysis eSight system (ACEA Biosciences). Target tumor cells were plated (1 × 10⁴ cells/well) and allowed to adhere for 24 h. Effector T cells were added at E:T ratios 10:1, 3:1, and 1:1, and the media was supplemented with Annexin V-CF594 (Biotium, Fremont, CA). Cell index (relative cell impedance) was monitored every 15 min for 5 days and normalized to the maximum cell index value immediately prior to effector-cell plating. Shaded lines reflect the mean of duplicate wells ± SD. Concurrent time lapse video monitoring was performed with acquisition of brightfield, green (GFP), and red (CF594) windows every hour.

CORL23 lung xenograft mouse model

All animal procedures were performed according to the approved University of Pennsylvania Institutional Animal Care and Use Committee (IACUC) protocol 804226. NSG (NOD/scid/γcnull) mice were purchased from the University of Pennsylvania Stem Cell Xenograft Core (SCXC) and housed in micro-isolator cages under sterile conditions in the SCXC AAALAC accredited animal facility at the University of Pennsylvania (temperature: 24 °C, humidity: 50–60%, dark/light cycle: 12 h/12 h). The CORL23 xenograft model was established by intravenous tail injections of 5 × 10⁵ CORL23-A3 or CORL23-A11 cells expressing CBR luciferase for bioluminescence imaging (BLI). On day 4 post tumor inoculation, mice were treated with 1 × 10⁷ mKRAS-TCR engineered CD8⁺ T cells (TCRA3V for CORL23-A3 Cohort, TCRA11V for CORL23-A11 Cohort). Control mice either received no treatment (Mock) or 1 × 10⁷ TCR KO CD8⁺ T cells. BLI was used to monitor tumor growth by intraperitoneal injections of 3 mg D-luciferin/150 ul PBS per mouse (~120 mg/kg) and imaged after 10 min using an IVIS Spectrum imaging system. Data were analyzed using Living Image Version 4.5.2 software (PerkinElmer).

Statistics and reproducibility

Statistical analysis of multiple comparisons was performed using one-way ANOVA with Tukey’s HST post-test, and comparisons between just two groups were performed using Student’s unpaired t-test. Significance of overall survival was determined via Kaplan–Meier analysis with log-rank (Mantel-Cox) analysis. Data were analyzed using Microsoft Excel (version 16.16.27), and all statistical analyses were performed with Graphpad Prism 7.0 (GraphPad), except analysis of tumor growth and survival curves, which were performed using linear mixed-effects model with Tukey’s HSD post-test using the Ime4 and the survival package in R. Error bars represent the SD, and p < 0.05 was considered statistically significant. *p < 0.05, **p < 0.01, ***p < 0.001, and ****p < 0.0001 unless otherwise indicated. ns denotes not significant. All data presented are representative of two or more independent experiments.

Reporting summary

Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability

MS raw data, DDA and PRM files supporting this study have been deposited in the PRIDE repository [https://www.ebi.ac.uk/pride/archive/projects/PXD024412]. TCRα and TCRβ sequencing data described in this study have been deposited in the NCBI GenBank nucleotide database under the following accession numbers: MZ355910, MZ355911, MZ355912, MZ355913, MZ355914, MZ355915, In Fig. 1b, HLA class I allelic frequencies are supported by data obtained from Allele Frequency Net Database [http://www.allelefrequencies.net] and mKRAS frequencies are supported by data obtained from International Cancer Genome Consortium Data Portal [https://dcc.icgc.org]. KRAS gene expression values reported in Supplementary Table 4 are supported by data obtained from European Bioinformatics Institute (EMBL-EBI) Expression Atlas Portal [https://www.ebi.ac.uk/gxa/]. The remaining data are available within the article, Supplementary information, or Source Data file. Source data are provided with this paper.

References

Downward, J. Targeting RAS signalling pathways in cancer therapy. Nat. Rev. Cancer 3, 11–22 (2003).
Article CAS PubMed Google Scholar
Karnoub, A. E. & Weinberg, R. A. Ras oncogenes: split personalities. Nat. Rev. Mol. Cell Biol. 9, 517–531 (2008).
Article CAS PubMed PubMed Central Google Scholar
Simanshu, D. K., Nissley, D. V. & McCormick, F. RAS proteins and their regulators in human disease. Cell 170, 17–33 (2017).
Article CAS PubMed PubMed Central Google Scholar
Cox, A. D., Fesik, S. W., Kimmelman, A. C., Luo, J. & Der, C. J. Drugging the undruggable RAS: mission possible? Nat. Rev. Drug Discov. 13, 828–851 (2014).
Article CAS PubMed PubMed Central Google Scholar
Canon, J. et al. The clinical KRAS(G12C) inhibitor AMG 510 drives anti-tumour immunity. Nature 575, 217–223 (2019).
Article ADS CAS PubMed Google Scholar
Hallin, J. et al. The KRASG12C inhibitor MRTX849 provides insight toward therapeutic susceptibility of KRAS-mutant cancers in mouse models and patients. Cancer Disco. 10, 54–71 (2020).
Article CAS Google Scholar
Yadav, M. et al. Predicting immunogenic tumour mutations by combining mass spectrometry and exome sequencing. Nature 515, 572–576 (2014).
Article ADS CAS PubMed Google Scholar
Bassani-Sternberg, M. et al. Direct identification of clinically relevant neoepitopes presented on native human melanoma tissue by mass spectrometry. Nat. Commun. 7, 13404 (2016).
Article ADS CAS PubMed PubMed Central Google Scholar
Ribas, A. & Wolchok, J. D. Cancer immunotherapy using checkpoint blockade. Science 359, 1350–1355 (2018).
Article ADS CAS PubMed PubMed Central Google Scholar
Cohen, C. J. et al. Isolation of neoantigen-specific T cells from tumor and peripheral lymphocytes. J. Clin. Invest. 125, 3981–3991 (2015).
Article PubMed PubMed Central Google Scholar
Carreno, B. M. et al. Cancer immunotherapy. A dendritic cell vaccine increases the breadth and diversity of melanoma neoantigen-specific T cells. Science 348, 803–808 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Strønen, E. et al. Targeting of cancer neoantigens with donor-derived T cell receptor repertoires. Science 2288, 1–11 (2016).
Google Scholar
Priestley, P. et al. Pan-cancer whole-genome analyses of metastatic solid tumours. Nature 575, 210–216 (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
ICGC/TCGA Pan-Cancer Analysis of Whole Genomes Consortium. Pan-cancer analysis of whole genomes. Nature 578, 82–93 (2020).
Article ADS CAS Google Scholar
Gjertsen, M. K. et al. Vaccination with mutant ras peptides and induction of T-cell responsiveness in pancreatic carcinoma patients carrying the corresponding RAS mutation. Lancet (Lond., Engl.) 346, 1399–1400 (1995).
Article CAS Google Scholar
Gjertsen, M. K., Bjorheim, J., Saeterdal, I., Myklebust, J. & Gaudernack, G. Cytotoxic CD4+ and CD8+ T lymphocytes, generated by mutant p21-ras (12Val) peptide vaccination of a patient, recognize 12Val-dependent nested epitopes present within the vaccine peptide and kill autologous tumour cells carrying this mutation. Int. J. Cancer 72, 784–790 (1997).
Article CAS PubMed Google Scholar
Abrams, S. I. et al. Generation of stable CD4+ and CD8+ T cell lines from patients immunized with ras oncogene-derived peptides reflecting codon 12 mutations. Cell. Immunol. 182, 137–151 (1997).
Article CAS PubMed Google Scholar
Chaft, J. E. et al. Phase II study of the GI-4000 KRAS vaccine after curative therapy in patients with stage I-III lung adenocarcinoma harboring a KRAS G12C, G12D, or G12V mutation. Clin. Lung Cancer 15, 405–410 (2014).
Article CAS PubMed Google Scholar
Cafri, G. et al. Memory T cells targeting oncogenic mutations detected in peripheral blood of epithelial cancer patients. Nat. Commun. 10, 449 (2019).
Article ADS CAS PubMed PubMed Central Google Scholar
Tran, E. et al. Immunogenicity of somatic mutations in human gastrointestinal cancers. Science 350, 1387–1390 (2015).
Article ADS CAS PubMed PubMed Central Google Scholar
Wang, Q. J. et al. Identification of T-cell receptors targeting KRAS-mutated human tumors. Cancer Immunol. Res. 4, 204–214 (2016).
Article CAS PubMed Google Scholar
Tran, E. et al. T-cell transfer therapy targeting mutant KRAS in cancer. N. Engl. J. Med. 375, 2255–2262 (2016).
Article CAS PubMed PubMed Central Google Scholar
Hundal, J. et al. pVAC-Seq: a genome-guided in silico approach to identifying tumor neoantigens. Genome Med. 8, 11 (2016).
Article PubMed PubMed Central CAS Google Scholar
Richman, L. P., Vonderheide, R. H. & Rech, A. J. Neoantigen dissimilarity to the self-proteome predicts immunogenicity and response to immune checkpoint blockade. Cell Syst. 9, 375–382.e4 (2019).
Article CAS PubMed PubMed Central Google Scholar
Sarkizova, S. et al. A large peptidome dataset improves HLA class I epitope prediction across most of the human population. Nat. Biotechnol. 38, 99–209 (2020). https://doi.org/10.1038/s41587-019-0322-9
Pearson, H. et al. MHC class I-associated peptides derive from selective regions of the human genome. J. Clin. Invest. 126, 4690–4701 (2016).
Article PubMed PubMed Central Google Scholar
Haen, S. P. & Rammensee, H.-G. The repertoire of human tumor-associated epitopes—identification and selection of antigens and their application in clinical trials. Curr. Opin. Immunol. 25, 277–283 (2013).
Article CAS PubMed Google Scholar
Neefjes, J., Jongsma, M. L. M., Paul, P. & Bakke, O. Towards a systems understanding of MHC class I and MHC class II antigen presentation. Nat. Rev. Immunol. 11, 823–836 (2011).
Article CAS PubMed Google Scholar
Harndahl, M. et al. Peptide-MHC class I stability is a better predictor than peptide affinity of CTL immunogenicity. Eur. J. Immunol. 42, 1405–1416 (2012).
Article CAS PubMed Google Scholar
Buchli, R. et al. Development and validation of a fluorescence polarization-based competitive peptide-binding assay for HLA-A*0201—a new tool for epitope discovery. Biochemistry 44, 12491–12507 (2005).
Article CAS PubMed Google Scholar
Abelin, J. G. et al. Mass spectrometry profiling of HLA-associated peptidomes in mono-allelic cells enables more accurate epitope prediction. Immunity 46, 315–326 (2017).
Article CAS PubMed PubMed Central Google Scholar
Peterson, A. C., Russell, J. D., Bailey, D. J., Westphall, M. S. & Coon, J. J. Parallel reaction monitoring for high resolution and high mass accuracy quantitative, targeted proteomics. Mol. Cell. Proteom. 11, 1475–1488 (2012).
Article CAS Google Scholar
Pino, L. K. et al. The Skyline ecosystem: informatics for quantitative mass spectrometry proteomics. Mass Spectrom. Rev. 39, 229–244 (2020).
Article ADS CAS PubMed Google Scholar
Bassani-Sternberg, M., Pletscher-Frankild, S., Jensen, L. J. & Mann, M. Mass spectrometry of human leukocyte antigen class I peptidomes reveals strong effects of protein abundance and turnover on antigen presentation. Mol. Cell. Proteom. 14, 658–673 (2015).
Article CAS Google Scholar
Trolle, T. et al. Automated benchmarking of peptide-MHC class I binding predictions. Bioinformatics 31, 2174–2181 (2015).
Article CAS PubMed PubMed Central Google Scholar
Wang, Q. et al. Direct detection and quantification of neoantigens. Cancer Immunol. Res. 7, 1748–1754 (2019).
Article CAS PubMed PubMed Central Google Scholar
Linette, G. P. et al. Immunological ignorance is an enabling feature of the oligo-clonal T cell response to melanoma neoantigens. Proc. Natl Acad. Sci. USA 116, 23662–23670 (2019).
Article CAS PubMed PubMed Central Google Scholar
Stadtmauer, E. A. et al. CRISPR-engineered T cells in patients with refractory cancer. Science 367, eaba7365 (2020). https://doi.org/10.1126/science.aba7365
Schober, K. et al. Orthotopic replacement of T-cell receptor α- and β-chains with preservation of near-physiological T-cell function. Nat. Biomed. Eng. 3, 974–984 (2019).
Article PubMed Google Scholar
Elias, J. E. & Gygi, S. P. Target-decoy search strategy for increased confidence in large-scale protein identifications by mass spectrometry. Nat. Methods 4, 207–214 (2007).
Article CAS PubMed Google Scholar
Jin, X. et al. A metastasis map of human cancer cell lines. Nature 588, 331–336 (2020).
Article ADS CAS PubMed PubMed Central Google Scholar
Mullard, A. Cracking KRAS. Nat. Rev. Drug Discov. 18, 887–891 (2019).
Article CAS PubMed Google Scholar
Oldfield, L. E., Connor, A. A. & Gallinger, S. Molecular events in the natural history of pancreatic cancer. Trends Cancer 3, 336–346 (2017).
Article CAS PubMed Google Scholar
Chan-Seng-Yue, M. et al. Transcription phenotypes of pancreatic cancer are driven by genomic events during tumor evolution. Nat. Genet. 52, 231–240 (2020). https://doi.org/10.1038/s41588-019-0566-9
Mishto, M. et al. An in silico-in vitro pipeline identifying an HLA-A*02:01+ KRAS G12V+ spliced epitope candidate for a broad tumor-immune response in cancer patients. Front. Immunol. 10, 2572 (2019).
Article CAS PubMed PubMed Central Google Scholar
Kubuschok, B. et al. Naturally occurring T-cell response against mutated p21 ras oncoprotein in pancreatic cancer. Clin. Cancer Res. 12, 1365–1372 (2006).
Article CAS PubMed Google Scholar
Ott, P. A. et al. An immunogenic personal neoantigen vaccine for patients with melanoma. Nature 547, 217–221 (2017).
Article ADS CAS PubMed PubMed Central Google Scholar
Sahin, U. et al. Personalized RNA mutanome vaccines mobilize poly-specific therapeutic immunity against cancer. Nature 547, 222–226 (2017).
Article ADS CAS PubMed Google Scholar
González-Galarza, F. F. et al. Allele frequency net 2015 update: new features for HLA epitopes, KIR and disease and HLA adverse drug reaction associations. Nucleic Acids Res 43, D784–D788 (2015).
Article PubMed CAS Google Scholar
Harndahl, M., Rasmussen, M., Roder, G. & Buus, S. Real-time, high-throughput measurements of peptide-MHC-I dissociation using a scintillation proximity assay. J. Immunol. Methods 374, 5–12 (2011).
Article CAS PubMed Google Scholar
Scholtalbers, J. et al. TCLP: an online cancer cell line catalogue integrating HLA type, predicted neo-epitopes, virus and gene expression. Genome Med. 7, 118 (2015).
Article PubMed PubMed Central CAS Google Scholar
Carreno, B. M. et al. Amino-terminal extended peptide single-chain trimers are potent synthetic agonists for memory human CD8+ T cells. J. Immunol. 188, 5839–5849 (2012).
Article CAS PubMed Google Scholar
Parry, R. V., Rumbley, C. A., Vandenberghe, L. H., June, C. H. & Riley, J. L. CD28 and inducible costimulatory protein Src homology 2 binding domains show distinct regulation of phosphatidylinositol 3-kinase, Bcl-xL, and IL-2 expression in primary human CD4 T lymphocytes. J. Immunol. 171, 166–174 (2003).
Article CAS PubMed Google Scholar
Yin, Y. et al. Checkpoint blockade reverses anergy in IL-13Rα2 humanized scFv-based CAR T cells to treat murine and canine gliomas. Mol. Ther. Oncolytics 11, 20–38 (2018).
Article CAS PubMed PubMed Central Google Scholar
MacLean, B. et al. Skyline: an open source document editor for creating and analyzing targeted proteomics experiments. Bioinformatics 26, 966–968 (2010).
Article CAS PubMed PubMed Central Google Scholar
Brademan, D. R., Riley, N. M., Kwiecien, N. W. & Coon, J. J. Interactive peptide spectral annotator: a versatile web-based tool for proteomic applications. Mol. Cell. Proteom. 18, S193–S201 (2019).
Article CAS Google Scholar
Carreno, B. M. et al. IL-12p70-producing patient DC vaccine elicits Tc1-polarized immunity. J. Clin. Invest. 123, 3383–3394 (2013).
Article CAS PubMed PubMed Central Google Scholar
Gottschalk, S. et al. Generating CTLs against the subdominant Epstein-Barr virus LMP1 antigen for the adoptive immunotherapy of EBV-associated malignancies. Blood 101, 1905–1912 (2003).
Article CAS PubMed Google Scholar
Janetzki, S. et al. Guidelines for the automated evaluation of Elispot assays. Nat. Protoc. 10, 1098–1115 (2015).
Article CAS PubMed Google Scholar
van Dongen, J. J. M. et al. Design and standardization of PCR primers and protocols for detection of clonal immunoglobulin and T-cell receptor gene recombinations in suspect lymphoproliferations: report of the BIOMED-2 concerted action BMH4-CT98-3936. Leukemia 17, 2257–2317 (2003).
Article PubMed Google Scholar
Rosenfeld, A. M. et al. Computational evaluation of B-cell clone sizes in bulk populations. Front. Immunol. 9, 1472 (2018).
Article PubMed PubMed Central CAS Google Scholar
Meng, W. et al. An atlas of B-cell clonal distribution in the human body. Nat. Biotechnol. 35, 879–884 (2017).
Article CAS PubMed PubMed Central Google Scholar
Bolotin, D. A. et al. MiXCR: software for comprehensive adaptive immunity profiling. Nat. Methods 12, 380–381 (2015).
Article CAS PubMed Google Scholar
Shugay, M. et al. VDJtools: unifying post-analysis of T cell receptor repertoires. PLOS Comput. Biol. 11, e1004503 (2015).
Article PubMed PubMed Central CAS Google Scholar

Download references

Acknowledgements

We thank the Human Immunology Core at the Univerisity of Pennsylvania for providing leukocytes for research, Drs Nina Luning Prak and Wenzhao Meng for assistance in TCR sequencing, Drs Rico Buchli, Thomas Osterbye, and Anette Stryhn for assistance in biochemical characterization of peptide epitopes, Dr Andrew Rech for bioinformatic analysis, Bereket Gebregziabher and Dr Nune Markosyan for statistical analysis, and Drs Carl June and Elizabeth Jaffee for advice and helpful discussions. This work was supported by NIH R01 CA204261 (to B.M.C. and G.P.L.), NIH P30 CA016520 (to R.H.V.), NIH CA196539 and CA232568 (to B.A.G.), SU2C/Lustgarten Foundation Pancreatic Cancer Collective (to R.H.V. and B.M.C.), the Parker Institute for Cancer Immunotherapy (to B.M.C., G.P.L., and R.H.V.), and Institute for Immunology at University of Pennsylvania (to A.S.B. and G.P.L.). Training grants provided support to A.S.B. (K12 CA076931) and T.B. (T32 CA009140).

Author information

These authors jointly supervised this work: Robert H. Vonderheide, Gerald P. Linette, Beatriz M. Carreno.

Authors and Affiliations

Division of Hematology-Oncology, Department of Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
Adham S. Bear, Mark H. O’Hara, Robert H. Vonderheide & Gerald P. Linette
Abramson Cancer Center, University of Pennsylvania, Philadelphia, PA, USA
Adham S. Bear, Lee P. Richman, Mark H. O’Hara & Robert H. Vonderheide
Center for Cellular Immunotherapies, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
Tatiana Blanchard, Chong Xu, Miren L. Baroja, Khatuna Gabunia, John Scholler, Avery D. Posey Jr., Anze Smole, Daniel J. Powell Jr., Robert H. Vonderheide, Gerald P. Linette & Beatriz M. Carreno
Department of Biochemistry and Biophysics, University of Pennsylvania, Philadelphia, PA, USA
Joseph Cesare
MSBioworks, Ann Arbor, MI, USA
Michael J. Ford
Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
Sarah McCuaig & Christina Costeas
Department of Systems Pharmacology and Translational Therapeutics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
Avery D. Posey Jr.
Corporal Michael J. Crescenz VA Medical Center, Philadelphia, PA, USA
Avery D. Posey Jr.
Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
Daniel J. Powell Jr. & Beatriz M. Carreno
Department of Biochemistry and Biophysics, Epigenetics Institute, University of Pennsylvania, Philadelphia, PA, USA
Benjamin A. Garcia
Parker Institute for Cancer Immunotherapy, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
Robert H. Vonderheide, Gerald P. Linette & Beatriz M. Carreno

Authors

Adham S. Bear
View author publications
You can also search for this author in PubMed Google Scholar
Tatiana Blanchard
View author publications
You can also search for this author in PubMed Google Scholar
Joseph Cesare
View author publications
You can also search for this author in PubMed Google Scholar
Michael J. Ford
View author publications
You can also search for this author in PubMed Google Scholar
Lee P. Richman
View author publications
You can also search for this author in PubMed Google Scholar
Chong Xu
View author publications
You can also search for this author in PubMed Google Scholar
Miren L. Baroja
View author publications
You can also search for this author in PubMed Google Scholar
Sarah McCuaig
View author publications
You can also search for this author in PubMed Google Scholar
Christina Costeas
View author publications
You can also search for this author in PubMed Google Scholar
Khatuna Gabunia
View author publications
You can also search for this author in PubMed Google Scholar
John Scholler
View author publications
You can also search for this author in PubMed Google Scholar
Avery D. Posey Jr.
View author publications
You can also search for this author in PubMed Google Scholar
Mark H. O’Hara
View author publications
You can also search for this author in PubMed Google Scholar
Anze Smole
View author publications
You can also search for this author in PubMed Google Scholar
Daniel J. Powell Jr.
View author publications
You can also search for this author in PubMed Google Scholar
Benjamin A. Garcia
View author publications
You can also search for this author in PubMed Google Scholar
Robert H. Vonderheide
View author publications
You can also search for this author in PubMed Google Scholar
Gerald P. Linette
View author publications
You can also search for this author in PubMed Google Scholar
Beatriz M. Carreno
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

Conceptualization: R.H.V., G.P.L., B.M.C.; methodology: A.S.B., T.B., B.M.C., G.P.L., J.C., M.J.F.; software: L.P.R.; validation: M.L.B., C.X., T.B.; formal analysis: B.A.G., M.J.F., B.M.C.; investigation: A.S.B., T.B., J.C., M.J.F., C.X., M.L.B., S.M.C., C.C., A.D.P., M.H.O.; resources: A.S., D.J.P., A.D.P.; data curation: A.S.B., B.M.C., L.P.R.; writing – original draft: A.S.B., G.P.L., B.M.C.; writing – review and editing: A.S.B., R.H.V., G.P.L., B.M.C.; visualization: A.S.B., J.C., R.H.V., G.P.L., B.M.C.; supervision: G.P.L., B.M.C., R.H.V., B.A.G.; funding: G.P.L., R.H.V., B.M.C., B.A.G.

Corresponding authors

Correspondence to Adham S. Bear or Beatriz M. Carreno.

Ethics declarations

Competing interests

R.H.V. reports having received consulting fees or honoraria from Medimmune and Verastem; and research funding from Fibrogen, Janssen, and Lilly. He is an inventor on a licensed patent relating to cancer cellular immunotherapy and cancer vaccines, and receives royalties from Children’s Hospital Boston for a licensed research-only monoclonal antibody. A.S.B., R.H.V., G.P.L., and B.M.C. have filed a patent application (PCT/US2020/014988) related to the targeting of KRAS for immunotherapy entitled “Compositions and Methods for Targeting Mutant RAS”. A.S., A.D.P., and D.J.P. have filed a patent application (PCT/US2019/026378) that incorporates discoveries related to the NFAT reporter system used in this study entitled “Methods and Compositions Comprising a Viral Vector for Expression of a Transgene and an Effector”. M.J.F. is a founder and an employee of MS Bioworks. All other authors have no competing interest to declare.

Additional information

Peer review information Nature Communications thanks Victor Engelhard, Angela Krackhardt, Juliane Liepe, and Benjamin Vincent for their contribution to the peer review of this work. Peer review reports are available.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Peer Review File

Descriptions of additional Supplementary Files

Supplementary Data 1

Supplementary Movie 1

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Bear, A.S., Blanchard, T., Cesare, J. et al. Biochemical and functional characterization of mutant KRAS epitopes validates this oncoprotein for immunological targeting. Nat Commun 12, 4365 (2021). https://doi.org/10.1038/s41467-021-24562-2

Download citation

Received: 23 April 2020
Accepted: 16 June 2021
Published: 16 July 2021
DOI: https://doi.org/10.1038/s41467-021-24562-2

This article is cited by

Challenges in developing personalized neoantigen cancer vaccines
- Peter D. Katsikis
- Ken J. Ishii
- Christopher Schliehe
Nature Reviews Immunology (2024)
Combinatorial strategies to target RAS-driven cancers
- Naiara Perurena
- Lisa Situ
- Karen Cichowski
Nature Reviews Cancer (2024)
Lymph-node-targeted, mKRAS-specific amphiphile vaccine in pancreatic and colorectal cancer: the phase 1 AMPLIFY-201 trial
- Shubham Pant
- Zev A. Wainberg
- Eileen M. O’Reilly
Nature Medicine (2024)
Targeting small GTPases: emerging grasps on previously untamable targets, pioneered by KRAS
- Guowei Yin
- Jing Huang
- Yijun Gao
Signal Transduction and Targeted Therapy (2023)
Conformational plasticity of RAS Q61 family of neoepitopes results in distinct features for targeted recognition
- Andrew C. McShan
- David Flores-Solis
- Nikolaos G. Sgourakis
Nature Communications (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Subjects

Abstract

Similar content being viewed by others

Introduction

Results

Identification and characterization of mKRAS neoepitopes

Validation of mKRAS antigen processing and presentation

mKRAS peptide immunogenicity and mKRAS-specific TCRs isolation

mKRAS-specific TCRs confer lytic activity against human tumor cells

Quantitation and recognition of mKRAS peptides presented by human tumor cells

In vivo anti-tumor activity of mKRAS-specific TCR+ CD8+ T cells

Discussion

Methods

Primary cells

Cell lines

Antibodies and flow cytometry

Computational prediction of mKRAS G12 neoepitopes

Peptides and biochemical peptide-MHC binding assays

HLA class I genotyping and phenotyping of cell lines

Lentiviral vector constructs and production

Monoallelic KRAS–tandem minigene construct (TMC) cell lines

Proteomic validation of mKRAS antigen processing and presentation

Sample preparation

Data-dependent acquisition (DDA)

Parallel reaction monitoring (PRM)

Proteomic quantitation of mKRAS peptides presented by human tumor cells

Sample preparation

Data-dependent acquisition (DDA)

Quantification of mKRAS epitopes by parallel reaction monitoring (PRM)

Generation of mKRAS-reactive T cell lines

IFN-γ ELISPOT assay

p-HLA multimer assay

Isolation and validation of KRAS-specific TCRs

TCRα/TCRβ sequencing

Data analysis

Jurkat reporter system to assess TCR antigen specificity and avidity

TCRs as probes for recognition of mKRAS epitopes

Expression of mKRAS-specific TCRs in primary CD8+ T cells

Primary T cell expansion

Cas9 RNP assembly and electroporation

Secondary T cell expansion

Cytotoxicity assays

51Cr-release assay

Real-time apoptotic cell death analysis

CORL23 lung xenograft mouse model

Statistics and reproducibility

Reporting summary

Data availability

References

Acknowledgements

Author information

Authors and Affiliations

Contributions

Corresponding authors

Ethics declarations

Competing interests

Additional information

Supplementary information

Source data

Rights and permissions

About this article

Cite this article

Share this article

This article is cited by

Comments

Search

Quick links

Expression of mKRAS-specific TCRs in primary CD8⁺ T cells

⁵¹Cr-release assay