Modeling and simulation of cancer immunoprevention vaccine

From Infection to Immunity: Understanding the Response to SARS-CoV2 Through In-Silico Modeling

Background

Immune system conditions of the patient is a key factor in COVID-19 infection survival. A growing number of studies have focused on immunological determinants to develop better biomarkers for therapies.

Aim

Studies of the insurgence of immunity is at the core of both SARS-CoV-2 vaccine development and therapies. This paper attempts to describe the insurgence (and the span) of immunity in COVID-19 at the population level by developing an in-silico model. We simulate the immune response to SARS-CoV-2 and analyze the impact of infecting viral load, affinity to the ACE2 receptor, and age in an artificially infected population on the course of the disease.

Methods

We use a stochastic agent-based immune simulation platform to construct a virtual cohort of infected individuals with age-dependent varying degrees of immune competence. We use a parameter set to reproduce known inter-patient variability and general epidemiological statistics.

Results

By assuming the viremia at day 30 of the infection to be the proxy for lethality, we reproduce in-silico several clinical observations and identify critical factors in the statistical evolution of the infection. In particular, we evidence the importance of the humoral response over the cytotoxic response and find that the antibody titers measured after day 25 from the infection are a prognostic factor for determining the clinical outcome of the infection. Our modeling framework uses COVID-19 infection to demonstrate the actionable effectiveness of modeling the immune response at individual and population levels. The model developed can explain and interpret observed patterns of infection and makes verifiable temporal predictions. Within the limitations imposed by the simulated environment, this work proposes quantitatively that the great variability observed in the patient outcomes in real life can be the mere result of subtle variability in the infecting viral load and immune competence in the population. In this work, we exemplify how computational modeling of immune response provides an important view to discuss hypothesis and design new experiments, in particular paving the way to further investigations about the duration of vaccine-elicited immunity especially in the view of the blundering effect of immunosenescence.

A mathematical model for the quantification of a patient's sensitivity to checkpoint inhibitors and long-term tumour burden

A large proportion of patients with cancer are unresponsive to treatment with immune checkpoint blockade and other immunotherapies. Here, we report a mathematical model of the time-course of tumour responses to immune-checkpoint inhibitors. The model takes into account intrinsic tumour-growth rates, the rates of immune activation and of tumour–immune-cell interactions, and the efficacy of immune-mediated tumour killing. For 124 patients, four cancer types and two immunotherapy agents, the model reliably described the immune responses and final tumour burden across all different cancers and drug combinations examined. In validation cohorts from four clinical trials of checkpoint inhibitors (with a total of 177 patients), the model accurately stratified the patients according to reduced or increased long-term tumour burden. We also provide model-derived quantitative measures of treatment sensitivity for specific drug–cancer combinations. The model can be used to predict responses to therapy and to quantify specific drug–cancer sensitivities in individual patients.

Computational modelling of modern cancer immunotherapy

Modern cancer immunotherapy has revolutionised oncology and carries the potential to radically change the approach to cancer treatment. However, numerous questions remain to be answered to understand immunotherapy response better and further improve the benefit for future cancer patients. Computational models are promising tools that can contribute to accelerated immunotherapy research by providing new clues and hypotheses that could be tested in future trials, based on preceding simulations in addition to the empirical rationale. In this topical review, we briefly summarise the history of cancer immunotherapy, including computational modelling of traditional cancer immunotherapy, and comprehensively review computational models of modern cancer immunotherapy, such as immune checkpoint inhibitors (as monotherapy and combination treatment), co-stimulatory agonistic antibodies, bispecific antibodies, and chimeric antigen receptor T cells. The modelling approaches are classified into one of the following categories: data-driven top-down vs mechanistic bottom-up, simplistic vs detailed, continuous vs discrete, and hybrid. Several common modelling approaches are summarised, such as pharmacokinetic/pharmacodynamic models, Lotka-Volterra models, evolutionary game theory models, quantitative systems pharmacology models, spatio-temporal models, agent-based models, and logic-based models. Pros and cons of each modelling approach are critically discussed, particularly with the focus on the potential for successful translation into immuno-oncology research and routine clinical practice. Specific attention is paid to calibration and validation of each model, which is a necessary prerequisite for any successful model, and at the same time, one of the main obstacles. Lastly, we provide guidelines and suggestions for the future development of the field.

Potential predictors of type-2 diabetes risk: machine learning, synthetic data and wearable health devices

BACKGROUND

The aim of a recent research project was the investigation of the mechanisms involved in the onset of type 2 diabetes in the absence of familiarity. This has led to the development of a computational model that recapitulates the aetiology of the disease and simulates the immunological and metabolic alterations linked to type-2 diabetes subjected to clinical, physiological, and behavioural features of prototypical human individuals.

RESULTS

We analysed the time course of 46,170 virtual subjects, experiencing different lifestyle conditions. We then set up a statistical model able to recapitulate the simulated outcomes.

CONCLUSIONS

