Modeling Pancreatic Cancer Dynamics with Immunotherapy

Bistability in a model of tumor-immune system interactions with an oncolytic viral therapy

A mathematical model of tumor-immune system interactions with an oncolytic virus therapy for which the immune system plays a twofold role against cancer cells is derived. The immune cells can kill cancer cells but can also eliminate viruses from the therapy. In addition, immune cells can either be stimulated to proliferate or be impaired to reduce their growth by tumor cells. It is shown that if the tumor killing rate by immune cells is above a critical value, the tumor can be eradicated for all sizes, where the critical killing rate depends on whether the immune system is immunosuppressive or proliferative. For a reduced tumor killing rate with an immunosuppressive immune system, that bistability exists in a large parameter space follows from our numerical bifurcation study. Depending on the tumor size, the tumor can either be eradicated or be reduced to a size less than its carrying capacity. However, reducing the viral killing rate by immune cells always increases the effectiveness of the viral therapy. This reduction may be achieved by manipulating certain genes of viruses via genetic engineering or by chemical modification of viral coat proteins to avoid detection by the immune cells.

MATHEMATICAL MODELING OF PANCREATIC CANCER TREATMENT WITH CANCER STEM CELLS

Of all cancers, pancreatic cancer has a significantly low rate of survival, mostly due to lack of early screening. Thus, once detected, pancreatic cancer is usually in later stages, reducing the likelihood of full recovery. The most common treatment strategy is chemotherapy, although several immunotherapeutic drugs show promising results in extending the patient’s lifespan. In this paper, we provide a validated mathematical model for the pancreatic cancer after fitting the parameter values, such as tumor growth rate, inverse carrying capacity, activation and decay rate of pancreatic stellate cells, with the use of the experimental data presented by Cioffi et al. cioffi2015inhibition For treatments with the chemotherapeutic drugs, Abraxane and Gemcitabine, and the immunotherapeutic drug, Anti-CD47, we modified the model accurately and compared the simulation results with the experimental data not only to model pancreatic cancer treatment correctly but also to move forward with other drug trials. Then, we include the cancer stem cells, which are known to initiate tumors and cause a relapse post-chemotherapy, per cancer stem cell hypothesis so that cancer progression can be assessed based on this phenomenon. In addition, we investigate optimal drug protocols. We find out that the most effective treatment is dual therapy due to extending survival time when compared to other drugs. Moreover, this study reveals that drug dose is more effectual than frequency of drug injection on account of different treatment scheduling with the same dose over a week. The model could be a starting point to investigate pancreatic cancer progression based on cancer stem cell hypothesis and shed light on novel drug discoveries.

Improving cancer treatments via dynamical biophysical models

Despite significant advances in oncological research, cancer nowadays remains one of the main causes of mortality and morbidity worldwide. New treatment techniques, as a rule, have limited efficacy, target only a narrow range of oncological diseases, and have limited availability to the general public due their high cost. An important goal in oncology is thus the modification of the types of antitumor therapy and their combinations, that are already introduced into clinical practice, with the goal of increasing the overall treatment efficacy. One option to achieve this goal is optimization of the schedules of drugs administration or performing other medical actions. Several factors complicate such tasks: the adverse effects of treatments on healthy cell populations, which must be kept tolerable; the emergence of drug resistance due to the intrinsic plasticity of heterogeneous cancer cell populations; the interplay between different types of therapies administered simultaneously. Mathematical modeling, in which a tumor and its microenvironment are considered as a single complex system, can address this complexity and can indicate potentially effective protocols, that would require experimental verification. In this review, we consider classical methods, current trends and future prospects in the field of mathematical modeling of tumor growth and treatment. In particular, methods of treatment optimization are discussed with several examples of specific problems related to different types of treatment.