Modeling the relationship between fluorodeoxyglucose uptake and tumor radioresistance as a function of the tumor microenvironment

Data Science Opportunities To Improve Radiotherapy Planning and Clinical Decision Making

Radiotherapy aims to achieve a high tumor control probability while minimizing damage to normal tissues. Personalizing radiotherapy treatments for individual patients, therefore, depends on integrating physical treatment planning with predictive models of tumor control and normal tissue complications. Predictive models could be improved using a wide range of rich data sources, including tumor and normal tissue genomics, radiomics, and dosiomics. Deep learning will drive improvements in classifying normal tissue tolerance, predicting intra-treatment tumor changes, tracking accumulated dose distributions, and quantifying the tumor response to radiotherapy based on imaging. Mechanistic patient-specific computer simulations ('digital twins') could also be used to guide adaptive radiotherapy. Overall, we are entering an era where improved modeling methods will allow the use of newly available data sources to better guide radiotherapy treatments.

A theoretical investigation of adequate range uncertainty margins in proton treatment planning to preserve tumor control probability

Background: Proton dose distributions are sensitive to range uncertainties, resulting in margins added to ensure adequate tumor control probability (TCP). We investigated the required margin and dose shape needed to ensure adequate TCP, for representative tumor cell distributions in the clinical target volume (CTV). Material and methods: A mechanistic tumor response model, validated for lung tumors, was used to estimate TCP. The tumor cell distribution ( ρ ) was assumed to decrease exponentially in the CTV with decay parameter λ toward the outer border ( xCTVmax ). It was investigated if a gradual dose fall-off could reduce the dose to normal tissues outside the CTV, while achieving adequate TCP. For various values of xCTVmax and λ, we derived adequate uniform dose margins ( m ), coupled to linear dose fall-off regions ( Δx, Δxnom=Δx-0.9 cm), that ensured TCP>TCPlimit, while delivering the least mean dose outside the CTV. To account for variabilities in patients and tumor types, variable probabilities ( p ) of finding tumor cells in the non-GTV part of the CTV for a given patient were also tested. Dose from a single beam or two opposing beams was simulated under the influence of a typical stopping power ratio uncertainty of 3.5%. Results: For large λ and xCTVmax, a dose distribution with a shallower dose fall-off ( Δx>0 ) was advantageous, and m could be smaller than xCTVmax. In the case of small xCTVmax values, however, a conventional dose distribution ( Δx=0 ) would generally perform better. For no CTV, m=0.4 cm in the case of two opposing beams, while it was 0.7 cm for a single beam, however, for two opposing beams Δx=1.2 cm ( Δxnom=0.3 cm), while it was zero for a single beam. Conclusion: The details of the underlying cancer cell distribution characteristics do impact the adequate dose arrangements, and for opposing beams a non-conventional dose distribution shape is often advantageous.

Modeling the Cellular Response of Lung Cancer to Radiation Therapy for a Broad Range of Fractionation Schedules

Purpose: To demonstrate that a mathematical model can be used to quantitatively understand tumor cellular dynamics during a course of radiotherapy and to predict the likelihood of local control as a function of dose and treatment fractions.Experimental Design: We model outcomes for early-stage, localized non-small cell lung cancer (NSCLC), by fitting a mechanistic, cellular dynamics-based tumor control probability that assumes a constant local supply of oxygen and glucose. In addition to standard radiobiological effects such as repair of sub-lethal damage and the impact of hypoxia, we also accounted for proliferation as well as radiosensitivity variability within the cell cycle. We applied the model to 36 published and two unpublished early-stage patient cohorts, totaling 2,701 patients.Results: Precise likelihood best-fit values were derived for the radiobiological parameters: α [0.305 Gy^-1; 95% confidence interval (CI), 0.120-0.365], the α/β ratio (2.80 Gy; 95% CI, 0.40-4.40), and the oxygen enhancement ratio (OER) value for intermediately hypoxic cells receiving glucose but not oxygen (1.70; 95% CI, 1.55-2.25). All fractionation groups are well fitted by a single dose-response curve with a high χ² P value, indicating consistency with the fitted model. The analysis was further validated with an additional 23 patient cohorts (n = 1,628). The model indicates that hypofractionation regimens overcome hypoxia (and cell-cycle radiosensitivity variations) by the sheer impact of high doses per fraction, whereas lower dose-per-fraction regimens allow for reoxygenation and corresponding sensitization, but lose effectiveness for prolonged treatments due to proliferation.Conclusions: This proposed mechanistic tumor-response model can accurately predict overtreatment or undertreatment for various treatment regimens. Clin Cancer Res; 23(18); 5469-79. ©2017 AACR.

A radiobiological model of radiotherapy response and its correlation with prognostic imaging variables

Radiobiological models of tumour control probability (TCP) can be personalized using imaging data. We propose an extension to a voxel-level radiobiological TCP model in order to describe patient-specific differences and intra-tumour heterogeneity. In the proposed model, tumour shrinkage is described by means of a novel kinetic Monte Carlo method for inter-voxel cell migration and tumour deformation. The model captures the spatiotemporal evolution of the tumour at the voxel level, and is designed to take imaging data as input. To test the performance of the model, three image-derived variables found to be predictive of outcome in the literature have been identified and calculated using the model's own parameters. Simulating multiple tumours with different initial conditions makes it possible to perform an in silico study of the correlation of these variables with the dose for 50% tumour control ([Formula: see text]) calculated by the model. We find that the three simulated variables correlate with the calculated [Formula: see text]. In addition, we find that different variables have different levels of sensitivity to the spatial distribution of hypoxia within the tumour, as well as to the dynamics of the migration mechanism. Finally, based on our results, we observe that an adequate combination of the variables may potentially result in higher predictive power.