Modeling Oncolytic Viral Therapy, Immune Checkpoint Inhibition, and the Complex Dynamics of Innate and Adaptive Immunity in Glioblastoma Treatment

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.

The role of immune cells in resistance to oncolytic viral therapy

Resistance to treatment poses a major challenge for cancer therapy, and oncoviral treatment encounters the issue of viral resistance as well. In this investigation, we introduce deterministic differential equation models to explore the effect of resistance on oncolytic viral therapy. Specifically, we classify tumor cells into resistant, sensitive, or infected with respect to oncolytic viruses for our analysis. Immune cells can eliminate both tumor cells and viruses. Our research shows that the introduction of immune cells into the tumor-virus interaction prevents all tumor cells from becoming resistant in the absence of conversion from resistance to sensitivity, given that the proliferation rate of immune cells exceeds their death rate. The inclusion of immune cells leads to an additional virus-free equilibrium when the immune cell recruitment rate is sufficiently high. The total tumor burden at this virus-free equilibrium is smaller than that at the virus-free and immune-free equilibrium. Therefore, immune cells are capable of reducing the tumor load under the condition of sufficient immune strength. Numerical investigations reveal that the virus transmission rate and parameters related to the immune response significantly impact treatment outcomes. However, monotherapy alone is insufficient for eradicating tumor cells, necessitating the implementation of additional therapies. Further numerical simulation shows that combination therapy with chimeric antigen receptor (CAR T-cell) therapy can enhance the success of treatment.

Mathematical model predicts tumor control patterns induced by fast and slow cytotoxic T lymphocyte killing mechanisms

Immunotherapy has dramatically transformed the cancer treatment landscape largely due to the efficacy of immune checkpoint inhibitors (ICIs). Although ICIs have shown promising results for many patients, the low response rates in many cancers highlight the ongoing challenges in cancer treatment. Cytotoxic T lymphocytes (CTLs) execute their cell-killing function via two distinct mechanisms: a fast-acting, perforin-mediated process and a slower, Fas ligand (FasL)-driven pathway. Evidence also suggests that the preferred killing mechanism of CTLs depends on the antigenicity of tumor cells. To determine the critical factors affecting responses to ICIs, we construct an ordinary differential equation model describing in vivo tumor-immune dynamics in the presence of active or blocked PD-1/PD-L1 immune checkpoint. Specifically, we identify important aspects of the tumor-immune landscape that affect tumor size and composition in the short and long term. We also generate a virtual cohort of mice with diverse tumor and immune attributes to simulate the outcomes of immune checkpoint blockade in a heterogeneous population. By identifying key tumor and immune characteristics associated with tumor elimination, dormancy, and escape, we predict which fraction of a population potentially responds well to ICIs and ways to enhance therapeutic outcomes with combination therapy.

Global stability and parameter analysis reinforce therapeutic targets of PD-L1-PD-1 and MDSCs for glioblastoma

Glioblastoma (GBM) is an aggressive primary brain cancer that currently has minimally effective treatments. Like other cancers, immunosuppression by the PD-L1-PD-1 immune checkpoint complex is a prominent axis by which glioma cells evade the immune system. Myeloid-derived suppressor cells (MDSCs), which are recruited to the glioma microenviroment, also contribute to the immunosuppressed GBM microenvironment by suppressing T cell functions. In this paper, we propose a GBM-specific tumor-immune ordinary differential equations model of glioma cells, T cells, and MDSCs to provide theoretical insights into the interactions between these cells. Equilibrium and stability analysis indicates that there are unique tumorous and tumor-free equilibria which are locally stable under certain conditions. Further, the tumor-free equilibrium is globally stable when T cell activation and the tumor kill rate by T cells overcome tumor growth, T cell inhibition by PD-L1-PD-1 and MDSCs, and the T cell death rate. Bifurcation analysis suggests that a treatment plan that includes surgical resection and therapeutics targeting immune suppression caused by the PD-L1-PD1 complex and MDSCs results in the system tending to the tumor-free equilibrium. Using a set of preclinical experimental data, we implement the approximate Bayesian computation (ABC) rejection method to construct probability density distributions that estimate model parameters. These distributions inform an appropriate search curve for global sensitivity analysis using the extended fourier amplitude sensitivity test. Sensitivity results combined with the ABC method suggest that parameter interaction is occurring between the drivers of tumor burden, which are the tumor growth rate and carrying capacity as well as the tumor kill rate by T cells, and the two modeled forms of immunosuppression, PD-L1-PD-1 immune checkpoint and MDSC suppression of T cells. Thus, treatment with an immune checkpoint inhibitor in combination with a therapeutic targeting the inhibitory mechanisms of MDSCs should be explored.

