The role of memory in non-genetic inheritance and its impact on cancer treatment resistance.

Intra-tumour heterogeneity is a leading cause of treatment failure and disease progression in cancer. While genetic mutations have long been accepted as a primary mechanism of generating this heterogeneity, the role of phenotypic plasticity is becoming increasingly apparent as a driver of intra-tumour heterogeneity. Consequently, understanding the role of this plasticity in treatment resistance and failure is a key component of improving cancer therapy. We develop a mathematical model of stochastic phenotype switching that tracks the evolution of drug-sensitive and drug-tolerant subpopulations to clarify the role of phenotype switching on population growth rates and tumour persistence. By including cytotoxic therapy in the model, we show that, depending on the strategy of the drug-tolerant subpopulation, stochastic phenotype switching can lead to either transient or permanent drug resistance. We study the role of phenotypic heterogeneity in a drug-resistant, genetically homogeneous population of non-small cell lung cancer cells to derive a rational treatment schedule that drives population extinction and avoids competitive release of the drug-tolerant sub-population. This model-informed therapeutic schedule results in increased treatment efficacy when compared against periodic therapy, and, most importantly, sustained tumour decay without the development of resistance.

In silico trials predict that combination strategies for enhancing vesicular stomatitis oncolytic virus are determined by tumor aggressivity.

BACKGROUND: Immunotherapies, driven by immune-mediated antitumorigenicity, offer the potential for significant improvements to the treatment of multiple cancer types. Identifying therapeutic strategies that bolster antitumor immunity while limiting immune suppression is critical to selecting treatment combinations and schedules that offer durable therapeutic benefits. Combination oncolytic virus (OV) therapy, wherein complementary OVs are administered in succession, offer such promise, yet their translation from preclinical studies to clinical implementation is a major challenge. Overcoming this obstacle requires answering fundamental questions about how to effectively design and tailor schedules to provide the most benefit to patients. METHODS: We developed a computational biology model of combined oncolytic vaccinia (an enhancer virus) and vesicular stomatitis virus (VSV) calibrated to and validated against multiple data sources. We then optimized protocols in a cohort of heterogeneous virtual individuals by leveraging this model and our previously established in silico clinical trial platform. RESULTS: Enhancer multiplicity was shown to have little to no impact on the average response to therapy. However, the duration of the VSV injection lag was found to be determinant for survival outcomes. Importantly, through treatment individualization, we found that optimal combination schedules are closely linked to tumor aggressivity. We predicted that patients with aggressively growing tumors required a single enhancer followed by a VSV injection 1 day later, whereas a small subset of patients with the slowest growing tumors needed multiple enhancers followed by a longer VSV delay of 15 days, suggesting that intrinsic tumor growth rates could inform the segregation of patients into clinical trials and ultimately determine patient survival. These results were validated in entirely new cohorts of virtual individuals with aggressive or non-aggressive subtypes. CONCLUSIONS: Based on our results, improved therapeutic schedules for combinations with enhancer OVs can be studied and implemented. Our results further underline the impact of interdisciplinary approaches to preclinical planning and the importance of computational approaches to drug discovery and development.

Characterizing Chemotherapy-Induced Neutropenia and Monocytopenia Through Mathematical Modelling.

In spite of the recent focus on the development of novel targeted drugs to treat cancer, cytotoxic chemotherapy remains the standard treatment for the vast majority of patients. Unfortunately, chemotherapy is associated with high hematopoietic toxicity that may limit its efficacy. We have previously established potential strategies to mitigate chemotherapy-induced neutropenia (a lack of circulating neutrophils) using a mechanistic model of granulopoiesis to predict the interactions defining the neutrophil response to chemotherapy and to define optimal strategies for concurrent chemotherapy/prophylactic granulocyte colony-stimulating factor (G-CSF). Here, we extend our analyses to include monocyte production by constructing and parameterizing a model of monocytopoiesis. Using data for neutrophil and monocyte concentrations during chemotherapy in a large cohort of childhood acute lymphoblastic leukemia patients, we leveraged our model to determine the relationship between the monocyte and neutrophil nadirs during cyclic chemotherapy. We show that monocytopenia precedes neutropenia by 3 days, and rationalize the use of G-CSF during chemotherapy by establishing that the onset of monocytopenia can be used as a clinical marker for G-CSF dosing post-chemotherapy. This work therefore has important clinical applications as a comprehensive approach to understanding the relationship between monocyte and neutrophils after cyclic chemotherapy with or without G-CSF support.

A mathematical model of viral oncology as an immuno-oncology instigator.

We develop and analyse a mathematical model of tumour-immune interaction that explicitly incorporates heterogeneity in tumour cell cycle duration by using a distributed delay differential equation. We derive a necessary and sufficient condition for local stability of the cancer-free equilibrium in which the amount of tumour-immune interaction completely characterizes disease progression. Consistent with the immunoediting hypothesis, we show that decreasing tumour-immune interaction leads to tumour expansion. Finally, by simulating the mathematical model, we show that the strength of tumour-immune interaction determines the long-term success or failure of viral therapy.

Determinants of combination GM-CSF immunotherapy and oncolytic virotherapy success identified through in silico treatment personalization.

Oncolytic virotherapies, including the modified herpes simplex virus talimogene laherparepvec (T-VEC), have shown great promise as potent instigators of anti-tumour immune effects. The OPTiM trial, in particular, demonstrated the superior anti-cancer effects of T-VEC as compared to systemic immunotherapy treatment using exogenous administration of granulocyte-macrophage colony-stimulating factor (GM-CSF). Theoretically, a combined approach leveraging exogenous cytokine immunotherapy and oncolytic virotherapy would elicit an even greater immune response and improve patient outcomes. However, regimen scheduling of combination immunostimulation and T-VEC therapy has yet to be established. Here, we calibrate a computational biology model of sensitive and resistant tumour cells and immune interactions for implementation into an in silico clinical trial to test and individualize combination immuno- and virotherapy. By personalizing and optimizing combination oncolytic virotherapy and immunostimulatory therapy, we show improved simulated patient outcomes for individuals with late-stage melanoma. More crucially, through evaluation of individualized regimens, we identified determinants of combination GM-CSF and T-VEC therapy that can be translated into clinically-actionable dosing strategies without further personalization. Our results serve as a proof-of-concept for interdisciplinary approaches to determining combination therapy, and suggest promising avenues of investigation towards tailored combination immunotherapy/oncolytic virotherapy.

Transit and lifespan in neutrophil production: implications for drug intervention.

A comparison of the transit compartment ordinary differential equation modelling approach to distributed and discrete delay differential equation models is studied by focusing on Quartino's extension to the Friberg transit compartment model of myelosuppression, widely relied upon in the pharmaceutical sciences to predict the neutrophil response after chemotherapy, and on a QSP delay differential equation model of granulopoiesis. An extension to the Quartino model is provided by considering a general number of transit compartments and introducing an extra parameter that allows for the decoupling of the maturation time from the production rate of cells. An overview of the well established linear chain technique, used to reformulate transit compartment models with constant transit rates as distributed delay differential equations (DDEs), is then given. A state-dependent time rescaling of the Quartino model is performed to apply the linear chain technique and rewrite the Quartino model as a distributed DDE, yielding a discrete DDE model in a certain parameter limit. Next, stability and bifurcation analyses are undertaken in an effort to situate such studies in a mathematical pharmacology context. We show that both the original Friberg and the Quartino extension models incorrectly define the mean maturation time, essentially treating the proliferative pool as an additional maturation compartment. This misspecification can have far reaching consequences on the development of future models of myelosuppression in PK/PD.