Mathematical Analysis of Virotherapy Treatment for Cancer

  • Renuka J Assistant Professor, Department of Mathematics, Women’s Christian College, Chennai, Tamil Nadu, India
  • Balaganesan P Department of Mathematics, AMET Deemed to be University, Chennai, Tamil Nadu, India.
  • Vijayakumar PN Research Scholar, AMET Deemed to be University, Chennai, Tamil Nadu, India.
Keywords: Non-linear Differential Equations, Cancer Cell Count, Count of Infected Cancer Cells, Count of Dead Cells, Death Rate of Infected Cancer Cells, Homotopy Perturbation Method


A virotherapy model for cancer has been identified and analysed in this
article. This mathematical model of cancer virotherapy was developed
by Friedman. It was constructed by means of non-linear differential
equations related to “the count of cancer cells, the count of infected
cancer cells, the count of dead cells, and the count of virus cells that
do not belong to tumour cells”. It also includes parameters for the rate
of viral infection of cancer cells, the proliferation rate of cancer cells,
the removal rate of debris of dead cells, the death rate of infected
cancer cells, and the count of virus particles that are not contained
in cancer cells. The Homotopy Perturbation Method was applied to
solve the non-linear differential equations and was analysed both
qualitatively and numerically. The values of the variation in parameters
were evaluated to investigate the impact of controls on the spread of
cancer cells. Through the numerical and graphical results, the optimal
control, which helps to significantly reduce the impact of cancer cells,
has been discussed. A significant agreement is noted with approximate
analytical results and numerical simulations.

How to cite this article:
Renuka J, Balaganesan P, Vijayakumar PN.
Mathematical Analysis of Virotherapy Treatment
for Cancer. Chettinad Health City Med J.



Cancer Research UK [Internet]. Worldwide cancer statistics; [cited 2022 Dec 10]. Available from:

National Cancer Institute [Internet]. Cancer statistics, technical report; [cited 2022 Dec 10]. Available from:

Mokhtari RB, Homayouni TS, Baluch N, Morgatskaya E, Kumar S, Das B, Yeger H. Combination therapy in combating cancer. Oncotarget. 2017;8(23):38022-43. [PubMed] [Google Scholar]

Khan KA, Kerbel RS. Improving immunotherapy outcomes with anti-angiogenic treatments and vice versa. Nat Rev Clin Oncol. 2018;15(5):310. [PubMed] [Google Scholar]

Teles FF, Lemos JM. Cancer therapy optimization based on multiple model adaptive control. Biomed Signal Process Control. 2019;48:255-64. [Google Scholar]

Babaei N, Salamci MU. Mixed therapy in cancer treatment for personalized drug administration using model reference adaptive control. Eur J Control. 2019;50:117-37. [Google Scholar]

Child JA, Morgan GJ, Davies FE, Owen RG, Bell SE, Hawkins K, Brown J, Drayson MT, Selby PJ; Medical Research Council Adult Leukaemia Working Party. High-dose chemotherapy with hematopoietic stem-cell rescue for multiple myeloma. N Engl J Med. 2003;348(19):1875-83. [PubMed] [Google Scholar]

Eftimie R, Bramson JL, Earn DJ. Interactions between the immune system and cancer: a brief review ofnon-spatial mathematical models. Bull Math Biol. 2011;73(1):2-32.[PubMed] [Google Scholar]

Chou CS, Friedman A. Introduction to mathematical biology: modeling, analysis, and simulations. Switzerland: Springer International Publishing; 2016. p. 153-5. [Google Scholar]

Pillis LG, Radunskaya A. A mathematical tumor model with immune resistance and drug therapy: an optimal control approach. J Theor Med. 2001;3(2):79-100. [Google Scholar]

Friedman A, Tian JP, Fulci G, Chiocca EA, Wang J. Glioma virotherapy: effects of innate immune suppression and increased viral replication capacity. Cancer Res. 2006;66(4):2314-9. [PubMed] [Google Scholar]