Andrea Hawkins-Daarud

Gene therapy enhances chemotherapy tolerance and efficacy in glioblastoma patients

BACKGROUND Temozolomide (TMZ) is one of the most potent chemotherapy agents for the treatment of glioblastoma. Unfortunately, almost half of glioblastoma tumors are TMZ resistant due to overexpression of methylguanine methyltransferase (MGMT(hi)). Coadministration of O6-benzylguanine (O6BG) can restore TMZ sensitivity, but causes off-target myelosuppression. Here, we conducted a prospective clinical trial to test whether gene therapy to confer O6BG resistance in hematopoietic stem cells (HSCs) improves chemotherapy tolerance and outcome. METHODS We enrolled 7 newly diagnosed glioblastoma patients with MGMT(hi) tumors. Patients received autologous gene-modified HSCs following single-agent carmustine administration. After hematopoietic recovery, patients underwent O6BG/TMZ chemotherapy in 28-day cycles. Serial blood samples and tumor images were collected throughout the study. Chemotherapy tolerance was determined by the observed myelosuppression and recovery following each cycle. Patient-specific biomathematical modeling of tumor growth was performed. Progression-free survival (PFS) and overall survival (OS) were also evaluated. RESULTS Gene therapy permitted a significant increase in the mean number of tolerated O6BG/TMZ cycles (4.4 cycles per patient, P < 0.05) compared with historical controls without gene therapy (n = 7 patients, 1.7 cycles per patient). One patient tolerated an unprecedented 9 cycles and demonstrated long-term PFS without additional therapy. Overall, we observed a median PFS of 9 (range 3.5-57+) months and OS of 20 (range 13-57+) months. Furthermore, biomathematical modeling revealed markedly delayed tumor growth at lower cumulative TMZ doses in study patients compared with patients that received standard TMZ regimens without O6BG. CONCLUSION These data support further development of chemoprotective gene therapy in combination with O6BG and TMZ for the treatment of glioblastoma and potentially other tumors with overexpression of MGMT. TRIAL REGISTRATION Clinicaltrials.gov NCT00669669. FUNDING R01CA114218, R01AI080326, R01HL098489, P30DK056465, K01DK076973, R01HL074162, R01CA164371, R01NS060752, U54CA143970.

SELECTION AND ASSESSMENT OF PHENOMENOLOGICAL MODELS OF TUMOR GROWTH

We address general approaches to the rational selection and validation of mathematical and computational models of tumor growth using methods of Bayesian inference. The model classes are derived from a general diffuse-interface, continuum mixture theory and focus on mass conservation of mixtures with up to four species. Synthetic data are generated using higher-order base models. We discuss general approaches to model calibration, validation, plausibility, and selection based on Bayesian-based methods, information theory, and maximum information entropy. We also address computational issues and provide numerical experiments based on Markov chain Monte Carlo algorithms and high performance computing implementations.

Bayesian calibration, validation, and uncertainty quantification of diffuse interface models of tumor growth

The idea that one can possibly develop computational models that predict the emergence, growth, or decline of tumors in living tissue is enormously intriguing as such predictions could revolutionize medicine and bring a new paradigm into the treatment and prevention of a class of the deadliest maladies affecting humankind. But at the heart of this subject is the notion of predictability itself, the ambiguity involved in selecting and implementing effective models, and the acquisition of relevant data, all factors that contribute to the difficulty of predicting such complex events as tumor growth with quantifiable uncertainty. In this work, we attempt to lay out a framework, based on Bayesian probability, for systematically addressing the questions of Validation, the process of investigating the accuracy with which a mathematical model is able to reproduce particular physical events, and Uncertainty quantification, developing measures of the degree of confidence with which a computer model predicts particular quantities of interest. For illustrative purposes, we exercise the process using virtual data for models of tumor growth based on diffuse-interface theories of mixtures utilizing virtual data.

A patient-specific computational model of hypoxia-modulated radiation resistance in glioblastoma using 18F-FMISO-PET

