Mathematical modeling of the immune system recognition to mammary carcinoma antigen
 Carlo Bianca^{1},
 Ferdinando Chiacchio^{2},
 Francesco Pappalardo^{3} and
 Marzio Pennisi^{4}Email author
https://doi.org/10.1186/1471210513S17S21
© Bianca et al.; licensee BioMed Central Ltd. 2012
Published: 13 December 2012
Abstract
The definition of artificial immunity, realized through vaccinations, is nowadays a practice widely developed in order to eliminate cancer disease. The present paper deals with an improved version of a mathematical model recently analyzed and related to the competition between immune system cells and mammary carcinoma cells under the action of a vaccine (Triplex). The model describes in detail both the humoral and cellular response of the immune system to the tumor associate antigen and the recognition process between B cells, T cells and antigen presenting cells. The control of the tumor cells growth occurs through the definition of different vaccine protocols. The performed numerical simulations of the model are in agreement with in vivo experiments on transgenic mice.
Background
The immune system (IS) is a complex system of organs, cells and molecules whose main task is to protect living organisms from external pathogens such as viruses and bacteria. Nevertheless the effectiveness of the IS in tumors disease is nowadays under discussion among biologists and physicians. As stated by the immunosurveillance theory [1, 2], biotechnology engineered naked mice (mice without immune system) show the developing of multiple variants of malignant tumors that are not usually visible in wild mice, thus suggesting that the immune system plays an important role also against tumors. Indeed most mutated malignant cells are recognized and eliminated by immune system mechanisms right after their birth, and tumors that usually arise are indeed poorly immunogenic tumors, originating from malignant cells which escape from immune surveillance. Some tumors are caused by exogenous factors (e.g., smoke for lung cancer), and the elimination of the exogenous cause would in theory prevent the risk of developing the tumor. However many other tumors are caused by endogenous factors and their developing cannot be easily predicted and controlled. Among human cancers, the mammary carcinoma represents a major cause of concerns in women, since it belongs to the class of endogenous cancers which escape immunosurveillance of the IS.
The risk of appearance of mammary carcinoma is usually estimated by analyzing the family history of cancer, and breast cancer screening in young women is highly recommended since the achievement of earlier diagnosis could greatly improve the outcomes of the treatment. Strong family history of cancer usually entitles higher risks of developing the tumor, thus suggesting that tumor hereditary is encoded into the DNA. Some gene tests such as the genetic screening for the BRCA genes [3] are nowadays possible and may determine the risk of cancer. Indeed the analysis of the genome of individuals will be useful to better estimate the risk of cancer.
Biologists and physicians are exploring novel immunopreventive treatments that can avoid the development of breast cancer in patients with high risks of malignant cell mutations. Among others, Lollini et al. [4] have developed a cellular vaccine, called Triplex, which is able to elicit complete prevention of mammary carcinogenesis in HER2/neu transgenic mice. Triplex combines three different elements (the tumor antigen with two adjuvants) that stimulate the immune system response with different actions [4]. Vaccine cells have been engineered to present and release the tumor associated antigen p185 (that is also the oncogene of the tumor) with the addition of Allogeneic MHCclass I molecules to easier recognizing by cytotoxic T cells. Moreover, thanks to transduction of interleukin12 genes, they release interleukin12 molecules that have a broad range of costimulatory functions in boosting the immune response against tumors.
It is worth stressing that differently from the vaccine for virus or bacteria, cancer vaccines require repeated administration for the the entire life of the host. This is due to the low antigenicity of the cancer cells, the capability to escape the immune system surveillance. Moreover present cancer immunoprevention research is unable to find better vaccines able to assure complete, longterm protection.
The repeated administration of the vaccine, realized with the aim to increase the antigenicity of the tumor associated antigens, maintains the immune system response ready against newborn cancer cells. However, even if vaccines are usually less toxic than standard drugs, uncontrolled administration of the vaccine can induce undesirable effects such as autoimmune diseases. Therefore the optimization of the vaccination protocol constitutes a fundamental and open problem.
In the in vivo experiments it is not usually possible to reach an optimum vaccination protocol that maximizes the efficacy of the tumor depletion on the one hand and minimizes the risk of side effects on the other hand, because of the large variability cases. Indeed vaccination protocols are usually determined heuristically basing on best practice and previous experience. Moreover the cost of in vivo experiments can be prohibitive.
In order to understand whether it was possible to gain complete prevention of mammary carcinogenesis with fewer injections, a (multi) agentbased model named SimTriplex [5] has been developed. It is worth noticing that SimTriplex has been also employed for other pathologies [6–10]. However agentbased models do not allow the development of asymptotic analysis of the competition and an easy investigation of parameters' space.
Different mathematical tools have been developed in order to model complex biological systems and among others, immune systemcancer competition. The most famous approach is the ODEbased model where the overall system is decomposed in different cell populations whose time evolution is depicted by solutions of a nonlinear ODE system (nonlinear terms take care of the interactions among two or more cell populations), see paper [11] for a review of ODE models available in the literature and [12] for a comparison between ODE models with and without delay.
Kinetic theory models have been also proposed for the immune systemcancer competition. These models consist in partial integrodifferential equations and allow both the modeling of proliferative/destructive events and the modeling of mutations occurring in the onset of tumor, [13]. Further modeling approaches for the immune systemcancer competition include cellular automata, agentbased models, see the recent expository paper [14].
Most of the mathematical models of the IS summarize the response of the immune system in a single population of cells, named effector cells, which perform the task of destroying cancer cells. This simplifying assumption allows to reduce the complexity of the dynamics of immune system but it neglects the recognition process that occurs among the different cells that constitutes the response of the IS to the tumor antigen.
The present paper is organized as follows: Section “The Triplex vaccine in vivo experiments” briefly deals with the phenomenological analysis of the biological system. Section “The ODEbased model” is devoted to the description of the ODEbased model. Section “Sensitivity analysis” introduces the sensitivity analysis technique. Section “Results and discussion” compares, by means of numerical simulations, the mathematical model with the computational model SimTriplex. Finally Section “Conclusions” concludes the paper with a critical analysis and research perspective of the model. For interested readers Additional File 1 presents a simplified version of the model by coupling the differential model with an algebraic model.
Materials and methods
The Triplex vaccine in vivo experiments
This section briefly deals with the in vivo experiment carried on BALBneuT neu virgin female mice groups which overexpress the activated rat HER2/neu oncogene. The description does not pretend to be exhaustive from the biological point of view but highlights the essentials of the experiments in order to motivate our study.
The Triplex vaccine has been obtained from a mammary carcinoma of a FVBneuN #202 (H2^{ q }) mouse, transgenic for the rat protooncogene cneu, and combines different stimuli:

