Novel radiation quality metrics accounting for proton energy spectra for RBE proton models.

BACKGROUND: For proton therapy, a relative biological effectiveness (RBE) of 1.1 is widely applied clinically. However, due to abundant evidence of variable RBE in vitro, and as suggested in studies of patient outcomes, RBE might increase by the end of the proton tracks, as described by several proposed variable RBE models. Typically, the dose averaged linear energy transfer ( LET d $\text{LET}_d$ ) has been used as a radiation quality metric (RQM) for these models. However, the optimal choice of RQM has not been fully explored. PURPOSE: This study aims to propose novel RQMs that effectively weight protons of different energies, and assess their predictive power for variable RBE in proton therapy. The overall objective is to identify an RQM that better describes the contribution of individual particles to the RBE of proton beams. METHODS: High-throughput experimental set-ups of in vitro cell survival studies for proton RBE determination are simulated utilizing the SHIELD-HIT12A Monte Carlo particle transport code. For every data point, the proton energy spectra are simulated, allowing the calculation of novel RQMs by applying different power levels to the spectra of LET or effective Q $Q$ ( Q eff $Q_\mathrm{eff}$ ) values. A phenomenological linear-quadratic-based RBE model is then applied to the in vitro data, using various RQMs as input variables, and the model performance is evaluated by root-mean-square-error (RMSE) for the logarithm of cell surviving fractions of each data point. RESULTS: Increasing the power level, that is, putting an even higher weight on higher LET particles when constructing the RQM is generally associated with an increased model performance, with dose averaged LET 3 $\text{LET}^3$ (i.e., dose averaged cubed LET, cLET d $\mathrm{cLET}_d$ ) resulting in a RMSE value 0.31, compared to 0.45 for a model based on (linearly weighted) LET d $\text{LET}_d$ , with similar trends also observed for track averaged and Q eff $Q_\mathrm{eff}$ -based RQMs. CONCLUSIONS: The results indicate that improved proton variable RBE models can be constructed assuming a non-linear RBE(LET) relationship for individual protons. If similar trends hold also for an in vitro-environment, variable RBE effects are likely better described by cLET d $\mathrm{cLET}_d$ or tracked averaged cubed LET ( cLET t $\mathrm{cLET}_t$ ), or corresponding Q eff $Q_\mathrm{eff}$ -based RQM, rather than linearly weighted LET d $\text{LET}_d$ or LET t $\text{LET}_t$ which is conventionally applied today.

Patterns of practice of image guided particle therapy for cranio-spinal irradiation: A site specific multi-institutional survey of European Particle Therapy Network.

PURPOSE: To investigate the current practice patterns in image-guided particle therapy (IGPT) for cranio-spinal irradiation (CSI). METHODS: A multi-institutional survey was distributed to European particle therapy centres to analyse all aspects of IGPT. Based on the survey results, a Delphi consensus analysis was developed to define minimum requirements and optimal workflow for clinical practice. The centres participating in the institutional survey were invited to join the Delphi process. RESULTS: Eleven centres participated in the survey. Imaging for treatment planning was rather similar among the centres with Computed Tomography (CT) being the main modality. For positioning verification, 2D IGPT was more commonly used than 3D IGPT. Two centres performed routinely imaging for plan adaptation, by the rest ad hoc. Eight centres participated in the Delphi consensus analysis. The full consensus was reached on the use of CT imaging without contrast for treatment planning and the role of magnetic resonance imaging (MRI) in target and organs-at-risk delineation. There was an agreement on the necessity to perform patient position verification and correction before each isocentre. The most important outcome was the clear need for standardization and harmonization of the workflow. CONCLUSION: There were differences in CSI IGPT clinical practice among the European particle therapy centres. Moreover, the optimal workflow as identified by experts was not yet reached. There is a strong need for consensus guidelines. The state-of-the-art imaging technology and protocols need to be implemented into clinical practice to improve the quality of IGPT for CSI.

Generation and evaluation of anatomy-preserving virtual CT for online adaptive proton therapy.

