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 patient's 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.

[1]  D. Silbergeld,et al.  Gene therapy enhances chemotherapy tolerance and efficacy in glioblastoma patients. , 2014, The Journal of clinical investigation.

[2]  Philip Hahnfeldt,et al.  Simple ODE models of tumor growth and anti-angiogenic or radiation treatment , 2001 .

[3]  Alan C. Evans,et al.  BrainWeb: Online Interface to a 3D MRI Simulated Brain Database , 1997 .

[4]  D. Gillespie,et al.  Hypoxia‐regulated protein expression, patient characteristics, and preoperative imaging as predictors of survival in adults with glioblastoma multiforme , 2008, Cancer.

[5]  R. Fisher THE WAVE OF ADVANCE OF ADVANTAGEOUS GENES , 1937 .

[6]  Alexander R A Anderson,et al.  Quantifying the Role of Angiogenesis in Malignant Progression of Gliomas: in Silico Modeling Integrates Imaging and Histology Nih Public Access Author Manuscript Introduction , 2011 .

[7]  J. Webster,et al.  Changes in tumor oxygen tension during radiation therapy. , 1969, Acta radiologica: therapy, physics, biology.

[8]  K Hendrickson,et al.  Predicting the efficacy of radiotherapy in individual glioblastoma patients in vivo: a mathematical modeling approach , 2010, Physics in medicine and biology.

[9]  E. Corey,et al.  Potentiation of angiostatic steroids by a synthetic inhibitor of arylsulfatase. , 1988, Laboratory investigation; a journal of technical methods and pathology.

[10]  Robert Jeraj,et al.  An imaging-based tumour growth and treatment response model: investigating the effect of tumour oxygenation on radiation therapy response , 2008, Physics in medicine and biology.

[11]  Laurent Capelle,et al.  Continuous growth of mean tumor diameter in a subset of grade II gliomas , 2003, Annals of neurology.

[12]  Peter Canoll,et al.  MRI-localized biopsies reveal subtype-specific differences in molecular and cellular composition at the margins of glioblastoma , 2014, Proceedings of the National Academy of Sciences.

[13]  Deric M. Park,et al.  Improving the radiosensitivity of radioresistant and hypoxic glioblastoma. , 2010, Future oncology.

[14]  F. O’Sullivan,et al.  Regional Hypoxia in Glioblastoma Multiforme Quantified with [18F]Fluoromisonidazole Positron Emission Tomography before Radiotherapy: Correlation with Time to Progression and Survival , 2008, Clinical Cancer Research.

[15]  Georgios Stamatakos,et al.  A Patient-specific in vivo Tumor and Normal Tissue Model for Prediction of the Response to Radiotherapy , 2007, Methods of Information in Medicine.

[16]  M. Bernstein,et al.  Advances in radiotherapy of brain tumors: radiobiology versus reality , 1998, Journal of Clinical Neuroscience.

[17]  C. Meltzer,et al.  Brain tumor volume measurement: comparison of manual and semiautomated methods. , 1999, Radiology.

[18]  J. Barnholtz-Sloan,et al.  CBTRUS statistical report: primary brain and central nervous system tumors diagnosed in the United States in 2007-2011. , 2012, Neuro-oncology.

[19]  S. Gabriel,et al.  Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. , 2010, Cancer cell.

[20]  Benjamin M Ellingson,et al.  Spatially quantifying microscopic tumor invasion and proliferation using a voxel‐wise solution to a glioma growth model and serial diffusion MRI , 2011, Magnetic resonance in medicine.

[21]  R. Guillevin,et al.  Simulation of anisotropic growth of low‐grade gliomas using diffusion tensor imaging , 2005, Magnetic resonance in medicine.

[22]  J. Murray,et al.  A quantitative model for differential motility of gliomas in grey and white matter , 2000, Cell proliferation.

[23]  V. Amberger-Murphy Hypoxia helps glioma to fight therapy. , 2009, Current cancer drug targets.

[24]  Soonmee Cha,et al.  Imaging Glioblastoma Multiforme , 2003, Cancer journal.

[25]  M. Graham,et al.  A modeling approach for quantifying tumor hypoxia with [F-18]fluoromisonidazole PET time-activity data. , 1995, Medical physics.

[26]  B. Liang Effects of hypoxia on drug resistance phenotype and genotype in human glioma cell lines , 1996, Journal of Neuro-Oncology.

[27]  Kristin R. Swanson,et al.  Response classification based on a minimal model of glioblastoma growth is prognostic for clinical outcomes and distinguishes progression from pseudoprogression. , 2013, Cancer research.

[28]  R K Sachs,et al.  The mechanistic basis of the linear-quadratic formalism. , 1998, Medical physics.

[29]  Minsun Kim,et al.  A hierarchical evolutionary algorithm for multiobjective optimization in IMRT. , 2010, Medical physics.

[30]  M. Westphal,et al.  Glioma cell adhesion and migration on human brain sections. , 1998, Anticancer research.

[31]  B. Scheithauer,et al.  The 2007 WHO classification of tumours of the central nervous system , 2007, Acta Neuropathologica.

[32]  V P Antipas,et al.  A four-dimensional computer simulation model of the in vivo response to radiotherapy of glioblastoma multiforme: studies on the effect of clonogenic cell density. , 2006, The British journal of radiology.

