An imaging-based tumour growth and treatment response model: investigating the effect of tumour oxygenation on radiation therapy response

A multiscale tumour simulation model employing cell-line-specific biological parameters and functional information derived from pre-therapy PET/CT imaging data was developed to investigate effects of different oxygenation levels on the response to radiation therapy. For each tumour voxel, stochastic simulations were performed to model cellular growth and therapeutic response. Model parameters were fitted to published preclinical experiments of head and neck squamous cell carcinoma (HNSCC). Using the obtained parameters, the model was applied to a human HNSCC case to investigate effects of different uniform and non-uniform oxygenation levels and results were compared for treatment efficacy. Simulations of the preclinical studies showed excellent agreement with published data and underlined the model's ability to quantitatively reproduce tumour behaviour within experimental uncertainties. When using a simplified transformation to derive non-uniform oxygenation levels from molecular imaging data, simulations of the clinical case showed heterogeneous tumour response and variability in radioresistance with decreasing oxygen levels. Once clinically validated, this model could be used to transform patient-specific data into voxel-based biological objectives for treatment planning and to investigate biologically optimized dose prescriptions.

[1]  D. Hanahan,et al.  The Hallmarks of Cancer , 2000, Cell.

[2]  A. Hanlon,et al.  Hypoxic regions exist in human prostate carcinoma. , 1999, Urology.

[3]  T K Lewellen,et al.  Quantifying regional hypoxia in human tumors with positron emission tomography of [18F]fluoromisonidazole: a pretherapy study of 37 patients. , 1996, International journal of radiation oncology, biology, physics.

[4]  R. Wasserman,et al.  A patient-specific in vivo tumor model. , 1996, Mathematical biosciences.

[5]  H Treuer,et al.  Reoxygenation of hypoxic cells by tumor shrinkage during irradiation. A computer simulation. , 1995, Strahlentherapie und Onkologie : Organ der Deutschen Rontgengesellschaft ... [et al].

[6]  A. van Waarde,et al.  Selectivity of 18F-FLT and 18F-FDG for differentiating tumor from inflammation in a rodent model. , 2004, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[7]  X Allen Li,et al.  BGRT: biologically guided radiation therapy-the future is fast approaching! , 2007, Medical physics.

[8]  F. Lohr,et al.  Measurement of the proliferative activity of three different sublines of the Dunning rat prostate tumor R3327. , 1993, Strahlentherapie und Onkologie : Organ der Deutschen Rontgengesellschaft ... [et al].

[9]  Olivier Clatz,et al.  Biocomputing: numerical simulation of glioblastoma growth using diffusion tensor imaging , 2008, Physics in medicine and biology.

[10]  B. Baguley Tumor stem cell niches: a new functional framework for the action of anticancer drugs. , 2006, Recent patents on anti-cancer drug discovery.

[11]  D. Mankoff,et al.  Positron-emission tomographic imaging of cancer: glucose metabolism and beyond. , 2001, Seminars in radiation oncology.

[12]  Michael J. Welch,et al.  In vivo assessment of tumor hypoxia in lung cancer with 60Cu-ATSM , 2003, European Journal of Nuclear Medicine and Molecular Imaging.

[13]  M. Dewhirst,et al.  Tumor hypoxia adversely affects the prognosis of carcinoma of the head and neck. , 1997, International journal of radiation oncology, biology, physics.

[14]  M Magnano,et al.  Effects of anatomical constraints on tumor growth. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[15]  Konstantina S. Nikita,et al.  A computer simulation of in vivo tumour growth and response to radiotherapy: New algorithms and parametric results , 2006, Comput. Biol. Medicine.

[16]  T L Phillips,et al.  Oxygen in human tumors: correlations between methods of measurement and response to therapy. Summary of a workshop held November 19-20, 1992, at the National Cancer Institute, Bethesda, Maryland. , 1993, Radiation research.

[17]  P. Harari,et al.  Modulation of radiation response after epidermal growth factor receptor blockade in squamous cell carcinomas: inhibition of damage repair, cell cycle kinetics, and tumor angiogenesis. , 2000, Clinical cancer research : an official journal of the American Association for Cancer Research.

[18]  J. Moore,et al.  Radiation response and cure rate of human colon adenocarcinoma spheroids of different size: the significance of hypoxia on tumor control modelling. , 2001, International journal of radiation oncology, biology, physics.