The resulting machine learning model adequately predicts the synthetic dataset and can, therefore, be used as a computationally-cheaper version of the detailed mathematical model, ready to be implemented on mobile devices to allow self-assessment by informed and aware individuals. The computational model used to generate the dataset of this work is available as a web-service at the following address: http://kraken.iac.rm.cnr.it/T2DM .

Moving forward through the in silico modeling of tuberculosis: a further step with UISS-TB

BACKGROUND

In 2018, about 10 million people were found infected by tuberculosis, with approximately 1.2 million deaths worldwide. Despite these numbers have been relatively stable in recent years, tuberculosis is still considered one of the top 10 deadliest diseases worldwide. Over the years, Mycobacterium tuberculosis has developed a form of resistance to first-line tuberculosis treatments, specifically to isoniazid, leading to multi-drug-resistant tuberculosis. In this context, the EU and Indian DBT funded project STriTuVaD-In Silico Trial for Tuberculosis Vaccine Development-is supporting the identification of new interventional strategies against tuberculosis thanks to the use of Universal Immune System Simulator (UISS), a computational framework capable of predicting the immunity induced by specific drugs such as therapeutic vaccines and antibiotics.

RESULTS

Here, we present how UISS accurately simulates tuberculosis dynamics and its interaction within the immune system, and how it predicts the efficacy of the combined action of isoniazid and RUTI vaccine in a specific digital population cohort. Specifically, we simulated two groups of 100 digital patients. The first group was treated with isoniazid only, while the second one was treated with the combination of RUTI vaccine and isoniazid, according to the dosage strategy described in the clinical trial design. UISS-TB shows to be in good agreement with clinical trial results suggesting that RUTI vaccine may favor a partial recover of infected lung tissue.

CONCLUSIONS

In silico trials innovations represent a powerful pipeline for the prediction of the effects of specific therapeutic strategies and related clinical outcomes. Here, we present a further step in UISS framework implementation. Specifically, we found that the simulated mechanism of action of RUTI and INH are in good alignment with the results coming from past clinical phase IIa trials.

Analysis of a breast cancer mathematical model by a new method to find an optimal protocol for HER2-positive cancer

Treatment of breast cancer (positive for HER2, i.e., ERBB2) is described by a mathematical model involving non-linear ordinary differential equations with a hidden hierarchy. To reveal the hierarchy of dynamical variables of the system being considered, we applied the singular perturbed vector field (SPVF) method, where a system of equations can be decomposed to fast and slow sub-systems with explicit small parameters. This new form of the model, which is called a singular perturbed system, enables us to apply a semi-analytical method called the method of directly defining inverse mapping (MDDiM), which is based on the homotopy analysis asymptotic method. We introduced the treatment protocol in explicit form, through an analytical function that describes the exact dose and intervals between treatments in a cyclical manner. In addition, a new algorithm for the optimal dosage that causes tumour shrinkage is presented in this study. Furthermore, we took the concept of protocol optimisation a step further and derived a differential equation that represents vaccination depending on tumour size and yields an optimal protocol of different doses at every time point. We introduced the treatment protocol in explicit form, through an analytical function that describes the exact dose and intervals between treatments in a cyclical manner. In addition, a new algorithm for finding the optimal dosage that causes tumour shrinkage is presented in this study. Additionally, we took the concept of protocol optimisation a step further and derived a differential equation that represents vaccination depending on tumour size and yields an optimal protocol of different doses at every time point.

Recent advances in theranostic polymeric nanoparticles for cancer treatment: A review

Nanotheranostics is fast-growing pharmaceutical technology for simultaneously monitoring drug release and its distribution, and to evaluate the real time therapeutic efficacy through a single nanoscale for treatment and diagnosis of deadly disease such as cancers. In recent two decades, biodegradable polymers have been discovered as important carriers to accommodate therapeutic and medical imaging agents to facilitate construction of multi-modal formulations. In this review, we summarize various multifunctional polymeric nano-sized formulations such as polymer-based super paramagnetic nanoparticles, ultrasound-triggered polymeric nanoparticles, polymeric nanoparticles bearing radionuclides, and fluorescent polymeric nano-sized formulations for purpose of theranostics. The use of such multi-modal nano-sized formulations for near future clinical trials can assist clinicians to predict therapeutic properties (for instance, depending upon the quantity of drug accumulated at the cancerous site) and observed the progress of tumor growth in patients, thus improving tailored medicines.

Computer Modeling of Clonal Dominance: Memory-Anti-Naïve and Its Curbing by Attrition

