Exercise:AntigenProcessing ans

• Q1: How many epitopes do you find with the initial filter?

• Q2: Do you agree with calling these peptides epitopes? Why?

Now refine your search by including only mass-spectrometry assays on HLA-DR restriction alleles. From the References tab, select the publication from the Journal of Proteome Research from 2017 by Wang et al.

• Q3: How many eluted ligands did they find in this study?

Under the Assays tab, select MHC ligand assays and export these results.

Open the exported file using Excel or your preferred program and inspect the different columns. We are interested in the epitope modification column (Modifications).

Note that there are two peptide sequences that contain post-translational modifications.

• Q4: What kind of post-translational modifications are present?

These modified peptides look very interesting. However, the MHC Motif Deconvolution software cannot currently deal with them, so exclude them from the dataset.

In addition, inspect the column containing the MHC Types present in the host.

• Q5: Which alleles are expressed in the patient used for this assay?

You should now have both the peptide dataset and the HLA typing information required to use MHCMotifDecon.

MHCMotifDecon

Go to the MHCMotifDecon server.

Enter the submission page and select MHC class II.

Paste the list of peptides into the input window. The input should contain only one column of peptide sequences.

Scroll down and select the HLA-DR alleles expressed by the cell line for which the immunopeptidomics assay was performed.

Press Submit.

Note: Be aware that HLA alleles are entered into this software using a particular naming format. For example:

HLA-DRB1*01:01 → DRB1_0101

• Q6: How many peptides were assigned to each of the alleles, and how many were assigned to the trash cluster?

• Q7: Can you see any differences between the trash-cluster peptides and the peptides assigned to a specific allele? Explain.

Go back to the previous submission and modify the %rank threshold.

• Q8: What is the effect of changing the %rank threshold?

We know that HLA-DRB1 is not the only HLA-DR protein expressed by humans. HLA-DRB3, HLA-DRB4, and HLA-DRB5 are located in close proximity to HLA-DRB1 on the chromosome and are often expressed together.

In particular, cell lines expressing HLA-DRB1*04:01 are strongly associated with the expression of a second allele called HLA-DRB4*01:01.

If you want to learn more about this topic, search for HLA-DR linkage disequilibrium.

Having this information, extend the HLA typing of the experiment by including:

DRB4_0101

Run MHCMotifDecon again with this additional allele.

• Q9: How many peptides are now associated with HLA-DRB4? Where did these peptides come from?

GibbsCluster

We will now analyse the same ligand dataset using an unsupervised approach.

Go to the GibbsCluster server.

Enter the submission page and paste the list of ligands into the input window.

Because GibbsCluster uses unsupervised learning, you do not need to specify the HLA alleles present in the sample.

Select the MHC class II parameters.

Change:

Number of iterations per sequence per temperature step = 100

and select:

Preference for hydrophobic AAs at P1

Click Submit.

• Q10: Does the number of clusters found by GibbsCluster correspond to the number of HLA alleles in your sample?

• Q11: Can you identify other differences between the GibbsCluster solution and the MHCMotifDecon solution?

---

You should now have experience with two different approaches for analysing complex immunopeptidomics datasets:

• MHCMotifDecon – supervised motif deconvolution using known HLA alleles and MHC-binding predictions.
• GibbsCluster – unsupervised clustering of peptide sequences without prior knowledge of the HLA alleles.

Consider how prior knowledge of the HLA type influences the interpretation of an immunopeptidomics dataset and how the results from the two approaches differ.

Done!