Glioblastoma multiforme (GBM) is a highly invasive primary brain tumour that has poor prognosis despite aggressive treatment. A hallmark of these tumours is diffuse invasion into the surrounding brain, necessitating a multi-modal treatment approach, including surgery, radiation and chemotherapy. We have previously demonstrated the ability of our model to predict radiographic response immediately following radiation therapy in individual GBM patients using a simplified geometry of the brain and theoretical radiation dose. Using only two pre-treatment magnetic resonance imaging scans, we calculate net rates of proliferation and invasion as well as radiation sensitivity for a patients disease. Here, we present the application of our clinically targeted modelling approach to a single glioblastoma patient as a demonstration of our method. We apply our model in the full three-dimensional architecture of the brain to quantify the effects of regional resistance to radiation owing to hypoxia in vivo determined by [18F]-fluoromisonidazole positron emission tomography (FMISO-PET) and the patient-specific three-dimensional radiation treatment plan. Incorporation of hypoxia into our model with FMISO-PET increases the model–data agreement by an order of magnitude. This improvement was robust to our definition of hypoxia or the degree of radiation resistance quantified with the FMISO-PET image and our computational model, respectively. This work demonstrates a useful application of patient-specific modelling in personalized medicine and how mathematical modelling has the potential to unify multi-modality imaging and radiation treatment planning.

From Patient-Specific Mathematical Neuro-Oncology to Precision Medicine

Gliomas are notoriously aggressive, malignant brain tumors that have variable response to treatment. These patients often have poor prognosis, informed primarily by histopathology. Mathematical neuro-oncology (MNO) is a young and burgeoning field that leverages mathematical models to predict and quantify response to therapies. These mathematical models can form the basis of modern “precision medicine” approaches to tailor therapy in a patient-specific manner. Patient-specific models (PSMs) can be used to overcome imaging limitations, improve prognostic predictions, stratify patients, and assess treatment response in silico. The information gleaned from such models can aid in the construction and efficacy of clinical trials and treatment protocols, accelerating the pace of clinical research in the war on cancer. This review focuses on the growing translation of PSM to clinical neuro-oncology. It will also provide a forward-looking view on a new era of patient-specific MNO.

Modeling Tumor-Associated Edema in Gliomas during Anti-Angiogenic Therapy and Its Impact on Imageable Tumor.

Glioblastoma, the most aggressive form of primary brain tumor, is predominantly assessed with gadolinium-enhanced T1-weighted (T1Gd) and T2-weighted magnetic resonance imaging (MRI). Pixel intensity enhancement on the T1Gd image is understood to correspond to the gadolinium contrast agent leaking from the tumor-induced neovasculature, while hyperintensity on the T2/FLAIR images corresponds with edema and infiltrated tumor cells. None of these modalities directly show tumor cells; rather, they capture abnormalities in the microenvironment caused by the presence of tumor cells. Thus, assessing disease response after treatments impacting the microenvironment remains challenging through the obscuring lens of MR imaging. Anti-angiogenic therapies have been used in the treatment of gliomas with spurious results ranging from no apparent response to significant imaging improvement with the potential for extremely diffuse patterns of tumor recurrence on imaging and autopsy. Anti-angiogenic treatment normalizes the vasculature, effectively decreasing vessel permeability and thus reducing tumor-induced edema, drastically altering T2-weighted MRI. We extend a previously developed mathematical model of glioma growth to explicitly incorporate edema formation allowing us to directly characterize and potentially predict the effects of anti-angiogenics on imageable tumor growth. A comparison of simulated glioma growth and imaging enhancement with and without bevacizumab supports the current understanding that anti-angiogenic treatment can serve as a surrogate for steroids and the clinically driven hypothesis that anti-angiogenic treatment may not have any significant effect on the growth dynamics of the overall tumor cell populations. However, the simulations do illustrate a potentially large impact on the level of edematous extracellular fluid, and thus on what would be imageable on T2/FLAIR MR. Additionally, by evaluating virtual tumors with varying growth kinetics, we see tumors with lower proliferation rates will have the most reduction in swelling from such treatments.

Analysis of glioblastoma tumor coverage by oncolytic virus-loaded neural stem cells using MRI-based tracking and histological reconstruction

