Abstract

In recent years, interest in the capability of virus particles as a treatment for cancer has increased. In this work, we present a mathematical model embodying the interaction between tumour cells and virus particles engineered to infect and destroy cancerous tissue. To quantify the effectiveness of oncolytic virotherapy, we conduct a local stability analysis and bifurcation analysis of our model. In the absence of tumour growth or viral decay, the model predicts that oncolytic virotherapy will successfully eliminate the tumour cell population for a large proportion of initial conditions. In comparison, for growing tumours and decaying viral particles there are no stable equilibria in the model; however, oscillations emerge for certain regions in our parameter space. We investigate how the period and amplitude of oscillations depend on tumour growth and viral decay. We find that higher tumour replication rates result in longer periods between oscillations and lower amplitudes for uninfected tumour cells. From our analysis, we conclude that oncolytic viruses can reduce growing tumours into a stable oscillatory state, but are insufficient to completely eradicate them. We propose that it is only with the addition of other anti-cancer agents that tumour eradication may be achieved by oncolytic virus.