Oncolytic virus oHSV2 combined with PD-1/PD-L1 inhibitors exert antitumor activity by mediating CD4 + T and CD8 + T cell infiltration in the lymphoma tumor microenvironment

A novel therapeutic regimen showed that the oncolytic type II herpes simplex virus (oHSV2) was able to prevent colorectal cancer growth, recurrence, and metastasis. However, no study has yet explored whether oHSV2 has an impact on the development of diffuse large B-cell lymphoma (DLBCL). We chose the clinical chemotherapeutic drug doxorubicin (DOX) as a positive control to evaluate the effect of oHSV2 infection on the apoptotic, invasive, and proliferative capacity of DLBCL cells. We next further explored the therapeutic efficacy of oncolytic virus oHSV2 or DOX in DLBCL tumor bearing BALB/c mice, and evaluated the infiltration of CD8 + T cells and CD4 + T cells in tumor tissues. A pathological approach was used to explore the effects of oHSV2 on various organs of tumor bearing mice, including the heart, liver, and kidney. Next, SU-DHL-4 cells were co-cultured with cytotoxic T lymphocytes (CTLs) to mimic the tumor immune microenvironment (TME), to explore the impact of oHSV2 on the immune environment at the cellular level, and then analyzed the relationship between oHSV2 and the PD-1/PD-L1 immune-checkpoint. Subsequently, we further validated the efficacy of combined oHSV2 and PD-L1 treatment on transplanted tumor growth in mice at the in vivo level. DLBCL cells were sensitive to the action of the oncolytic virus oHSV2, and the decline in their proliferative activity showed a time-and dose-dependent manner. oHSV2 and DOX intervention preeminently increased the cell apoptosis, restrained cell proliferation and invasion, with the greatest changes occurring in response to oHSV2 infection. oHSV2 application effectively improved the immune status of the tumor microenvironment, favoring the invasion of CD8 + T and CD4 + T cells, thereby enhancing their antitumor effects. Besides, oHSV2 treatment has a safety profile in the organs of tumor bearing mice and indeed inhibits the PD-1/PD-L1 immune checkpoint in DLBCL. Interestingly, the combination of oHSV2 and PD-L1 antibodies results in more profound killing of DLBCL cells than oHSV2 infection alone, with a significant increase in the proportion of CD4 + T cells and CD8 + T cells. The antitumor effect was the best after combining oHSV2 and PD-L1 antibodies, suggesting that the combination therapy of oHSV2 and PD-L1 would have a better prospect for clinical application.

Nonlinear dynamics of estrogen receptor-positive breast cancer integrating experimental data: A novel spatial modeling approach

Oncology research has focused extensively on estrogen hormones and their function in breast cancer proliferation. Mathematical modeling is essential for the analysis and simulation of breast cancers. This research presents a novel approach to examine the therapeutic and inhibitory effects of hormone and estrogen therapies on the onset of breast cancer. Our proposed mathematical model comprises a nonlinear coupled system of partial differential equations, capturing intricate interactions among estrogen, cytotoxic T lymphocytes, dormant cancer cells, and active cancer cells. The model's parameters are meticulously estimated through experimental studies, and we conduct a comprehensive global sensitivity analysis to assess the uncertainty of these parameter values. Remarkably, our findings underscore the pivotal role of hormone therapy in curtailing breast tumor growth by blocking estrogen's influence on cancer cells. Beyond this crucial insight, our proposed model offers an integrated framework to delve into the complexity of tumor progression and immune response under hormone therapy. We employ diverse experimental datasets encompassing gene expression profiles, spatial tumor morphology, and cellular interactions. Integrating multidimensional experimental data with mathematical models enhances our understanding of breast cancer dynamics and paves the way for personalized treatment strategies. Our study advances our comprehension of estrogen receptor-positive breast cancer and exemplifies a transformative approach that merges experimental data with cutting-edge mathematical modeling. This framework promises to illuminate the complexities of cancer progression and therapy, with broad implications for oncology.

Agent-Based Modelling Reveals the Role of the Tumor Microenvironment on the Short-Term Success of Combination Temozolomide/Immune Checkpoint Blockade to Treat Glioblastoma