Experimental and computational studies have revealed that T-cell cross-reactivity is a widespread phenomenon that can either be advantageous or detrimental to the host. In particular, detrimental effects can occur whenever the clonal dominance of memory cells is not justified by their infection-clearing capacity. Using an agent-based model of the immune system, we recently predicted the "memory anti-naïve" phenomenon, which occurs when the secondary challenge is similar but not identical to the primary stimulation. In this case, the pre-existing memory cells formed during the primary infection may be rapidly deployed in spite of their low affinity and can actually prevent a potentially higher affinity naïve response from emerging, resulting in impaired viral clearance. This finding allowed us to propose a mechanistic explanation for the concept of "antigenic sin" originally described in the context of the humoral response. However, the fact that antigenic sin is a relatively rare occurrence suggests the existence of evolutionary mechanisms that can mitigate the effect of the memory anti-naïve phenomenon. In this study we use computer modeling to further elucidate clonal dominance and the memory anti-naïve phenomenon, and to investigate a possible mitigating factor called attrition. Attrition has been described in the experimental and computational literature as a combination of competition for space and apoptosis of lymphocytes via type-I interferon in the early stages of a viral infection. This study systematically explores the relationship between clonal dominance and the mechanism of attrition. Our results suggest that attrition can indeed mitigate the memory anti-naïve effect by enabling the emergence of a diverse, higher affinity naïve response against the secondary challenge. In conclusion, modeling attrition allows us to shed light on the nature of clonal interaction and dominance.

Quantitative Mechanistic Modeling in Support of Pharmacological Therapeutics Development in Immuno-Oncology

Following the approval, in recent years, of the first immune checkpoint inhibitor, there has been an explosion in the development of immuno-modulating pharmacological modalities for the treatment of various cancers. From the discovery phase to late-stage clinical testing and regulatory approval, challenges in the development of immuno-oncology (IO) drugs are multi-fold and complex. In the preclinical setting, the multiplicity of potential drug targets around immune checkpoints, the growing list of immuno-modulatory molecular and cellular forces in the tumor microenvironment-with additional opportunities for IO drug targets, the emergence of exploratory biomarkers, and the unleashed potential of modality combinations all have necessitated the development of quantitative, mechanistically-oriented systems models which incorporate key biology and patho-physiology aspects of immuno-oncology and the pharmacokinetics of IO-modulating agents. In the clinical setting, the qualification of surrogate biomarkers predictive of IO treatment efficacy or outcome, and the corresponding optimization of IO trial design have become major challenges. This mini-review focuses on the evolution and state-of-the-art of quantitative systems models describing the tumor vs. immune system interplay, and their merging with quantitative pharmacology models of IO-modulating agents, as companion tools to support the addressing of these challenges.

Computational Approaches in Theranostics: Mining and Predicting Cancer Data

The ability to understand the complexity of cancer-related data has been prompted by the applications of (1) computer and data sciences, including data mining, predictive analytics, machine learning, and artificial intelligence, and (2) advances in imaging technology and probe development. Computational modelling and simulation are systematic and cost-effective tools able to identify important temporal/spatial patterns (and relationships), characterize distinct molecular features of cancer states, and address other relevant aspects, including tumor detection and heterogeneity, progression and metastasis, and drug resistance. These approaches have provided invaluable insights for improving the experimental design of therapeutic delivery systems and for increasing the translational value of the results obtained from early and preclinical studies. The big question is: Could cancer theranostics be determined and controlled in silico? This review describes the recent progress in the development of computational models and methods used to facilitate research on the molecular basis of cancer and on the respective diagnosis and optimized treatment, with particular emphasis on the design and optimization of theranostic systems. The current role of computational approaches is providing innovative, incremental, and complementary data-driven solutions for the prediction, simplification, and characterization of cancer and intrinsic mechanisms, and to promote new data-intensive, accurate diagnostics and therapeutics.

Modeling the Effect of High Calorie Diet on the Interplay between Adipose Tissue, Inflammation, and Diabetes

BACKGROUND

Type 2 diabetes (T2D) is a chronic metabolic disease potentially leading to serious widespread tissue damage. Human organism develops T2D when the glucose-insulin control is broken for reasons that are not fully understood but have been demonstrated to be linked to the emergence of a chronic inflammation. Indeed such low-level chronic inflammation affects the pancreatic production of insulin and triggers the development of insulin resistance, eventually leading to an impaired control of the blood glucose concentration. On the contrary, it is well-known that obesity and inflammation are strongly correlated.

AIM

In this study, we investigate in silico the effect of overfeeding on the adipose tissue and the consequent set up of an inflammatory state. We model the emergence of the inflammation as the result of adipose mass increase which, in turn, is a direct consequence of a prolonged excess of high calorie intake.

RESULTS

The model reproduces the fat accumulation due to excessive caloric intake observed in two clinical studies. Moreover, while showing consistent weight gains over long periods of time, it reveals a drift of the macrophage population toward the proinflammatory phenotype, thus confirming its association with fatness.

Review of the systems biology of the immune system using agent-based models