BACKGROUND: Daily CTs generated by CBCT correction are required for daily replanning in online-adaptive proton therapy (APT) to effectively deal with inter-fractional changes. Out of the currently available methods, the suitability of a daily CT generation method for proton dose calculation also depends on the anatomical site. PURPOSE: We propose an anatomy-preserving virtual CT (APvCT) method as a hybrid method of CBCT correction, which is especially suitable for large anatomy deformations. The accuracy of the hybrid method was assessed by comparison with the corrected CBCT (cCBCT) and virtual CT (vCT) methods in the context of online APT. METHODS: Seventy-one daily CBCTs of four prostate cancer patients treated with intensity modulated proton therapy (IMPT) were converted to daily CTs using cCBCT, vCT, and the newly proposed APvCT method. In APvCT, planning CT (pCT) were mapped to CBCT geometry using deformable image registration with boundary conditions on controlling regions of interest (ROIs) created with deep learning segmentation on cCBCT. The relative frequency distribution (RFD) of HU, mass density and stopping power ratio (SPR) values were assessed and compared with the pCT. The ROIs in the APvCT and vCT were compared with cCBCT in terms of Dice similarity coefficient (DSC) and mean distance-to-agreement (mDTA). For each patient, a robustly optimized IMPT plan was created on the pCT and subsequent daily adaptive plans on daily CTs. For dose distribution comparison on the same anatomy, the daily adaptive plans on cCBCT and vCT were recalculated on the corresponding APvCT. The dose distributions were compared in terms of isodose volumes and 3D global gamma-index passing rate (GPR) at γ(2%, 2 mm) criterion. RESULTS: For all patients, no noticeable difference in RFDs was observed amongst APvCT, vCT, and pCT except in cCBCT, which showed a noticeable difference. The minimum DSC value was 0.96 and 0.39 for contours in APvCT and vCT respectively. The average value of mDTA for APvCT was 0.01 cm for clinical target volume and ≤0.01 cm for organs at risk, which increased to 0.18 cm and ≤0.52 cm for vCT. The mean GPR value was 90.9%, 64.5%, and 67.0% for APvCT versus cCBCT, vCT versus cCBCT, and APvCT versus vCT, respectively. When recalculated on APvCT, the adaptive cCBCT and vCT plans resulted in mean GPRs of 89.5 ± 5.1% and 65.9 ± 19.1%, respectively. The mean DSC values for 80.0%, 90.0%, 95.0%, 98.0%, and 100.0% isodose volumes were 0.97, 0.97, 0.97, 0.95, and 0.91 for recalculated cCBCT plans, and 0.89, 0.88, 0.87, 0.85, and 0.81 for recalculated vCT plans. Hausdorff distance for the 100.0% isodose volume in some cases of recalculated cCBCT plans on APvCT exceeded 1.00 cm. CONCLUSIONS: APvCT contours showed good agreement with reference contours of cCBCT which indicates anatomy preservation in APvCT. A vCT with erroneous anatomy can result in an incorrect adaptive plan. Further, slightly lower values of GPR between the APvCT and cCBCT-based adaptive plans can be explained by the difference in the cCBCT's SPR RFD from the pCT.

Impact of the Mean Cochlear Biologically Effective Dose on Hearing Preservation After Stereotactic Radiosurgery for Vestibular Schwannoma: A Retrospective Longitudinal Analysis.