Glioblastoma is the most common and deadly primary brain tumor in adults. All glioblastoma patients receiving standard-of-care surgery-radiotherapy-chemotherapy (i.e., temozolomide (TMZ)) recur, with an average survival time of only 15 months. New approaches to the treatment of glioblastoma, including immune checkpoint blockade and oncolytic viruses, offer the possibility of improving glioblastoma outcomes and have as such been under intense study. Unfortunately, these treatment modalities have thus far failed to achieve approval. Recently, in an attempt to bolster efficacy and improve patient outcomes, regimens combining chemotherapy and immune checkpoint inhibitors have been tested in trials. Unfortunately, these efforts have not resulted in significant increases to patient survival. To better understand the various factors impacting treatment outcomes of combined TMZ and immune checkpoint blockade, we developed a systems-level, computational model that describes the interplay between glioblastoma, immune, and stromal cells with this combination treatment. Initializing our model to spatial resection patient samples labeled using imaging mass cytometry, our model's predictions show how the localization of glioblastoma cells, influence therapeutic success. We further validated these predictions in samples of brain metastases from patients given they generally respond better to checkpoint blockade compared with primary glioblastoma. Ultimately, our model provides novel insights into the mechanisms of therapeutic success of immune checkpoint inhibitors in brain tumors and delineates strategies to translate combination immunotherapy regimens more effectively into the clinic. SIGNIFICANCE STATEMENT: Extending survival times for glioblastoma patients remains a critical challenge. Although immunotherapies in combination with chemotherapy hold promise, clinical trials have not shown much success. Here, systems models calibrated to and validated against patient samples can improve preclinical and clinical studies by shedding light on the factors distinguishing responses/failures. By initializing our model with imaging mass cytometry visualization of patient samples, we elucidate how factors such as localization of glioblastoma cells and CD8+ T cell infiltration impact treatment outcomes.

Modeling of oncolytic viruses in a heterogeneous cell population to predict spread into non-cancerous cells

New cancer treatment modalities that limit patient discomfort need to be developed. One possible new therapy is the use of oncolytic (cancer-killing) viruses. It is only recently that our ability to manipulate viral genomes has allowed us to consider deliberately infecting cancer patients with viruses. One key consideration is to ensure that the virus exclusively targets cancer cells and does not harm nearby non-cancerous cells. Here, we use a mathematical model of viral infection to determine the characteristics a virus would need to have in order to eradicate a tumor, but leave non-cancerous cells untouched. We conclude that the virus must differ in its ability to infect the two different cell types, with the infection rate of non-cancerous cells needing to be less than one hundredth of the infection rate of cancer cells. Differences in viral production rate or infectious cell death rate alone are not sufficient to protect non-cancerous cells.

Global stability and parameter analysis reinforce therapeutic targets of PD-L1-PD-1 and MDSCs for glioblastoma

Glioblastoma (GBM) is an aggressive primary brain cancer that currently has minimally effective treatments. Like other cancers, immunosuppression by the PD-L1-PD-1 immune checkpoint complex is a prominent axis by which glioma cells evade the immune system. Myeloid-derived suppressor cells (MDSCs), which are recruited to the glioma microenviroment, also contribute to the immunosuppressed GBM microenvironment by suppressing T cell functions. In this paper, we propose a GBM-specific tumor-immune ordinary differential equations model of glioma cells, T cells, and MDSCs to provide theoretical insights into the interactions between these cells. Equilibrium and stability analysis indicates that there are unique tumorous and tumor-free equilibria which are locally stable under certain conditions. Further, the tumor-free equilibrium is globally stable when T cell activation and the tumor kill rate by T cells overcome tumor growth, T cell inhibition by PD-L1-PD-1 and MDSCs, and the T cell death rate. Bifurcation analysis suggests that a treatment plan that includes surgical resection and therapeutics targeting immune suppression caused by the PD-L1-PD1 complex and MDSCs results in the system tending to the tumor-free equilibrium. Using a set of preclinical experimental data, we implement the Approximate Bayesian Computation (ABC) rejection method to construct probability density distributions that estimate model parameters. These distributions inform an appropriate search curve for global sensitivity analysis using the extended Fourier Amplitude Sensitivity Test (eFAST). Sensitivity results combined with the ABC method suggest that parameter interaction is occurring between the drivers of tumor burden, which are the tumor growth rate and carrying capacity as well as the tumor kill rate by T cells, and the two modeled forms of immunosuppression, PD-L1-PD-1 immune checkpoint and MDSC suppression of T cells. Thus, treatment with an immune checkpoint inhibitor in combination with a therapeutic targeting the inhibitory mechanisms of MDSCs should be explored.

