<?xml version="1.0"?>
<feed xmlns="http://www.w3.org/2005/Atom" xml:lang="en">
	<id>https://teaching.healthtech.dtu.dk/22145/index.php?action=history&amp;feed=atom&amp;title=Protein_Drug_Deimmunization</id>
	<title>Protein Drug Deimmunization - Revision history</title>
	<link rel="self" type="application/atom+xml" href="https://teaching.healthtech.dtu.dk/22145/index.php?action=history&amp;feed=atom&amp;title=Protein_Drug_Deimmunization"/>
	<link rel="alternate" type="text/html" href="https://teaching.healthtech.dtu.dk/22145/index.php?title=Protein_Drug_Deimmunization&amp;action=history"/>
	<updated>2026-05-02T22:05:19Z</updated>
	<subtitle>Revision history for this page on the wiki</subtitle>
	<generator>MediaWiki 1.41.0</generator>
	<entry>
		<id>https://teaching.healthtech.dtu.dk/22145/index.php?title=Protein_Drug_Deimmunization&amp;diff=23&amp;oldid=prev</id>
		<title>WikiSysop: Created page with &quot;== Protein Drug Deimmunization - Exercise ==  ----  Over the past few decades, biologics have gone from being experimental treatment options to standard care. They are very targeted and useful molecules in health care, but they come with their own immunological challenges. If a biologic is recognized as foreign by the immune system, Anti-Drug Antibodies can be formed that neutralize the biologic and render it useless. Deimmunization provides a potential solution to this...&quot;</title>
		<link rel="alternate" type="text/html" href="https://teaching.healthtech.dtu.dk/22145/index.php?title=Protein_Drug_Deimmunization&amp;diff=23&amp;oldid=prev"/>
		<updated>2024-03-20T14:43:37Z</updated>

		<summary type="html">&lt;p&gt;Created page with &amp;quot;== Protein Drug Deimmunization - Exercise ==  ----  Over the past few decades, biologics have gone from being experimental treatment options to standard care. They are very targeted and useful molecules in health care, but they come with their own immunological challenges. If a biologic is recognized as foreign by the immune system, Anti-Drug Antibodies can be formed that neutralize the biologic and render it useless. Deimmunization provides a potential solution to this...&amp;quot;&lt;/p&gt;