BACKGROUND AND OBJECTIVES: Stereotactic radiosurgery (SRS) is a useful alternative for small- to medium-sized vestibular schwannoma. To evaluate whether biologically effective dose (BED Gy2.47 ), calculated for mean (BED Gy2.47 mean) and maximal (BED Gy2.47 max) cochlear dose, is relevant for hearing preservation. METHODS: This is a retrospective longitudinal single-center study. Were analyzed 213 patients with useful baseline hearing. Risk of hearing decline was assessed for Gardner-Robertson classes and pure tone average (PTA) loss. The mean follow-up period was 39 months (median 36, 6-84). RESULTS: Hearing decline (Gardner-Robertson class) 3 years after SRS was associated with higher cochlear BED Gy2.47 mean (odds ratio [OR] 1.39, P = .009). Moreover, BED Gy2.47 mean was more relevant as compared with BED Gy2.47 max (OR 1.13, P = .04). Risk of PTA loss (continuous outcome, follow-up minus baseline) was significantly corelated with BED Gy2.47 mean at 24 (beta coefficient 1.55, P = .002) and 36 (beta coefficient 2.01, P = .004) months after SRS. Risk of PTA loss (>20 dB vs ≤) was associated with higher BED Gy2.47 mean at 6 (OR 1.36, P = .002), 12 (OR 1.36, P = .007), and 36 (OR 1.37, P = .02) months. Risk of hearing decline at 36 months for the BED Gy2.47 mean of 7-8, 10, and 12 Gy 2.47 was 28%, 57%, and 85%, respectively. CONCLUSION: Cochlear BED Gy2.47 mean is relevant for hearing decline after SRS and more relevant as compared with BED Gy2.47 max. Three years after SRS, this was sustained for all hearing decline evaluation modalities. Our data suggest the BED Gy2.47 mean cut-off of ≤8 Gy 2.47 for better hearing preservation rates .

Clinically Driven Alpha/Beta Ratios for Melanoma Brain Metastases and Investigation of Biologically Effective Dose as a Predictor for Local Control After Radiosurgery: A Proof of Concept in a Retrospective Longitudinal Series of 274 Consecutive Lesions.

BACKGROUND AND OBJECTIVES: Brain metastases (BM) develop in nearly half of the patients with advanced melanoma. The aim of this retrospective historical cohort study was to analyze radiological response of melanoma BM to single-fraction Gamma Knife radiosurgery (GKRS), in relation to biologically effective dose (BED) for various alpha/beta ratios. METHODS: Included in the study were 274 lesions. Primary outcome was local control (LC). Mean marginal dose was 21.6 Gy (median 22, range 15-25). Biologically effective dose was calculated for an alpha/beta ratio of 3 (Gy 3 ), 5 (Gy 10 ), 10 (Gy 10 ), and 15 (Gy 15 ). RESULTS: Receiver operating characteristic value for LC and BED was 85% (most statistically significant odds ratio 1.14 for BED Gy 15 , P = .006), while for LC and physical dose was 79% ( P = .02). When comparing equality of 2 receiver operating characteristic areas, this was statistically significant ( P = .02 and .03). Fractional polynomial regression revealed BED (Gy 10 and Gy 15 ) as statistically significant ( P = .05) with BED of more than 63 Gy 10 or 49 Gy 15 as relevant, also for higher probability of quick decrease in volume first month after GKRS and lower probability of radiation necrosis. Shorter irradiation time was associated with better LC ( P = .001), particularly less than 40 minutes (LC below 90%, P = .05). CONCLUSION: BED Gy 10 and particularly Gy 15 were more statistically significant than physical dose for LC after GKRS for radioresistant melanoma BM. Irradiation time (per lesion) longer than 40 minutes was predictive for lower rates of LC. Such results need to be validated in larger cohorts.

Range shifter contribution to neutron exposure of patients undergoing proton pencil beam scanning.

BACKGROUND: Superficial targets require the use of the lowest energies within the available energy range in proton pencil-beam scanning (PBS) technique. However, the lower efficiency of the energy selection system at these energies and the requirement of a greater number of layers may represent disadvantages for this approach. The alternative is to use a range shifter (RS) at nozzle exit. However, one of the concerns of using this beamline element is that it becomes an additional source of neutrons that could irradiate organs situated far from the target. PURPOSE: The purpose of this study is to assess the increase in neutron dose due to the RS in proton PBS technique. Additionally, an analytical model for the neutron production is tested. METHODS: Two clinical plans, designed to achieve identical target coverage, were created for an anthropomorphic phantom. These plans consisted of a lateral field delivering an absorbed dose of 60 Gy (RBE) to the target. One of the plans employed the RS. The MCNP code was used to simulate the plans, evaluating the distribution of neutron dose equivalent (Hn) and the equivalent dose in organ. In the plan with the RS plan, neutron production from both the patient and the RS were assessed separately. Hn values were also fitted versus the distance to field edge using a Gaussian function. RESULTS: Hn per prescription dose, in the plan using the RS, ranged between 1.4 and 3.7 mSv/Gy at the field edge, whereas doses at 40 cm from the edge ranged from 9.9 to 32 µSv/Gy. These values are 1.2 to 10 times higher compared to those obtained without the RS. Both this factor and the contribution of neutrons originating from the RS increases with the distance from field edge. A triple-Gaussian function was able to reproduce the equivalent dose in organs within a factor of 2, although underestimating the values. CONCLUSIONS: The dose deposited in the patient by the neutrons originating from the RS predominantly affects areas away from the target (beyond approximately 25 cm from field edge), resulting in a neutron dose equivalent of the order of mSv. This indicates an overall low neutron contribution from the use of RS in PBS.