[19]  Brian W. Pogue,et al.  Estimation of Oxygen Distribution in RIF-1 Tumors by Diffusion Model-Based Interpretation of Pimonidazole Hypoxia and Eppendorf Measurements , 2001, Radiation research.

[20]  D. Brizel,et al.  Prognostic value of tumor oxygenation in 397 head and neck tumors after primary radiation therapy. An international multi-center study. , 2005, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[21]  Jelena Pjesivac-Grbovic,et al.  A multiscale model for avascular tumor growth. , 2005, Biophysical journal.

[22]  R. Hill Tumor progression: Potential role of unstable genomic changes , 1990, Cancer and Metastasis Reviews.

[23]  T. Alper,et al.  Cellular radiobiology. , 1960, Annual review of nuclear science.

[24]  P Vaupel,et al.  Association between tumor hypoxia and malignant progression in advanced cancer of the uterine cervix. , 1996, Cancer research.

[25]  Y. Fujibayashi,et al.  Delineation of hypoxia in canine myocardium using PET and copper(II)-diacetyl-bis(N(4)-methylthiosemicarbazone). , 2002, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[26]  S. Schnell,et al.  A multiscale mathematical model of cancer, and its use in analyzing irradiation therapies , 2006, Theoretical Biology and Medical Modelling.

[27]  Nikolaos K. Uzunoglu,et al.  A spatiotemporal, patient individualized simulation model of solid tumor response to chemotherapy in vivo: the paradigm of glioblastoma multiforme treated by temozolomide , 2006, IEEE Transactions on Biomedical Engineering.

[28]  Christian P Karger,et al.  Single-cell-based computer simulation of the oxygen-dependent tumour response to irradiation , 2007, Physics in medicine and biology.

[29]  Y. Fujibayashi,et al.  Tumor uptake of copper-diacetyl-bis(N(4)-methylthiosemicarbazone): effect of changes in tissue oxygenation. , 2001, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[30]  M. Weil,et al.  An equation for the oxygen hemoglobin dissociation curve. , 1973, Journal of applied physiology.

[31]  Peter Peschke,et al.  Modeling and Computer Simulations of Tumor Growth and Tumor Response to Radiotherapy , 2004, Radiation research.

[32]  P. Monaghan,et al.  Ten human carcinoma cell lines derived from squamous carcinomas of the head and neck. , 1981, British journal of cancer.

[33]  M Kardar,et al.  Dynamics of tumor growth and combination of anti-angiogenic and cytotoxic therapies , 2007, Physics in medicine and biology.

[34]  P. Vaupel,et al.  Tumor hypoxia: definitions and current clinical, biologic, and molecular aspects. , 2001, Journal of the National Cancer Institute.

[35]  B. Wouters,et al.  Cells at intermediate oxygen levels can be more important than the "hypoxic fraction" in determining tumor response to fractionated radiotherapy. , 1997, Radiation research.

[36]  M. Martel,et al.  Fractionated regimens for stereotactic radiotherapy of recurrent tumors in the brain. , 1991, International journal of radiation oncology, biology, physics.

[37]  J. F. De Los Santos,et al.  Anemia correction in malignancy management: threat or opportunity? , 2007, Gynecologic oncology.

[38]  TIKVAH ALPER,et al.  Role of Oxygen in Modifying the Radiosensitivity of E. Coli B. , 1956, Nature.

[39]  J. Nilsson Accurate description of heterogeneous tumors for biologically optimized radiation therapy , 2005 .

[40]  Otto Muzik,et al.  Imaging proliferation in vivo with [F-18]FLT and positron emission tomography , 1998, Nature Medicine.

[41]  S Mutic,et al.  A novel approach to overcome hypoxic tumor resistance: Cu-ATSM-guided intensity-modulated radiation therapy. , 2001, International journal of radiation oncology, biology, physics.

[42]  G. Dalrymple,et al.  Initial studies with a line of radioresistant rat tumor cells. , 1974, Radiation research.

[43]  O. S. Nielsen,et al.  Hypoxia in human soft tissue sarcomas: Adverse impact on survival and no association with p53 mutations , 2001, British Journal of Cancer.

[44]  P Vaupel,et al.  Oxygenation of human tumors: evaluation of tissue oxygen distribution in breast cancers by computerized O2 tension measurements. , 1991, Cancer research.