The immune system is an inherent protection system in vertebrate animals including human beings that exhibit properties such as self-organisation, self-adaptation, learning, and recognition. It interacts with the other allied systems such as the gut and lymph nodes. There is a need for immune system modelling to know about its complex internal mechanism, to understand how it maintains the homoeostasis, and how it interacts with the other systems. There are two types of modelling techniques used for the simulation of features of the immune system: equation-based modelling (EBM) and agent-based modelling. Owing to certain shortcomings of the EBM, agent-based modelling techniques are being widely used. This technique provides various predictions for disease causes and treatments; it also helps in hypothesis verification. This study presents a review of agent-based modelling of the immune system and its interactions with the gut and lymph nodes. The authors also review the modelling of immune system interactions during tuberculosis and cancer. In addition, they also outline the future research directions for the immune system simulation through agent-based techniques such as the effects of stress on the immune system, evolution of the immune system, and identification of the parameters for a healthy immune system.

Computational modeling of immune system of the fish for a more effective vaccination in aquaculture

Motivation

A computational model equipped with the main immunological features of the sea bass (Dicentrarchus labrax L.) immune system was used to predict more effective vaccination in fish. The performance of the model was evaluated by using the results of two in vivo vaccinations trials against L. anguillarum and P. damselae.

Results

Tests were performed to select the appropriate doses of vaccine and infectious bacteria to set up the model. Simulation outputs were compared with the specific antibody production and the expression of BcR and TcR gene transcripts in spleen. The model has shown a good ability to be used in sea bass and could be implemented for different routes of vaccine administration even with more than two pathogens. The model confirms the suitability of in silico methods to optimize vaccine doses and the immune response to them. This model could be applied to other species to optimize the design of new vaccination treatments of fish in aquaculture.

Availability and implementation

The method is available at http://www.iac.cnr.it/∼filippo/c-immsim/.

Contact

nromano@unitus.it.

Supplementary information

Supplementary data are available at Bioinformatics online.

Addressing current challenges in cancer immunotherapy with mathematical and computational modelling

The goal of cancer immunotherapy is to boost a patient's immune response to a tumour. Yet, the design of an effective immunotherapy is complicated by various factors, including a potentially immunosuppressive tumour microenvironment, immune-modulating effects of conventional treatments and therapy-related toxicities. These complexities can be incorporated into mathematical and computational models of cancer immunotherapy that can then be used to aid in rational therapy design. In this review, we survey modelling approaches under the umbrella of the major challenges facing immunotherapy development, which encompass tumour classification, optimal treatment scheduling and combination therapy design. Although overlapping, each challenge has presented unique opportunities for modellers to make contributions using analytical and numerical analysis of model outcomes, as well as optimization algorithms. We discuss several examples of models that have grown in complexity as more biological information has become available, showcasing how model development is a dynamic process interlinked with the rapid advances in tumour-immune biology. We conclude the review with recommendations for modellers both with respect to methodology and biological direction that might help keep modellers at the forefront of cancer immunotherapy development.

The pulse vaccination effects in mammary carcinoma

Induction of antitumor immunity by vaccination is one of the major current immunotherapy strategies. We present a mathematical model of the competition between immune cells and mammary carcinogenesis under the effect of Triplex vaccine. The model describes both humoral and cell-mediated immune responses against cancer cells. The control of the cancer cells growth occurs through the application of the pulse vaccination. Here we determine the relationship between the strength of the vaccine and the time required to eradicate cancer cells, and we present some simulations to illustrate our theoretical results, namely, the total cancer cells depletion, which is influenced by competition occurs among the immune and cancer cells.

Mathematical Modelling and Analysis of the Tumor Treatment Regimens with Pulsed Immunotherapy and Chemotherapy

To begin with, in this paper, single immunotherapy, single chemotherapy, and mixed treatment are discussed, and sufficient conditions under which tumor cells will be eliminated ultimately are obtained. We analyze the impacts of the least effective concentration and the half-life of the drug on therapeutic results and then find that increasing the least effective concentration or extending the half-life of the drug can achieve better therapeutic effects. In addition, since most types of tumors are resistant to common chemotherapy drugs, we consider the impact of drug resistance on therapeutic results and propose a new mathematical model to explain the cause of the chemotherapeutic failure using single drug. Based on this, in the end, we explore the therapeutic effects of two-drug combination chemotherapy, as well as mixed immunotherapy with combination chemotherapy. Numerical simulations indicate that combination chemotherapy is very effective in controlling tumor growth. In comparison, mixed immunotherapy with combination chemotherapy can achieve a better treatment effect.

Mathematical modeling for novel cancer drug discovery and development

INTRODUCTION

Mathematical modeling enables: the in silico classification of cancers, the prediction of disease outcomes, optimization of therapy, identification of promising drug targets and prediction of resistance to anticancer drugs. In silico pre-screened drug targets can be validated by a small number of carefully selected experiments.

AREAS COVERED

This review discusses the basics of mathematical modeling in cancer drug discovery and development. The topics include in silico discovery of novel molecular drug targets, optimization of immunotherapies, personalized medicine and guiding preclinical and clinical trials. Breast cancer has been used to demonstrate the applications of mathematical modeling in cancer diagnostics, the identification of high-risk population, cancer screening strategies, prediction of tumor growth and guiding cancer treatment.

