1
|
Kumar P, Lacroix M, Dupré P, Arslan J, Fenou L, Orsetti B, Le Cam L, Racoceanu D, Radulescu O. Deciphering oxygen distribution and hypoxia profiles in the tumor microenvironment: a data-driven mechanistic modeling approach. Phys Med Biol 2024; 69:125023. [PMID: 38815610 DOI: 10.1088/1361-6560/ad524a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2024] [Accepted: 05/30/2024] [Indexed: 06/01/2024]
Abstract
Objective. The distribution of hypoxia within tissues plays a critical role in tumor diagnosis and prognosis. Recognizing the significance of tumor oxygenation and hypoxia gradients, we introduce mathematical frameworks grounded in mechanistic modeling approaches for their quantitative assessment within a tumor microenvironment. By utilizing known blood vasculature, we aim to predict hypoxia levels across different tumor types.Approach. Our approach offers a computational method to measure and predict hypoxia using known blood vasculature. By formulating a reaction-diffusion model for oxygen distribution, we derive the corresponding hypoxia profile.Main results. The framework successfully replicates observed inter- and intra-tumor heterogeneity in experimentally obtained hypoxia profiles across various tumor types (breast, ovarian, pancreatic). Additionally, we propose a data-driven method to deduce partial differential equation models with spatially dependent parameters, which allows us to comprehend the variability of hypoxia profiles within tissues. The versatility of our framework lies in capturing diverse and dynamic behaviors of tumor oxygenation, as well as categorizing states of vascularization based on the dynamics of oxygen molecules, as identified by the model parameters.Significance. The proposed data-informed mechanistic method quantitatively assesses hypoxia in the tumor microenvironment by integrating diverse histopathological data and making predictions across different types of data. The framework provides valuable insights from both modeling and biological perspectives, advancing our comprehension of spatio-temporal dynamics of tumor oxygenation.
Collapse
Affiliation(s)
- P Kumar
- Laboratory of Pathogens and Host Immunity, University of Montpellier, CNRS, INSERM, Montpellier, France
- Sorbonne Université, CNRS, INSERM, AP-HP, Inria, Paris Brain Institute (ICM), Paris, France
| | - M Lacroix
- Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, University of Montpellier, Institut régional du Cancer de Montpellier (ICM), Montpellier, France
- Equipe labélisée Ligue Contre le Cancer, Paris, France
| | - P Dupré
- Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, University of Montpellier, Institut régional du Cancer de Montpellier (ICM), Montpellier, France
- Equipe labélisée Ligue Contre le Cancer, Paris, France
| | - J Arslan
- Sorbonne Université, CNRS, INSERM, AP-HP, Inria, Paris Brain Institute (ICM), Paris, France
- Centre for Eye Research Australia, University of Melbourne, Royal Victorian Eye & Ear Hospital, East Melbourne, Australia
| | - L Fenou
- Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, University of Montpellier, Institut régional du Cancer de Montpellier (ICM), Montpellier, France
| | - B Orsetti
- Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, University of Montpellier, Institut régional du Cancer de Montpellier (ICM), Montpellier, France
| | - L Le Cam
- Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, University of Montpellier, Institut régional du Cancer de Montpellier (ICM), Montpellier, France
- Equipe labélisée Ligue Contre le Cancer, Paris, France
| | - D Racoceanu
- Sorbonne Université, CNRS, INSERM, AP-HP, Inria, Paris Brain Institute (ICM), Paris, France
| | - O Radulescu
- Laboratory of Pathogens and Host Immunity, University of Montpellier, CNRS, INSERM, Montpellier, France
| |
Collapse
|
2
|
Taylor E. A simple mathematical model of cyclic hypoxia and its impact on hypofractionated radiotherapy. Med Phys 2023; 50:1893-1904. [PMID: 36594511 DOI: 10.1002/mp.16200] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Revised: 08/29/2022] [Accepted: 12/16/2022] [Indexed: 01/04/2023] Open
Abstract
PURPOSE There is evidence that the population of cells that experience fluctuating oxygen levels ("acute," or, "cyclic" hypoxia) are more radioresistant than chronically hypoxic ones and hence, this population may determine radiotherapy (RT) response, in particular for hypofractionated RT, where reoxygenation may not be as prominent. A considerable effort has been devoted to examining the impact of hypoxia on hypofractionated RT; however, much less attention has been paid to cyclic hypoxia specifically and the role its kinetics may play in determining the efficacy of these treatments. Here, a simple mathematical model of cyclic hypoxia and fractionation effects was worked out to quantify this. METHODS Cancer clonogen survival fraction was estimated using the linear quadratic model, modified to account for oxygen enhancement effects. An analytic approximation for oxygen transport away from a random network of capillaries with fluctuating oxygen levels was used to model inter-fraction tissue oxygen kinetics. The resulting survival fraction formula was used to derive an expression for the iso-survival biologically effective dose (BED), BEDiso-SF . These were computed for some common extra-cranial hypofractionated RT regimens. RESULTS Using relevant literature parameter values, inter-fraction fluctuations in oxygenation were found to result in an added 1-2 logs of clonogen survival fraction in going from five fractions to one for the same nominal BED (i.e., excluding the effects of oxygen levels on radiosensitivity). BEDiso-SF 's for most ultra-hypofractionated (five or fewer fractions) regimens in a given tumor site are similar in magnitude, suggesting iso-efficacy for common fractionation schedules. CONCLUSIONS Although significant, the loss of cell-killing with increasing hypofractionation is not nearly as large as previous estimates based on the assumption of complete reoxygenation between fractions. Most ultra-hypofractionated regimens currently in place offer sufficiently high doses to counter this loss of cell killing, although care should be taken in implementing single-fraction regimens.
Collapse
Affiliation(s)
- Edward Taylor
- Princess Margaret Cancer Centre, University Health Network, Toronto, Canada
- Department of Radiation Oncology, University of Toronto, Toronto, Canada
| |
Collapse
|
3
|
Abnormal morphology biases hematocrit distribution in tumor vasculature and contributes to heterogeneity in tissue oxygenation. Proc Natl Acad Sci U S A 2020; 117:27811-27819. [PMID: 33109723 PMCID: PMC7668105 DOI: 10.1073/pnas.2007770117] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023] Open
Abstract
Oxygen heterogeneity in solid tumors is recognized as a limiting factor for therapeutic efficacy. This heterogeneity arises from the abnormal tumor vascular structure. We investigate the role that anomalies in red blood cell transport plays in establishing oxygen heterogeneity in tumor tissue. We introduce a metric to characterize tumor vasculature (mean vessel length-to-diameter ratio, λ) and demonstrate how it predicts tissue-oxygen heterogeneity. We also report an increase in λ following treatment with the antiangiogenic agent DC101. Together, we propose λ as an effective way of monitoring the action of antiangiogenic agents and a proxy measure of oxygen heterogeneity in tumor tissue. Unraveling the causal relationship between tumor vascular structure and tissue oxygenation will pave the way for new personalized therapeutic approaches. Oxygen heterogeneity in solid tumors is recognized as a limiting factor for therapeutic efficacy. This heterogeneity arises from the abnormal vascular structure of the tumor, but the precise mechanisms linking abnormal structure and compromised oxygen transport are only partially understood. In this paper, we investigate the role that red blood cell (RBC) transport plays in establishing oxygen heterogeneity in tumor tissue. We focus on heterogeneity driven by network effects, which are challenging to observe experimentally due to the reduced fields of view typically considered. Motivated by our findings of abnormal vascular patterns linked to deviations from current RBC transport theory, we calculated average vessel lengths L¯ and diameters d¯ from tumor allografts of three cancer cell lines and observed a substantial reduction in the ratio λ=L¯/d¯ compared to physiological conditions. Mathematical modeling reveals that small values of the ratio λ (i.e., λ<6) can bias hematocrit distribution in tumor vascular networks and drive heterogeneous oxygenation of tumor tissue. Finally, we show an increase in the value of λ in tumor vascular networks following treatment with the antiangiogenic cancer agent DC101. Based on our findings, we propose λ as an effective way of monitoring the efficacy of antiangiogenic agents and as a proxy measure of perfusion and oxygenation in tumor tissue undergoing antiangiogenic treatment.
Collapse
|
4
|
Kuperman VY, Lubich LM. Effect of reoxygenation on hypofractionated radiotherapy of prostate cancer. Med Phys 2020; 47:5383-5391. [DOI: 10.1002/mp.14343] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2020] [Revised: 06/05/2020] [Accepted: 06/09/2020] [Indexed: 11/07/2022] Open
|
5
|
Rühle A, Grosu AL, Wiedenmann N, Mix M, Stoian R, Niedermann G, Baltas D, Werner M, Weber WA, Kayser G, Nicolay NH. Hypoxia dynamics on FMISO-PET in combination with PD-1/PD-L1 expression has an impact on the clinical outcome of patients with Head-and-neck Squamous Cell Carcinoma undergoing Chemoradiation. Am J Cancer Res 2020; 10:9395-9406. [PMID: 32802199 PMCID: PMC7415814 DOI: 10.7150/thno.48392] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2020] [Accepted: 06/29/2020] [Indexed: 12/24/2022] Open
Abstract
Tumor-associated hypoxia influences the radiation response of head-and-neck cancer (HNSCC) patients, and a lack of early hypoxia resolution during treatment considerably deteriorates outcomes. As the detrimental effects of hypoxia are partly related to the induction of an immunosuppressive microenvironment, we investigated the interaction between tumor hypoxia dynamics and the PD-1/PD-L1 axis in HNSCC patients undergoing chemoradiation and its relevance for patient outcomes in a prospective trial. Methods: 49 patients treated with definitive chemoradiation for locally advanced HNSCC were enrolled in this trial and received longitudinal hypoxia PET imaging using fluorine-18 misonidazole ([18F]FMISO) at weeks 0, 2 and 5 during treatment. Pre-therapeutic tumor biopsies were immunohistochemically analyzed regarding the PD-1/PD-L1 expression both on immune cells and on tumor cells, and potential correlations between the PD-1/PD-L1 axis and tumor hypoxia dynamics during chemoradiation were assessed using Spearman's rank correlations. Hypoxia dynamics during treatment were quantified by subtracting the standardized uptake value (SUV) index at baseline from the SUV values at weeks 2 or 5, whereby SUV index was defined as ratio of maximum tumor [18F]FMISO SUV to mean SUV in the contralateral sternocleidomastoid muscle (i.e. tumor-to-muscle ratio). The impact of the PD-1/PD-L1 expression alone and in combination with persistent tumor hypoxia on locoregional control (LRC), progression-free survival (PFS) and overall survival (OS) was examined using log-rank tests and Cox proportional hazards models. Results: Neither PD-L1 nor PD-1 expression levels on tumor-infiltrating immune cells influenced LRC (HR = 0.734; p = 0.480 for PD-L1, HR = 0.991; p = 0.989 for PD-1), PFS (HR = 0.813; p = 0.597 for PD-L1, HR = 0.796; p = 0.713 for PD-1) or OS (HR = 0.698; p = 0.405 for PD-L1, HR = 0.315; p = 0.265 for PD-1). However, patients with no hypoxia resolution between weeks 0 and 2 and PD-L1 expression on tumor cells, quantified by a tumor proportional score (TPS) of at least 1%, showed significantly worse LRC (HR = 3.374, p = 0.022) and a trend towards reduced PFS (HR = 2.752, p = 0.052). In the multivariate Cox regression analysis, the combination of absent tumor hypoxia resolution and high tumoral PD-L1 expression remained a significant prognosticator for impaired LRC (HR = 3.374, p = 0.022). On the other side, tumoral PD-L1 expression did not compromise the outcomes of patients whose tumor-associated hypoxia declined between week 0 and 2 during chemoradiation (LRC: HR = 1.186, p = 0.772, PFS: HR = 0.846, p = 0.766). Conclusion: In this exploratory analysis, we showed for the first time that patients with both persistent tumor-associated hypoxia during treatment and PD-L1 expression on tumor cells exhibited a worse outcome, while the tumor cells' PD-L1 expression did not influence the outcomes of patients with early tumor hypoxia resolution. While the results have to be validated in an independent cohort, these findings form a foundation to investigate the combination of hypoxic modification and immune checkpoint inhibitors for the unfavorable subgroup, moving forward towards personalized radiation oncology treatment.
Collapse
|
6
|
Manem VS, Lambie M, Smith I, Smirnov P, Kofia V, Freeman M, Koritzinsky M, Abazeed ME, Haibe-Kains B, Bratman SV. Modeling Cellular Response in Large-Scale Radiogenomic Databases to Advance Precision Radiotherapy. Cancer Res 2019; 79:6227-6237. [PMID: 31558563 DOI: 10.1158/0008-5472.can-19-0179] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2019] [Revised: 07/03/2019] [Accepted: 09/17/2019] [Indexed: 12/22/2022]
Abstract
Radiotherapy is integral to the care of a majority of patients with cancer. Despite differences in tumor responses to radiation (radioresponse), dose prescriptions are not currently tailored to individual patients. Recent large-scale cancer cell line databases hold the promise of unravelling the complex molecular arrangements underlying cellular response to radiation, which is critical for novel predictive biomarker discovery. Here, we present RadioGx, a computational platform for integrative analyses of radioresponse using radiogenomic databases. We fit the dose-response data within RadioGx to the linear-quadratic model. The imputed survival across a range of dose levels (AUC) was a robust radioresponse indicator that correlated with biological processes known to underpin the cellular response to radiation. Using AUC as a metric for further investigations, we found that radiation sensitivity was significantly associated with disruptive mutations in genes related to nonhomologous end joining. Next, by simulating the effects of different oxygen levels, we identified putative genes that may influence radioresponse specifically under hypoxic conditions. Furthermore, using transcriptomic data, we found evidence for tissue-specific determinants of radioresponse, suggesting that tumor type could influence the validity of putative predictive biomarkers of radioresponse. Finally, integrating radioresponse with drug response data, we found that drug classes impacting the cytoskeleton, DNA replication, and mitosis display similar therapeutic effects to ionizing radiation on cancer cell lines. In summary, RadioGx provides a unique computational toolbox for hypothesis generation to advance preclinical research for radiation oncology and precision medicine. SIGNIFICANCE: The RadioGx computational platform enables integrative analyses of cellular response to radiation with drug responses and genome-wide molecular data. GRAPHICAL ABSTRACT: http://cancerres.aacrjournals.org/content/canres/79/24/6227/F1.large.jpg.See related commentary by Spratt and Speers, p. 6076.
Collapse
Affiliation(s)
- Venkata Sk Manem
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada.,Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
| | - Meghan Lambie
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada.,Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada
| | - Ian Smith
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada.,Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada.,Vector Institute, Toronto, Ontario, Canada
| | - Petr Smirnov
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada.,Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada.,Vector Institute, Toronto, Ontario, Canada
| | - Victor Kofia
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
| | - Mark Freeman
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada
| | - Marianne Koritzinsky
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada.,Department of Radiation Oncology, University of Toronto, Toronto, Ontario, Canada.,Institute of Medical Sciences, University of Toronto, Toronto, Ontario, Canada
| | - Mohamed E Abazeed
- Department of Translational Hematology Oncology Research, Cleveland, Ohio.,Department of Radiation Oncology, Cleveland Clinic, Cleveland, Ohio
| | - Benjamin Haibe-Kains
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada. .,Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada.,Vector Institute, Toronto, Ontario, Canada.,Department of Computer Science, University of Toronto, Toronto, Ontario, Canada.,Ontario Institute of Cancer Research, Toronto, Ontario, Canada
| | - Scott V Bratman
- Princess Margaret Cancer Centre, University Health Network, Toronto, Ontario, Canada. .,Department of Medical Biophysics, University of Toronto, Toronto, Ontario, Canada.,Department of Radiation Oncology, University of Toronto, Toronto, Ontario, Canada
| |
Collapse
|
7
|
Forster JC, Marcu LG, Bezak E. Approaches to combat hypoxia in cancer therapy and the potential for in silico models in their evaluation. Phys Med 2019; 64:145-156. [PMID: 31515013 DOI: 10.1016/j.ejmp.2019.07.006] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/28/2019] [Revised: 06/17/2019] [Accepted: 07/09/2019] [Indexed: 02/07/2023] Open
Abstract
AIM The negative impact of tumour hypoxia on cancer treatment outcome has been long-known, yet there has been little success combating it. This paper investigates the potential role of in silico modelling to help test emerging hypoxia-targeting treatments in cancer therapy. METHODS A Medline search was undertaken on the current landscape of in silico models that simulate cancer therapy and evaluate their ability to test hypoxia-targeting treatments. Techniques and treatments to combat tumour hypoxia and their current challenges are also presented. RESULTS Hypoxia-targeting treatments include tumour reoxygenation, hypoxic cell radiosensitization with nitroimidazoles, hypoxia-activated prodrugs and molecular targeting. Their main challenges are toxicity and not achieving adequate delivery to hypoxic regions of the tumour. There is promising research toward combining two or more of these techniques. Different types of in silico therapy models have been developed ranging from temporal to spatial and from stochastic to deterministic models. Numerous models have compared the effectiveness of different radiotherapy fractionation schedules for controlling hypoxic tumours. Similarly, models could help identify and optimize new treatments for overcoming hypoxia that utilize novel hypoxia-targeting technology. CONCLUSION Current therapy models should attempt to incorporate more sophisticated modelling of tumour angiogenesis/vasculature and vessel perfusion in order to become more useful for testing hypoxia-targeting treatments, which typically rely upon the tumour vasculature for delivery of additional oxygen, (pro)drugs and nanoparticles.
Collapse
Affiliation(s)
- Jake C Forster
- SA Medical Imaging, Department of Nuclear Medicine, The Queen Elizabeth Hospital, Woodville South, SA 5011, Australia; Department of Physics, University of Adelaide, North Terrace, Adelaide SA 5005, Australia
| | - Loredana G Marcu
- Faculty of Science, University of Oradea, Oradea 410087, Romania; Cancer Research Institute and School of Health Sciences, University of South Australia, Adelaide SA 5001, Australia.
| | - Eva Bezak
- Department of Physics, University of Adelaide, North Terrace, Adelaide SA 5005, Australia; Cancer Research Institute and School of Health Sciences, University of South Australia, Adelaide SA 5001, Australia
| |
Collapse
|
8
|
Grogan JA, Markelc B, Connor AJ, Muschel RJ, Pitt-Francis JM, Maini PK, Byrne HM. Predicting the Influence of Microvascular Structure On Tumor Response to Radiotherapy. IEEE Trans Biomed Eng 2017; 64:504-511. [PMID: 27623567 DOI: 10.1109/tbme.2016.2606563] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2024]
Abstract
OBJECTIVE The purpose of this study is to investigate how theoretical predictions of tumor response to radiotherapy (RT) depend on the morphology and spatial representation of the microvascular network. METHODS A hybrid multiscale model, which couples a cellular automaton model of tumor growth with a model for oxygen transport from blood vessels, is used to predict the viable fraction of cells following one week of simulated RT. Both artificial and biologically derived three-dimensional (3-D) vessel networks of well vascularized tumors are considered and predictions compared with 2-D descriptions. RESULTS For literature-derived values of the cellular oxygen consumption rate there is little difference in predicted viable fraction when 3-D network representations of biological or artificial vessel networks are employed. Different 2-D representations are shown to either over- or under-estimate viable fractions relative to the 3-D cases, with predictions based on point-wise descriptions shown to have greater sensitivity to vessel network morphology. CONCLUSION The predicted RT response is relatively insensitive to the morphology of the microvessel network when 3-D representations are adopted, however, sensitivity is greater in certain 2-D representations. SIGNIFICANCE By using realistic 3-D vessel network geometries this study shows that real and artificial network descriptions and assumptions of spatially uniform oxygen distributions lead to similar RT response predictions in relatively small tissue volumes. This suggests that either a more detailed description of oxygen transport in the microvasculature is required or that the oxygen enhancement ratio used in the well known linear-quadratic RT response model is relatively insensitive to microvascular structure.
Collapse
Affiliation(s)
- James A Grogan
- Mathematical Institute, University of Oxford, Oxford, U.K
| | - Bostjan Markelc
- Department of OncologyCRUK/MRC Oxford Institute for Radiation Oncology, University of Oxford
| | | | - Ruth J Muschel
- Department of OncologyCRUK/MRC Oxford Institute for Radiation Oncology, University of Oxford
| | | | | | | |
Collapse
|
9
|
Lindblom E, Dasu A, Beskow C, Toma-Dasu I. High brachytherapy doses can counteract hypoxia in cervical cancer—a modelling study. Phys Med Biol 2016; 62:560-572. [DOI: 10.1088/1361-6560/aa520f] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
10
|
Warren DR, Partridge M. The role of necrosis, acute hypoxia and chronic hypoxia in 18F-FMISO PET image contrast: a computational modelling study. Phys Med Biol 2016; 61:8596-8624. [PMID: 27880734 PMCID: PMC5717515 DOI: 10.1088/1361-6560/61/24/8596] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2016] [Revised: 09/14/2016] [Accepted: 10/26/2016] [Indexed: 12/22/2022]
Abstract
Positron emission tomography (PET) using 18F-fluoromisonidazole (FMISO) is a promising technique for imaging tumour hypoxia, and a potential target for radiotherapy dose-painting. However, the relationship between FMISO uptake and oxygen partial pressure ([Formula: see text]) is yet to be quantified fully. Tissue oxygenation varies over distances much smaller than clinical PET resolution (<100 μm versus ∼4 mm), and cyclic variations in tumour perfusion have been observed on timescales shorter than typical FMISO PET studies (∼20 min versus a few hours). Furthermore, tracer uptake may be decreased in voxels containing some degree of necrosis. This work develops a computational model of FMISO uptake in millimetre-scale tumour regions. Coupled partial differential equations govern the evolution of oxygen and FMISO distributions, and a dynamic vascular source map represents temporal variations in perfusion. Local FMISO binding capacity is modulated by the necrotic fraction. Outputs include spatiotemporal maps of [Formula: see text] and tracer accumulation, enabling calculation of tissue-to-blood ratios (TBRs) and time-activity curves (TACs) as a function of mean tissue oxygenation. The model is characterised using experimental data, finding half-maximal FMISO binding at local [Formula: see text] of 1.4 mmHg (95% CI: 0.3-2.6 mmHg) and half-maximal necrosis at 1.2 mmHg (0.1-4.9 mmHg). Simulations predict a non-linear non-monotonic relationship between FMISO activity (4 hr post-injection) and mean tissue [Formula: see text] : tracer uptake rises sharply from negligible levels in avascular tissue, peaking at ∼5 mmHg and declining towards blood activity in well-oxygenated conditions. Greater temporal variation in perfusion increases peak TBRs (range 2.20-5.27) as a result of smaller predicted necrotic fraction, rather than fundamental differences in FMISO accumulation under acute hypoxia. Identical late FMISO uptake can occur in regions with differing [Formula: see text] and necrotic fraction, but simulated TACs indicate that additional early-phase information may allow discrimination of hypoxic and necrotic signals. We conclude that a robust approach to FMISO interpretation (and dose-painting prescription) is likely to be based on dynamic PET analysis.
Collapse
Affiliation(s)
- Daniel R Warren
- CRUK/MRC Oxford Institute for Radiation Oncology, University of Oxford, Old Road Campus Research Building, Roosevelt Drive, Oxford OX3 7DQ, UK
| | - Mike Partridge
- CRUK/MRC Oxford Institute for Radiation Oncology, University of Oxford, Old Road Campus Research Building, Roosevelt Drive, Oxford OX3 7DQ, UK
| |
Collapse
|
11
|
Moghaddasi L, Bezak E, Harriss-Phillips W. Monte-Carlo model development for evaluation of current clinical target volume definition for heterogeneous and hypoxic glioblastoma. Phys Med Biol 2016; 61:3407-26. [PMID: 27046324 DOI: 10.1088/0031-9155/61/9/3407] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
Clinical target volume (CTV) determination may be complex and subjective. In this work a microscopic-scale tumour model was developed to evaluate current CTV practices in glioblastoma multiforme (GBM) external radiotherapy. Previously, a Geant4 cell-based dosimetry model was developed to calculate the dose deposited in individual GBM cells. Microscopic extension probability (MEP) models were then developed using Matlab-2012a. The results of the cell-based dosimetry model and MEP models were combined to calculate survival fractions (SF) for CTV margins of 2.0 and 2.5 cm. In the current work, oxygenation and heterogeneous radiosensitivity profiles were incorporated into the GBM model. The genetic heterogeneity was modelled using a range of α/β values (linear-quadratic model parameters) associated with different GBM cell lines. These values were distributed among the cells randomly, taken from a Gaussian-weighted sample of α/β values. Cellular oxygen pressure was distributed randomly taken from a sample weighted to profiles obtained from literature. Three types of GBM models were analysed: homogeneous-normoxic, heterogeneous-normoxic, and heterogeneous-hypoxic. The SF in different regions of the tumour model and the effect of the CTV margin extension from 2.0-2.5 cm on SFs were investigated for three MEP models. The SF within the beam was increased by up to three and two orders of magnitude following incorporation of heterogeneous radiosensitivities and hypoxia, respectively, in the GBM model. However, the total SF was shown to be overdominated by the presence of tumour cells in the penumbra region and to a lesser extent by genetic heterogeneity and hypoxia. CTV extension by 0.5 cm reduced the SF by a maximum of 78.6 ± 3.3%, 78.5 ± 3.3%, and 77.7 ± 3.1% for homogeneous and heterogeneous-normoxic, and heterogeneous hypoxic GBMs, respectively. Monte-Carlo model was developed to quantitatively evaluate SF for genetically heterogeneous and hypoxic GBM with two CTV margins and three MEP distributions. The results suggest that photon therapy may not provide cure for hypoxic and genetically heterogeneous GBM. However, the extension of the CTV margin by 0.5 cm could be beneficial to delay the recurrence time for this tumour type due to significant increase in tumour cell irradiation.
Collapse
Affiliation(s)
- L Moghaddasi
- Department of Medical Physics, Royal Adelaide Hospital, Adelaide, SA, Australia. School of Chemistry & Physics, University of Adelaide, Adelaide, SA, Australia
| | | | | |
Collapse
|
12
|
Lindblom E, Dasu A, Toma-Dasu I. Optimal fractionation in radiotherapy for non-small cell lung cancer--a modelling approach. Acta Oncol 2015. [PMID: 26217986 DOI: 10.3109/0284186x.2015.1061207] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
BACKGROUND Conventionally fractionated radiotherapy (CFRT) has proven ineffective in treating non-small cell lung cancer while more promising results have been obtained with stereotactic body radiotherapy (SBRT). Hypoxic tumours, however, might present a challenge to extremely hypofractionated schedules due to the decreased possibility for inter-fraction fast reoxygenation. A potentially successful compromise might be found in schedules employing several fractions of varying fractional doses. In this modelling study, a wide range of fractionation schedules from single-fraction treatments to heterogeneous, multifraction schedules taking into account repair, repopulation, reoxygenation and radiosensitivity of the tumour cells, has been explored with respect to the probability of controlling lung tumours. MATERIAL AND METHODS The response to radiation of tumours with heterogeneous spatial and temporal oxygenation was simulated including the effects of accelerated repopulation and intra-fraction repair. Various treatments with respect to time, dose and fractionation were considered and the outcome was estimated as Poisson-based tumour control probability for local control. RESULTS For well oxygenated tumours, heterogeneous fractionation could increase local control while hypoxic tumours are not efficiently targeted by such treatments despite reoxygenation. For hypofractionated treatments employing large doses per fraction, a synergistic effect was observed between intra-fraction repair and inter-fraction fast reoxygenation of the hypoxic cells as demonstrated by a reduction in D50 from 53.3 Gy for 2 fractions to 52.7 Gy for 5 fractions. CONCLUSIONS For well oxygenated tumours, heterogeneous fractionation schedules could increase local control rates substantially compared to CFRT. For hypoxic tumours, SBRT-like hypofractionated schedules might be optimal despite the increased risk of intra-fraction repair due to a synergistic effect with inter-fraction reoxygenation.
Collapse
Affiliation(s)
- Emely Lindblom
- a Medical Radiation Physics, Department of Physics , Stockholm University , Stockholm , Sweden
| | - Alexandru Dasu
- b Department of Radiation Physics and Department of Medical and Health Sciences , Linköping University , Linköping , Sweden
| | - Iuliana Toma-Dasu
- a Medical Radiation Physics, Department of Physics , Stockholm University , Stockholm , Sweden
- c Medical Radiation Physics, Department of Oncology and Pathology , Karolinska Institutet , Stockholm , Sweden
| |
Collapse
|
13
|
Espinoza I, Peschke P, Karger CP. A voxel-based multiscale model to simulate the radiation response of hypoxic tumors. Med Phys 2015; 42:90-102. [PMID: 25563250 DOI: 10.1118/1.4903298] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE In radiotherapy, it is important to predict the response of tumors to irradiation prior to the treatment. This is especially important for hypoxic tumors, which are known to be highly radioresistant. Mathematical modeling based on the dose distribution, biological parameters, and medical images may help to improve this prediction and to optimize the treatment plan. METHODS A voxel-based multiscale tumor response model for simulating the radiation response of hypoxic tumors was developed. It considers viable and dead tumor cells, capillary and normal cells, as well as the most relevant biological processes such as (i) proliferation of tumor cells, (ii) hypoxia-induced angiogenesis, (iii) spatial exchange of cells leading to tumor growth, (iv) oxygen-dependent cell survival after irradiation, (v) resorption of dead cells, and (vi) spatial exchange of cells leading to tumor shrinkage. Oxygenation is described on a microscopic scale using a previously published tumor oxygenation model, which calculates the oxygen distribution for each voxel using the vascular fraction as the most important input parameter. To demonstrate the capabilities of the model, the dependence of the oxygen distribution on tumor growth and radiation-induced shrinkage is investigated. In addition, the impact of three different reoxygenation processes is compared and tumor control probability (TCP) curves for a squamous cells carcinoma of the head and neck (HNSSC) are simulated under normoxic and hypoxic conditions. RESULTS The model describes the spatiotemporal behavior of the tumor on three different scales: (i) on the macroscopic scale, it describes tumor growth and shrinkage during radiation treatment, (ii) on a mesoscopic scale, it provides the cell density and vascular fraction for each voxel, and (iii) on the microscopic scale, the oxygen distribution may be obtained in terms of oxygen histograms. With increasing tumor size, the simulated tumors develop a hypoxic core. Within the model, tumor shrinkage was found to be significantly more important for reoxygenation than angiogenesis or decreased oxygen consumption due to an increased fraction of dead cells. In the studied HNSSC-case, the TCD50 values (dose at 50% TCP) decreased from 71.0 Gy under hypoxic to 53.6 Gy under the oxic condition. CONCLUSIONS The results obtained with the developed multiscale model are in accordance with expectations based on radiobiological principles and clinical experience. As the model is voxel-based, radiological imaging methods may help to provide the required 3D-characterization of the tumor prior to irradiation. For clinical application, the model has to be further validated with experimental and clinical data. If this is achieved, the model may be used to optimize fractionation schedules and dose distributions for the treatment of hypoxic tumors.
Collapse
Affiliation(s)
- I Espinoza
- Institute of Physics, Pontificia Universidad Católica de Chile, Santiago 7820436, Chile and Department of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
| | - P Peschke
- Clinical Cooperation Unit Molecular Radiooncology, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
| | - C P Karger
- Department of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
| |
Collapse
|
14
|
Antonovic L, Dasu A, Furusawa Y, Toma-Dasu I. Relative clinical effectiveness of carbon ion radiotherapy: theoretical modelling for H&N tumours. JOURNAL OF RADIATION RESEARCH 2015; 56:639-645. [PMID: 25858182 PMCID: PMC4497389 DOI: 10.1093/jrr/rrv016] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 10/31/2014] [Revised: 02/16/2015] [Accepted: 03/06/2015] [Indexed: 06/04/2023]
Abstract
Comparison of the efficiency of photon and carbon ion radiotherapy (RT) administered with the same number of fractions might be of limited clinical interest, since a wide range of fractionation patterns are used clinically today. Due to advanced photon treatment techniques, hypofractionation is becoming increasingly accepted for prostate and lung tumours, whereas patients with head and neck tumours still benefit from hyperfractionated treatments. In general, the number of fractions is considerably lower in carbon ion RT. A clinically relevant comparison would be between fractionation schedules that are optimal within each treatment modality category. In this in silico study, the relative clinical effectiveness (RCE) of carbon ions was investigated for human salivary gland tumours, assuming various radiation sensitivities related to their oxygenation. The results indicate that, for hypoxic tumours in the absence of reoxygenation, the RCE (defined as the ratio of D(50) for photons to carbon ions) ranges from 3.5 to 5.7, corresponding to carbon ion treatments given in 36 and 3 fractions, respectively, and 30 fractions for photons. Assuming that interfraction local oxygenation changes take place, results for RCE are lower than that for an oxic tumour if only a few fractions of carbon ions are used. If the carbon ion treatment is given in more than 12 fractions, the RCE is larger for the hypoxic than for the well-oxygenated tumour. In conclusion, this study showed that in silico modelling enables the study of a wide range of factors in the clinical considerations and could be an important step towards individualisation of RT treatments.
Collapse
Affiliation(s)
- Laura Antonovic
- Medical Radiation Physics, Department of Physics, Stockholm University, Stockholm, Sweden
| | - Alexandru Dasu
- Department of Radiation Physics and Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
| | - Yoshiya Furusawa
- Next Generation Medical Physics Research Program and International Open Laboratories, National Institute of Radiological Sciences, Chiba 263-8555, Japan
| | - Iuliana Toma-Dasu
- Medical Radiation Physics, Department of Physics, Stockholm University, Stockholm, Sweden Medical Radiation Physics, Department of Oncology and Pathology, Karolinska Institutet, Stockholm, Sweden
| |
Collapse
|
15
|
Manem VSK, Dhawan A, Kohandel M, Sivaloganathan S. Efficacy of dose escalation on TCP, recurrence and second cancer risks: a mathematical study. Br J Radiol 2014; 87:20140377. [PMID: 25210783 DOI: 10.1259/bjr.20140377] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
OBJECTIVE We investigated the effects of conventional and hypofractionation protocols by modelling tumour control probability (TCP) and tumour recurrence time, and examined their impact on second cancer risks. The main objectives of this study include the following: (a) incorporate tumour recurrence time and second cancer risks into the TCP framework and analyse the effects of variable doses and (b) investigate an efficient protocol to reduce the risk of a secondary malignancy while maximizing disease-free survival and tumour control. METHODS A generalized mathematical formalism was developed that incorporated recurrence and second cancer risk models into the TCP dynamics. RESULTS Our results suggest that TCP and relapse time are almost identical for conventional and hypofractionated regimens; however, second cancer risks resulting from hypofractionation were reduced by 22% when compared with the second cancer risk associated with a conventional protocol. The hypofractionated regimen appears to be sensitive to dose escalation and the corresponding impact on tumour recurrence time and reduction in second cancer risks. The reduction in second cancer risks is approximately 20% when the dose is increased from 60 to 72 Gy in a hypofractionated protocol. CONCLUSION Our results suggest that hypofractionation may be a more efficient regimen in the context of TCP, relapse time and second cancer risks. Overall, our study demonstrates the importance of including a second cancer risk model in designing an efficient radiation regimen. ADVANCES IN KNOWLEDGE The impact of various fractionation protocols on TCP and relapse in conjunction with second cancer risks is an important clinical question that is as yet unexplored.
Collapse
Affiliation(s)
- V S K Manem
- 1 Department of Applied Mathematics, University of Waterloo, Waterloo, ON, Canada
| | | | | | | |
Collapse
|
16
|
Antonovic L, Lindblom E, Dasu A, Bassler N, Furusawa Y, Toma-Dasu I. Clinical oxygen enhancement ratio of tumors in carbon ion radiotherapy: the influence of local oxygenation changes. JOURNAL OF RADIATION RESEARCH 2014; 55:902-11. [PMID: 24728013 PMCID: PMC4240637 DOI: 10.1093/jrr/rru020] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/09/2013] [Revised: 02/25/2014] [Accepted: 03/07/2014] [Indexed: 05/28/2023]
Abstract
The effect of carbon ion radiotherapy on hypoxic tumors has recently been questioned because of low linear energy transfer (LET) values in the spread-out Bragg peak (SOBP). The aim of this study was to investigate the role of hypoxia and local oxygenation changes (LOCs) in fractionated carbon ion radiotherapy. Three-dimensional tumors with hypoxic subvolumes were simulated assuming interfraction LOCs. Different fractionations were applied using a clinically relevant treatment plan with a known LET distribution. The surviving fraction was calculated, taking oxygen tension, dose and LET into account, using the repairable-conditionally repairable (RCR) damage model with parameters for human salivary gland tumor cells. The clinical oxygen enhancement ratio (OER) was defined as the ratio of doses required for a tumor control probability of 50% for hypoxic and well-oxygenated tumors. The resulting OER was well above unity for all fractionations. For the hypoxic tumor, the tumor control probability was considerably higher if LOCs were assumed, rather than static oxygenation. The beneficial effect of LOCs increased with the number of fractions. However, for very low fraction doses, the improvement related to LOCs did not compensate for the increase in total dose required for tumor control. In conclusion, our results suggest that hypoxia can influence the outcome of carbon ion radiotherapy because of the non-negligible oxygen effect at the low LETs in the SOBP. However, if LOCs occur, a relatively high level of tumor control probability is achievable with a large range of fractionation schedules for tumors with hypoxic subvolumes, but both hyperfractionation and hypofractionation should be pursued with caution.
Collapse
Affiliation(s)
- Laura Antonovic
- Department of Physics, Stockholm University, Stockholm, Sweden
| | - Emely Lindblom
- Department of Physics, Stockholm University, Stockholm, Sweden
| | - Alexandru Dasu
- Department of Radiation Physics and Department of Medical and Health Sciences, Linköping University, Linköping, Sweden
| | - Niels Bassler
- Department of Physics and Astronomy, Aarhus University, Aarhus, Denmark and Department of Experimental Clinical Oncology, Aarhus University Hospital, Aarhus, Denmark
| | - Yoshiya Furusawa
- Next Generation Medical Physics Research Program and International Open Laboratories, National Institute of Radiological Sciences, Chiba 263-8555, Japan
| | - Iuliana Toma-Dasu
- Department of Physics, Stockholm University, Stockholm, Sweden Department of Oncology and Pathology, Karolinska Institute, Stockholm, Sweden
| |
Collapse
|
17
|
Lindblom E, Dasu A, Lax I, Toma-Dasu I. Survival and tumour control probability in tumours with heterogeneous oxygenation: a comparison between the linear-quadratic and the universal survival curve models for high doses. Acta Oncol 2014; 53:1035-40. [PMID: 24957551 DOI: 10.3109/0284186x.2014.925582] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
BACKGROUND The validity of the linear-quadratic (LQ) model at high doses has been questioned due to a decreasing agreement between predicted survival and experimental cell survival data. A frequently proposed alternative is the universal survival curve (USC) model, thought to provide a better fit in the high-dose region. The comparison between the predictions of the models has mostly been performed for uniform populations of cells with respect to sensitivity to radiation. This study aimed to compare the two models in terms of cell survival and tumour control probability (TCP) for cell populations with mixed sensitivities related to their oxygenation. METHODS The study was performed in two parts. For the first part, cell survival curves were calculated with both models assuming various homogeneous populations of cells irradiated with uniform doses. For the second part, a realistic three-dimensional (3D) model of complex tumour oxygenation was used to study the impact of the differences in cell survival on the modelled TCP. Cellular response was assessed with the LQ and USC models at voxel level and a Poisson TCP model at tumour level. RESULTS For hypoxic tumours, the disputed continuous bend of the LQ survival curve was counteracted by the increased radioresistance of the hypoxic cells and the survival curves started to diverge only at much higher doses than for oxic tumours. This was also reflected by the TCP curves for hypoxic tumours for which the difference in D50 values for the LQ and USC models was reduced from 5.4 to 0.2 Gy for 1 and 3 fractions, respectively, in a tumour with only 1.1% hypoxia and from 9.5 to 0.4 Gy in a tumour with 11.1% hypoxia. CONCLUSIONS For a large range of fractional doses including hypofractionated schemes, the difference in predicted survival and TCP between the LQ and USC models for tumours with heterogeneous oxygenation was found to be negligible.
Collapse
Affiliation(s)
- Emely Lindblom
- Medical Radiation Physics, Department of Physics, Stockholm University , Stockholm , Sweden
| | | | | | | |
Collapse
|
18
|
Lindblom E, Antonovic L, Dasu A, Lax I, Wersäll P, Toma-Dasu I. Treatment fractionation for stereotactic radiotherapy of lung tumours: a modelling study of the influence of chronic and acute hypoxia on tumour control probability. Radiat Oncol 2014; 9:149. [PMID: 24974778 PMCID: PMC4091751 DOI: 10.1186/1748-717x-9-149] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2014] [Accepted: 06/19/2014] [Indexed: 11/10/2022] Open
Abstract
Background Stereotactic body radiotherapy (SBRT) for non-small-cell lung cancer (NSCLC) has led to promising local control and overall survival for fractionation schemes with increasingly high fractional doses. A point has however been reached where the number of fractions used might be too low to allow efficient local inter-fraction reoxygenation of the hypoxic cells residing in the tumour. It was therefore the purpose of this study to investigate the impact of hypoxia and extreme hypofractionation on the tumour control probability (TCP) from SBRT. Methods A three-dimensional model of tumour oxygenation able to simulate oxygenation changes on the microscale was used. The TCP was determined for clinically relevant SBRT fractionation schedules of 1, 3 and 5 fractions assuming either static tumour oxygenation or that the oxygenation changes locally between fractions due to fast reoxygenation of acute hypoxia without an overall reduction in chronic hypoxia. Results For the schedules applying three or five fractions the doses required to achieve satisfying levels of TCP were considerably lower when local oxygenation changes were assumed compared to the case of static oxygenation; a decrease in D50 of 17.7 Gy was observed for a five-fractions schedule applied to a 20% hypoxic tumour when fast reoxygenation was modelled. Assuming local oxygenation changes, the total doses required for a tumor control probability of 50% were of similar size for one, three and five fractions. Conclusions Although attractive from a practical point of view, extreme hypofractionation using just one single fraction may result in impaired local control of hypoxic tumours, as it eliminates the possibility for any kind of reoxygenation.
Collapse
Affiliation(s)
- Emely Lindblom
- Medical Radiation Physics, Department of Physics, Stockholm University, Stockholm, Sweden.
| | | | | | | | | | | |
Collapse
|
19
|
Dhawan A, Kohandel M, Hill R, Sivaloganathan S. Tumour control probability in cancer stem cells hypothesis. PLoS One 2014; 9:e96093. [PMID: 24811314 PMCID: PMC4014481 DOI: 10.1371/journal.pone.0096093] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2013] [Accepted: 04/02/2014] [Indexed: 12/02/2022] Open
Abstract
The tumour control probability (TCP) is a formalism derived to compare various treatment regimens of radiation therapy, defined as the probability that given a prescribed dose of radiation, a tumour has been eradicated or controlled. In the traditional view of cancer, all cells share the ability to divide without limit and thus have the potential to generate a malignant tumour. However, an emerging notion is that only a sub-population of cells, the so-called cancer stem cells (CSCs), are responsible for the initiation and maintenance of the tumour. A key implication of the CSC hypothesis is that these cells must be eradicated to achieve cures, thus we define TCPS as the probability of eradicating CSCs for a given dose of radiation. A cell surface protein expression profile, such as CD44high/CD24low for breast cancer or CD133 for glioma, is often used as a biomarker to monitor CSCs enrichment. However, it is increasingly recognized that not all cells bearing this expression profile are necessarily CSCs, and in particular early generations of progenitor cells may share the same phenotype. Thus, due to the lack of a perfect biomarker for CSCs, we also define a novel measurable TCPCD+, that is the probability of eliminating or controlling biomarker positive cells. Based on these definitions, we use stochastic methods and numerical simulations parameterized for the case of gliomas, to compare the theoretical TCPS and the measurable TCPCD+. We also use the measurable TCP to compare the effect of various radiation protocols.
Collapse
Affiliation(s)
- Andrew Dhawan
- Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada
- School of Medicine, Queen's University, Kingston, Ontario, Canada
| | - Mohammad Kohandel
- Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada
- Center for Mathematical Medicine, Fields Institute for Research in Mathematical Sciences, Toronto, Ontario, Canada
- * E-mail:
| | - Richard Hill
- Department of Medical Biophysics, Princess Margaret Hospital, Ontario Cancer Institute, University of Toronto, Toronto, Ontario, Canada
| | - Sivabal Sivaloganathan
- Department of Applied Mathematics, University of Waterloo, Waterloo, Ontario, Canada
- Center for Mathematical Medicine, Fields Institute for Research in Mathematical Sciences, Toronto, Ontario, Canada
| |
Collapse
|
20
|
Espinoza I, Peschke P, Karger CP. A model to simulate the oxygen distribution in hypoxic tumors for different vascular architectures. Med Phys 2014; 40:081703. [PMID: 23927300 DOI: 10.1118/1.4812431] [Citation(s) in RCA: 24] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022] Open
Abstract
PURPOSE As hypoxic cells are more resistant to photon radiation, it is desirable to obtain information about the oxygen distribution in tumors prior to the radiation treatment. Noninvasive techniques are currently not able to provide reliable oxygenation maps with sufficient spatial resolution; therefore mathematical models may help to simulate microvascular architectures and the resulting oxygen distributions in the surrounding tissue. Here, the authors present a new computer model, which uses the vascular fraction of tumor voxels, in principle measurable noninvasively in vivo, as input parameter for simulating realistic PO2 histograms in tumors, assuming certain 3D vascular architectures. METHODS Oxygen distributions were calculated by solving a reaction-diffusion equation in a reference volume using the particle strength exchange method. Different types of vessel architectures as well as different degrees of vascular heterogeneities are considered. Two types of acute hypoxia (ischemic and hypoxemic) occurring additionally to diffusion-limited (chronic) hypoxia were implemented as well. RESULTS No statistically significant differences were observed when comparing 2D- and 3D-vessel architectures (p>0.79 in all cases) and highly heterogeneously distributed linear vessels show good agreement, when comparing with published experimental intervessel distance distributions and PO2 histograms. It could be shown that, if information about additional acute hypoxia is available, its contribution to the hypoxic fraction (HF) can be simulated as well. Increases of 128% and 168% in the HF were obtained when representative cases of ischemic and hypoxemic acute hypoxia, respectively, were considered in the simulations. CONCLUSIONS The presented model is able to simulate realistic microscopic oxygen distributions in tumors assuming reasonable vessel architectures and using the vascular fraction as macroscopic input parameter. The model may be used to generate PO2 histograms, which are needed as input in models predicting the radiation response of hypoxic tumors.
Collapse
Affiliation(s)
- Ignacio Espinoza
- Department of Medical Physics in Radiation Oncology, German Cancer Research Center (DKFZ), Heidelberg 69120, Germany
| | | | | |
Collapse
|
21
|
Hedman M, Björk-Eriksson T, Brodin O, Toma-Dasu I. Predictive value of modelled tumour control probability based on individual measurements of in vitro radiosensitivity and potential doubling time. Br J Radiol 2013; 86:20130015. [PMID: 23479396 DOI: 10.1259/bjr.20130015] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
OBJECTIVE The aim of this study was to compare patient-specific radiobiological parameters with population averages in predicting the clinical outcome after radiotherapy (RT) using a tumour control probability (TCP) model based on the biological effective dose (BED). METHODS A previously published study of 46 head and neck carcinomas with individually identified radiobiological parameters, radiosensitivity and potential doubling time (Tpot), and known tumour size was investigated. These patients had all been treated with external beam RT, and the majority had also received brachytherapy. The TCP for each individual based on the BED using patient-specific radiobiological parameters was compared with the TCP based on the BED using average radiobiological parameters (α=0.3 Gy(-1), Tpot=3 days). RESULTS 43 patients remained in the final analysis. There was only a weak trend for increasing local tumour control with increasing BED in both groups. However, when the TCP was calculated, the use of patient-specific parameters was better for identifying local control correctly. The sensitivity and specificity for tumour-specific parameters were 63% and 80%, respectively. The corresponding values for population-based averages were 0% and 91%, respectively. The positive predictive value was 92% when tumour-specific parameters were used compared with 0% for population-based averages. A receiver operating characteristic curve confirmed the superiority of patient-specific parameters over population averages in predicting local control. CONCLUSION Individual radiobiological parameters are better than population-derived averages when used in a mathematical model to predict TCP after curative RT in head and neck carcinomas. ADVANCES IN KNOWLEDGE TCP based on individual radiobiological parameters is better than TCP based on population-based averages for identifying local control correctly.
Collapse
Affiliation(s)
- M Hedman
- Department of Oncology and Pathology, Karolinska Institutet, Stockholm, Sweden.
| | | | | | | |
Collapse
|
22
|
Harriss-Phillips WM, Bezak E, Yeoh EK. Altered fractionation outcomes for hypoxic head and neck cancer using the HYP-RT Monte Carlo model. Br J Radiol 2013; 86:20120443. [PMID: 23392195 DOI: 10.1259/bjr.20120443] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022] Open
Abstract
OBJECTIVE Altered fractionation radiotherapy is simulated on a set of virtual tumours to assess the total doses required for tumour control compared with clinical head and neck data and the doses required to control hypoxic vs well-oxygenated tumours with different radiobiological properties. METHODS The HYP-RT model is utilised to explore the impact of tumour oxygenation and the onset times of accelerated repopulation (AR) and reoxygenation (ROx) during radiotherapy. A biological effective dose analysis is used to rank the schedules based on their relative normal tissue toxicities. RESULTS Altering the onset times of AR and ROx has a large impact on the doses required to achieve tumour control. Immediate onset of ROx and 2-week onset time of AR produce results closely predicting average human outcomes in terms of the total prescription doses in clinical trials. Modifying oxygen enhancement ratio curves based on dose/fraction significantly reduces the dose (5-10 Gy) required for tumour control for hyperfractionated schedules. HYP-RT predicts 10×1.1 Gy per week to be most beneficial, whereas the conventional schedule is predicted as beneficial for early toxicity but has average-poor late toxicity. CONCLUSION HYP-RT predicts that altered radiotherapy schedules increase the therapeutic ratio and may be used to make predictions about the prescription doses required to achieve tumour control for tumours with different oxygenation levels and treatment responses. ADVANCES IN KNOWLEDGE Oxic and hypoxic tumours have large differences in total radiation dose requirements, affected by AR and ROx onset times by up to 15-25 Gy for the same fractionation schedule.
Collapse
Affiliation(s)
- W M Harriss-Phillips
- Department of Medical Physics, Royal Adelaide Hospital Cancer Centre, South Australia, Australia.
| | | | | |
Collapse
|
23
|
Lagerlöf JH, Kindblom J, Cortez E, Pietras K, Bernhardt P. Image-based 3D modeling study of the influence of vessel density and blood hemoglobin concentration on tumor oxygenation and response to irradiation. Med Phys 2013; 40:024101. [DOI: 10.1118/1.4773886] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
|
24
|
In silico modelling of treatment-induced tumour cell kill: developments and advances. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2012; 2012:960256. [PMID: 22852024 PMCID: PMC3407630 DOI: 10.1155/2012/960256] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/14/2012] [Revised: 05/10/2012] [Accepted: 05/14/2012] [Indexed: 12/04/2022]
Abstract
Mathematical and stochastic computer (in silico) models of tumour growth and treatment response of the past and current eras are presented, outlining the aims of the models, model methodology, the key parameters used to describe the tumour system, and treatment modality applied, as well as reported outcomes from simulations. Fractionated radiotherapy, chemotherapy, and combined therapies are reviewed, providing a comprehensive overview of the modelling literature for current modellers and radiobiologists to ignite the interest of other computational scientists and health professionals of the ever evolving and clinically relevant field of tumour modelling.
Collapse
|
25
|
Modeling the spatial distribution of chronic tumor hypoxia: implications for experimental and clinical studies. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2012; 2012:410602. [PMID: 22400049 PMCID: PMC3287099 DOI: 10.1155/2012/410602] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/30/2011] [Accepted: 10/17/2011] [Indexed: 11/18/2022]
Abstract
Tumor oxygenation status is considered one of the important prognostic markers in cancer since it strongly influences the response of cancer cells to various treatments; in particular, to radiation therapy. Thus, a proper and accurate assessment of tumor oxygen distribution before the treatment may highly affect the outcome of the treatment. The heterogeneous nature of tumor hypoxia, mainly influenced by the complex tumor microenvironment, often makes its quantification very difficult. The usual methods used to measure tumor hypoxia are biomarkers and the polarographic needle electrode. Although these techniques may provide an acceptable assessment of hypoxia, they are invasive and may not always give a spatial distribution of hypoxia, which is very useful for treatment planning. An alternative method to quantify the tumor hypoxia is to use theoretical simulations with the knowledge of tumor vasculature. The purpose of this paper is to model tumor hypoxia using a known spatial distribution of tumor vasculature obtained from image data, to analyze the accuracy of polarographic needle electrode measurements in quantifying hypoxia, to quantify the optimum number of measurements required to satisfactorily evaluate the tumor oxygenation status, and to study the effects of hypoxia on radiation response. Our results indicate that the model successfully generated an accurate oxygenation map for tumor cross-sections with known vascular distribution. The method developed here provides a way to estimate tumor hypoxia and provides guidance in planning accurate and effective therapeutic strategies and invasive estimation techniques. Our results agree with the previous findings that the needle electrode technique gives a good estimate of tumor hypoxia if the sampling is done in a uniform way with 5-6 tracks of 20–30 measurements each. Moreover, the analysis indicates that the accurate measurement of oxygen profile can be very useful in determining right radiation doses to the patients.
Collapse
|
26
|
Tuckwell W, Bezak E, Yeoh E, Marcu L. Efficient Monte Carlo modelling of individual tumour cell propagation for hypoxic head and neck cancer. Phys Med Biol 2008; 53:4489-507. [PMID: 18677039 DOI: 10.1088/0031-9155/53/17/002] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Abstract
A Monte Carlo tumour model has been developed to simulate tumour cell propagation for head and neck squamous cell carcinoma. The model aims to eventually provide a radiobiological tool for radiation oncology clinicians to plan patient treatment schedules based on properties of the individual tumour. The inclusion of an oxygen distribution amongst the tumour cells enables the model to incorporate hypoxia and other associated parameters, which affect tumour growth. The object oriented program FORTRAN 95 has been used to create the model algorithm, with Monte Carlo methods being employed to randomly assign many of the cell parameters from probability distributions. Hypoxia has been implemented through random assignment of partial oxygen pressure values to individual cells during tumour growth, based on in vivo Eppendorf probe experimental data. The accumulation of up to 10 million virtual tumour cells in 15 min of computer running time has been achieved. The stem cell percentage and the degree of hypoxia are the parameters which most influence the final tumour growth rate. For a tumour with a doubling time of 40 days, the final stem cell percentage is approximately 1% of the total cell population. The effect of hypoxia on the tumour growth rate is significant. Using a hypoxia induced cell quiescence limit which affects 50% of cells with and oxygen levels less than 1 mm Hg, the tumour doubling time increases to over 200 days and the time of tumour growth for a clinically detectable tumour (10(9) cells) increases from 3 to 8 years. A biologically plausible Monte Carlo model of hypoxic head and neck squamous cell carcinoma tumour growth has been developed for real time assessment of the effects of multiple biological parameters which impact upon the response of the individual patient to fractionated radiotherapy.
Collapse
Affiliation(s)
- W Tuckwell
- School of Chemistry and Physics, University of Adelaide, South Australia, Australia.
| | | | | | | |
Collapse
|
27
|
Kapanen M, Tenhunen M, Parkkinen R, Sipilä P, Järvinen H. The influence of output measurement time interval and tolerance on treatment dose deviation in photon external beam radiotherapy. Phys Med Biol 2006; 51:4857-67. [PMID: 16985275 DOI: 10.1088/0031-9155/51/19/009] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Abstract
A method was developed to estimate the spread of treatment doses due to a given output measurement, time interval and action level in external photon beam radiotherapy. The method exploits empirical or simulated data of both output stability and the reproducibility of the absolute output measurements and relative constancy checks. The spread of an output and the probability that an output is outside the predefined tolerance limits were estimated for different quality control procedures of the output. The use of internationally recommended output measurement intervals of 3, 1 and 0.25 months together with an action level of +/-2% results in the out-of-tolerance probability up to 49%, 26% and 3%, respectively, while the output is kept within up to +/-5.5%, +/-4.0% and +/-2.0%, respectively. It was also found that the action level has more prominent relative influence than the measurement interval. The results suggest that similar treatment quality can be achieved by doubling the output measurement time interval and lowering the action level from +/-2% to +/-1.5% and the constancy check action level from +/-3% to +/-2%.
Collapse
Affiliation(s)
- Mika Kapanen
- Radiation and Nuclear Safety Authority, STUK, Laippatie 4, 00881 Helsinki, Finland.
| | | | | | | | | |
Collapse
|