In preclinical studies, neural stem cell (NSC)-based delivery of oncolytic virus has shown great promise in the treatment of malignant glioma. Ensuring the success of this therapy will require critical evaluation of the spatial distribution of virus after NSC transplantation. In this study, the patient-derived GBM43 human glioma line was established in the brain of athymic nude mice, followed by the administration of NSCs loaded with conditionally replicating oncolytic adenovirus (NSC-CRAd-S-pk7). We determined the tumor coverage potential of oncolytic adenovirus by examining NSC distribution using magnetic resonance (MR) imaging and by three-dimensional reconstruction from ex vivo tissue specimens. We demonstrate that unmodified NSCs and NSC-CRAd-S-pk7 exhibit a similar distribution pattern with most prominent localization occurring at the tumor margins. We were further able to visualize the accumulation of these cells at tumor sites via T2-weighted MR imaging as well as the spread of viral particles using immunofluorescence. Our analyses reveal that a single administration of oncolytic virus-loaded NSCs allows for up to 31% coverage of intracranial tumors. Such results provide valuable insights into the therapeutic potential of this novel viral delivery platform.

Tumor cells in search for glutamate: an alternative explanation for increased invasiveness of IDH1 mutant gliomas

We thank the authors for a very thoughtful letter and agree that there are a number of different mechanisms through which isocitrate dehydrogenase (IDH) mutation, the downstream 2-hydroxyglutarate (2HG), can lead to a number of different state changes within tumor cells. The acidification of the tumor microenvironment was solely an interpretation of our Fig. 1. Schematic overview of the proposed effects of mutated IDH1 on cellular metabolism. a-KG is reduced to 2-HG, which when exported out of the cell could lead to acidification of tumor microenvironment. This may promote local invasive growth of tumor cells. We postulate the hypothesis that, because of depletion of cytosolic a-KG, glutamate is imported via EAAT2 and converted to a-KG by GDH (thick arrows). In this way glutamate could act as a chemotactic source that also promotes invasive tumor cell growth. Letters to the editor

In silico analysis suggests differential response to bevacizumab and radiation combination therapy in newly diagnosed glioblastoma.

Recently, two phase III studies of bevacizumab, an anti-angiogenic, for newly diagnosed glioblastoma (GBM) patients were released. While they were unable to statistically significantly demonstrate that bevacizumab in combination with other therapies increases the overall survival of GBM patients, there remains a question of potential benefits for subpopulations of patients. We use a mathematical model of GBM growth to investigate differential benefits of combining surgical resection, radiation and bevacizumab across observed tumour growth kinetics. The differential hypoxic burden after gross total resection (GTR) was assessed along with the change in radiation cell kill from bevacizumab-induced tissue re-normalization when starting therapy for tumours at different diagnostic sizes. Depending on the tumour size at the time of treatment, our model predicted that GTR would remove a variable portion of the hypoxic burden ranging from 11% to 99.99%. Further, our model predicted that the combination of bevacizumab with radiation resulted in an additional cell kill ranging from 2.6×107 to 1.1×1010 cells. By considering the outcomes given individual tumour kinetics, our results indicate that the subpopulation of patients who would receive the greatest benefit from bevacizumab and radiation combination therapy are those with large, aggressive tumours and who are not eligible for GTR.

Distinct Phenotypic Clusters of Glioblastoma Growth and Response Kinetics Predict Survival

PURPOSE Despite the intra- and intertumoral heterogeneity seen in glioblastoma multiforme (GBM), there is little definitive data on the underlying cause of the differences in patient survivals. Serial imaging assessment of tumor growth allows quantification of tumor growth kinetics (TGK) measured in terms of changes in the velocity of radial expansion seen on imaging. Because a systematic study of this entire TGK phenotype-growth before treatment and during each treatment to recurrence -has never been coordinately studied in GBMs, we sought to identify whether patients cluster into discrete groups on the basis of their TGK. PATIENTS AND METHODS From our multi-institutional database, we identified 48 patients who underwent maximally safe resection followed by radiotherapy with imaging follow-up through the time of recurrence. The patients were then clustered into two groups through a k-means algorithm taking as input only the TGK before and during treatment. RESULTS There was a significant survival difference between the clusters ( P = .003). Paradoxically, patients among the long-lived cluster had significantly larger tumors at diagnosis ( P = .027) and faster growth before treatment ( P = .003) but demonstrated a better response to adjuvant chemotherapy ( P = .048). A predictive model was built to identify which cluster patients would likely fall into on the basis of information that would be available to clinicians immediately after radiotherapy (accuracy, 90.3%). CONCLUSION Dichotomizing the heterogeneity of GBMs into two populations-one faster growing yet more responsive with increased survival and one slower growing yet less responsive with shorter survival-suggests that many patients who receive standard-of-care treatments may get better benefit from select alternative treatments.