Oncolytic viruses are of growing interest to cancer researchers and clinicians, due to their selectivity for tumor cells over healthy cells and their immunostimulatory properties. The immune response to an oncolytic virus plays a critical role in treatment efficacy. However, uncertainty remains regarding the circumstances under which the immune system either assists in eliminating tumor cells or inhibits treatment via rapid viral clearance, leading to the cessation of the immune response. In this work, we develop an ordinary differential equation model of treatment for a lethal brain tumor, glioblastoma, using an oncolytic Herpes Simplex Virus. 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 (OVT), and including a mechanism of adaptive immune suppression via the PD-1/PD-L1 checkpoint pathway. We focus on the tradeoff between viral clearance by innate immune cells and the innate immune cell-mediated recruitment of antiviral and antitumor adaptive immune cells. Our model suggests that when a tumor is treated with OVT alone, the innate immune cells' ability to clear the virus quickly after administration has a much larger impact on the treatment outcome than the adaptive immune cells' antitumor activity. Even in a highly antigenic tumor with a strong innate immune response, the faster recruitment of antitumor adaptive immune cells is not sufficient to offset the rapid viral clearance. This motivates our subsequent incorporation of an immunotherapy that inhibits the PD-1/PD-L1 checkpoint pathway by blocking PD-1, which we combine with OVT within the model. The combination therapy is most effective for a highly antigenic tumor or for intermediate levels of innate immune localization. Extreme levels of innate immune cell activity either clear the virus too quickly or fail to activate a sufficiently strong adaptive response, yielding ineffective combination therapy of GBM. Hence, we show that the innate and adaptive immune interactions significantly influence treatment response and that combining OVT with an immune checkpoint inhibitor expands the range of immune conditions that allow for tumor size reduction or clearance.

We use topological data analysis and machine learning to study a seminal model of collective motion in biology [M. R. D’Orsogna et al., Phys. Rev. Lett. 96, 104302 (2006)]. This model describes agents interacting nonlinearly via attractive-repulsive social forces and gives rise to collective behaviors such as flocking and milling. To classify the emergent collective motion in a large library of numerical simulations and to recover model parameters from the simulation data, we apply machine learning techniques to two different types of input. First, we input time series of order parameters traditionally used in studies of collective motion. Second, we input measures based on topology that summarize the time-varying persistent homology of simulation data over multiple scales. This topological approach does not require prior knowledge of the expected patterns. For both unsupervised and supervised machine learning methods, the topological approach outperforms the one that is based on traditional order parameters.

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.

Tumor recurrence in glioblastoma multiforme (GBM) is often attributed to acquired resistance to the standard chemotherapeutic agent, temozolomide (TMZ). Promoter methylation of the DNA repair gene MGMT (O⁶-methylguanine–DNA methyltransferase) has been associated with sensitivity to TMZ, whereas increased expression of MGMT has been associated with TMZ resistance. Clinical studies have observed a downward shift in MGMT methylation percentage from primary to recurrent stage tumors; however, the evolutionary processes that drive this shift and more generally the emergence and growth of TMZ-resistant tumor subpopulations are still poorly understood. Here, we develop a mathematical model, parameterized using clinical and experimental data, to investigate the role of MGMT methylation in TMZ resistance during the standard treatment regimen for GBM—surgery, chemotherapy, and radiation. We first found that the observed downward shift in MGMT promoter methylation status between detection and recurrence cannot be explained solely by evolutionary selection. Next, our model suggests that TMZ has an inhibitory effect on maintenance methylation of MGMT after cell division. Finally, incorporating this inhibitory effect, we study the optimal number of TMZ doses per adjuvant cycle for patients with GBM with high and low levels of MGMT methylation at diagnosis.

Neuroendocrine neoplasms (NEN) encompass a diverse set of malignancies with limited precision therapy options. Recently, therapies targeting DLL3 have shown clinical efficacy in aggressive NENs, including small cell lung cancers and neuroendocrine prostate cancers. Given the continued development and expansion of DLL3-targeted therapies, we sought to characterize the expression of DLL3 and identify its clinical and molecular correlates across diverse neuroendocrine and non-neuroendocrine cancers. Here, we interrogated paired DNA and RNA-sequencing from 1,589 NENs across 29 sites, as well as 203,252 tumors across 47 cancer types. We found that high transcriptomic levels of DLL3 correlated with more aggressive histologic and mutational patterns in NENs, with adverse survival outcomes being reflected in NENs originating from the lung, pancreas, stomach, and small bowel. The heterogeneity in DLL3 expression across NENs was largely explained by site of origin, with lung, prostate, and bladder NENs exhibiting relatively high levels of DLL3, whereas gastroenteropancreatic NENs displayed relatively low expression levels. Although the therapeutic targeting of DLL3 may be less applicable for gastroenteropancreatic NENs, we did find an upregulation of alternative targets such as SEZ6, CELSR3, and SSTR2 in these settings. Lastly, expanding our investigation into non-neuroendocrine cancers, we detected an enrichment of DLL3 in both low-grade and high-grade gliomas, Merkel cell carcinomas, medulloblastomas, and melanomas, with such enrichment being associated with prolonged overall survival in gliomas, but worse overall survival in melanomas. Altogether, we demonstrate that DLL3 represents an attractive target for subsets of neuroendocrine and non-neuroendocrine cancers and uncover opportunities for future therapeutic strategies. Significance: DLL3-targeted therapies have recently shown robust clinical efficacy in aggressive neuroendocrine cancers, positioning them to fulfill a great unmet need in these settings. Here, we examine the clinical and biological correlates of DLL3 expression in both neuroendocrine and non-neuroendocrine cancers. Our findings may stimulate the development and application of DLL3-targeted therapies, as well as other precision therapies, in neuroendocrine cancers and beyond.