Analysis of virotherapy in solid tumor invasion

Oscillations in a Spatial Oncolytic Virus Model

Virotherapy treatment is a new and promising target therapy that selectively attacks cancer cells without harming normal cells. Mathematical models of oncolytic viruses have shown predator-prey like oscillatory patterns as result of an underlying Hopf bifurcation. In a spatial context, these oscillations can lead to different spatio-temporal phenomena such as hollow-ring patterns, target patterns, and dispersed patterns. In this paper we continue the systematic analysis of these spatial oscillations and discuss their relevance in the clinical context. We consider a bifurcation analysis of a spatially explicit reaction-diffusion model to find the above mentioned spatio-temporal virus infection patterns. The desired pattern for tumor eradication is the hollow ring pattern and we find exact conditions for its occurrence. Moreover, we derive the minimal speed of travelling invasion waves for the cancer and for the oncolytic virus. Our numerical simulations in 2-D reveal complex spatial interactions of the virus infection and a new phenomenon of a periodic peak splitting. An effect that we cannot explain with our current methods.

Mathematical Modeling of Oncolytic Virus Therapy Reveals Role of the Immune Response

Oncolytic adenoviruses (OAds) present a promising path for cancer treatment due to their selectivity in infecting and lysing tumor cells and their ability to stimulate the immune response. In this study, we use an ordinary differential equation (ODE) model of tumor growth inhibited by oncolytic virus activity to parameterize previous research on the effect of genetically re-engineered OAds in A549 lung cancer tumors in murine models. We find that the data are best fit by a model that accounts for an immune response, and that the immune response provides a mechanism for elimination of the tumor. We also find that parameter estimates for the most effective OAds share characteristics, most notably a high infection rate and low viral clearance rate, that might be potential reasons for these viruses' efficacy in delaying tumor growth. Further studies observing E1A and P19 recombined viruses in different tumor environments may further illuminate the extent of the effects of these genetic modifications.

Dynamic analysis of an age structure model for oncolytic virus therapy

Cancer is recognized as one of the serious diseases threatening human health. Oncolytic therapy is a safe and effective new cancer treatment method. Considering the limited ability of uninfected tumor cells to infect and the age of infected tumor cells have a significant effect on oncolytic therapy, an age-structured model of oncolytic therapy involving Holling-Ⅱ functional response is proposed to investigate the theoretical significance of oncolytic therapy. First, the existence and uniqueness of the solution is obtained. Furthermore, the stability of the system is confirmed. Then, the local stability and global stability of infection-free homeostasis are studied. The uniform persistence and local stability of the infected state are studied. The global stability of the infected state is proved by constructing the Lyapunov function. Finally, the theoretical results are verified by numerical simulation. The results show that when the tumor cells are at the appropriate age, injection of the right amount of oncolytic virus can achieve the purpose of tumor treatment.

Nonlocal multiscale modelling of tumour-oncolytic viruses interactions within a heterogeneous fibrous/non-fibrous extracellular matrix

In this study we investigate computationally tumour-oncolytic virus (OV) interactions that take place within a heterogeneous extracellular matrix (ECM). The ECM is viewed as a mixture of two constitutive phases, namely a fibre phase and a non-fibre phase. The multiscale mathematical model presented here focuses on the nonlocal cell-cell and cell-ECM interactions, and how these interactions might be impacted by the infection of cancer cells with the OV. At macroscale we track the kinetics of cancer cells, virus particles and the ECM. At microscale we track (i) the degradation of ECM by matrix degrading enzymes (MDEs) produced by cancer cells, which further influences the movement of tumour boundary; (ii) the re-arrangement of the microfibres that influences the re-arrangement of macrofibres (i.e., fibres at macroscale). With the help of this new multiscale model, we investigate two questions: (i) whether the infected cancer cell fluxes are the result of local or non-local advection in response to ECM density; and (ii) what is the effect of ECM fibres on the the spatial spread of oncolytic viruses and the outcome of oncolytic virotherapy.

Investigating Macrophages Plasticity Following Tumour-Immune Interactions During Oncolytic Therapies