EXPERT OPINION

Mathematical models are the key components of the toolkit used in the fight against cancer. The combinatorial complexity of new drugs discovery is enormous, making systematic drug discovery, by experimentation, alone difficult if not impossible. The biggest challenges include seamless integration of growing data, information and knowledge, and making them available for a multiplicity of analyses. Mathematical models are essential for bringing cancer drug discovery into the era of Omics, Big Data and personalized medicine.

Current trends in cancer vaccines--a bioinformatics perspective

Cancer vaccine development is in the process of becoming reality in future, due to successful phase II/III clinical trials. However, there are still problems due to the specificity of tumor antigens and weakness of tumor associated antigens in eliciting an effective immune response. Computational models to assess the vaccine efficacy have helped to improve and understand what is necessary for personalized treatment. Further research is needed to elucidate the mechanisms of activation of antigen specific cytotoxic T lymphocytes, decreased TREG number functionality and antigen cascade, so that overall improvement in vaccine efficacy and disease free survival can be attained. T cell epitomic based in sillico approaches might be very effective for the design and development of novel cancer vaccines.

A study of tumour growth based on stoichiometric principles: a continuous model and its discrete analogue

In this paper, we consider a continuous mathematically tractable model and its discrete analogue for the tumour growth. The model formulation is based on stoichiometric principles considering tumour-immune cell interactions in potassium (K (+))-limited environment. Our both continuous and discrete models illustrate 'cancer immunoediting' as a dynamic process having all three phases namely elimination, equilibrium and escape. The stoichiometric principles introduced into the model allow us to study its dynamics with the variation in the total potassium in the surrounding of the tumour region. It is found that an increase in the total potassium may help the patient fight the disease for a longer period of time. This result seems to be in line with the protective role of the potassium against the risk of pancreatic cancer as has been reported by Bravi et al. [Dietary intake of selected micronutrients and risk of pancreatic cancer: An Italian case-control study, Ann. Oncol. 22 (2011), pp. 202-206].

Agent based modeling of Treg-Teff cross regulation in relapsing-remitting multiple sclerosis

Background

Multiple sclerosis (MS) is a disease of central nervous system that causes the removal of fatty myelin sheath from axons of the brain and spinal cord. Autoimmunity plays an important role in this pathology outcome and body's own immune system attacks on the myelin sheath causing the damage. The etiology of the disease is partially understood and the response to treatment cannot easily be predicted.

Results

We presented the results obtained using 8 genetically predisposed randomly chosen individuals reproducing both the absence and presence of malfunctions of the Teff-Treg cross-balancing mechanisms at a local level. For simulating the absence of a local malfunction we supposed that both Teff and Treg populations had similar maximum duplication rates. Results presented here suggest that presence of a genetic predisposition is not always a sufficient condition for developing the disease. Other conditions such as a breakdown of the mechanisms that regulate and allow peripheral tolerance should be involved.

Conclusions

The presented model allows to capture the essential dynamics of relapsing-remitting MS despite its simplicity. It gave useful insights that support the hypothesis of a breakdown of Teff-Treg cross balancing mechanisms.

VaccImm: simulating peptide vaccination in cancer therapy

Background

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.

Results

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.

Conclusions

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.

Cancer vaccines: state of the art of the computational modeling approaches

Cancer vaccines are a real application of the extensive knowledge of immunology to the field of oncology. Tumors are dynamic complex systems in which several entities, events, and conditions interact among them resulting in growth, invasion, and metastases. The immune system includes many cells and molecules that cooperatively act to protect the host organism from foreign agents. Interactions between the immune system and the tumor mass include a huge number of biological factors. Testing of some cancer vaccine features, such as the best conditions for vaccine administration or the identification of candidate antigenic stimuli, can be very difficult or even impossible only through experiments with biological models simply because a high number of variables need to be considered at the same time. This is where computational models, and, to this extent, immunoinformatics, can prove handy as they have shown to be able to reproduce enough biological complexity to be of use in suggesting new experiments. Indeed, computational models can be used in addition to biological models. We now experience that biologists and medical doctors are progressively convinced that modeling can be of great help in understanding experimental results and planning new experiments. This will boost this research in the future.

Mathematical modeling of the immune system recognition to mammary carcinoma antigen

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.

Optimal vaccination schedule search using genetic algorithm over MPI technology

BACKGROUND

Immunological strategies that achieve the prevention of tumor growth are based on the presumption that the immune system, if triggered before tumor onset, could be able to defend from specific cancers. In supporting this assertion, in the last decade active immunization approaches prevented some virus-related cancers in humans. An immunopreventive cell vaccine for the non-virus-related human breast cancer has been recently developed. This vaccine, called Triplex, targets the HER-2-neu oncogene in HER-2/neu transgenic mice and has shown to almost completely prevent HER-2/neu-driven mammary carcinogenesis when administered with an intensive and life-long schedule.

METHODS