The p185neu oncoantigen;

The H2^{ q }MHC molecules (allogeneic for H2^{ d }BALBneuT mice);

The interleukin12 (vaccine cells are engineered with the genes coding for murine IL12).
The experiment starts at the sixth week of age, where BALBneuT mice start the vaccination protocol. Mice are divided in different groups, one for control untreated group, and one for each protocol tested. All vaccine protocols that have been tested are built upon the same 4week cycle which consists in twiceweekly intra peritoneum vaccinations (Tuesday and Friday) for the first 2 weeks followed by 2 weeks of rest.
The Prophylactic, lifelong Chronic vaccination protocol of cancerprone HER2/neu transgenic mice with cells expressing HER2, allogeneic MHC antigens and IL12 demonstrated able to completely prevent the onset of mammary carcinoma. The Early vaccination protocol (which counts only three 4week cycles at the beginning of the experiment) produces a significant delay in the onset of tumors, but all mice eventually succumbe to mammary carcinoma. Other tested protocols demonstrated much less effective, with little or no gain in efficacy when compared to untreated control mice.
It is worth stressing that maximal prevention against mammary carcinoma required all the three vaccine components (HER2/neu, allogeneic MHC antigens, and IL12) and was due to the induction of both cellular and humoral immune responses. Although cellular and humoral immune responses are taken into account in the vaccine administration, the relative importance of antibody subclasses for successful cancer prevention indicates that humoral immune responses is more important than cellular responses driven by cytotoxic T cells [4].
It should be therefore clear that any vaccination protocol should be started early enough to avoid carcinoma in situ formation. On the other hand it should be advisable to minimize the number of administrations in order to both maintain complete efficacy and reduce the risk of any undesirable effect. In order to help biologists in finding better vaccination protocols, a (multi) agent model named SimTriplex has been developed in [5]. The model has been inspired by the work of Celada and Seiden [22] and uses an approach that models ab initio the interaction and diffusion kinetics of each relevant biological entity. SimTriplex has been tuned with the in vivo experiments and demonstrated able to coherently reproduce the behaviors of the entities involved in the in vivo immunoprevention experiment. In addition the use of SimTriplex as a predictive tool yielded to encouraging results [23].
The ODEbased model
The competition between immune system cells and cancer cells reminds the well known PredatorPrey (PP) model described by LotkaVolterra equations. There is a population of prey, represented by the cancer cells, with an infinite set of food resources (nutrients coming from the host blood) and, differently from the classical PP, multiple populations of predators (the effectors cells) cooperate through cellcell and cellmolecule interactions to neutralize the prey. Differently from the classical predatorprey models, predator survival does not depends on the number of prey, since predator populations exist normally and, in absence of the prey, their number oscillate around given equilibrium levels.
If cancer cells are able to escape immunosurveillance, the cancer takes over, tends to compete with the healthy cells for nutrients, and could be able to kill the host. However the immune system response can be helped in recognizing harmful cells by an external agent represented in our case by a vaccine. The induced immune response is the result of a complex network of interactions between IS cells which mainly depends from cells receptors. IS cells which present, through their receptors, specific tumor antigens can trigger a complex process whose final result is the eradication of the tumor. Specific interactions, which involve cell receptors, cannot be described with an ODE population model, see [5], however if we assume that the vaccine cells activate IS cells at given ratios we can model the subsequent immune response of activated cells.
The network of organs, cells and molecules involved in the immune system is very large. In the model of the present paper we include only the entities recognized as fundamental for the biological process. We also assume that the IScancer competition occurs only in one hybrid organ which includes all the physical involved compartments (peritoneum, mammary gland, lymph nodes and so on).
Model variables
Variable  Description  Short name 