The Impact of Heterogeneous Cell Density in Hypoxic Tumors Treated with Radiotherapy.

Hypoxia is frequently found in solid tumors and is known to increase the resistance to several kinds of treatment modalities including radiation therapy. Besides, the treatment response is also largely determined by the total number of clonogenic cells, i.e., cells with unlimited proliferative capacity. Depending on the duration of hypoxia, the rate of proliferation and hence also the clonogen density could be expected to differ in hypoxic compartments. The combination at the microscale between heterogeneous tumor oxygenation and clonogen density could therefore be crucial with respect to the outcome of a radiotherapy treatment. In this study it was investigated the impact of heterogeneous clonogen density on the outcome of stereotactic radiotherapy treatments of hypoxic tumors. A recently developed three-dimensional model for tissue vasculature and oxygenation was used to create realistic in silico tumors with heterogeneous oxygenation. Stereotactic radiotherapy treatments were simulated, and cell survival was calculated on a voxel-level accounting for the oxygenation. For a tumor with a diameter of 1 cm and a baseline clonogenic density of 107/cm3 for the normoxic subvolume, when the relative density for the hypoxic cells drops by a factor of 10 the tumor control probability (TCP) decreases by about 10% when relatively small hypoxic volumes and few fractions are considered; longer treatments tend to level out the results. With increasing size of the hypoxic subvolume, the TCP decreased overall as expected, and the difference in TCP between a homogeneous and a heterogeneous distribution of cells increased. The results demonstrate a delicate interplay between the heterogeneous distribution of tumor oxygenation and clonogenic cells that could significantly impact on the treatment outcome of radiotherapy.

Hypoxia dose painting in SBRT - the virtual clinical trial approach.

BACKGROUND: Treating hypoxic tumours remains a challenge in radiotherapy as hypoxia leads to enhanced tumour aggressiveness and resistance to radiation. As escalating the doses is rarely feasible within the healthy tissue constraints, dose-painting strategies have been explored. Consensus about the best of care for hypoxic tumours has however not been reached because, among other reasons, the limits of current functional in-vivo imaging systems in resolving the details and dynamics of oxygen transport in tissue. Computational modelling of the tumour microenvironment enables the design and conduction of virtual clinical trials by providing relationships between biological features and treatment outcomes. This study presents a framework for assessing the therapeutic influence of the individual characteristics of the vasculature and the resulting oxygenation of hypoxic tumours in a virtual clinical trial on dose painting in stereotactic body radiotherapy (SBRT) circumventing the limitations of the imaging systems. MATERIAL AND METHODS: The homogeneous doses required to overcome hypoxia in simulated SBRT treatments of 1, 3 or 5 fractions were calculated for tumours with heterogeneous oxygenation derived from virtual vascular networks. The tumour control probability (TCP) was calculated for different scenarios for oxygenation dynamics resulting on cellular reoxygenation. RESULTS: A three-fractions SBRT treatment delivering 41.9 Gy (SD 2.8) and 26.5 Gy (SD 0.1) achieved only 21% (SD 12) and 48% (SD 17) control in the hypoxic and normoxic subvolumes, respectively whereas fast reoxygenation improved the control by 30% to 50%. TCP values for the individual tumours with similar characteristics, however, might differ substantially, highlighting the crucial role of the magnitude and time evolution of hypoxia at the microscale. CONCLUSION: The results show that local microvascular heterogeneities may affect the predicted outcome in the hypoxic core despite escalated doses, emphasizing the role of theoretical modelling in understanding of and accounting for the dominant factors of the tumour microenvironment.

