VaccImm: simulating peptide vaccination in cancer therapy
- Joachim von Eichborn†1,
- Anna Lena Woelke†1, 2,
- Filippo Castiglione3 and
- Robert Preissner1Email author
© von Eichborn et al.; licensee BioMed Central Ltd. 2013
Received: 27 August 2012
Accepted: 19 March 2013
Published: 15 April 2013
Despite progress in conventional cancer therapies, cancer is still one of the leading causes of death in industrial nations. Therefore, an urgent need of progress in fighting cancer remains. A promising alternative to conventional methods is immune therapy. This relies on the fact that low-immunogenic tumours can be eradicated if an immune response against them is induced. Peptide vaccination is carried out by injecting tumour peptides into a patient to trigger a specific immune response against the tumour in its entirety. However, peptide vaccination is a highly complicated treatment and currently many factors like the optimal number of epitopes are not known precisely. Therefore, it is necessary to evaluate how certain parameters influence the therapy.
We present the VaccImm Server that allows users to simulate peptide vaccination in cancer therapy. It uses an agent-based model that simulates peptide vaccination by explicitly modelling the involved cells (immune system and cancer) as well as molecules (antibodies, antigens and semiochemicals). As a new feature, our model uses real amino acid sequences to represent molecular binding sites of relevant immune cells. The model is used to generate detailed statistics of the population sizes and states of the single cell types over time. This makes the VaccImm web server well suited to examine the parameter space of peptide vaccination in silico. VaccImm is publicly available without registration on the web at http://bioinformatics.charite.de/vaccimm; all major browsers are supported.
The VaccImm Server provides a convenient way to analyze properties of peptide vaccination in cancer therapy. Using the server, we could gain interesting insights into peptide vaccination that reveal the complex and patient-specific nature of peptide vaccination.
KeywordsSystems biology Immunoinformatics Cancer Modelling Proteins Protein interaction
To mount an immune response against tumours is the objective of immunotherapy and constitutes a promising alternative to conventional methods . However, the applicability of this approach remains limited so far, because of largely undefined treatment parameters. For instance, the success of immunotherapy depends strongly on the cancer epitope sequence, the individual major histocompatibility complex (MHC) genotype of the individual, and on the overall immune response dynamics. To explore this parameter space systematically, we have developed an agent-based model simulating a certain line of immunotherapy called peptide vaccination . Our model VaccImm has the ability to take the MHC genotype, the amino acid epitope sequences and the spatial cell dynamics into account.
Here, we present a server that allows users to run these simulations online. Several parameters, like the cancer antigens and the MHC genotype of the virtual individual, can be selected for the simulation.
Agent-based models are well suited to model the immune system. They consist of independent agents, each represented with individual properties and behavioural rules. For example, every single cell is localized within a given simulated volume, has a specific type, clonotype (in case of lymphocytes) and developmental stage. This allows agent-based models to simulate the immune system, in which the different cells and cell types may be distributed very inhomogenously .
Back in 1992, Celada and Seiden published a cellular automaton called ImmSim to model cellular interactions in the immune system . In this model, the interactions of immune cells were rather simplistically modelled by bit string complementarity. Nevertheless, this automaton was able to reproduce basic immune properties like clonal expansion of B-cells and T-cells after stimulation as well as more advanced immune system traits like the competition between cross reacting clones . Since 1992, several other rule based models have been developed in order to simulate the immune system [6, 7]. The original ImmSim has been upgraded and extended in the meantime, leading to different forks like C-ImmSim  and the mice-specific SimTriplex .
C-ImmSim is a refined version of the original ImmSim model that was ported to ANSI C language. However, the interaction between immune cells is still modelled as bit strings in C-ImmSim, thus having no direct translation to biologically meaningful amino acid sequences. Therefore, Rapin et al. developed an extension of C-ImmSim simulating the branch of bacterial infection [10, 11]. As immune reactions against cancer cells follow different rules, they left out this part of C-ImmSim.
In a completely independent approach, we have created VaccImm - another extension of C-ImmSim - that now simulates immune reactions against cancer using amino acid sequences and knowledge-based interaction potentials to predict cell interaction. The parameters used in VaccImm were carefully examined elsewhere .
VaccImm models a three-dimensional Cartesian lattice which contains the simulated cells (cancer cells, helper T-cells, cytotoxic T-cells, B-cells, dendritic cells, macrophages) and molecules (antigens, antibodies, interleukin-2 and a danger signal that acts as a general activator of macrophages). These cells and molecules can interact and move within the lattice according to their behavioural rules. The model simulates an adaptive immune response against cancer. The thymus selection is modelled implicitly. All T-cells are checked for reactivity against MHC and self-peptides before being introduced in the simulation. Only cells having sufficiently high reactivity against the own MHC complexes (positive selection) and sufficiently low reactivity against self-peptides (negative selection) will be able to enter the simulation.
As one input parameter for the model, a certain antigenic sequence is given. Peptides for injection originating from this protein sequence are predicted for their binding capability to MHC I or MHC II using acknowledged prediction algorithms (consensus  for MHC I, smm_align  for MHC II) from the Immune Epitope Database .
To assess the binding probability of a receptor-ligand pair in amino acid dependent manner, knowledge-based interaction potentials are a convenient solution. We developed different interaction potentials for B-cell receptors and T-cell receptors to be used in VaccImm and we have shown that they are able to clearly distinguish between random complexes and experimentally observed ones (details are described elsewhere ).
The simulation is carried out in discrete time steps, each corresponding to eight hours of real life, which corresponds to one cell division cycle. In each time step the cells can interact, move and generally follow their behavioural properties with respect to their environment.
Injection of peptides, either emulsified in adjuvant or not. Adjuvant is modelled using the general danger signal mentioned before.
Antigen presenting cells take up the peptides. In the case of macrophages and dendritic cells, phagocytosis happens in an unspecific manner whereas B-cells take up only antigens they recognize. Whether antigens are recognized by B-cells or not is computed using the B-cell interaction potentials.
MHC I and II processing and presentation of peptides by antigen presenting cells.
If a peptide presented on MHC I or MHC II is recognized by cytotoxic T-cells (TC) or helper T-cells (TH), respectively, T-cells get activated, duplicate and create memory cells. Whether recognition takes place or not is assessed using the T-cell interaction potentials.
Both humoral and cytotoxic immune response is initiated. Activated TCs kill the cells they recognize, thereby eliminating the tumour. Recognition is again computed using the T-cell interaction potential. B-cells that are stimulated by THs duplicate into memory cells and plasma cells. These plasma cells produce antibodies that clear the antigen.
Results and discussion
To start a new simulation, the user can decide, whether urological cancer (kidney or prostate) or another cancer type should be simulated. If one of the former types is chosen, the user can select the antigen from a list of proteins that are over-expressed in this tumour type . Alternatively, a protein-sequence can be entered.
In the next step, the MHC-genotype of the virtual individual can be selected from a list of all human MHC-alleles. With this information, the strongest MHC-binding peptides are successively computed using prediction tools from the Immune Epitope Database . These are the ones used for vaccination; the user can then set how many different peptides are to be injected.
Finally, some general parameters of the simulation can be changed, i.e., the simulated volume, the duration of the simulation and the initial tumour size.
The simulation results consist of the detailed dynamics of immune and cancer cell populations over time. They are presented as a number of interactive plots that enable the user to hide and show individual data lines in the plot and to get the exact cell population size at any given time. Each plot is annotated with a description that states what it is showing. Furthermore, the plots can be downloaded as image files and the raw simulation output data is accessible.
In the next sections, we will analyze how strong the influences of the above-mentioned parameters are on the outcome of the simulation.
Dependence on the type of antigen
Clinical studies of peptide vaccination in cancer treatment usually use peptides from several different cancer proteins in the hope that the immune system will react to at least one of them . The problem of finding an immunogenic epitope able to induce an effective immune response in vivo is one of the critical steps in immune therapy.
Dependence on the MHC-genotype
The peptide binding characteristics of MHC alleles differ significantly and susceptibility to several diseases has been associated with the MHC-genotype [19, 20]. Therefore, we wanted to investigate the influence of the MHC-genotype on the success of peptide vaccination therapy.
This demonstrates, once more, the complex nature of the immune response against tumours. Injecting peptides from the same antigen might induce a strong immune response having a certain MHC-allele combination, while the response could be completely absent when having a different allele combination, as can be seen in Figure 2. It is worth emphasizing that in VaccImm, the injected peptides are chosen in accordance with the MHC alleles. Therefore, a possible explanation for the different responses to the same antigen presented on different MHC-alleles could be that the peptide sequences binding to one MHC are different from the peptides binding to other MHCs, which results in a difference in immunogenicity of these different peptides.
This experiment is an indication that peptide vaccination remains a personalized treatment because injecting the peptides derived from the same antigen might have a completely different outcome depending on the individual MHC-alleles present.
Dependence on the number of injected peptides
When a possible cancer target is found by expression analysis, the respective peptide presented on MHC I or MHC II is often unknown. Consequently, prediction algorithms are usually used to find cancer epitopes presented on MHCs. Likewise, we have used this kind of algorithm, bearing in mind that they are of limited accuracy  and the highest ranked peptide could turn out not to be the one eventually presented by APCs. Therefore, we investigated whether injecting more peptides from the same antigen increases the chance of successful treatment.
It seems that a strong responsiveness of TCs can complement a weaker TH responsiveness and vice versa; for example, injecting one MHC II and nine MHC I peptides or one MHC I and nine MHC II peptides induces an equally strong immune response whereas injecting one peptide for both MHC types does not. This observation highlights the interconnection of the immune cells. TCs need the help of THs to become activated, while their reactivity depends on the immunogenicity and binding properties of MHC I peptides, but also on the reactivity of THs. In some cases, the reactivity of TCs and THs may complement each other, although the actual mechanisms behind are more complex than simple dose-dependent reactions.
It should be mentioned that this dependency on the number of injected peptides is not observed for all antigens and MHC-allele combinations. Some antigens fail to induce an immune response regardless of the number of injected peptides whereas others lead to tumour eradication even when only one peptide for each MHC type is injected (data not shown).
The observation above suggests using as many peptides as possible in planning a clinical study. However, this is not practically because every injected peptide bears the risk of inducing an autoimmune disease. What must be carefully taken into account is which peptides are likely to induce a strong favoured response against the tumour but a low or absent response against self proteins.
Dependence on the initial tumour size
Cancer dormancy describes the phenomenon that small tumours can be kept in check without treatment for a longer period of time ; in contrast, larger tumours often grow exponentially. In general, the chosen therapy depends on the tumour size. For example, for a prostate tumour of a considerable size, radical prostatectomy is often the treatment of choice, whereas smaller tumours are frequently kept under active surveillance for a long time . It is therefore of special interest to analyze how the initial tumour size influences the success rate of the immune therapy.
The results in Figure 5 indicate that if the immune system is activated sufficiently against a tumour, eradication is only a matter of time. Yet, the problem with cancer cells is that they are usually composed of an ever growing diversity of cells with different expression patterns; a fact that is not considered in our model so far. Thus, a larger tumour has more time to mutate or to change the expression pattern and hence a cure with peptide vaccination is less likely to have a positive effect.
Our future plan is to include the diversity of cells, since without this parameter we cannot deduce from the results of this experiment whether the initial tumour size is an indication of success for peptide vaccination therapy.
We have presented VaccImm, a user-friendly server to simulate the effect of peptide vaccination in cancer therapy. We believe this tool is very useful for analyzing parameters of peptide vaccination in cancer therapy. It is the first web server that can model cancer immunotherapy based on cancer epitope sequences and MHC genotypes. All sequence-related parameters can be selected by the user, along with parameters concerning size and duration of the simulation. The output shows detailed cell population statistics for all modelled cell types. Their population sizes and activation states can be studied over time.
As stated in the introduction, VaccImm can be very useful to explore the parameter space of peptide vaccination. For instance, we have performed some studies analyzing the influence of several parameters on the success of peptide vaccination. The results presented above revealed the complex nature of peptide vaccination as the influences of the single parameters strongly depend upon each other. The selection of the antigen is crucial for peptide vaccination, because the antigen must be able to evoke an immune response. However, the quality of a potential antigen in turn depends on the MHC-genotype of the patient. This is because the antigens are processed and subsequently different MHCs may present different parts of the antigens. VaccImm accounts for this by predicting the MHC binding peptides based on the MHC-genotype. It could also be shown, that a successful therapy gets more likely the more different peptides are injected, still bearing the risk of an autoimmune reaction. One feature that is currently not included in our model is mutation of cancer cells. Therefore, the initial size of the tumour has no qualitative effect on the outcome of the simulation. If mutation of tumour cells was possible, one would expect a higher probability for a cell in a large tumour to mutate in such a way that it is not or less affected by the evoked immune response. Thereby the probability of a successful treatment can be expected to become smaller with increasing tumour size. Due to the modular structure of our model, it is easily possible to add such behaviour in a future version.
Availability and requirements
Project name: VaccImm
Project home page: http://bioinformatics.charite.de/vaccimm/
Operating systems: Platform independent
Programming language: C and PHP
License: VaccImm is available free of charge without registration
Any restrictions to use by non-academics: None
Major histocompatibility complex
JvE was funded by the international research training group “Genomics and Systems Biology of Molecular Networks” (GRK1360 of the DFG). JvE and AW were funded by the research training group “Computational Systems Biology” (GRK1772/1 of the DFG). The project was supported by the MedSys programme (BMBF).
- Palucka K, Ueno H, Banchereau J: Recent developments in cancer vaccines. J Immunol. 2011, 186: 1325-1331. 10.4049/jimmunol.0902539.PubMed CentralView ArticlePubMedGoogle Scholar
- Woelke AL, Von Eichborn J, Murgueitio MS, Worth CL, Castiglione F, Preissner R: Development of immune-specific interaction potentials and their application in the multi-agent-system VaccImm. PLoS One. 2011, 6: e23257-10.1371/journal.pone.0023257.PubMed CentralView ArticlePubMedGoogle Scholar
- Louzoun Y, Solomon S, Atlan H, Cohen IR: Modeling complexity in biology. Physica A: Stat Mech and its Appl. 2001, 297: 242-252. 10.1016/S0378-4371(01)00201-1.View ArticleGoogle Scholar
- Celada F, Seiden PE: A computer model of cellular interactions in the immune system. Immunol Today. 1992, 13: 56-62. 10.1016/0167-5699(92)90135-T.View ArticlePubMedGoogle Scholar
- Cheng Y, Ghersi D, Calcagno C, Selin LK, Puzone R, Celada F: A discrete computer model of the immune system reveals competitive interactions between the humoral and cellular branch and between cross-reacting memory and naïve responses. Vaccine. 2009, 27: 833-845. 10.1016/j.vaccine.2008.11.109.PubMed CentralView ArticlePubMedGoogle Scholar
- Folcik VA, An GC, Orosz CG: The Basic Immune Simulator: an agent-based model to study the interactions between innate and adaptive immunity. Theor Biol Med Model. 2007, 4: 39-10.1186/1742-4682-4-39.PubMed CentralView ArticlePubMedGoogle Scholar
- Mallet DG, De Pillis LG: A cellular automata model of tumor-immune system interactions. J Theor Biol. 2006, 239: 334-350. 10.1016/j.jtbi.2005.08.002.View ArticlePubMedGoogle Scholar
- Castiglione F, Toschi F, Bernaschi M, Succi S, Benedetti R, Falini B, Liso A: Computational modeling of the immune response to tumor antigens. J Theor Biol. 2005, 237: 390-400. 10.1016/j.jtbi.2005.04.024.View ArticlePubMedGoogle Scholar
- Pappalardo F, Lollini P-L, Castiglione F, Motta S: Modeling and simulation of cancer immunoprevention vaccine. Bioinformatics. 2005, 21: 2891-2897. 10.1093/bioinformatics/bti426.View ArticlePubMedGoogle Scholar
- Rapin N, Lund O, Castiglione F: Immune System Simulation Online. Bioinformatics. 2011, 27: 2013-2014. 10.1093/bioinformatics/btr335.View ArticlePubMedGoogle Scholar
- Rapin N, Lund O, Bernaschi M, Castiglione F: Computational immunology meets bioinformatics: the use of prediction tools for molecular binding in the simulation of the immune system. PLoS One. 2010, 5: e9862-10.1371/journal.pone.0009862.PubMed CentralView ArticlePubMedGoogle Scholar
- Moutaftsi M, Peters B, Pasquetto V, Tscharke DC, Sidney J, Bui H-H, Grey H, Sette A: A consensus epitope prediction approach identifies the breadth of murine T(CD8+)-cell responses to vaccinia virus. Nat Biotechnol. 2006, 24: 817-819. 10.1038/nbt1215.View ArticlePubMedGoogle Scholar
- Nielsen M, Lundegaard C, Lund O: Prediction of MHC class II binding affinity using SMM-align, a novel stabilization matrix alignment method. BMC Bioinforma. 2007, 8: 238-10.1186/1471-2105-8-238.View ArticleGoogle Scholar
- Peters B, Sidney J, Bourne P, Bui H-H, Buus S, Doh G, Fleri W, Kronenberg M, Kubo R, Lund O, Nemazee D, Ponomarenko JV, Sathiamurthy M, Schoenberger S, Stewart S, Surko P, Way S, Wilson S, Sette A: The immune epitope database and analysis resource: from vision to blueprint. PLoS Biol. 2005, 3: e91-10.1371/journal.pbio.0030091.PubMed CentralView ArticlePubMedGoogle Scholar
- Jones J, Otu H, Spentzos D, Kolia S, Inan M, Beecken WD, Fellbaum C, Gu X, Joseph M, Pantuck AJ, Jonas D, Libermann TA: Gene signatures of progression and metastasis in renal cell cancer. Clin Cancer Res. 2005, 11: 5730-5739. 10.1158/1078-0432.CCR-04-2225.View ArticlePubMedGoogle Scholar
- Vita R, Zarebski L, Greenbaum JA, Emami H, Hoof I, Salimi N, Damle R, Sette A, Peters B: The immune epitope database 2.0. Nucleic Acids Res. 2010, 38: D854-D862. 10.1093/nar/gkp1004.PubMed CentralView ArticlePubMedGoogle Scholar
- Stevanovic S: Identification of tumour-associated T-cell epitopes for vaccine development. Nat Rev Cancer. 2002, 2: 514-520. 10.1038/nrc841.View ArticlePubMedGoogle Scholar
- Nagy G, Emody L, Pál T: Strategies for the development of vaccines conferring broad-spectrum protection. Int J Med Microbiol. 2008, 298: 379-395. 10.1016/j.ijmm.2008.01.012.View ArticlePubMedGoogle Scholar
- Ridgway WM, Fathman CG: The association of MHC with autoimmune diseases: understanding the pathogenesis of autoimmune diabetes. Clin Immunol Immunopathol. 1998, 86: 3-10. 10.1006/clin.1997.4449.View ArticlePubMedGoogle Scholar
- Riley EM: The role of MHC- and non-MHC-associated genes in determining the human immune response to malaria antigens. Parasitology. 1996, 112 (Suppl): S39-S51.PubMedGoogle Scholar
- Nielsen M, Lund O, Buus S, Lundegaard C: MHC class II epitope predictive algorithms. Immunology. 2010, 130: 319-328. 10.1111/j.1365-2567.2010.03268.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Uhr JW, Scheuermann RH, Street NE, Vitetta ES: Cancer dormancy: opportunities for new therapeutic approaches. Nat Med. 1997, 3: 505-509. 10.1038/nm0597-505.View ArticlePubMedGoogle Scholar
- Rayala HJ, Richie JP: Radical prostatectomy reigns supreme. Oncology (Williston Park). 2009, 23: 863-867.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.