Over the last few years, oncolytic virus therapy has been recognised as a promising approach in cancer treatment, due to the potential of these viruses to induce systemic anti-tumour immunity and selectively killing tumour cells. However, the effectiveness of these viruses depends significantly on their interactions with the host immune responses, both innate (e.g., macrophages, which accumulate in high numbers inside solid tumours) and adaptive (e.g., \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {CD8}^{+}$$\end{document}CD8+ T cells). In this article, we consider a mathematical approach to investigate the possible outcomes of the complex interactions between two extreme types of macrophages (M1 and M2 cells), effector \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {CD8}^{+}$$\end{document}CD8+ T cells and an oncolytic Vesicular Stomatitis Virus (VSV), on the growth/elimination of B16F10 melanoma. We discuss, in terms of VSV, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {CD8}^{+}$$\end{document}CD8+ and macrophages levels, two different types of immune responses which could ensure tumour control and eventual elimination. We show that both innate and adaptive anti-tumour immune responses, as well as the oncolytic virus, could be very important in delaying tumour relapse and eventually eliminating the tumour. Overall this study supports the use mathematical modelling to increase our understanding of the complex immune interaction following oncolytic virotherapies. However, the complexity of the model combined with a lack of sufficient data for model parametrisation has an impact on the possibility of making quantitative predictions.

Mathematical Analysis of a Mathematical Model of Chemovirotherapy: Effect of Drug Infusion Method

A mathematical model for the treatment of cancer using chemovirotherapy is developed with the aim of determining the efficacy of three drug infusion methods: constant, single bolus, and periodic treatments. The model is in the form of ODEs and is further extended into DDEs to account for delays as a result of the infection of tumor cells by the virus and chemotherapeutic drug responses. Analysis of the model is carried out for each of the three drug infusion methods. Analytic solutions are determined where possible and stability analysis of both steady state solutions for the ODEs and DDEs is presented. The results indicate that constant and periodic drug infusion methods are more efficient compared to a single bolus injection. Numerical simulations show that with a large virus burst size, irrespective of the drug infusion method, chemovirotherapy is highly effective compared to either treatments. The simulations further show that both delays increase the period within which a tumor can be cleared from body tissue.

A mathematical approach to effects of CTLs on cancer virotherapy in the second injection of virus

This paper proposes a planar delay differential equation for cancer virotherapy. The model simulates the situation in which an oncolytic virus is injected for the second time, and the immune system suppresses the viral infection with a time delay. Our purpose is to provide theoretical conditions so that the therapy can be continued successfully. With the help of the characteristic equation, we examine the singularities and their local stability. Hopf bifurcation is also investigated around the endemic singularity. It is shown that there is a sequence of Hopf bifurcations, but the Hopf cycles do not persist continuously between the two sequential bifurcations. Finally, we see that virotherapy can be conducted successfully by controlling the delay parameter.

Modelling the spatiotemporal dynamics of chemovirotherapy cancer treatment

Chemovirotherapy is a combination therapy with chemotherapy and oncolytic viruses. It is gaining more interest and attracting more attention in the clinical setting due to its effective therapy and potential synergistic interactions against cancer. In this paper, we develop and analyse a mathematical model in the form of parabolic non-linear partial differential equations to investigate the spatiotemporal dynamics of tumour cells under chemovirotherapy treatment. The proposed model consists of uninfected and infected tumour cells, a free virus, and a chemotherapeutic drug. The analysis of the model is carried out for both the temporal and spatiotemporal cases. Travelling wave solutions to the spatiotemporal model are used to determine the minimum wave speed of tumour invasion. A sensitivity analysis is performed on the model parameters to establish the key parameters that promote cancer remission during chemovirotherapy treatment. Model analysis of the temporal model suggests that virus burst size and virus infection rate determine the success of the virotherapy treatment, whereas travelling wave solutions to the spatiotemporal model show that tumour diffusivity and growth rate are critical during chemovirotherapy. Simulation results reveal that chemovirotherapy is more effective and a good alternative to either chemotherapy or virotherapy, which is in agreement with the recent experimental studies.

Synthesis of minoxidil conjugates and their evaluation as HL-60 differentiation agents

Activation of minoxidil (MNX) with N,N'-carbonyldiimidazole and coupling with natural polyamines (PAs) and commercially available aliphatic or aromatic amines provided a series of new conjugates which were evaluated for their ability to induce differentiation to HL-60 acute myeloid leukemia cancer cells, using a modified NBTZ reduction test. Although neither MNX nor 4,4'-methylenedianiline (MDA) or 2,7-diaminofluorene (DAF), alone or in combination, had any effect, the MNX-spermine (SPM) conjugate (11) and the conjugates 7 and 8 of MNX with MDA and DAF exhibited a differentiation-inducing effect at a concentration of 10 μM without being toxic on proliferating human peripheral blood mononuclear cells.