[45]  E. T. Gawlinski,et al.  A Cellular Automaton Model of Early Tumor Growth and Invasion: The Effects of Native Tissue Vascularity and Increased Anaerobic Tumor Metabolism , 2001 .

[46]  Mark A Mintun,et al.  Assessing tumor hypoxia in cervical cancer by positron emission tomography with 60Cu-ATSM: relationship to therapeutic response-a preliminary report. , 2003, International journal of radiation oncology, biology, physics.

[47]  Konstantina S. Nikita,et al.  Modeling tumor growth and irradiation response in vitro-a combination of high-performance computing and Web-based technologies including VRML visualization , 2001, IEEE Transactions on Information Technology in Biomedicine.

[48]  P. Harari,et al.  Modulation of radiation response and tumor-induced angiogenesis after epidermal growth factor receptor inhibition by ZD1839 (Iressa). , 2002, Cancer research.

[49]  E. Hall,et al.  Radiobiology for the radiologist , 1973 .

[50]  C. Ling,et al.  Radiation-induced apoptosis: effects of cell age and dose fractionation. , 1995, Cancer research.

[51]  H. Frieboes,et al.  Predictive oncology: A review of multidisciplinary, multiscale in silico modeling linking phenotype, morphology and growth , 2007, NeuroImage.

[52]  M. Dorie,et al.  Tumor oxygenation and reoxygenation during radiation therapy: their importance in predicting tumor response. , 1986, International journal of radiation oncology, biology, physics.

[53]  N. Johnson,et al.  DNA damage-inducible genes as biomarkers for exposures to environmental agents. , 1997, Environmental health perspectives.

[54]  J F Fowler,et al.  Biological Factors Influencing Optimum Fractionation in Radiation Therapy , 2001, Acta oncologica.

[55]  T Ido,et al.  Intratumoral distribution of fluorine-18-fluorodeoxyglucose in vivo: high accumulation in macrophages and granulation tissues studied by microautoradiography. , 1992, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[56]  Morand Piert,et al.  Hypoxia-specific tumor imaging with 18F-fluoroazomycin arabinoside. , 2005, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[57]  G. Hanks,et al.  Single-photon emission computed tomography and positron-emission tomography assays for tissue oxygenation. , 2001, Seminars in radiation oncology.

[58]  Y Yonekura,et al.  Copper-62-ATSM: a new hypoxia imaging agent with high membrane permeability and low redox potential. , 1997, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[59]  G. Steel,et al.  Basic Clinical Radiobiology , 1997 .

[60]  H D Thames,et al.  Does variation in the in vitro cellular radiosensitivity explain the shallow clinical dose-control curve for malignant melanoma? , 1990, International journal of radiation biology.

[61]  J H Hendry,et al.  A realistic closed-form radiobiological model of clinical tumor-control data incorporating intertumor heterogeneity. , 1998, International journal of radiation oncology, biology, physics.

[62]  H. Withers,et al.  Four R's of radiotherapy , 1975 .

[63]  R. Pötter,et al.  Intratumoral pO2-measurements as predictive assay in the treatment of carcinoma of the uterine cervix. , 1999, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[64]  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.

[65]  J. Overgaard,et al.  Pretreatment oxygenation predicts radiation response in advanced squamous cell carcinoma of the head and neck. , 1996, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[66]  G Powathil,et al.  Mathematical modeling of brain tumors: effects of radiotherapy and chemotherapy , 2007, Physics in medicine and biology.

[67]  Georgios S Stamatakos,et al.  A four-dimensional simulation model of tumour response to radiotherapy in vivo: parametric validation considering radiosensitivity, genetic profile and fractionation. , 2004, Journal of theoretical biology.

[68]  M J Welch,et al.  Evaluation of 64Cu-ATSM in vitro and in vivo in a hypoxic tumor model. , 1999, Journal of nuclear medicine : official publication, Society of Nuclear Medicine.

[69]  E. Rofstad,et al.  Radiobiological and immunohistochemical assessment of hypoxia in human melanoma xenografts: acute and chronic hypoxia in individual tumours. , 1999, International journal of radiation biology.

[70]  J Denekamp,et al.  Theoretical simulation of oxygen tension measurement in tissues using a microelectrode: I. The response function of the electrode. , 2001, Physiological measurement.