Proton beam radiation therapy for vestibular schwannomas-tumor control and hearing preservation rates: a systematic review and meta-analysis.

OBJECTIVE: Proton beam therapy is considered, by some authors, as having the advantage of delivering dose distributions more conformal to target compared with stereotactic radiosurgery (SRS). Here, we performed a systematic review and meta-analysis of proton beam for VSs, evaluating tumor control and cranial nerve preservation rates, particularly with regard to facial and hearing preservation. METHODS: We reviewed, using the Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) articles published between 1968 and September 30, 2022. We retained 8 studies reporting 587 patients. RESULTS: Overall rate of tumor control (both stability and decrease in volume) was 95.4% (range 93.5-97.2%, p heterogeneity= 0.77, p<0.001). Overall rate of tumor progression was 4.6% (range 2.8-6.5%, p heterogeneity < 0.77, p<0.001). Overall rate of trigeminal nerve preservation (absence of numbness) was 95.6% (range 93.5-97.7%, I2 = 11.44%, p heterogeneity= 0.34, p<0.001). Overall rate of facial nerve preservation was 93.7% (range 89.6-97.7%, I2 = 76.27%, p heterogeneity<0.001, p<0.001). Overall rate of hearing preservation was 40.6% (range 29.4-51.8%, I2 = 43.36%, p heterogeneity= 0.1, p<0.001). CONCLUSION: Proton beam therapy for VSs achieves high tumor control rates, as high as 95.4%. Facial rate preservation overall rates are 93%, which is lower compared to the most SRS series. Compared with most currently reported SRS techniques, proton beam radiation therapy for VSs does not offer an advantage for facial and hearing preservation compared to most of the currently reported SRS series.

The Relevance of Biologically Effective Dose for Hearing Preservation After Stereotactic Radiosurgery for Vestibular Schwannomas: A Retrospective Longitudinal Study.