Assessment of Clinical Response to V937 Oncolytic Virus After Intravenous or Intratumoral Administration Using Physiologically-Based Modeling

Oncolytic viruses (OVs) represent a potential therapeutic strategy in cancer treatment. However, there is currently a lack of comprehensive quantitative models characterizing clinical OV kinetics and distribution to the tumor. In this work, we present a mechanistic modeling framework for V937 OV, after intratumoral (i.t.) or intravascular (i.v.) administration in patients with cancer. A minimal physiologically-based pharmacokinetic model was built to characterize biodistribution of OVs in humans. Viral dynamics was incorporated at the i.t. cellular level and linked to tumor response, enabling the characterization of a direct OV killing triggered by the death of infected tumor cells and an indirect killing induced by the immune response. The model provided an adequate description of changes in V937 mRNA levels and tumor size obtained from phase I/II clinical trials after V937 administration. The model showed prominent role of viral clearance from systemic circulation and infectivity in addition to known tumor aggressiveness on clinical response. After i.v. administration, i.t. exposure of V937 was predicted to be several orders of magnitude lower compared with i.t. administration. These differences could be overcome if there is high virus infectivity and/or replication. Unfortunately, the latter process could not be identified at the current clinical setting. This work provides insights on selecting optimal OV considering replication rate and infectivity.

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.

Global asymptotic stability of a delay differential equation model for SARS-CoV-2 virus infection mediated by ACE2 receptor protein

The COVID-19 pandemic has brought a serious threat to human life safety worldwide. SARS-CoV-2 virus mainly binds to the target cell surface receptor ACE2 (Angiotensin-converting enzyme 2 ) through the S protein expressed on the surface of the virus, resulting in infection of target cells. During this infection process, the target cell ACE2 receptor plays a very important mediating role. In this paper, a delay differential equation model containing the mediated effect of target cell receptor is established according to the mechanism of SARS-CoV-2 virus invasion of target cells, and the global stability of the infection-free equilibrium and the infected equilibrium of the model is obtained by using the basic reproduction number  ℛ 0  and constructing the appropriate Lyapunov functional. The expression of the basic reproduction number  ℛ 0  intuitively gives the dependence on the expression ratio of the target cell surface ACE2 receptor, which is helpful for the understanding of the mechanism of SARS-CoV-2 virus infection.

Exploring the Interactions of Oncolytic Viral Therapy and Immunotherapy of Anti-CTLA-4 for Malignant Melanoma Mice Model

Oncolytic ability to direct target and lyse tumor cells makes oncolytic virus therapy (OVT) a promising approach to treating cancer. Despite its therapeutic potential to stimulate anti-tumor immune responses, it also has immunosuppressive effects. The efficacy of OVTs as monotherapies can be enhanced by appropriate adjuvant therapy such as anti-CTLA-4. In this paper, we propose a mathematical model to explore the interactions of combined therapy of oncolytic viruses and a checkpoint inhibitor, anti-CTLA-4. The model incorporates both the susceptible and infected tumor populations, natural killer cell population, virus population, tumor-specific immune populations, virus-specific immune populations, tumor suppressive cytokine IFN-g, and the effect of immune checkpoint inhibitor CTLA-4. In particular, we distinguish the tumor-specific immune abilities of CD8⁺ T, NK cells, and CD4⁺ T cells and describe the destructive ability of cytokine on tumor cells as well as the inhibitory capacity of CTLA-4 on various components. Our model is validated through the experimental results. We also investigate various dosing strategies to improve treatment outcomes. Our study reveals that tumor killing rate by cytokines, cytokine decay rate, and tumor growth rate play important roles on both the OVT monotherapy and the combination therapy. Moreover, parameters related to CD8⁺ T cell killing have a large impact on treatment outcomes with OVT alone, whereas parameters associated with IFN-g strongly influence treatment responses for the combined therapy. We also found that virus killing by NK cells may halt the desired spread of OVs and enhance the probability of tumor escape during the treatment. Our study reveals that it is the activation of host anti-tumor immune system responses rather than its direct destruction of the tumor cells plays a major biological function of the combined therapy.

Oncolytic viral kinetics mechanistic modeling of Talimogene Laherparepvec (T-VEC) a first-in-class oncolytic viral therapy in patients with advanced melanoma