&lt;p&gt;&lt;b&gt;New page&lt;/b&gt;&lt;/p&gt;&lt;div&gt;== Protein Drug Deimmunization - Exercise ==&lt;br /&gt;
&lt;br /&gt;
----&lt;br /&gt;
&lt;br /&gt;
Over the past few decades, biologics have gone from being experimental treatment options to standard care. They&lt;br /&gt;
are very targeted and useful molecules in health care, but they come with their own immunological challenges.&lt;br /&gt;
If a biologic is recognized as foreign by the immune system, Anti-Drug Antibodies can be formed that neutralize&lt;br /&gt;
the biologic and render it useless. Deimmunization provides a potential solution to this challenge.&lt;br /&gt;
&lt;br /&gt;
The general strategy of deimmunization is to target T-Cell epitopes with mutation, thereby reducing immunogenicity of the biologic. Epitopes can be identified experimentally, but in silico prediction makes the process more&lt;br /&gt;
economically tractable.&lt;br /&gt;
&lt;br /&gt;
Here we will go through a short exercise on in silico deimmunization of human Erythropoietin (EPO). EPO&lt;br /&gt;
is a protein that signals for the production of red blood cells in response to hypoxia. It is used in treatment of&lt;br /&gt;
anemia, but is perhaps better known for its use in blood doping. We will use the [http://www.IEDB.org IEDB] to get epitope data on&lt;br /&gt;
EPO, [https://www.uniprot.org/ Uniprot] to get its sequence, and [http://www.IEDB.org IEDB] again for tools on CD4 epitope prediction and deimmunization.&lt;br /&gt;
&lt;br /&gt;
The goal of today’s exercise is to get you acquainted with existing CD4 epitope prediction- and deimmunization-&lt;br /&gt;
tools. You will handle and process data and question some technical aspects of the tools and results.&lt;br /&gt;
&lt;br /&gt;
== Data Gathering ==&lt;br /&gt;
&lt;br /&gt;
Start by going to the [http://www.IEDB.org IEDB homepage] and find T Cell epitope results for human erythropoietin. If nothing comes up, try reformulating your query. It should result in 19 epitopes from a single reference. In this publication, EPO was tested for immunogenicity by testing 15-mers; overlapping by 10-mers, covering the entire sequence. Under the antigen tab, navigate to the immunome browser to get an overview of the results.&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q1: Looking at the T-cell epitope response frequency profile of EPO, in what region/regions of the protein is immunogenicity mostly found?&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q2: Again, looking at the response frequency profile, an irregular pattern is observed for the first 27 amino acids. In the table of results below, it becomes apparent that no epitope assays were performed on this region. Why do you think that is? Hint: EPO is an extracellular protein&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
Export the immunome browser results to a .csv file.&lt;br /&gt;
&lt;br /&gt;
The answer to question 2 can be found on [https://www.uniprot.org/ Uniprot], where we will also get the sequence for human EPO. Go to [https://www.uniprot.org/ Uniprot] and search for Erythropoietin. Select the human entry and download the .fasta file with the EPO sequence (click &amp;#039;format&amp;#039; tab at top of page and select fasta). On the Uniprot page for human EPO, you can learn that EPO is a glycoprotein.&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q3: Can you find any other information on the EPO Uniprot site that could explain the irregular response frequency pattern in the N-terminal of EPO(as seen in Question 2)?&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
This N-terminal sequence will interfere with our prediction so remove the first 27 residues in your EPO fasta and save it as e.g. EPO-processed.fasta&lt;br /&gt;
&lt;br /&gt;
== CD4 Epitope Prediction ==&lt;br /&gt;
&lt;br /&gt;
To deimmunize the EPO protein we will use two prediction tools available on the IEDB: &amp;quot;CD4 T-cell immunogenicity prediction&amp;quot; and &amp;quot;Deimmunization&amp;quot;. On the [http://www.IEDB.org IEDB home page], under &amp;quot;Analysis Resource&amp;quot; select &amp;quot;T Cell Epitope Prediction&amp;quot;. You should see links to the two tools with a paragraph explaining how they work. Read this carefully!&lt;br /&gt;
&lt;br /&gt;
To gain confidence in the validity of the predictors, let us plot together the predicted immunogenicity and the experimentally measured response frequency.&lt;br /&gt;
&lt;br /&gt;
Using the &amp;quot;CD4 T-cell immunogenicity prediction&amp;quot; tool, upload your processed EPO fasta file and run the prediction with the following parameters, &amp;quot;Prediction method&amp;quot;: &amp;quot;IEDB recommended (combined)&amp;quot; and &amp;quot;Select maximum percentile rank threshold&amp;quot;:&amp;quot;Show all peptides&amp;quot;. This will start a prediction on 15-mers, overlapping by 10-mers, from the EPO protein and output results from two prediction methods(and their combined score). This will take a few seconds and once done, save the resulting .csv file.&lt;br /&gt;
&lt;br /&gt;
Now we have two .csv files to compare, one with experimental data and the other with predictions. Import these with your program of choice (Python, R, Excel,...) so that you can manipulate, compare and plot the data.&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q4: Compare the &amp;quot;Mapped Postion&amp;quot; column in the experimental result table to the &amp;quot;Start&amp;quot; and &amp;quot;End&amp;quot; columns in the prediction result table. Why are they different? Create a new &amp;quot;Start_mod&amp;quot; column in the prediction table so that it aligns with the first value in the &amp;quot;Mapped Position&amp;quot; column.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q5: The two tables have a different number of peptides. Why do you think that is? Hint: look at newly created &amp;quot;Start_mod&amp;quot; column in the prediction result table. &amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
Remove the row with peptides that are not shared between the tables.&lt;br /&gt;
&lt;br /&gt;
Now, the prediction scores and the response frequency have different units, so to make them comparable, we must process them. In the prediction score table, the &amp;quot;Combined score&amp;quot; column shows the weighed score of two models. It has the unit of rank, which can be interpreted as &amp;quot;where does the prediction score rank among prediction scores of random peptides&amp;quot; i.e. it is a percentage of random peptides with lower prediction score. That means that a high rank indicates low predicted immunogenicity, which is the reverse of the response frequency. Create a new column (e.g. &amp;quot;Combined score mod&amp;quot;) where the &amp;quot;Combined score&amp;quot; become comparable to response frequency.&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q6: What mapping function would you use to make the &amp;quot;Combined score&amp;quot; fall in the same range(and directionality) as response frequency?&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
Now create a line plot and a scatter plot comparing &amp;quot;Combined score mod&amp;quot; with the &amp;quot;Response Freq.&amp;quot; columns.&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q7: Does the predicted immunogenicity align with the measured epitopes?&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q8: What is the Pearson correlation coefficient between &amp;quot;Combined Score Mod&amp;quot; and &amp;quot;Response Freq.&amp;quot;?&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
== Deimmunization ==&lt;br /&gt;
&lt;br /&gt;
Now that we have some confidence in the prediction tools, let us deimmunize EPO.&lt;br /&gt;
&lt;br /&gt;
Open the &amp;quot;Deimmunization&amp;quot; tool on the IEDB and input your processed EPO .fasta file. Run the program with default parameters. This will predict and present the most immunogenic peptides in the input sequence. Select the top 5 peptides, enter a job name and your e-mail to receive your results when they are ready. Start running the job.&lt;br /&gt;
&lt;br /&gt;
This takes quite a lot longer than performing a CD4 epitope prediction of the full EPO sequence, despite the same CD4 epitope prediction tools being used.&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q9: Why does a deimmunization prediction for just a few epitopes take longer than a CD4 epitope prediction for the full EPO sequence?&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
This takes quite a while so I have uploaded the results of the [http://www.cbs.dtu.dk/courses/27685.imm/exercise_deimmune/deimmunization_EPO.csv deimmunization]. Download the result and inspect. Sort the table based on drop in immunogenicity to find the most influential mutations. Have a look at the suggested amino acid substitutions.&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q10: Using the BLOSUM substitution matrix as a guide, do you expect these top suggested mutations to be structurally/functionally conservative?&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
Select an influential mutation (perhaps even several) and create a new EPO .fasta file where you substitute the wild type variant with the less immunogenic mutant. Input this fasta into the CD4 epitope prediction tool and run the prediction with the same parameters as in the &amp;quot;CD4 Epitope Prediction&amp;quot; section. Download the results when ready.&lt;br /&gt;
&lt;br /&gt;
Now, after going through the same processing steps as in the &amp;quot;CD4 Epitope Prediction&amp;quot; section, plot the predicted immunogenicity for your mutant together with the predicted immunogenicity of the wild type.&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q11: Do you notice a drop in predicted immunogenicity?&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
&lt;br /&gt;
You now have tools to improve biologics!&lt;br /&gt;
&lt;br /&gt;
Done!&lt;br /&gt;
&lt;br /&gt;
&amp;lt;!--&lt;br /&gt;
== Answers ==&lt;br /&gt;
&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q1: The highest response frequency is in the region 127-141. Lesser immunogenicity in regions 160-175, 70-85 and 99-110.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q2: EPO is an extracellular protein that must be excreted from the cell. The first 27 amino acids are a signal peptide that codes for the excretion.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q3: Under the section &amp;quot;PTM/Processing&amp;quot;it is shown that residues 1-27 belong to a signal peptide.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q4: This comes back to the 27 residue signal peptide. The immunome browser starts counting from the whole protein sequence, but we input the EPO sequence without the signal peptide to the predictor.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q5: The prediction algorithm processes 15-mers, overlapping by 10-mers, and so it cuts the input sequence into fragments of this length until it reaches the end. If the length of the input is not divisible by 5 the whole sequence will not be covered by this, so an extra peptide must be included. The researchers who performed the immunoassay opted to skip the second to last peptide, since it is likely redundant.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q6: &amp;quot;Combined score mod&amp;quot;=1-&amp;quot;Combined score&amp;quot;/100&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q7: Roughly, yes. All measured epitopes have corresponding peaks in the prediction scores&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q8: The Pearson correlation coefficient is 0.696&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q9: A brute force approach entails prediction for 19*15 mutant peptide variants for each peptide(19 site variants, 15 sites in an epitope). This is further complicated by the fact that the mutation should not introduce a novel epitope in neighbouring positions. This makes for p*19*15*3 where p is the number of peptides to deimmunize and 3 accounts for neighbour mutations. EPO contains ~30 overlapping 15-mers which means that ~30 times as many predictions are involved in the deimmunization computation, per epitope. The challenge grows combinatorially with each added mutation. Applying a structutral conservation restriction, with e.g. BLOSUM would limit the mutants we would have to predict and therefore speed up the computation.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q10: No, these mutations are not conservative. This tool does not take structural aspects into consideration.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
*&amp;#039;&amp;#039;&amp;#039;Q11: Yes! There should be a drop in predicted immunogenicity. It might be slight, but you could introduce more mutations if you feel up for it.&amp;#039;&amp;#039;&amp;#039;&lt;br /&gt;
--&amp;gt;&lt;/div&gt;</summary>
		<author><name>WikiSysop</name></author>
	</entry>
</feed>