To better understand the preventive efficacy of the Triplex vaccine in reduced schedules we employed a computational approach. The computer model developed allowed us to test in silico specific vaccination schedules in the quest for optimality. Specifically here we present a parallel genetic algorithm able to suggest optimal vaccination schedule.

RESULTS & CONCLUSIONS

The enormous complexity of combinatorial space to be explored makes this approach the only possible one. The suggested schedule was then tested in vivo, giving good results. Finally, biologically relevant outcomes of optimization are presented.

Relaxation estimation of RMSD in molecular dynamics immunosimulations

Molecular dynamics simulations have to be sufficiently long to draw reliable conclusions. However, no method exists to prove that a simulation has converged. We suggest the method of "lagged RMSD-analysis" as a tool to judge if an MD simulation has not yet run long enough. The analysis is based on RMSD values between pairs of configurations separated by variable time intervals Δt. Unless RMSD(Δt) has reached a stationary shape, the simulation has not yet converged.

Fine-tuning anti-tumor immunotherapies via stochastic simulations

BACKGROUND

Anti-tumor therapies aim at reducing to zero the number of tumor cells in a host within their end or, at least, aim at leaving the patient with a sufficiently small number of tumor cells so that the residual tumor can be eradicated by the immune system. Besides severe side-effects, a key problem of such therapies is finding a suitable scheduling of their administration to the patients. In this paper we study the effect of varying therapy-related parameters on the final outcome of the interplay between a tumor and the immune system.

RESULTS

This work generalizes our previous study on hybrid models of such an interplay where interleukins are modeled as a continuous variable, and the tumor and the immune system as a discrete-state continuous-time stochastic process. The hybrid model we use is obtained by modifying the corresponding deterministic model, originally proposed by Kirschner and Panetta. We consider Adoptive Cellular Immunotherapies and Interleukin-based therapies, as well as their combination. By asymptotic and transitory analyses of the corresponding deterministic model we find conditions guaranteeing tumor eradication, and we tune the parameters of the hybrid model accordingly. We then perform stochastic simulations of the hybrid model under various therapeutic settings: constant, piece-wise constant or impulsive infusion and daily or weekly delivery schedules.

CONCLUSIONS

Results suggest that, in some cases, the delivery schedule may deeply impact on the therapy-induced tumor eradication time. Indeed, our model suggests that Interleukin-based therapies may not be effective for every patient, and that the piece-wise constant is the most effective delivery to stimulate the immune-response. For Adoptive Cellular Immunotherapies a metronomic delivery seems more effective, as it happens for other anti-angiogenesis therapies and chemotherapies, and the impulsive delivery seems more effective than the piece-wise constant. The expected synergistic effects have been observed when the therapies are combined.

Combining cellular automata and Lattice Boltzmann method to model multiscale avascular tumor growth coupled with nutrient diffusion and immune competition

In the last decades the Lattice Boltzmann method (LB) has been successfully used to simulate a variety of processes. The LB model describes the microscopic processes occurring at the cellular level and the macroscopic processes occurring at the continuum level with a unique function, the probability distribution function. Recently, it has been tried to couple deterministic approaches with probabilistic cellular automata (probabilistic CA) methods with the aim to model temporal evolution of tumor growths and three dimensional spatial evolution, obtaining hybrid methodologies. Despite the good results attained by CA-PDE methods, there is one important issue which has not been completely solved: the intrinsic stochastic nature of the interactions at the interface between cellular (microscopic) and continuum (macroscopic) level. CA methods are able to cope with the stochastic phenomena because of their probabilistic nature, while PDE methods are fully deterministic. Even if the coupling is mathematically correct, there could be important statistical effects that could be missed by the PDE approach. For such a reason, to be able to develop and manage a model that takes into account all these three level of complexity (cellular, molecular and continuum), we believe that PDE should be replaced with a statistic and stochastic model based on the numerical discretization of the Boltzmann equation: The Lattice Boltzmann (LB) method. In this work we introduce a new hybrid method to simulate tumor growth and immune system, by applying Cellular Automata Lattice Boltzmann (CA-LB) approach.

Emerging vaccine informatics

Vaccine informatics is an emerging research area that focuses on development and applications of bioinformatics methods that can be used to facilitate every aspect of the preclinical, clinical, and postlicensure vaccine enterprises. Many immunoinformatics algorithms and resources have been developed to predict T- and B-cell immune epitopes for epitope vaccine development and protective immunity analysis. Vaccine protein candidates are predictable in silico from genome sequences using reverse vaccinology. Systematic transcriptomics and proteomics gene expression analyses facilitate rational vaccine design and identification of gene responses that are correlates of protection in vivo. Mathematical simulations have been used to model host-pathogen interactions and improve vaccine production and vaccination protocols. Computational methods have also been used for development of immunization registries or immunization information systems, assessment of vaccine safety and efficacy, and immunization modeling. Computational literature mining and databases effectively process, mine, and store large amounts of vaccine literature and data. Vaccine Ontology (VO) has been initiated to integrate various vaccine data and support automated reasoning.