Talimogene Laherparepvec (T-VEC) is a first-in-class oncolytic virotherapy approved for the treatment of unresectable melanoma recurrent after initial surgery. Biodistribution data from a phase II study was used to develop a viral kinetic mechanistic model describing the interaction between cytokines such as granulocyte-macrophage colony-stimulating factor (GM-CSF), the immune system, and T-VEC treatment. Our analysis found that (1) the viral infection rate has a great influence on T-VEC treatment efficacy; (2) an increase in T-VEC dose of 10² plaque-forming units/ml 21 days and beyond after the initial dose of T-VEC resulted in an ~12% increase in response; and (3) at the systemic level, the ratio of resting innate immune cells to the death rate of innate immune impact T-VEC treatment efficacy. This analysis clarifies under which condition the immune system either assists in eliminating tumor cells or inhibits T-VEC treatment efficacy, which is critical to both efficiently design future oncolytic agents and understand cancer development.

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.

Stochastic Analysis of Nonlinear Cancer Disease Model through Virotherapy and Computational Methods

Cancer is a common term for many diseases that can affect anybody. A worldwide leading cause of death is cancer, according to the World Health Organization (WHO) report. In 2020, ten million people died from cancer. This model identifies the interaction of cancer cells, viral therapy, and immune response. In this model, the cell population has four parts, namely uninfected cells (x), infected cells (y), virus-free cells (v), and immune cells (z). This study presents the analysis of the stochastic cancer virotherapy model in the cell population dynamics. The model results have restored the properties of the biological problem, such as dynamical consistency, positivity, and boundedness, which are the considerable requirements of the models in these fields. The existing computational methods, such as the Euler Maruyama, Stochastic Euler, and Stochastic Runge Kutta, fail to restore the abovementioned properties. The proposed stochastic nonstandard finite difference method is efficient, cost-effective, and accommodates all the desired feasible properties. The existing standard stochastic methods converge conditionally or diverge in the long run. The solution by the nonstandard finite difference method is stable and convergent over all time steps.

A validated mathematical model of FGFR3-mediated tumor growth reveals pathways to harness the benefits of combination targeted therapy and immunotherapy in bladder cancer

Bladder cancer is a common malignancy with over 80,000 estimated new cases and nearly 18,000 deaths per year in the United States alone. Therapeutic options for metastatic bladder cancer had not evolved much for nearly four decades, until recently, when five immune checkpoint inhibitors were approved by the U.S. Food and Drug Administration (FDA). Despite the activity of these drugs in some patients, the objective response rate for each is less than 25%. At the same time, fibroblast growth factor receptors (FGFRs) have been attractive drug targets for a variety of cancers, and in 2019 the FDA approved the first therapy targeted against FGFR3 for bladder cancer. Given the excitement around these new receptor tyrosine kinase and immune checkpoint targeted strategies, and the challenges they each may face on their own, emerging data suggest that combining these treatment options could lead to improved therapeutic outcomes. In this paper, we develop a mathematical model for FGFR3-mediated tumor growth and use it to investigate the impact of the combined administration of a small molecule inhibitor of FGFR3 and a monoclonal antibody against the PD-1/PD-L1 immune checkpoint. The model is carefully calibrated and validated with experimental data before survival benefits, and dosing schedules are explored. Predictions of the model suggest that FGFR3 mutation reduces the effectiveness of anti-PD-L1 therapy, that there are regions of parameter space where each monotherapy can outperform the other, and that pretreatment with anti-PD-L1 therapy always results in greater tumor reduction even when anti-FGFR3 therapy is the more effective monotherapy.

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.

An Agent-Based Model of Combination Oncolytic Viral Therapy and Anti-PD-1 Immunotherapy Reveals the Importance of Spatial Location When Treating Glioblastoma

Simple Summary

A combination of oncolytic viral therapy and immunotherapy provides an alternative option to the standard of care for treating the lethal brain tumor glioblastoma (GBM). Although this combination therapy shows promise, there are many unknown questions regarding how the tumor landscape and spatial dosing strategies impact the effectiveness of the treatment. Our study aims to shed light on these questions using a novel spatially explicit computational model of GBM response to treatment. Our results suggest that oncolytic viral dosing in the location of highest tumor cell density leads to substantial tumor size reduction over viral dosing in the center of the tumor. These results can help to inform future clinical trials and more effective treatment strategies for oncolytic viral therapy in GBM patients.

Abstract

