The role of viral infectivity in oncolytic virotherapy outcomes: A mathematical study

Finding Hopf bifurcation islands and identifying thresholds for success or failure in oncolytic viral therapy

We model interactions between cancer cells and viruses during oncolytic viral therapy. One of our primary goals is to identify parameter regions that yield treatment failure or success. We show that the tumor size under therapy at a particular time is less than the size without therapy. Our analysis demonstrates two thresholds for the horizontal transmission rate: a "failure threshold" below which treatment fails, and a "success threshold" above which infection prevalence reaches 100% and the tumor shrinks to its smallest size. Moreover, we explain how changes in the virulence of the virus alter the success threshold and the minimum tumor size. Our study suggests that the optimal virulence of an oncolytic virus depends on the timescale of virus dynamics. We identify a threshold for the virulence of the virus and show how this threshold depends on the timescale of virus dynamics. Our results suggest that when the timescale of virus dynamics is fast, administering a more virulent virus leads to a greater reduction in the tumor size. Conversely, when the viral timescale is slow, higher virulence can induce oscillations with high amplitude in the tumor size. Furthermore, we introduce the concept of a "Hopf bifurcation Island" in the parameter space, an idea that has applications far beyond the results of this paper and is applicable to many mathematical models. We elucidate what a Hopf bifurcation Island is, and we prove that small Islands can imply very slowly growing oscillatory solutions.

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.

Optimal strategies of oncolytic virus-bortezomib therapy via the apoptotic, necroptotic, and oncolysis signaling network

Bortezomib and oncolytic virotherapy are two emerging targeted cancer therapies. Bortezomib, a proteasome inhibitor, disrupts protein degradation in cells, leading to the accumulation of unfolded proteins that induce apoptosis. On the other hand, virotherapy uses genetically modified oncolytic viruses (OVs) to infect cancer cells, trigger cell lysis, and activate anti-tumor response. Despite progress in cancer treatment, identifying administration protocols for therapeutic agents remains a significant concern, aiming to strike a balance between efficacy, minimizing toxicity, and administrative costs. In this work, optimal control theory was employed to design a cost-effective and efficient co-administration protocols for bortezomib and OVs that could significantly diminish the population of cancer cells via the cell death program with the NF$ \kappa $B-BAX-RIP1 signaling network. Both linear and quadratic control strategies were explored to obtain practical treatment approaches by adapting necroptosis protocols to efficient cell death programs. Our findings demonstrated that a combination therapy commencing with the administration of OVs followed by bortezomib infusions yields an effective tumor-killing outcome. These results could provide valuable guidance for the development of clinical administration protocols in cancer treatment.

Agent-Based and Continuum Models for Spatial Dynamics of Infection by Oncolytic Viruses

The use of oncolytic viruses as cancer treatment has received considerable attention in recent years, however the spatial dynamics of this viral infection is still poorly understood. We present here a stochastic agent-based model describing infected and uninfected cells for solid tumours, which interact with viruses in the absence of an immune response. Two kinds of movement, namely undirected random and pressure-driven movements, are considered: the continuum limit of the models is derived and a systematic comparison between the systems of partial differential equations and the individual-based model, in one and two dimensions, is carried out. In the case of undirected movement, a good agreement between agent-based simulations and the numerical and well-known analytical results for the continuum model is possible. For pressure-driven motion, instead, we observe a wide parameter range in which the infection of the agents remains confined to the center of the tumour, even though the continuum model shows traveling waves of infection; outcomes appear to be more sensitive to stochasticity and uninfected regions appear harder to invade, giving rise to irregular, unpredictable growth patterns. Our results show that the presence of spatial constraints in tumours' microenvironments limiting free expansion has a very significant impact on virotherapy. Outcomes for these tumours suggest a notable increase in variability. All these aspects can have important effects when designing individually tailored therapies where virotherapy is included.

A combination therapy of oncolytic viruses and chimeric antigen receptor T cells: a mathematical model proof-of-concept

Combining chimeric antigen receptor T (CAR-T) cells with oncolytic viruses (OVs) has recently emerged as a promising treatment approach in preclinical studies that aim to alleviate some of the barriers faced by CAR-T cell therapy. In this study, we address by means of mathematical modeling the main question of whether a single dose or multiple sequential doses of CAR-T cells during the OVs therapy can have a synergetic effect on tumor reduction. To that end, we propose an ordinary differential equations-based model with virus-induced synergism to investigate potential effects of different regimes that could result in efficacious combination therapy against tumor cell populations. Model simulations show that, while the treatment with a single dose of CAR-T cells is inadequate to eliminate all tumor cells, combining the same dose with a single dose of OVs can successfully eliminate the tumor in the absence of virus-induced synergism. However, in the presence of virus-induced synergism, the same combination therapy fails to eliminate the tumor. Furthermore, it is shown that if the intensity of virus-induced synergy and/or virus oncolytic potency is high, then the induced CAR-T cell response can inhibit virus oncolysis. Additionally, the simulations show a more robust synergistic effect on tumor cell reduction when OVs and CAR-T cells are administered simultaneously compared to the combination treatment where CAR-T cells are administered first or after OV injection. Our findings suggest that the combination therapy of CAR-T cells and OVs seems unlikely to be effective if the virus-induced synergistic effects are included when genetically engineering oncolytic viral vectors.