Criticality of timing for anti-HIV therapy initiation

The time of initiation of antiretroviral therapy in HIV-1 infected patients has a determinant effect on the viral dynamics. The question is, how far can the therapy be delayed? Is sooner always better? We resort to clinical data and to microsimulations to forecast the dynamics of the viral load at therapy interruption after prolonged antiretroviral treatment. A computational model previously evaluated, produces results that are statistically adherent to clinical data. In addition, it allows a finer grain analysis of the impact of the therapy initiation point to the disease course. We find a swift increase of the viral density as a function of the time of initiation of the therapy measured when the therapy is stopped. In particular there is a critical time delay with respect to the infection instant beyond which the therapy does not affect the viral rebound. Initiation of the treatment is beneficial because it can down-regulate the immune activation, hence limiting viral replication and spread.

Theoretical modeling techniques and their impact on tumor immunology

Currently, cancer is one of the leading causes of death in industrial nations. While conventional cancer treatment usually results in the patient suffering from severe side effects, immunotherapy is a promising alternative. Nevertheless, some questions remain unanswered with regard to using immunotherapy to treat cancer hindering it from being widely established. To help rectify this deficit in knowledge, experimental data, accumulated from a huge number of different studies, can be integrated into theoretical models of the tumor-immune system interaction. Many complex mechanisms in immunology and oncology cannot be measured in experiments, but can be analyzed by mathematical simulations. Using theoretical modeling techniques, general principles of tumor-immune system interactions can be explored and clinical treatment schedules optimized to lower both tumor burden and side effects. In this paper, we aim to explain the main mathematical and computational modeling techniques used in tumor immunology to experimental researchers and clinicians. In addition, we review relevant published work and provide an overview of its impact to the field.

Modeling the competition between lung metastases and the immune system using agents

BACKGROUND

The Triplex cell vaccine is a cancer cellular vaccine that can prevent almost completely the mammary tumor onset in HER-2/neu transgenic mice. In a translational perspective, the activity of the Triplex vaccine was also investigated against lung metastases showing that the vaccine is an effective treatment also for the cure of metastases. A future human application of the Triplex vaccine should take into account several aspects of biological behavior of the involved entities to improve the efficacy of therapeutic treatment and to try to predict, for example, the outcomes of longer experiments in order to move faster towards clinical phase I trials. To help to address this problem, MetastaSim, a hybrid Agent Based - ODE model for the simulation of the vaccine-elicited immune system response against lung metastases in mice is presented. The model is used as in silico wet-lab. As a first application MetastaSim is used to find protocols capable of maximizing the total number of prevented metastases, minimizing the number of vaccine administrations.

RESULTS

The model shows that it is possible to obtain "in silico" a 45% reduction in the number of vaccinations. The analysis of the results further suggests that any optimal protocol for preventing lung metastases formation should be composed by an initial massive vaccine dosage followed by few vaccine recalls.

CONCLUSIONS

Such a reduction may represent an important result from the point of view of translational medicine to humans, since a downsizing of the number of vaccinations is usually advisable in order to minimize undesirable effects. The suggested vaccination strategy also represents a notable outcome. Even if this strategy is commonly used for many infectious diseases such as tetanus and hepatitis-B, it can be in fact considered as a relevant result in the field of cancer-vaccines immunotherapy. These results can be then used and verified in future "in vivo" experiments, and their outcome can be used to further improve and refine the model.

In silico modeling and in vivo efficacy of cancer-preventive vaccinations

Cancer vaccine feasibility would benefit from reducing the number and duration of vaccinations without diminishing efficacy. However, the duration of in vivo studies and the huge number of possible variations in vaccination protocols have discouraged their optimization. In this study, we employed an established mouse model of preventive vaccination using HER-2/neu transgenic mice (BALB-neuT) to validate in silico-designed protocols that reduce the number of vaccinations and optimize efficacy. With biological training, the in silico model captured the overall in vivo behavior and highlighted certain critical issues. First, although vaccinations could be reduced in number without sacrificing efficacy, the intensity of early vaccinations was a key determinant of long-term tumor prevention needed for predictive utility in the model. Second, after vaccinations ended, older mice exhibited more rapid tumor onset and sharper decline in antibody levels than young mice, emphasizing immune aging as a key variable in models of vaccine protocols for elderly individuals. Long-term studies confirmed predictions of in silico modeling in which an immune plateau phase, once reached, could be maintained with a reduced number of vaccinations. Furthermore, that rapid priming in young mice is required for long-term antitumor protection, and that the accuracy of mathematical modeling of early immune responses is critical. Finally, that the design and modeling of cancer vaccines and vaccination protocols must take into account the progressive aging of the immune system, by striving to boost immune responses in elderly hosts. Our results show that an integrated in vivo-in silico approach could improve both mathematical and biological models of cancer immunoprevention.

ImmunoGrid: towards agent-based simulations of the human immune system at a natural scale