Oncolytic viral therapies and immunotherapies are of growing clinical interest due to their selectivity for tumor cells over healthy cells and their immunostimulatory properties. These treatment modalities provide promising alternatives to the standard of care, particularly for cancers with poor prognoses, such as the lethal brain tumor glioblastoma (GBM). However, uncertainty remains regarding optimal dosing strategies, including how the spatial location of viral doses impacts therapeutic efficacy and tumor landscape characteristics that are most conducive to producing an effective immune response. We develop a three-dimensional agent-based model (ABM) of GBM undergoing treatment with a combination of an oncolytic Herpes Simplex Virus and an anti-PD-1 immunotherapy. We use a mechanistic approach to model the interactions between distinct populations of immune cells, incorporating both innate and adaptive immune responses to oncolytic viral therapy and including a mechanism of adaptive immune suppression via the PD-1/PD-L1 checkpoint pathway. We utilize the spatially explicit nature of the ABM to determine optimal viral dosing in both the temporal and spatial contexts. After proposing an adaptive viral dosing strategy that chooses to dose sites at the location of highest tumor cell density, we find that, in most cases, this adaptive strategy produces a more effective treatment outcome than repeatedly dosing in the center of the tumor.

Computation of the conditions for anti-angiogenesis and gene therapy synergistic effects: Sensitivity analysis and robustness of target solutions

Both anti-angiogenesis and gene therapy involve complex processes depending on non-point parameters belonging to a space of values. To successfully overcome the challenges involved in their therapeutic approaches, there is a need to analyze the sensitivity of these parameters. In this paper, a new mathematical model that combines immune system stimulations, inflammatory processes associated with tumor development, and gene therapy aimed at enhancing the efficacy of both treatments are explored. Using the global sensitivity methods of Sobol and Morris, the most important parameters are estimated. Estimation of the sensitivity variance revealed a strong interdependence between the parameters. Also, determinations of the conditions for effective therapy lead to a target of reducing the cancer cell numbers by at least 50%. This opened the way for delimiting the parameter spaces making it possible to reach the treatment target in addition to enhancing the estimation of the minimum time of remission. The combination of therapies and sensitivity analysis have demonstrated the robustness of therapy success.

Natural Killer Cells Recruitment in Oncolytic Virotherapy: A Mathematical Model

In this paper, we investigate how natural killer (NK) cell recruitment to the tumor microenvironment (TME) affects oncolytic virotherapy. NK cells play a major role against viral infections. They are, however, known to induce early viral clearance of oncolytic viruses, which hinders the overall efficacy of oncolytic virotherapy. Here, we formulate and analyze a simple mathematical model of the dynamics of the tumor, OV and NK cells using currently available preclinical information. The aim of this study is to characterize conditions under which the synergistic balance between OV-induced NK responses and required viral cytopathicity may or may not result in a successful treatment. In this study, we found that NK cell recruitment to the TME must take place neither too early nor too late in the course of OV infection so that treatment will be successful. NK cell responses are most influential at either early (partly because of rapid response of NK cells to viral infections or antigens) or later (partly because of antitumoral ability of NK cells) stages of oncolytic virotherapy. The model also predicts that: (a) an NK cell response augments oncolytic virotherapy only if viral cytopathicity is weak; (b) the recruitment of NK cells modulates tumor growth; and (c) the depletion of activated NK cells within the TME enhances the probability of tumor escape in oncolytic virotherapy. Taken together, our model results demonstrate that OV infection is crucial, not just to cytoreduce tumor burden, but also to induce the stronger NK cell response necessary to achieve complete or at least partial tumor remission. Furthermore, our modeling framework supports combination therapies involving NK cells and OV which are currently used in oncolytic immunovirotherapy to treat several cancer types.

The Current State of Oncolytic Herpes Simplex Virus for Glioblastoma Treatment

Glioblastoma (GBM) is a lethal primary malignant brain tumor with no current effective treatments. The recent emergence of immuno-virotherapy and FDA approval of T-VEC have generated a great expectation towards oncolytic herpes simplex viruses (oHSVs) as a promising treatment option for GBM. Since the generation and testing of the first genetically engineered oHSV in glioma in the early 1990s, oHSV-based therapies have shown a long way of great progress in terms of anti-GBM efficacy and safety, both preclinically and clinically. Here, we revisit the literature to understand the recent advancement of oHSV in the treatment of GBM. In addition, we discuss current obstacles to oHSV-based therapies and possible strategies to overcome these pitfalls.