x _{1}  Number of injected vaccine cells  VC 
x _{2}  Number of P185 tumor associated antigens  TAA 
x _{3}  Number of activated B cells  B 
x _{4}  Number of activated T helper cells  TH 
x _{5}  Number of interleukin 12 molecules  IL12 
x _{6}  Number of interleukin 2 molecules  IL2 
x _{7}  Number of released antibodies  AB 
x _{8}  Number of cancer cells  CC 
x _{9}  Number of activated cytotoxic cells  TC 
x _{0}  Number of activated antigen presenting cells  APC 
Model parameters
The model contains 44 parameters which have a specific biological meaning. These parameters, assumed as constants, modify the rate of variations of the populations due to natural death, interactions with other population and release of new quantities. Accordingly: parameters referring to natural death of entities are denoted by μ_{ i }, where i identifies the population under consideration; parameters referring to the interaction between populations i and j are identified by α_{ i j }; finally parameters referring to releasing processes are identified by γ_{ i j }, where i refers to the released entity and j to the releasing entity.
Vaccine cells
P185 tumor associated antigens
T helper cells
The key role of T helper cells is to stimulate both the humoral and cellular branches of the immune response by direct receptor binding or the releasing of specific cytokines that boost the immune response, such as interleukin 2. T helper cells are activated by specialized APC such as dendritic cells, macrophages or presenting B cells which present major histocompatibility class II/peptide complexes.
Plasma B cells
Plasma B lymphocytes may absolve to multiple functions in building the immune response chain against pathogens. In a first stage they can act as specialized antigen presenting cells, by recognizing pathogens through their specialized "Yshaped" receptors, and can then present peptidic sequences to T helper cells. As a consequence of a successful interaction with T helper cells, they can be stimulated to differentiate into plasma B cells, which release antibodies with the same receptors shape, or B memory cells, which readily act against new appearance of previously encountered pathogens. Since there is no in vivo experimental evidence of B memory cells, as also suggested by the need of a chronic vaccination to achieve complete protection against tumor onset, we decided to do not include for now memory B cells into the model [4].
Interleukin 12
Interleukin 2
Antibodies
Cancer cells
Cytotoxic T cells
Antigen presenting cells
With the term antigen presenting cells we indicate a class of different types of cells, such as dendritic cells, macrophages, but also B cells, whose focal mission is to recognize, capture, and process antigens in order to present small antigenic sequences named peptides in conjunction with MHC class molecules to both cytotoxic and helper T cells.
Model parameters
Param  Description  Value(estimate)  Ref 