[33]  Loredana G. Marcu,et al.  In Silico Modelling of Tumour Margin Diffusion and Infiltration: Review of Current Status , 2012, Comput. Math. Methods Medicine.

[34]  Peter T. Nelson,et al.  Hypoxia Is Important in the Biology and Aggression of Human Glial Brain Tumors , 2004, Clinical Cancer Research.

[35]  D. Wilkins,et al.  Fitting the linear–quadratic model to detailed data sets for different dose ranges , 2006, Physics in medicine and biology.

[36]  Gargi Chakraborty,et al.  Quantitative metrics of net proliferation and invasion link biological aggressiveness assessed by MRI with hypoxia assessed by FMISO-PET in newly diagnosed glioblastomas. , 2009, Cancer research.

[37]  D. Olsen,et al.  Optimization of tumour control probability in hypoxic tumours by radiation dose redistribution: a modelling study , 2007, Physics in medicine and biology.

[38]  Peter Canoll,et al.  Magnetic Resonance Imaging Characteristics of Glioblastoma Multiforme: Implications for Understanding Glioma Ontogeny , 2010, Neurosurgery.

[39]  Mark A. Brown,et al.  MRI: Basic Principles and Applications , 1995 .

[40]  V. P. Collins,et al.  Intratumor heterogeneity in human glioblastoma reflects cancer evolutionary dynamics , 2013, Proceedings of the National Academy of Sciences.

[41]  D. Silbergeld,et al.  Assessment of brain tumor cell motility in vivo and in vitro. , 1995, Journal of neurosurgery.

[42]  R. Mirimanoff,et al.  Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. , 2005, The New England journal of medicine.

[43]  J. Knisely,et al.  Importance of hypoxia in the biology and treatment of brain tumors. , 2002, Neuroimaging clinics of North America.

[44]  Kristin R. Swanson,et al.  The Evolution of Mathematical Modeling of Glioma Proliferation and Invasion , 2007, Journal of neuropathology and experimental neurology.

[45]  Paul E Kinahan,et al.  Applying a patient-specific bio-mathematical model of glioma growth to develop virtual [18F]-FMISO-PET images. , 2011, Mathematical medicine and biology : a journal of the IMA.

[46]  Kristin R. Swanson,et al.  Complementary but Distinct Roles for MRI and 18F-Fluoromisonidazole PET in the Assessment of Human Glioblastomas , 2008, Journal of Nuclear Medicine.

[47]  G. Holmes Epilepsy surgery in children: When, why, and how , 2002, Neurology.

[48]  F Barkhof,et al.  Interobserver variability in the radiological assessment of response to chemotherapy in glioma , 2003, Neurology.

[49]  T W Griffin,et al.  Imaging of hypoxia in human tumors with [F-18]fluoromisonidazole. , 1992, International journal of radiation oncology, biology, physics.

[50]  David J. Carlson,et al.  Effects of oxygen on intrinsic radiation sensitivity: A test of the relationship between aerobic and hypoxic linear-quadratic (LQ) model parametersa). , 2006, Medical physics.

[51]  J. Chapman,et al.  The detection and measurement of hypoxic cells in solid tumors , 1984, Cancer.

[52]  R K Sachs,et al.  A convenient extension of the linear-quadratic model to include redistribution and reoxygenation. , 1995, International journal of radiation oncology, biology, physics.

[53]  N Uzunoglu,et al.  Critical Parameters Determining Standard Radiotherapy Treatment Outcome for Glioblastoma Multiforme: A Computer Simulation , 2008, The open biomedical engineering journal.

[54]  J. Barnholtz-Sloan,et al.  CBTRUS statistical report: Primary brain and central nervous system tumors diagnosed in the United States in 2006-2010. , 2013, Neuro-oncology.

[55]  Karl J. Friston,et al.  Statistical parametric mapping (SPM) , 2008, Scholarpedia.

[56]  R. Jensen,et al.  Brain tumor hypoxia: tumorigenesis, angiogenesis, imaging, pseudoprogression, and as a therapeutic target , 2009, Journal of Neuro-Oncology.

[57]  R. Dale,et al.  Mathematical models of tumour and normal tissue response. , 1999, Acta oncologica.

[58]  Santosh Kesari,et al.  Malignant gliomas in adults. , 2008, The New England journal of medicine.

[59]  J. Hendry,et al.  Radiobiology for the Radiologist , 1979, British Journal of Cancer.

[60]  H. Withers,et al.  Biologically effective dose distribution based on the linear quadratic model and its clinical relevance. , 1995, International journal of radiation oncology, biology, physics.

[61]  Vladimir A Semenenko,et al.  Effects of oxygen on intrinsic radiation sensitivity: A test of the relationship between aerobic and hypoxic linear-quadratic (LQ) model parameters. , 2006, Medical physics.

[62]  R K Sachs,et al.  The linear-quadratic model and most other common radiobiological models result in similar predictions of time-dose relationships. , 1998, Radiation research.

[63]  J. Murray,et al.  Virtual brain tumours (gliomas) enhance the reality of medical imaging and highlight inadequacies of current therapy , 2002, British Journal of Cancer.

[64]  James B. Smathers,et al.  Biologically effective dose distribution based on the linear quadratic model and its clinical relevance , 1995 .