Modeling treatment-dependent glioma growth including a dormant tumor cell subpopulation

Understanding and leveraging phenotypic plasticity during metastasis formation

Cancer metastasis is the process of detrimental systemic spread and the primary cause of cancer-related fatalities. Successful metastasis formation requires tumor cells to be proliferative and invasive; however, cells cannot be effective at both tasks simultaneously. Tumor cells compensate for this trade-off by changing their phenotype during metastasis formation through phenotypic plasticity. Given the changing selection pressures and competitive interactions that tumor cells face, it is poorly understood how plasticity shapes the process of metastasis formation. Here, we develop an ecology-inspired mathematical model with phenotypic plasticity and resource competition between phenotypes to address this knowledge gap. We find that phenotypically plastic tumor cell populations attain a stable phenotype equilibrium that maintains tumor cell heterogeneity. Considering treatment types inspired by chemo- and immunotherapy, we highlight that plasticity can protect tumors against interventions. Turning this strength into a weakness, we corroborate current clinical practices to use plasticity as a target for adjuvant therapy. We present a parsimonious view of tumor plasticity-driven metastasis that is quantitative and experimentally testable, and thus potentially improving the mechanistic understanding of metastasis at the cell population level, and its treatment consequences.

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.

Agent-based computational modeling of glioblastoma predicts that stromal density is central to oncolytic virus efficacy

Oncolytic viruses (OVs) are emerging cancer immunotherapy. Despite notable successes in the treatment of some tumors, OV therapy for central nervous system cancers has failed to show efficacy. We used an ex vivo tumor model developed from human glioblastoma tissue to evaluate the infiltration of herpes simplex OV rQNestin (oHSV-1) into glioblastoma tumors. We next leveraged our data to develop a computational, model of glioblastoma dynamics that accounts for cellular interactions within the tumor. Using our computational model, we found that low stromal density was highly predictive of oHSV-1 therapeutic success, suggesting that the efficacy of oHSV-1 in glioblastoma may be determined by stromal-to-tumor cell regional density. We validated these findings in heterogenous patient samples from brain metastatic adenocarcinoma. Our integrated modeling strategy can be applied to suggest mechanisms of therapeutic responses for central nervous system cancers and to facilitate the successful translation of OVs into the clinic.

Glioblastoma Recurrence and the Role of O6-Methylguanine-DNA Methyltransferase Promoter Methylation

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.

From cells to tissue: How cell scale heterogeneity impacts glioblastoma growth and treatment response

Glioblastomas are aggressive primary brain tumors known for their inter- and intratumor heterogeneity. This disease is uniformly fatal, with intratumor heterogeneity the major reason for treatment failure and recurrence. Just like the nature vs nurture debate, heterogeneity can arise from intrinsic or environmental influences. Whilst it is impossible to clinically separate observed behavior of cells from their environmental context, using a mathematical framework combined with multiscale data gives us insight into the relative roles of variation from different sources. To better understand the implications of intratumor heterogeneity on therapeutic outcomes, we created a hybrid agent-based mathematical model that captures both the overall tumor kinetics and the individual cellular behavior. We track single cells as agents, cell density on a coarser scale, and growth factor diffusion and dynamics on a finer scale over time and space. Our model parameters were fit utilizing serial MRI imaging and cell tracking data from ex vivo tissue slices acquired from a growth-factor driven glioblastoma murine model. When fitting our model to serial imaging only, there was a spectrum of equally-good parameter fits corresponding to a wide range of phenotypic behaviors. When fitting our model using imaging and cell scale data, we determined that environmental heterogeneity alone is insufficient to match the single cell data, and intrinsic heterogeneity is required to fully capture the migration behavior. The wide spectrum of in silico tumors also had a wide variety of responses to an application of an anti-proliferative treatment. Recurrent tumors were generally less proliferative than pre-treatment tumors as measured via the model simulations and validated from human GBM patient histology. Further, we found that all tumors continued to grow with an anti-migratory treatment alone, but the anti-proliferative/anti-migratory combination generally showed improvement over an anti-proliferative treatment alone. Together our results emphasize the need to better understand the underlying phenotypes and tumor heterogeneity present in a tumor when designing therapeutic regimens.

Integrated cancer tissue engineering models for precision medicine

Tumors are not merely cancerous cells that undergo mindless proliferation. Rather, they are highly organized and interconnected organ systems. Tumor cells reside in complex microenvironments in which they are subjected to a variety of physical and chemical stimuli that influence cell behavior and ultimately the progression and maintenance of the tumor. As cancer bioengineers, it is our responsibility to create physiologic models that enable accurate understanding of the multi-dimensional structure, organization, and complex relationships in diverse tumor microenvironments. Such models can greatly expedite clinical discovery and translation by closely replicating the physiological conditions while maintaining high tunability and control of extrinsic factors. In this review, we discuss the current models that target key aspects of the tumor microenvironment and their role in cancer progression. In order to address sources of experimental variation and model limitations, we also make recommendations for methods to improve overall physiologic reproducibility, experimental repeatability, and rigor within the field. Improvements can be made through an enhanced emphasis on mathematical modeling, standardized in vitro model characterization, transparent reporting of methodologies, and designing experiments with physiological metrics. Taken together these considerations will enhance the relevance of in vitro tumor models, biological understanding, and accelerate treatment exploration ultimately leading to improved clinical outcomes. Moreover, the development of robust, user-friendly models that integrate important stimuli will allow for the in-depth study of tumors as they undergo progression from non-transformed primary cells to metastatic disease and facilitate translation to a wide variety of biological and clinical studies.