BACKGROUND: Stereotactic radiosurgery has become a common treatment approach for small-to-medium size vestibular schwannomas. OBJECTIVE: To evaluate relationship between time (beam-on and treatment) and risk of hearing decline after stereotactic radiosurgery for vestibular schwannomas in patients with Gardner-Robertson (GR) baseline classes I and II. METHODS: This retrospective longitudinal single-center study included 213 patients with GR I and II treated between June 2010 and December 2019. Risk of passing from GR classes I and II (coded 0) to other classes III, IV, and V (coded 1) and the increase in pure tone average (continuous outcome) were evaluated using a mixed-effect regression model. Biologically effective dose (BED) was further assessed for an alpha/beta ratio of 2.47 (Gy 2.47 ). RESULTS: Binary outcome analysis revealed sex, dose rate, integral dose, time [beam-on time odds ratio 1.03, P = .03, 95% CI 1.00-1.06; treatment time ( P = .02) and BED ( P = .001) as relevant. Fitted multivariable model included the sex, dose rate, and BED. Pure tone average analysis revealed age, integral dose received by tumor, isocenter number, time (beam-on time odds ratio 0.20, P = .001, 95% CI 0.083-0.33) and BED ( P = .005) as relevant. CONCLUSION: Our analysis showed that risk of hearing decline was associated with male sex, higher radiation dose rate (cutoff 2.5 Gy/minute), higher integral dose received by the tumor, higher beam-on time ≥20 minutes, and lower BED. A BED between 55 and 61 was considered as optimal for hearing preservation.

Modeling RBE with other quantities than LET significantly improves prediction of in vitro cell survival for proton therapy.

BACKGROUND: For proton therapy, a relative biological effectiveness (RBE) of 1.1 has broadly been applied clinically. However, as unexpected toxicities have been observed by the end of the proton tracks, variable RBE models have been proposed. Typically, the dose-averaged linear energy transfer (LETd ) has been used as an input variable for these models but the way the LETd was defined, calculated, or determined was not always consistent, potentially impacting the corresponding RBE value. PURPOSE: This study compares consistently calculated LETd with other quantities as input variables for a phenomenological RBE model and attempts to determine which quantity that can best predicts proton RBE. The comparison was performed within the frame of introducing a new model for the proton RBE. METHODS: High-throughput experimental setups of in vitro cell survival studies for proton RBE determination are simulated using the SHIELD-HIT12A Monte Carlo particle transport code. Together with LET, z ∗ 2 / ß 2 $z^{*2}/\beta ^2$ , here called effective Q (Qeff ), and Q are scored. Each quantity is calculated using the dose and track averaging methods, because the scoring includes all hadronic particles, all protons or only primaries. A phenomenological linear-quadratic-based RBE model is subsequently applied to the in vitro data with the various beam quality descriptors used as input variables and the goodness of fit is determined and compared using a bootstrapping approach. Both linear and nonlinear fit functions were tested. RESULTS: Versions of Qeff and Q outperform LET with a statistically significant margin, with the best nonlinear and linear fit having a relative root mean square error (RMSE) for RBE2Gy ± one standard error of 1.55 ± 0.04 (Qeff, t, primary ) and 2.84 ± 0.07 (Qeff, d, primary ), respectively. For comparison, the corresponding best nonlinear and linear fits for LETd, all protons had a relative RMSE of 2.07 ± 0.06 and 3.39 ± 0.08, respectively. Applying Welch's t-test for comparing the calculated RMSE of RBE2Gy resulted in two-tailed p-values of <0.002 for all Q and Qeff quantities compared to LETd, all protons . CONCLUSIONS: The study shows that Q or Qeff could be better RBE descriptors that dose averaged LET.

Photon and Proton Dose Painting Based on Oxygen Distribution - Feasibility Study and Tumour Control Probability Assessment.

Solid tumours may present hypoxic sub-regions of increased radioresistance. Hypoxia quantification requires of clinically implementable, non-invasive and reproducible techniques as positron emission tomography (PET). PET-based dose painting strategies aiming at targeting those sub-regions may be limited by the resolution gap between the PET imaging resolution and the smaller scale at which hypoxia occurs. The ultimate benefit of the usage of dose painting may be reached if the planned dose distribution can be performed and delivered consistently. This study aimed at assessing the feasibility of two PET-based dose painting strategies using two beam qualities (photon or proton beams) in terms of tumour control probability (TCP), accounting for underlying oxygen distribution at sub-millimetre scale.A tumour oxygenation model at submillimetre scale was created consisting of three regions with different oxygen partial pressure distributions, being hypoxia decreasing from core to periphery. A published relationship between uptake and oxygen partial pressure was used and a PET image of the tumour was simulated. The fundamental effects that limit the PET camera resolution were considered by processing the uptake distribution with a Gaussian 3D filter and re-binning to a PET image voxel size of 2 mm. Prescription doses to overcome tumour hypoxia were calculated based on the processed images, and planned using robust optimisation.Normal tissue complication probabilities and TCPs after the delivery of the planned doses were calculated for the nominal plan and the lowest bounds of the dose volume histograms resulting from the robust scenarios planned, taking into account the underlying oxygenation at submillimetre scale. Results were presented for the two beam qualities and the two dose painting strategies: by contours (DPBC) and by using a voxel grouping-based approach (DPBOX).In the studied case, DPBOX outperforms DPBC with respect to TCP regardless the beam quality, although both dose painting strategy plans demonstrated robust target coverage.

Perfusion-Limited Hypoxia Determines the Outcome of Radiation Therapy of Hypoxic Tumours.

Despite advancements in functional imaging, the resolution of modern techniques is still limited with respect to the tumour microenvironment. Radiotherapy strategies to counteract e.g., tumour hypoxia based on functional imaging therefore carry an inherent uncertainty that could compromise the outcome of the treatment. It was the aim of this study to investigate the impact of variations in the radiosensitivity of hypoxic tumours in small regions in comparison to the resolution of current imaging techniques on the probability of obtaining tumour control. A novel in silico model of three-dimensional tumour vasculature and oxygenation was used to model three tumours with different combinations of diffusion-limited, perfusion-limited and anaemic hypoxia. Specifically, cells in the transition region from a tumour core with diffusion-limited hypoxia to the well-oxygenated tumour rim were considered with respect to their differential radiosensitivity depending on the character of the hypoxia. The results showed that if the cells in the transition region were under perfusion-limited hypoxia, the tumour control probability was substantially lower in comparison to the case when the cells were anaemic (or under diffusion-limited hypoxia). This study therefore demonstrates the importance of differentiating between different forms of hypoxia on a scale currently unattainable to functional imaging techniques, lending support to the use and importance of radiobiological modelling of the cellular radiosensitivity and response at microscale.

Cell Type-Specific Patterns in the Accumulation of DNA Damage Following Multifractional Radiation Exposure.

Predicting the risk of second malignant neoplasms is complicated by uncertainties regarding the shape of the dose-response relationship at high doses. Limited understanding of the competitive relationship between cell killing and the accumulation of DNA lesions at high doses, as well as the effects of other modulatory factors unique to radiation exposure during radiotherapy, such as dose heterogeneity across normal tissue and dose fractionation, contribute to these uncertainties. The aim of this study was to analyze the impact of fractionated irradiations on two cell systems, focusing on the endpoints relevant for cancer induction. To simulate the heterogeneous dose distribution across normal tissue during radiotherapy, exponentially growing VH10 fibroblasts and AHH-1 lymphoblasts were irradiated with 9 and 12 fractions (VH10) and 10 fractions (AHH-1) at 0.25, 0.5, 1, or 2 Gy per fraction. The effects on cell growth, cell survival, radiosensitivity and the accumulation of residual DNA damage lesions were analyzed as functions of dose per fraction and the total absorbed dose. Residual γH2AX foci and other DNA damage markers (micronuclei, nuclear buds, and giant nuclei) were accumulated at high doses in both cell types, but in a cell type-dependent manner. The competitive relationship between cell killing and the accumulation of carcinogenic DNA damage following multifractional radiation exposure is cell type-specific.

Towards the virtual tumor for optimizing radiotherapy treatments of hypoxic tumors: A novel model of heterogeneous tissue vasculature and oxygenation.

PURPOSE: Tumor oxygenation is one of the key features influencing the response of cells to radiation and chemo therapies. This study presents a novel in silico tumor model simulating realistic 3D microvascular structures and related oxygenation maps, featuring regions with different levels and typologies of hypoxia (chronic, acute and anemic). Such model, if integrated into a treatment planning system, could allow evaluations and comparisons of various scenarios when deciding the therapy to administer. METHODS AND MATERIALS: Spherical tumors between 0.6 and 1.5 cm in diameter encompassed uniformly by vascular trees generated starting from pseudo-fractal principles were simulated with a voxel resolution of 10 µm. The approach ensures a continuous transition from a well-perfused rim to a core with poor vascularization. The oxygen diffusion equation in the tumor is solved by a finite difference method. Several quantities, such as the fractal dimension (FD), the microvascular density (MVD) and the hypoxic fraction (HF) were assessed and compared. RESULTS: Different tumors with various degrees of chronic hypoxia were simulated by varying the tumor size and the number of bifurcations in the vascular networks. The simulations showed that for the case of chronically hypoxic tumors, in well-oxygenated volumes FD = 2.53 ± 0.07, MVD = 3460 ± 2180 vessels/mm3 and HF = 4.0 ± 3.4%, while in hypoxic volumes FD = 2.34 ± 0.09, MVD = 365 ± 156 vessels/mm3, HF = 49.8 ± 18.3%. The superimposition of acute or anemic hypoxia accentuated the oxygen deprivation in the core of the volumes. CONCLUSIONS: Tumors varying in diameter and extension of their vasculature were simulated, showing features that define two distinctive subvolumes in terms of oxygenation. The model could be regarded as a testbed for simulations of key radiobiological features governing the tumor response to radio- and chemotherapy and thus for treatment outcome simulations.

CTV Delineation for High-Grade Gliomas: Is There Agreement With Tumor Cell Invasion Models?

Purpose: High-grade glioma (HGG) is a common form of malignant primary brain cancer with poor prognosis. The diffusive nature of HGGs implies that tumor cell invasion of normal tissue extends several centimeters away from the visible gross tumor volume (GTV). The standard methodology for clinical volume target (CTV) delineation is to apply a 2- to 3-cm margin around the GTV. However, tumor recurrence is extremely frequent. The purpose of this paper was to introduce a framework and computational model for the prediction of normal tissue HGG cell invasion and to investigate the agreement of the conventional CTV delineation with respect to the predicted tumor invasion. Methods and Materials: A model for HGG cell diffusion and proliferation was implemented and used to assess the tumor invasion patterns for 112 cases of HGGs. Normal brain structures and tissues as well as the GTVs visible on diagnostic images were delineated using automated methods. The volumes encompassed by different tumor cell concentration isolines calculated using the model for invasion were compared with the conventionally delineated CTVs, and the differences were analyzed. The 3-dimensional-Hausdorff distance between the CTV and the volumes encompassed by various isolines was also calculated. Results: In 50% of cases, the CTV failed to encompass regions containing tumor cell concentrations of 614 cells/mm³ or greater. In 84% of cases, the lowest cell concentration completely encompassed by the CTV was ≥1 cell/mm³. In the remaining 16%, the CTV overextended into normal tissue. The Hausdorff distance was on average comparable to the CTV margin. Conclusions: The standard methodology for CTV delineation appears to be inconsistent with HGG invasion patterns in terms of size and shape. Tumor invasion modeling could therefore be useful in assisting in the CTV delineation for HGGs.

Determining Out-of-Field Doses and Second Cancer Risk From Proton Therapy in Young Patients-An Overview.

Proton therapy has the potential to provide survival and tumor control outcomes comparable and frequently superior to photon therapy. This has led to a significant concern in the medical physics community on the risk for the induction of second cancers in all patients and especially in younger patients, as they are considered more radiosensitive than adults and have an even longer expected lifetime after treatment. Thus, our purpose is to present an overview of the research carried out on the evaluation of out-of-field doses linked to second cancer induction and the prediction of this risk. Most investigations consisted of Monte Carlo simulations in passive beam facilities for clinical scenarios. These works established that equivalent doses in organs could be up to 200 mSv or 900 mSv for a brain or a craniospinal treatment, respectively. The major contribution to this dose comes from the secondary neutrons produced in the beam line elements. Few works focused on scanned-beam facilities, but available data show that, for these facilities, equivalent doses could be between 2 and 50 times lower. Patient age is a relevant factor in the dose level, especially for younger patients (by means of the size of the body) and, in addition, in the predicted risk by models (due to the age dependence of the radiosensitivity). For risks, the sex of the patient also plays an important role, as female patients show higher sensitivity to radiation. Thus, predicted risks of craniospinal irradiation can range from 8% for a 15-year-old male patient to 58% for a 2-year-old female patient, using a risk model from a radiological protection field. These values must be taken with caution due to uncertainties in risk models, and then dosimetric evaluation of stray radiation becomes mandatory in order to complement epidemiological studies and be able to model appropriate dose-response functions for this dose range. In this sense, analytical models represent a useful tool and some models have been implemented to be used for young patients. Research carried out so far confirmed that proton beam therapy reduces the out-of-field doses and second cancer risk. However, further investigations may be required in scanned-beam delivery systems.

Radiobiological Evaluation of Combined Gamma Knife Radiosurgery and Hyperthermia for Pediatric Neuro-Oncology.

Combining radiotherapy (RT) with hyperthermia (HT) has been proven effective in the treatment of a wide range of tumours, but the combination of externally delivered, focused heat and stereotactic radiosurgery has never been investigated. We explore the potential of such treatment enhancement via radiobiological modelling, specifically via the linear-quadratic (LQ) model adapted to thermoradiotherapy through modulating the radiosensitivity of temperature-dependent parameters. We extend this well-established model by incorporating oxygenation effects. To illustrate the methodology, we present a clinically relevant application in pediatric oncology, which is novel in two ways. First, it deals with medulloblastoma, the most common malignant brain tumour in children, a type of brain tumour not previously reported in the literature of thermoradiotherapy studies. Second, it makes use of the Gamma Knife for the radiotherapy part, thereby being the first of its kind in this context. Quantitative metrics like the biologically effective dose (BED) and the tumour control probability (TCP) are used to assess the efficacy of the combined plan.