μ _{1}  VC (VC) death rate  ln(2) = 9  In vivo 
α _{19}  VC killing rate by TC cells (TC) (killing)  0.001  Estimated 
α _{17}  VC killing rate by AB (AB) molecules (killing)  0.001  Estimated 
q  No. of cancer cells to inject at every vaccine administration  50  SimTriplex 
γ _{21}  released TAA (TAA) rate by killed VC  3  Estimated 
γ _{28}  released TAA rate by killed CC (CC)  3  Estimated 
μ _{2}  TAA natural degradation rate  ln(2) = 9  In vivo 
α _{20}  Binding rate between TAA and APC cells  0.0005  Estimated 
α _{27}  Binding rate between TAA and AB (IC formation)  0.00001  Estimated 
γ _{34}  plasma B cells (B) activation rate by TH cells (TH)  0.05  Estimated 
α _{36}  B stimulation rate by IL2 (IL 2)  0.0035  Estimated 
s _{3}  B duplication stimulation threshold due to IL2  400  Estimated 
μ _{3}  B cells natural death rate (half life)  ln(2) = 15  [32] 
γ _{40}  TH cells (TH) activation rate by APC (APC) cells  0.15  Estimated 
α _{46}  TH cells stimulation rate by IL2 (IL 2) (duplication)  0.009  Estimated 
s _{1}  duplication stimulation threshold due to IL2  1000  Estimated 
α _{45}  TH cells cells stimulation rate by (IL 12) IL12 (duplication)  0.009  Estimated 
s _{2}  duplication stimulation threshold due to IL12  1000  Estimated 
μ _{4}  TH cells natural death rate (half life)  ln(2) = 15  Estimated 
γ _{51}  IL12 molecules release rate by VC  10  SimTriplex 
α _{54}  absorbed IL12 rate by TH cells for mitotic signals  0.00009  Estimated 
α _{59}  absorbed IL12 rate by TC cells for mitotic signals  0.001  Estimated 
μ _{5}  IL12 molecules natural degradation rate  ln(2) = 9  [33] 
γ _{64}  IL2 release rate by TH  5  Estimated 
α _{63}  absorbed IL2 rate by B cells for mitotic signals  0.0001  Estimated 
α _{69}  absorbed IL2 rate by TC cells for mitotic signals  0.0001  Estimated 
μ _{6}  IL2 molecules natural degradation rate  ln(2) = 3  Estimated 
γ _{73}  Released AB molecules rate by B cells  3  SimTriplex 
α _{78}  AB  CC binding rate  0.0001  Estimated 
α _{71}  AB  VC binding rate  0.001  Estimated 
α _{72}  AB  TAA binding rate (IC formation)  = α_{72}   
μ _{7}  AB natural degradation rate  ln(2) = 7  [34] 
c _{ max }  CC (CC) growth saturation threshold  10^{7}  Estimated 
k  CC duplication rate  0.0226  SimTriplex 
p  No. of newborn CC due to transgenic nature of mice  3  SimTriplex 
α _{88}  CC death rate due to other IS entities  0.0000001176  Estimated 
α _{89}  CC killing rate by TC cells  0.00004  Estimated 
α _{87}  CC killing rate by AB  0.00004  Estimated 
γ _{91}  TC cells activation rate by VC  0.2  Estimated 
α _{96}  TC cells duplication rate due to IL2  0.05  Estimated 
s _{96}  duplication stimulation threshold thanks to IL2  400  Estimated 
μ _{9}  TC cells natural death rate  ln(2) = 21  [35] 
γ _{02}  APC (APC) activation rate due to TAA (TAA)  0.07  Estimated 
μ _{0}  APC natural death rate  ln(2)/15  [36] 
During the in vivo experiment, biological dynamics is observed in time slices that are not smaller than eight hours. For this reason, we set the simulation time step equal to (Δ(t) = 8 hrs). This biological motivation also determined the SimTriplex timestep. The choice of the physical timestep allows to compare the results of the two models. Both models are supposed to simulate the dynamics of entities inside a volume of 1μl, which corresponds to a small portion of mammary gland of mice.
Sensitivity analysis
In order to understand which parameter may be considered fundamental in this process, it is significant to investigate the sensitivity of the model to the alteration of the parameters. Choosing a parameter in a suitable range while retaining fixed the others, represents the classical way to do sensitivity analysis. This methodology clearly owns limitations i.e., results are strongly bounded to the values of fixed parameters, and different sets of values for the fixed parameters may entitle completely different results.
Partial rank correlated coefficients (PRCC) [25] is a statistical approach used to bypass the above mentioned limitations. It works by calculating the partial correlation on ranktransformed data between input (model parameters) and output (entities behaviors). Such a technique does not depend on the values of fixed parameters and permits to vary all the parameters at the same time, allowing to study the influence of input parameters on the model outcomes. Nevertheless the methodology can be in principle easily applied and used with any kind of continuous or discrete model.
The methodology we used to perform sensitivity analysis (LHSPRCC) is briefly described as follows. The interested reader can found more information about the methodology in [26]. Parameters space is initially sampled using a MonteCarlo technique. In this case we use a technique named LatinHypercubeSampling (LHS) [27]. The technique divides the random parameter distributions into N (where N represents the chosen sample size) equal probability intervals that are then sampled. The choice for N should be at least k + 1, where k is the number of parameters varied, but usually much larger to ensure accuracy. In our trials we set N = 1000.
After sampling an LHS matrix X of sampled parameters is built. Each row represents an unique set of variables for the model sampled without replacement.
The model is then solved for each row of X, and the model output values are stored into an output matrix Y. Each matrix is then ranktransformed (X_{ R } and Y_{ R }). X and Y can be used to calculate the Pearson correlation coefficient. X_{ R } and Y_{ R } can be used to calculate the Spearman or rank correlation coefficient (RCC) and the partial rank correlation coefficient (PRCC).
Results and discussion
One of the first problems in modeling the process was to determine how to translate the biological concept of death in mathematical/computational terms. When developing the SimTriplex model, it has been decided to stop the simulation and, therefore, to consider a mouse as dead if the total number of cancer cells reached 10^{5} cells. Over such a threshold the formation of carcinoma in situ can be considered an inevitable circumstance. Since in vivo experiments demonstrated that the vaccine progressively loses its efficacy when such an event occurs [21], this threshold represents a point of no return that halves between survival and death.
Carcinoma in situ formation entitles a lot of different processes such as formation of physical barriers around the tumor mass and vascularization processes that are not described at this stage, even because this goes beyond the scope of the model. This means that both the ODEbased model and the SimTriplex models cannot be considered accurate in describing the in vivo experiment if the cancer cells threshold is overcome. Therefore the numerical simulations presented here refer to interactions where the number of cells do not go beyond this threshold.
As previously stated, the success of failure of a treatment has been determined mainly by the survival rates of the mice involved in the experiment. Even if some measurements were made during the in vivo experiment, it was not possible to keep track of the time evolution of the involved entities. Such measurements are not possible in vivo experiments, or can be achieved just partially in vitro for multiple reasons, i.e. it is not possible to do the measure too frequently due to wetlab requirements, it is not possible to take the measure at present time with current technology, or simply because the measure entitles the need to kill the host.
One of SimTriplex main features is represented by the possibility of simulating different individuals. Tuning of free parameters has been executed in order to reproduce the same population survival curves for the vaccination protocols tested in vivo [5]. Moreover, during its tuning phase, SimTriplex entities behaviors have been accurately checked by biologists in order to verify that they were qualitatively in line with both biologists assumptions and last immunological knowledge. The use of SimTriplex as a predictive tool, in conjunction with various optimization techniques [28–30], to find better vaccination protocols showed indeed that it represented a good approximation of the in vivo experiment [23], and therefore can be used to substitute missing in vivo data.
Bearing all the above in mind, we initially checked that the mathematical model mice survivals for all the tested vaccine protocols were in tune with mean survivals showed in the in vivo experiment, obtaining a good agreement between the two experiments. For the missing in vivo data, mainly represented by entities timebehaviors, we compared ODE behaviors obtained numerically with the ones obtained by SimTriplex, highlighting similarities and differences.
We would note here that, in order to compare the results, we looked in SimTriplex for "mean virtual mouse", i.e. a mouse whose death occurs near the middle of the Kaplan Meier curves for the tested protocols.
The ODE model demonstrated able to reproduce the available in vivo experimental data, in particular the in silico mice survivals for all vaccine protocols tested were in good agreement with mean survivals showed in in vivo experiment.
Since the biological behavior of the involved entities may change in a consistent manner even from mouse to mouse, we mainly focused in qualitatively analyzing cancer cells behaviors and the response times of the principal outcomes of immune response, i.e. antibodies and cytotoxic T cells behaviors for the Chronic, Early and Untreated protocols.
The above behavior does not happen in the Early case (see Figure 4, center panel), where the vaccination protocol is only able to delay the development of the cancer, and the threshold on the number of cancer cells that entitles high risks of carcinoma in situ is reached at around at 44 weeks of age in SimTriplex, and at 47 weeks in the ODE model. In in vivo experiments the middle of the Kaplan Meyer survival curve for the early protocol is reached approximately at 52 weeks of age [4], with carcinoma in situ formation between 5 to 9 weeks earlier.
Moreover the cytotoxic T cells peaks observed in SimTriplex for both the Chronic and Early protocols (see Figure 5, left and center panels) are a lot higher than those showed by the ODE model. In in vivo experiments it was observed that antibodies covered a major role in eradicating the tumor, whereas cytotoxic activity was estimated to be of secondary importance [23]. So from this point of view the ODE model may indeed be more accurate in describing this aspect of the immune response.
We used PRCC to analize the effects of the most important input parameters which influence more the behavior of Cancer Cells. We plotted for these entities the PRCCs over the entire time course of the experiment to how the parameters sensitivity varies as the process behavior advances. The analysis has been executed by supposing that the administration of the vaccine follows the Chronic protocol. In this way it is possible to study which mechanisms mainly drive the immune response against cancer cells and which parameters should be tweaked in vivo in order to obtain a strong immune response with the minimal effort. To this end we kept constant the parameter related to the quantity of injected vaccine cells (q) and the parameters related to the tumor growth (k, p, c_{ max }).
Conclusions
The mathematical model proposed in this paper is based on nonlinear ordinary differential equations. The model simulates the competition between the immune system and the mammary carcinoma under the action of an external force field (the vaccine). Three different protocols of the vaccine have been taken into account: Untreated, Early, and Chronic. The biological role of vaccine cells, cancer cells, tumor associated antigens, plasma B cells, thymus cytotoxic lymphocytes, thymus helper lymphocytes, antibodies, interleukins 2 and 12, and antigen presenting cells has been taken into account.
Numerical simulations of the model have been performed for different vaccination protocols and results were compared with a previously developed multiagent model, called SimTriplex. For the tested vaccination protocols, the ODEbased model is able to qualitatively reproduce the time evolution not only for the number of cancer cells, but also for antibodies and cytotoxic T cells, main outcomes of humoral and cell mediated immune responses. From a quantitative point of view the mathematical model showed, respectively, a weaker and a stronger immune response of cytotoxic T cells and antibodies with respect to the SimTripex model, showing indeed better agreement with the in vivo observations and speculations.
The sensitivity analysis gave two major results. First it confirmed the major role of humoral immune response also observed in in vivo experiments [23], then showed that during later stages of the experiment antigens loose their role of activating the immune response and in some cases may negatively influence the immune response. It is then possible to conclude that a reduction of the intensity of vaccine administrations in later stages, when the immune response is already set, is advisable. This has been also highlighted in [8], where an ABM model developed to illustrate the effects of the same vaccine in cancer immunotherapy, suggested to apply the golden standard vaccination procedure (initial boost followed by sparse recalls) also to cancer vaccines.
These results are certainly useful to research activity in immunology addressed to improve the efficacy of the treatment and to modulate the activation of the immune system in order to prevent side effects such as autoimmune diseases. Of course, different choices of initial conditions and of the parameters may modify the competition dynamics.
We plan to investigate the optimal protocol using mathematical tecniques which are currently under investigation. Results will be pubblished in due course.
Declarations
Acknowledgements
The first author was partially supported by the FIRB projectRBID08PP3JMetodi matematici e relativi strumenti per la modellizzazione e la simulazione della formazione di tumori, competizione con il sistema immunitario, e conseguenti suggerimenti terapeutici. The research was also partially supported by the INDAMGNFMYoung Researchers Projectprot.n.43Mathematical Modelling for the CancerImmune System Competition Elicited by a Vaccine and by University of Catania under PRA grant.
The authors would like to thank Prof. P.L. Lollini and coworkers for supplying data and Prof. S. Motta for useful discussions.
This article has been published as part of BMC Bioinformatics Volume 13 Supplement 17, 2012: Eleventh International Conference on Bioinformatics (InCoB2012): Bioinformatics. The full contents of the supplement are available online at http://www.biomedcentral.com/bmcbioinformatics/supplements/13/S17.
Authors’ Affiliations
References
 Burnet M: Cancer: a biological approach. III. Viruses associated with neoplastic conditions. IV. Practical applications. Br Med J. 1957, 1: 841847. 10.1136/bmj.1.5023.841.PubMed CentralView ArticlePubMedGoogle Scholar
 Thomas L: Reactions to homologous tissue antigens in relation to hypersensitivity. Cellular and Humoral Aspects of the Hypersensitive States. HoebersHarper. Edited by: Lawrence HS. 1959, New York, 529532.Google Scholar
 Berry DA, Iversen , Gudbjartsson DF, Hiller EH, Garber JE, Peshkin BN, Lerman C, Watson P, Lynch HT, Hilsenbeck SG, Rubinstein WS, Hughes KS, Parmigiani G: RCAPRO validation, sensitivity of genetic testing of BRCA1/BRCA2, and prevalence of other breast cancer susceptibility genes. J Clin Oncol. 2002, 20 (11): 270112. 10.1200/JCO.2002.05.121.View ArticlePubMedGoogle Scholar
 De Giovanni C, Nicoletti G, Landuzzi L, Astolfi A, Croci S, Comes A, Ferrini S, Meazza R, Iezzi M, Di Carlo E, Musiani P, Cavallo F, Nanni P, Lollini PL: Immunoprevention of HER2/neu transgenic mammary carcinoma through an interleukin 12engineered allogeneic cell vaccine. Cancer Res. 2004, 64 (11): 40014009. 10.1158/00085472.CAN032984.View ArticlePubMedGoogle Scholar
 Pappalardo F, Castiglione F, Lollini PL, Motta S: Modelling and Simulation of Cancer Immunoprevention vaccine. Bioinformatics. 2005, 21 (12): 28912897. 10.1093/bioinformatics/bti426.View ArticlePubMedGoogle Scholar
 HallingBrown M, Pappalardo F, Rapin N, Zhang P, Alemani D, Emerson A, Castiglione F, Doroux P, Pennisi M, Miotto O, Churchill D, Rossi E, Moss DS, Sansom CE, Bernaschi M, Lefranc MP, Brunak S, Motta S, Lollini PL, Basford KE, Brusic V, Shepherd AJ: ImmunoGrid: Towards Agentbased Simulations of the Human Immune System at a Natural Scale. Philosophical Transactions A. 2010, 368 (1920): 27992815. 10.1098/rsta.2010.0067.View ArticleGoogle Scholar
 Pappalardo F, HallingBrown MD, Rapin N, Zhang P, Alemani D, Emerson A, Paci P, Duroux P, Pennisi M, Palladini A, Miotto O, Churchill D, Rossi E, Shepherd AJ, Moss DS, Castiglione F, Bernaschi M, Lefranc MP, Brunak S, Motta S, Lollini PL, Basford KE, Brusic V: ImmunoGrid, an integrative environment for largescale simulation of the immune system for vaccine discovery, design, and optimization. Briefings in Bioinformatics. 2009, 10 (3): 330340.View ArticlePubMedGoogle Scholar
 Pennisi M, Pappalardo F, Palladini A, Nicoletti G, Nanni P, Lollini PL, Motta S: Modeling the competition between lung metastases and the immune system using agents. BMC Bioinformatics. 2010, 11 (Suppl 7): S1310.1186/1471210511S7S13.PubMed CentralView ArticlePubMedGoogle Scholar
 Pennisi M, Pappalardo F, Motta S: Agent based modeling of lung metastasisimmune system competition. Lecture Notes in Computer Science. 2009, 5666: 13. 10.1007/9783642032462_1.View ArticleGoogle Scholar
 Pappalardo F, Forero IM, Pennisi M, Palazon A, Melero I, Motta S: Simb16: Modeling induced immune system response against B16melanoma. PLoS ONE. 2011, 6 (10): art. no. e26523Google Scholar
 Eftimie R, Bramson JL, Earn DJD: Interactions between the immune system and cancer: A brief review of nonspatial mathematical models. Bull Math Biol. 2011, 73: 232. 10.1007/s1153801095263.View ArticlePubMedGoogle Scholar
 Baker CTH, Bocharov GA, Paul CAH: Mathematical modelling of the interleukin2 Tcell system:A comparative study of approaches based on ordinary and delay differential equations. Journal of Theoretical Medicine. 1997, 2: 117128.View ArticleGoogle Scholar
 Bianca C: Mathematical modeling for keloid formation triggered by virus: Malignant effects and immune system competition. Math Models Methods Appl Sci. 2011, 21: 389419. 10.1142/S021820251100509X.View ArticleGoogle Scholar
 Bianca C, Pennisi M: Immune system modeling by topdown and bottomup approaches. International Mathematical Forum. 2012, 7 (3): 109128.Google Scholar
 Bianca C, Pennisi M: The triplex vaccine effects in mammary carcinoma: A nonlinear model in tune with SimTriplex. Nonlinear Analysis: Real World Applications. 2012, 13: 19131940. 10.1016/j.nonrwa.2011.12.019.View ArticleGoogle Scholar
 Bianca C, Pennisi M, Motta S, Ragusa MA: Immune System Network and Cancer Vaccine. AIP Conference Proceedings. 2011, 1389: 945948.View ArticleGoogle Scholar
 Bianca C: On the modelling of space dynamics in the kinetic theory for active particles. Math Comput Modelling. 2010, 51: 7283. 10.1016/j.mcm.2009.08.044.View ArticleGoogle Scholar
 Bellouquid A, Bianca C: Modelling aggregationfragmentation phenomena from kinetic to macroscopic scales. Math Comput Modelling. 2010, 52: 802813. 10.1016/j.mcm.2010.05.010.View ArticleGoogle Scholar
 Bianca C: Kinetic theory for active particles modelling coupled to Gaussian thermostats. Applied Mathematical Sciences. 2012, 6: 651660.Google Scholar
 Bianca C: An existence and uniqueness theorem for the Cauchy problem for thermostattedKTAP models. Int Journal of Math Analysis. 2012, 6: 813824.Google Scholar
 Nanni P, Nicoletti G, Palladini A, Croci S, Murgo A, Antognoli A, Landuzzi L, Fabbi M, Ferrini S, Musiani P, Iezzi M, De Giovanni C, Lollini PL: Antimetastatic activity of a preventive cancer vaccine. Cancer Res. 2007, 67 (22): 1103711044. 10.1158/00085472.CAN072499. Erratum in: Cancer Res. 2007 Dec 15;67(24).12034View ArticlePubMedGoogle Scholar
 Seiden PE, Celada F: A Model for Simulating Cognate Recognition and Response in the Immune System. J Ther Biol. 1992, 158 (3): 329357. 10.1016/S00225193(05)807374.View ArticleGoogle Scholar
 Palladini A, Nicoletti G, Pappalardo F, Murgo A, Grosso V, Stivani V, Ianzano ML, Antognoli A, Croci S, Landuzzi L, De Giovanni C, Nanni P, Motta S, Lollini PL: In silico modeling and in vivo efficacy of cancer preventive vaccinations. Cancer Research. 2010, 70 (20): 77557763. 10.1158/00085472.CAN100701.View ArticlePubMedGoogle Scholar
 Nanni P, Nicoletti G, De Giovanni C, Landuzzi L, Di Carlo E, Cavallo F, Pupa SM, Rossi I, Colombo MP, Ricci C, Astolfi A, Musiani P, Forni G, Lollini PL: Combined allogeneic tumor cell vaccination and systemic interleukin 12 prevents mammary carcinogenesis in HER2/neu transgenic mice. J Exp Med. 2001, 194 (9): 11951205. 10.1084/jem.194.9.1195.PubMed CentralView ArticlePubMedGoogle Scholar
 Saltelli A: Sensitivity Analysis in Practice: A Guide to Assessing Scientific Models. 2008, Hoboken, NJ: WileyGoogle Scholar
 Marino S, Hogue IB, Ray CJ, Kirschner DE: A methodology for performing global uncertainty and sensitivity analysis in systems biology. Journal of Theoretical Biology. 2008, 254: 178196. 10.1016/j.jtbi.2008.04.011.PubMed CentralView ArticlePubMedGoogle Scholar
 Mckay M, Beckman R, Conover W: Comparison of 3 methods for selecting values of input variables in the analysis of output from a computer code. Technometrics. 1979, 21: 239245.Google Scholar
 Pappalardo F, Mastriani E, Lollini PL, Motta S: Genetic Algorithm against Cancer. Lecture Notes in Computer Science. 2006, 3849: 223228. 10.1007/11676935_27.View ArticleGoogle Scholar
 Pappalardo F, Pennisi M, Castiglione F, Motta S: Vaccine protocols optimization: in silico experiences. Biotechnology Advances. 2010, 28: 8293. 10.1016/j.biotechadv.2009.10.001.View ArticlePubMedGoogle Scholar
 Pennisi M, Catanuto R, Pappalardo F, Motta S: Optimal vaccination schedules using Simulated Annealing. Bioinformatics. 2008, 24 (15): 17401742. 10.1093/bioinformatics/btn260.View ArticlePubMedGoogle Scholar
 Abbas AK, Litchman AH, Pilllai S: Cellular and molecular immunology. 2011, Elsevier, 7Google Scholar
 Mattioli CA, Tomasi TB: The lifespan of IgA plasma Cells From the Mouse intestine. J Exp Med. 1973, 138 (2): 452460. 10.1084/jem.138.2.452.PubMed CentralView ArticlePubMedGoogle Scholar
 Remick DG, Friedland JS: Cytokines in Health and Disease. 1997, Marcel Dekker, New York, 2nd Review and ExpandGoogle Scholar
 Peppard JV, Orlans E: The biological halflives of four rat immunoglobulin isotypes. Immunology. 1980, 40: 683686.PubMed CentralPubMedGoogle Scholar
 Pardoll D: T cells take aim at cancer. PNAS. 2002, 99 (25): 1584015842. 10.1073/pnas.262669499.PubMed CentralView ArticlePubMedGoogle Scholar
 Zehn D, Cohen CJ, Reiter Y, Walden P: Extended presentation of specific MHCpeptide complexes by mature dendritic cells compared to other types of antigenpresenting cells. European Journal of Immunology. 2004, 34 (6): 15511560. 10.1002/eji.200324355.View ArticlePubMedGoogle Scholar
Copyright
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.