The ultimate aim of the EU-funded ImmunoGrid project is to develop a natural-scale model of the human immune system-that is, one that reflects both the diversity and the relative proportions of the molecules and cells that comprise it-together with the grid infrastructure necessary to apply this model to specific applications in the field of immunology. These objectives present the ImmunoGrid Consortium with formidable challenges in terms of complexity of the immune system, our partial understanding about how the immune system works, the lack of reliable data and the scale of computational resources required. In this paper, we explain the key challenges and the approaches adopted to overcome them. We also consider wider implications for the present ambitious plans to develop natural-scale, integrated models of the human body that can make contributions to personalized health care, such as the European Virtual Physiological Human initiative. Finally, we ask a key question: How long will it take us to resolve these challenges and when can we expect to have fully functional models that will deliver health-care benefits in the form of personalized care solutions and improved disease prevention?

Interactions between the immune system and cancer: a brief review of non-spatial mathematical models

We briefly review spatially homogeneous mechanistic mathematical models describing the interactions between a malignant tumor and the immune system. We begin with the simplest (single equation) models for tumor growth and proceed to consider greater immunological detail (and correspondingly more equations) in steps. This approach allows us to clarify the necessity for expanding the complexity of models in order to capture the biological mechanisms we wish to understand. We conclude by discussing some unsolved problems in the mathematical modeling of cancer-immune system interactions.

A novel paradigm for cell and molecule interaction ontology: from the CMM model to IMGT-ONTOLOGY

Background

Biology is moving fast toward the virtuous circle of other disciplines: from data to quantitative modeling and back to data. Models are usually developed by mathematicians, physicists, and computer scientists to translate qualitative or semi-quantitative biological knowledge into a quantitative approach. To eliminate semantic confusion between biology and other disciplines, it is necessary to have a list of the most important and frequently used concepts coherently defined.

Results

We propose a novel paradigm for generating new concepts for an ontology, starting from model rather than developing a database. We apply that approach to generate concepts for cell and molecule interaction starting from an agent based model. This effort provides a solid infrastructure that is useful to overcome the semantic ambiguities that arise between biologists and mathematicians, physicists, and computer scientists, when they interact in a multidisciplinary field.

Conclusions

This effort represents the first attempt at linking molecule ontology with cell ontology, in IMGT-ONTOLOGY, the well established ontology in immunogenetics and immunoinformatics, and a paradigm for life science biology. With the increasing use of models in biology and medicine, the need to link different levels, from molecules to cells to tissues and organs, is increasingly important.

Immune control of HIV-1 infection after therapy interruption: immediate versus deferred antiretroviral therapy

BACKGROUND

The optimal stage for initiating antiretroviral therapies in HIV-1 bearing patients is still a matter of debate.

METHODS

We present computer simulations of HIV-1 infection aimed at identifying the pro et contra of immediate as compared to deferred Highly Active Antiretroviral Therapy (HAART).

RESULTS

Our simulations highlight that a prompt specific CD8+ cytotoxic T lymphocytes response is detected when therapy is delayed. Compared to very early initiation of HAART, in deferred treated patients CD8+ T cells manage to mediate the decline of viremia in a shorter time and, at interruption of therapy, the virus experiences a stronger immune pressure. We also observe, however, that the immunological effects of the therapy fade with time in both therapeutic regimens. Thus, within one year from discontinuation, viral burden recovers to the value at which it would level off in the absence of therapy.In summary, simulations show that immediate therapy does not prolong the disease-free period and does not confer a survival benefit when compared to treatment started during the chronic infection phase.

CONCLUSION

Our conclusion is that, since there is no therapy to date that guarantees life-long protection, deferral of therapy should be preferred in order to minimize the risk of adverse effects, the occurrence of drug resistances and the costs of treatment.

ImmunoGrid, an integrative environment for large-scale simulation of the immune system for vaccine discovery, design and optimization

Vaccine research is a combinatorial science requiring computational analysis of vaccine components, formulations and optimization. We have developed a framework that combines computational tools for the study of immune function and vaccine development. This framework, named ImmunoGrid combines conceptual models of the immune system, models of antigen processing and presentation, system-level models of the immune system, Grid computing, and database technology to facilitate discovery, formulation and optimization of vaccines. ImmunoGrid modules share common conceptual models and ontologies. The ImmunoGrid portal offers access to educational simulators where previously defined cases can be displayed, and to research simulators that allow the development of new, or tuning of existing, computational models. The portal is accessible at <igrid-ext.cryst.bbk.ac.uk/immunogrid>.

Agent-based modeling of host-pathogen systems: The successes and challenges

Agent-based models have been employed to describe numerous processes in immunology. Simulations based on these types of models have been used to enhance our understanding of immunology and disease pathology. We review various agent-based models relevant to host–pathogen systems and discuss their contributions to our understanding of biological processes. We then point out some limitations and challenges of agent-based models and encourage efforts towards reproducibility and model validation.