Rachel L.
Norman
a,
Rajinder
Singh
a,
Frederick W.
Muskett
bc,
Emma L.
Parrott
a,
Alessandro
Rufini
a,
James I.
Langridge
d,
Franscois
Runau
a,
Ashley
Dennison
a,
Jacqui A.
Shaw
a,
Elena
Piletska
ef,
Francesco
Canfarotta
e,
Leong L.
Ng
g,
Sergey
Piletsky
ef and
Donald J. L.
Jones
*ag
aLeicester Cancer Research Centre, Leicester Royal Infirmary, University of Leicester, Leicester, LE1 5WW, UK. E-mail: Djlj1@le.ac.uk
bDepartment of Molecular and Cell Biology, University of Leicester, LE1 7RH Leicester, UK
cLeicester Institute of Structural and Chemical Biology, University of Leicester, LE1 7RH Leicester, UK
dWaters Corporation, Wilmslow, SK9 4AX, UK
eMIP Diagnostics, The Exchange Building, Colworth Park, MK44 1LQ, Bedford, UK
fSchool of Chemistry, University of Leicester, University Road, Leicester, LE1 7RH, UK
gDepartment of Cardiovascular Sciences, University of Leicester and National Institute for Health Research Leicester Biomedical Research Centre, Glenfield Hospital, Leicester, LE1 7RH, UK
First published on 29th November 2021
Cancer is a disease of cellular evolution where single base changes in the genetic code can have significant impact on the translation of proteins and their activity. Thus, in cancer research there is significant interest in methods that can determine mutations and identify the significant binding sites (epitopes) of antibodies to proteins in order to develop novel therapies. Nano molecularly imprinted polymers (nanoMIPs) provide an alternative to antibodies as reagents capable of specifically capturing target molecules depending on their structure. In this study, we used nanoMIPs to capture KRAS, a critical oncogene, to identify mutations which when present are indicative of oncological progress. Herein, coupling nanoMIPs (capture) and liquid chromatography-mass spectrometry (detection), LC-MS has allowed us to investigate mutational assignment and epitope discovery. Specifically, we have shown epitope discovery by generating nanoMIPs to a recombinant KRAS protein and identifying three regions of the protein which have been previously assigned as epitopes using much more time-consuming protocols. The mutation status of the released tryptic peptide was identified by LC-MS following capture of the conserved region of KRAS using nanoMIPS, which were tryptically digested, thus releasing the sequence of a non-conserved (mutated) region. This approach was tested in cell lines where we showed the effective genotyping of a KRAS cell line and in the plasma of cancer patients, thus demonstrating its ability to diagnose precisely the mutational status of a patient. This work provides a clear line-of-sight for the use of nanoMIPs to its translation from research into diagnostic and clinical utility.
KRAS (Kirsten RAt Sarcoma) is a member of the RAS family of genes6 and represents a classic oncogene that has clinical utility.7 The three RAS proto-oncogenes encode four proteins: NRAS, HRAS and the two splice variants KRAS4A and KRAS4B. The four proteins have high sequence homology, mainly differing at C terminal residues 165-188/9, which is also known as the hypervariable region. Ras proteins act as molecular switches by alternating between an active GTP-bound state and an inactive GDP-bound conformation. Ras proteins constitute major players in the MAPK pathway and are activated in response to ligand binding to Epidermal Growth Factor Receptor (EGFR).8 Ligand binding to EGFR leads to dimerisation of the receptor, changing its conformation so that growth factor receptor-bound protein 2 (GRB2) can bind. GRB2 recruits son of sevenless (SOS), a guanine nucleotide exchange factor (GEF) which dissociates guanine nucleotides from KRAS, allowing the more abundant GTP to bind and activate KRAS.9,10 Active KRAS binds to and enables activation of the BRAF kinase, which, in turn, stimulates cell division via the MEK-ERK signalling cascade.
At the clinical level, mutation of KRAS occurs in nearly 30% of all cancers and thus provides a clear target for diagnostics and targeted therapies.11 The most common mutation in KRAS is at codon 12, which in the wild type (WT) protein is occupied by a glycine (G) residue. When mutations occur, the glycine can alter to a number of alternatives including cysteine, aspartic acid, valine, alanine and serine. Other mutations occur at position 13, 61 and 146.12 Clinically, KRAS mutation status is widely used for colorectal and lung cancer but has utility in pancreatic and thyroid cancers.13 A number of methods can be used including pyrosequencing, Sanger sequencing, high-resolution melting analysis, single-strand conformation polymorphism analysis and allele-specific PCR methods that target mutation hotspots in KRAS exons 2 and 3 are currently employed in clinical testing of KRAS mutations.13 Attempts have been made to use immune capture prior to MS with reasonable levels of success in tissues but not plasma, preventing their application in a non-invasive test.14–16
Epitope discovery seeks to identify the region of a protein that binds most strongly to an antibody. There are many approaches, all of which are quite technically challenging. Modern methods for epitope discovery include Cryo-Electron Microscopy (CryoEM),17 Mass Spectrometry (e.g. Hydrogen–Deuterium Exchange (HDX)18 and MALDI-MS using crosslinking approaches19), and Nuclear Magnetic Resonance (NMR)20 to determine sites of interaction. Other biological methods include site-directed mutagenesis21 or ELISA inhibitor competition assays.22 All these methods provide relatively high resolution determination of epitopes but are all technically challenging, and very time consuming. Most importantly, they all require an antibody, the production of which can take months.
Molecularly imprinted polymers (MIPs) are synthetic binders, able to capture specific molecules by means of their target-specific binding sites. They are formed by polymerising a mixture of monomers in the presence of a target molecule which acts as a template.23 Once the template has been removed, it leaves a cavity in the polymer which retains its configuration and has the correct shape and functional group orientation to bind the target molecule with specific affinity.24 As they have memory of the size shape and functional groups, they can effectively be reused. They have shown to be highly selective in a number of applications for small molecules, peptides and proteins.25–27 They have had many different applications, as sensors,28 purification tools29,30 and in protein production.31 More recently, MIPs have been used in diagnostic and sensor applications providing a superior alternative to antibodies thanks to their robustness and reproducible synthesis that does not require the use of animals.26,32,33 Moreover, due to their high affinity, they allow for preconcentration of the analyte in question which provides a distinct advantage in scenarios where the analyte is expected to be at very low concentrations within a complex matrix.
In this paper, we describe the synthesis and application of MIPs as a diagnostic tool to recognise the C-terminus of the KRAS protein, which is conserved in both the WT and mutated protein. We use a synthetic mimic that contains the C-terminus of KRAS to generate the nanoMIPs. It is at the N-terminus that the clinically relevant mutations exist. The combination of MIP-assisted capture and mass spectrometry was used to distinguish between WT and mutated forms of KRAS, and we successfully tested this approach in three distinct biological systems demonstrating the utility of this approach for the diagnosis of cancer. We were also able to apply the nanoMIPs as a tool for epitope discovery, by identifying the three main epitopes of KRAS.
Fig. 1B shows that the recovery of the bridged peptide from plasma by the nanoMIPs was unaffected by the complex matrix which is demonstrated by comparing within the same figure the recovery from water. Moreover, we were able to demonstrate the nanoMIP re-usability by eluting the peptide and then repeating the experiment with another aliquot of plasma. Fig. 1C shows the pattern of intensity between 1st contact, elution followed by a further 2nd contact is very similar.
Once the initial evaluation took place, the nanoMIPs were tested for their efficacy to phenotype cancer cells by detecting KRAS protein in cell culture. Two lung adenocarcinoma cancer cell lines were selected, namely A549 cells which are homozygous for the G12S mutation and H1650 cells which contains WT KRAS. Cells were lysed and the lysate was then incubated with MIP-coated glass beads for 2 hours and the captured proteins tryptically digested. Resultant tryptic peptides were measured using LC-MS/MS (SRM). SRM chromatograms targeting WT KRAS signal in H1650 cells showed several confirmatory transitions for the analyte (Fig. 2A) which eluted with the same retention time as the stable isotope labelled (SIL) standard (Fig. 2B). Conversely in the A549 cells, a homozygous cell line containing the G12S mutation presented several confirmatory transitions for the G12S mutated KRAS protein (Fig. 2C) which co-eluted with the isotopically labelled standard for G12S (Fig. 2D). These data indicate that nanoMIPs are able to select and identify mutations from in vitro biological models which would enable investigation of these models and, with appropriate calibration curves, allow the identification and quantitation of mutant KRAS in clinical specimens.
Notwithstanding our successful detection of KRAS bridged peptides in spiked plasma samples, we reasoned that measurement of endogenous KRAS in such milieu is likely to be hampered due to the inherent low copy number of circulating KRAS proteins.14 Thus, we explored whether we could use nanoMIPs to effectively capture and enrich circulating KRAS. Plasma samples from eight treatment naïve patients with pancreatic cancer were obtained from Prof Ashley Dennison and Franscois Runau (University of Leicester). Plasma samples from 25 patients with non-small cell lung cancer (NSCLC) were obtained from Prof Jacqui Shaw (University of Leicester). The genotype of these samples was unknown, but about 90% of pancreatic cancer patients and about a 20% of NSCLC patients harbor KRAS mutations.11 Calibration lines were established for 6 of the main mutations and other RAS forms.34 For every sample, we were able to detect and measure WT KRAS (Fig. 3), and to the best of our knowledge, this is the first occasion that the gene product of KRAS has been measured in cancer patient plasma using LC-SRM.
Fig. 3 Quantification of KRAS WT 6–16 peptide in pancreatic cancer (blue) and NSCLC (red) plasma samples. The figure shows a spread of KRAS measurements for each cancer. |
However, mutated KRAS forms would inevitably be present in extremely low amounts possibly orders of magnitude less than the WT form. By using our integrated approach, we have detected in 4 pancreatic cancer patients a mutated form (G12C, G12S, G12A and G12D) but at the very limit of detection where a certain degree of doubt would exist in mutation status assignment (Fig. 4). This finding together with the very low concentrations of WT KRAS indicate that the detection of mutated KRAS would require sensitivity beyond the currently achievable limit of detection.
To further validate our approach, we carried out western blot for KRAS using an anti-panKRAS antibody (this antibody captures KRAS without discerning the isoform of the KRAS protein) in seven selected patient samples (5 NSCLC (lanes 2–5, 7), one pancreatic cancer (lane 6) and one healthy control (lane 8) where we had previously identified the WT KRAS by LC-SRM. We also examined two lysates from colorectal cancer cell lines (HCT-N or HT-29) that act as positive controls (lanes 9&10). KRAS was highly expressed in the two positive controls (see G12C) whilst the KRAS protein could only be detected in three out of seven plasma samples, thus confirming our previous finding that MIPs with LC-SRM detect KRAS proteins in plasma samples with better sensitivity than western blot (Fig. 5).
Finally, we examined the ability of nanoMIPs to discover epitopes for the whole KRAS protein. The workflow is shown in Fig. 6. We produced recombinant KRAS protein using a KRAS cDNA sequence supplied as a plasmid expressed in a DH10B TonA E. coli host (all details described in ESI†). The purified recombinant KRAS protein was used as a template to generate nanoMIPs using a modified solution imprinting process.23 The protein was then tryptically digested whilst still bound to the nanoMIP. Sections of the protein that were not directly protected by the MIP binding site and underwent digestion were washed away. Only the bound peptides of KRAS that were within the binding site of the nanoMIPs would represent epitopes. The epitope peptides were eluted from the nanoMIPs using hot water (95 °C) and identified using LC-MS/MS. In total, three peptides were identified by this epitope discovery process (see Table 1). The experiment was repeated independently three times and the same three peptides emerged. When the recombinant protein was tryptically digested in the absence of MIPs, a set of tryptic peptides resulted, many of which were stronger in intensity than the three putative epitope peptides. This result indicates that the peptides identified as epitopes are not merely the most abundant tryptic peptides, but they reflect the direct and specific interaction between nanoMIPs and KRAS.
Fig. 6 Identification of peptide sequences exposed on protein surfaces using molecular imprinting and mass spectrometry. |
Epitope | Peptide MW | Sequence | Residues | R t (min) | Ave. intens. | Peptides listed in the immune epitope database35–37 |
---|---|---|---|---|---|---|
1 | 955.5935 | (K)LVVVGAGGVGK(S) | 6–16 | 29.4 | 22603 | LVVVGAGGVGKSA/LVVVGAGGVGKSALT |
2 | 1312.678 | (R)SYGIPFIETSAK(T) | 136–147 | 36.5 | 41018 | IPFIETSAKTRQGVDD |
3 | 1383.69 | (R)QGVDDAFYTLVR(E) | 150–161 | 37.1 | 5680 | IPFIETSAKTRQGVDD |
To further support this result, NMR spectroscopy was used to map the binding site, or sites, of the MIPs on KRAS. 2D 15N-TROSY spectra of 15N labelled KRAS were recorded in the absence and presence of MIPs. The effects on the spectra are reported as a chemical shift perturbation analysis (reviewed by Williamson 201338) and three candidate epitopes were mapped in-silico onto the crystal structure of KRAS4B. The nanoMIPs produced small but significant chemical shift changes in the KRAS spectrum. These were mapped onto the crystal structure as shown in Fig. 7. Orange represents histidine residues and yellow regions are areas where there is no data. The chemical shifts are represented by the white and red areas where white indicates no change and darker red indicates larger chemical shift changes. The NMR data showed that there were multiple binding sites rather than a single site, corroborating the mass spectrometry analysis. The red region, indicating a greater chemical shift at the bottom right-hand corner of Fig. 7D, could match epitope three thus corroborating the MS data and suggesting that the peptides identified by mass spectrometry represent binding sites on the protein.
The three potential epitopes identified using nanoMIPs and mass spectrometry also overlap with confirmed antibody epitopes. Most commercially available antibodies are raised to the N-terminal region of KRAS, which includes the 6–16 peptide which was identified as a binding site by the nanoMIPs and mass spectrometry (Epitope 1). In addition, Santa Cruz Biotechnology, Inc. have developed an antibody raised against a region spanning amino acids 54–189, which includes epitopes two and three. In addition, the three epitopes identified using nanoMIPs and mass spectrometry were all either completely or partially covered by peptides listed in the Immune Epitope Database, as shown in Table 1. Overall, these data collectively show the ability of nanoMIPs and MS to identify protein epitopes.
The synthesis of the nanoMIPs to the whole protein was carried out in order to discover KRAS epitopes. The experiment was repeated three times and demonstrated reproducible identification of epitopes which were substantiated by NMR. The three epitopes were consistent with those found in the Epitope databases37 as shown in Table 1. A study by Bergseng et al.35 into coeliac disease-associated human leukocyte antigen (HLA) molecules identified peptides eluted from HLA molecules. Epitope 1 was listed as being identified by MaxQuant in peptides eluted from HLA molecules DQ2.5, DQ2.2 and DQ7.5. Epitopes 2 and 3 were partially identified in the immunopeptidome eluted from major histocompatibility complex (MHC) class II molecules in A20 mouse lymphoma cells.36 This suggests that the peptides identified using MIPs and mass spectrometry are all part of known immune epitopes.
This work demonstrates that the data obtained using this novel method – recently patented – is analogous to current methods to some extent, and could therefore be a feasible alternative. The method works on the assumption that the polymer protects the region of the protein that it binds to from proteolytic digestion. However, tryptic peptides may be longer than the region enveloped by the polymer. To develop this method further in future, combinations of different proteolytic enzymes could be investigated to obtain complete epitope sequences. This approach provides the possibility of identifying regions of the protein surface which have not been demonstrated to be antigenic in vivo, but which may offer improved affinity.
In conclusion we show that the use of nanoMIPs epitope discovery coupled with MS analysis allows epitope identification which, in contrast to the use of antibodies with NMR, is significantly quicker (mainly through the avoidance of antibody production which can take months), does not require in vivo experiments and could potentially be automated into a high-throughput screening approach. We also demonstrate mutation detection in clinical samples, which has great potential in the deployment of precision medicine.
Following polymerisation, dynamic light scattering (DLS – Zetasizer, Malvern Instruments, UK) was used to assess the average hydrodynamic diameter. The nanoMIPs imprinted against the peptide epitope showed a hydrodynamic diameter of 215.9 nm (±5 (n = 6)), whilst nanoMIPs imprinted against a recombinant KRAS protein showed an average particle size of 185.6 nm (±5 (n = 6)).
50 μL of each plasma sample was denatured with 350 μL 8 M urea, adjusted to pH 7.8. Denatured plasma was added to 350 mg KRAS4A peptide MIPs, immobilised on glass beads, on a 0.45 μm Spin X centrifugal filter and incubated for 2 hours on a sample rotator. Filters were centrifuged at 8000 rpm for 10 minutes at 4 °C. The filtrate was discarded then the beads were washed with 200 μL PBS three times, followed by 200 μL ammonium bicarbonate, pH 7.8, three times, centrifuging for 5 minutes at 8000 rpm and discarding the filtrate each time. Beads were then washed with 200 μL 10% acetonitrile/10% formic acid twice and the filtrate was collected and stored at −80 °C. Filters were transferred to clean tubes then 200 μL ammonium bicarbonate with 15 mM DTT were added and samples were incubated for 30 minutes at 60 °C. 20 mM IAA was added and beads were incubated for 30 minutes in the dark at room temperature. 16 μg trypsin was added and samples were incubated at 37 °C overnight.
Samples were centrifuged for 10 minutes at 8000 rpm and 20 °C. 1% formic acid was added to the filtrate then samples were incubated for 20 minutes at 37 °C before centrifuging for 10 minutes at 14000g. The supernatant was desalted via solid phase extraction which employed empore SPE cartridges. After lyophilisation, samples were reconstituted in 0.1% formic acid/3% acetonitrile with 25 fmol μL−1 stable isotope-labelled standards in a matrix of 1 mg mL−1 BSA. BSA was added to prevent non-specific binding to the sample vials because the calibration line in solvent displayed reduced peak intensity for many of the peptides. Non-specific binding was eradicated by the addition of digested human plasma as a sample matrix in the standards for the calibration line. In the MIP on-bead digestion samples, most of the patient plasma was removed during the washing steps. A human plasma background would interfere with the analysis as it would contain endogenous RAS protein thus BSA was used as an alternative sample matrix to coat the vials, reducing non-specific binding of the analytes.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1nr03180e |
This journal is © The Royal Society of Chemistry 2021 |