Radiobiological Meta-Analysis of the Response of Prostate Cancer to Different Fractionations: Evaluation of the Linear-Quadratic Response at Large Doses and the Effect of Risk and ADT.

The purpose of this work was to investigate the response of prostate cancer to different radiotherapy schedules, including hypofractionation, to evaluate potential departures from the linear-quadratic (LQ) response, to obtain the best-fitting parameters for low-(LR), intermediate-(IR), and high-risk (HR) prostate cancer and to investigate the effect of ADT on the radiobiological response. We constructed a dataset of the dose-response containing 87 entries/16,536 patients (35/5181 LR, 32/8146 IR, 20/3209 HR), with doses per fraction ranging from 1.8 to 10 Gy. These data were fit to tumour control probability models based on the LQ model, linear-quadratic-linear (LQL) model, and a modification of the LQ (LQmod) model accounting for increasing radiosensitivity at large doses. Fits were performed with the maximum likelihood expectation methodology, and the Akaike information criterion (AIC) was used to compare the models. The AIC showed that the LQ model was superior to the LQL and LQmod models for all risks, except for IR, where the LQL model outperformed the other models. The analysis showed a low α/ß for all risks: 2.0 Gy for LR (95% confidence interval: 1.7-2.3), 3.4 Gy for IR (3.0-4.0), and 2.8 Gy for HR (1.4-4.2). The best fits did not show proliferation for LR and showed moderate proliferation for IR/HR. The addition of ADT was consistent with a suppression of proliferation. In conclusion, the LQ model described the response of prostate cancer better than the alternative models. Only for IR, the LQL model outperformed the LQ model, pointing out a possible saturation of radiation damage with increasing dose. This study confirmed a low α/ß for all risks.

Dosimetric Evaluation of Cardiac Structures on Left Breast Cancer Radiotherapy: Impact of Movement, Dose Calculation Algorithm and Treatment Technique.

Background: Breast cancer is the most frequently diagnosed and leading cause of cancer-related deaths among females. The treatment of breast cancer with radiotherapy, albeit effective, has been shown to be toxic to the heart, resulting in an elevated risk of cardiovascular disease and associated fatalities. Methods: In this study, we evaluated the impact of respiratory movement, treatment plans and dose calculation algorithm on the dose delivered to the heart and its substructures during left breast radiotherapy over a cohort of 10 patients. We did this through three image sets, four different treatment plans and the employment of three algorithms on the same treatment plan. The dose parameters were then employed to estimate the impact on the 9-year excess cumulative risk for acute cardiac events by applying the model proposed by Darby. Results: The left ventricle was the structure most irradiated. Due to the lack of four-dimensional computed tomography (4DCT), we used a set of images called phase-average CT that correspond to the average of the images from the respiratory cycle (exhale, exhale 50%, inhale, inhale 50%). When considering these images, nearly 10% of the heart received more than 5 Gy and doses were on average 27% higher when compared to free breathing images. Deep inspiration breath-hold plans reduced cardiac dose for nine out of 10 patients and reduced mean heart dose in about 50% when compared to reference plans. We also found that the implementation of deep inspiration breath-hold would reduce the relative lifetime risk of ischemic heart disease to 10%, in comparison to 21% from the reference plan. Conclusion: Our findings illustrate the importance of a more accurate determination of the dose and its consideration in cardiologists' consultation, a factor often overlooked during clinical examination. They also motivate the evaluation of the dose to the heart substructures to derive new heart dose constraints, and a more mindful and individualized clinical practice depending on the treatment employed.

Correction factors for A1SL ionization chamber dosimetry in TomoTherapy: machine-specific, plan-class, and clinical fields.

PURPOSE: Recently, an international working group on nonstandard fields presented a new formalism for ionization chamber reference dosimetry of small and nonstandard fields [Alfonso et al., Med. Phys. 35, 5179-5186 (2008)] which has been adopted by AAPM TG-148. This work presents an experimental determination of the correction factors for reference dosimetry with an Exradin A1SL thimble ionization chamber in a TomoTherapy unit, focusing on: (i) machine-specific reference field, (ii) plan-class-specific reference field, and (iii) two clinical treatments. METHODS: Ionization chamber measurements were performed in the TomoTherapy unit for intermediate (machine-specific and plan-class-specific) calibration fields, based on the reference conditions defined by AAPM TG-148, and two clinical treatments (lung and head-and-neck). Alanine reference dosimetry was employed to determine absorbed dose to water at the point of interest for the fields under investigation. The corresponding chamber correction factors were calculated from alanine to ionization chamber measurements ratios. RESULTS: Two different methods of determining the beam quality correction factor k(Q,Q(0) ) for the A1SL ionization chamber in this TomoTherapy unit, where reference conditions for conventional beam quality determination cannot be met, result in consistent values. The observed values of overall correction factors obtained for intermediate and clinical fields are consistently around 0.98 with a typical expanded relative uncertainty of 2% (k = 2), which when considered make such correction factors compatible with unity. However, all of them are systematically lower than unity, which is shown to be significant when a hypothesis test assuming a t-student distribution is performed (p=1.8×10(-2)). Correction factors k(Q(clin),Q(pcsr) ) (f(clin),f(pcsr) ) and k(Q(clin),Q(msr) ) (f(clin),f(msr) ), which are needed for the computation of field factors for relative dosimetry of clinical beams, have been found to be very close to unity for two clinical treatments. CONCLUSIONS: The results indicate that the helical field deliveries in this study (including two clinical fields) do not introduce changes on the ion chamber correction factors for dosimetry. For those two specific clinical cases, ratios of chamber readings accurately represent field output factors. The values observed here for intermediate calibration fields are in agreement with previously published data based on alanine dosimetry but differ from values recently reported obtained via radiochromic dosimetry.

Evaluation of indirect damage and damage saturation effects in dose-response curves of hypofractionated radiotherapy of early-stage NSCLC and brain metastases.

BACKGROUND AND PURPOSE: To investigate the possible contribution of indirect damage and damage saturation to tumour control obtained with SBRT/SRS treatments for early-stage NSCLC and brain metastases. METHODS AND MATERIALS: We have constructed a dataset of early-stage NSCLC and brain metastases dose-response. These data were fitted to models based on the linear-quadratic (LQ), the linear-quadratic-linear (LQL), and phenomenological modifications of the LQ-model to account for indirect cell damage. We use the Akaike-Information-Criterion formalism to compare performance, and studied the stability of the results with changes in fitting parameters and perturbations on dose/TCP values. RESULTS: In NSCLC, a modified LQ-model with a beta-term increasing with dose yields the best-fits for α/ß = 10 Gy. Only the inclusion of very fast accelerated proliferation or low α/ß values can eliminate such superiority. In brain, the LQL model yields the best-fits, and the ranking is not affected by variations of fitting parameters or dose/TCP perturbations. CONCLUSIONS: For α/ß = 10 Gy, a modified LQ-model with a beta-term increasing with dose provides better fits to NSCLC dose-response curves. For brain metastases, the LQL provides the best fit. This might be interpreted as a hint of indirect damage in NSCLC, and damage saturation in brain metastases. The results for NSCLC are strongly dependent on the value of α/ß and may require further investigation, while those for brain seem to be clearly significant. Our results can assist in the design of improved radiotherapy for NSCLC and brain metastases, aiming at avoiding over/under-treatment.

Development of a Compartmental Pharmacokinetic Model for Molecular Radiotherapy with 131I-CLR1404.

Pharmacokinetic modeling of the radiopharmaceuticals used in molecular radiotherapy is an important step towards accurate radiation dosimetry of such therapies. In this paper, we present a pharmacokinetic model for CLR1404, a phospholipid ether analog that, labeled with 124I/131I, has emerged as a promising theranostic agent. We follow a systematic approach for the model construction based on a decoupling process applied to previously published experimental data, and using the goodness-of-fit, Sobol's sensitivity analysis, and the Akaike Information Criterion to construct the optimal form of the model, investigate potential simplifications, and study factor prioritization. This methodology was applied to previously published experimental human time-activity curves for 9 organs. The resulting model consists of 17 compartments involved in the CLR1404 metabolism. Activity dynamics in most tissues are well described by a blood contribution plus a two-compartment system, describing fast and slow uptakes. The model can fit both clinical and pre-clinical kinetic data of 124I/131I. In addition, we have investigated how simple fits (exponential and biexponential) differ from the complete model. Such fits, despite providing a less accurate description of time-activity curves, may be a viable alternative when limited data is available in a practical case.

Quantification of internal dosimetry in PET patients II: Individualized Monte Carlo-based dosimetry for [18F]fluorocholine PET.

PURPOSE: To obtain individualized internal doses with a Monte Carlo (MC) method in patients undergoing diagnostic [18F]FCH-PET studies and to compare such doses with the MIRD method calculations. METHODS: A patient cohort of 17 males were imaged after intravenous administration of a mean [18F]FCH activity of 244.3 MBq. The resulting PET/CT images were processed in order to generate individualized input source and geometry files for dose computation with the MC tool GATE. The resulting dose estimates were studied and compared to the MIRD method with two different computational phantoms. Mass correction of the S-factors was applied when possible. Potential sources of uncertainty were closely examined: the effect of partial body images, urinary bladder emptying, and biokinetic modeling. RESULTS: Large differences in doses between our methodology and the MIRD method were found, generally in the range ±25%, and up to ±120% for some cases. The mass scaling showed improvements, especially for non-walled and high-uptake tissues. Simulations of the urinary bladder emptying showed negligible effects on doses to other organs, with the exception of the prostate. Dosimetry based on partial PET/CT images (excluding the legs) resulted in an overestimation of mean doses to bone, skin, and remaining tissues, and minor differences in other organs/tissues. Estimated uncertainties associated with the biokinetics of FCH introduce variations of cumulated activities in the range of ±10% in the high-uptake organs. CONCLUSIONS: The MC methodology allows for a higher degree of dosimetry individualization than the MIRD methodology, which in some cases leads to important differences in dose values. Dosimetry of FCH-PET based on a single partial PET study seems viable due to the particular biokinetics of FCH, even though some correction factors may need to be applied to estimate mean skin/bone doses.

A kinetic model of continuous radiation damage to populations of cells: comparison to the LQ model and application to molecular radiotherapy.

The linear-quadratic (LQ) model to describe the survival of irradiated cells may be the most frequently used biomathematical model in radiotherapy. There has been an intense debate on the mechanistic origin of the LQ model. An interesting approach is that of obtaining LQ-like behavior from kinetic models, systems of differential equations that model the induction and repair of damage. Development of such kinetic models is particularly interesting for application to continuous dose rate therapies, such as molecular radiotherapy or brachytherapy. In this work, we present a simple kinetic model that describes the kinetics of populations of tumor cells, rather than lethal/sub-lethal lesions, which may be especially useful for application to continuous dose rate therapies, as in molecular radiotherapy. The multi-compartment model consists of a set of three differential equations. The model incorporates in an easy way different cross-interacting compartments of cells forming a tumor, and may be of especial interest for studying dynamics of treated tumors. In the fast dose delivery limit, the model can be analytically solved, obtaining a simple closed-form expression. Fitting of several surviving curves with both this solution and the LQ model shows that they produce similar fits, despite being functionally different. We have also investigated the operation of the model in the continuous dose rate scenario, firstly by fitting pre-clinical data of tumor response to 131I-CLR1404 therapy, and secondly by showing how damage repair and proliferation rates can cause a treatment to achieve control or not. Kinetic models like the one presented in this work may be of special interest when modeling response to molecular radiotherapy.

Simulation of hypoxia PET-tracer uptake in tumours: Dependence of clinical uptake-values on transport parameters and arterial input function.

Poor radiotherapy outcome is in many cases related to hypoxia, due to the increased radioresistance of hypoxic tumour cells. Positron emission tomography may be used to non-invasively assess the oxygenation status of the tumour using hypoxia-specific radiotracers. Quantification and interpretation of these images remains challenging, since radiotracer binding and oxygen tension are not uniquely related. Computer simulation is a useful tool to improve the understanding of tracer dynamics and its relation to clinical uptake parameters currently used to quantify hypoxia. In this study, a model for simulating oxygen and radiotracer distribution in tumours was implemented to analyse the impact of physiological transport parameters and of the arterial input function (AIF) on: oxygenation histograms, time-activity curves, tracer binding and clinical uptake-values (tissue-to-blood ratio, TBR, and a composed hypoxia-perfusion metric, FHP). Results were obtained for parallel and orthogonal vessel architectures and for vascular fractions (VFs) of 1% and 3%. The most sensitive parameters were the AIF and the maximum binding rate (Kmax). TBR allowed discriminating VF for different AIF, and FHP for different Kmax, but neither TBR nor FHP were unbiased in all cases. Biases may especially occur in the comparison of TBR- or FHP-values between different tumours, where the relation between measured and actual AIF may vary. Thus, these parameters represent only surrogates rather than absolute measurements of hypoxia in tumours.

Quantification of internal dosimetry in PET patients: individualized Monte Carlo vs generic phantom-based calculations.

PURPOSE: The purpose of this work is to calculate individualized dose distributions in patients undergoing 18 F-FDG PET/CT studies through a methodology based on full Monte Carlo (MC) simulations and PET/CT patient images, and to compare such values with those obtained by employing nonindividualized phantom-based methods. METHODS: We developed a MC-based methodology for individualized internal dose calculations, which relies on CT images (for organ segmentation and dose deposition), PET images (for organ segmentation and distributions of activities), and a biokinetic model (which works with information provided by PET and CT images) to obtain cumulated activities. The software vGATE version 8.1. was employed to carry out the Monte Carlo calculations. We also calculated deposited doses with nonindividualized phantom-based methods (Cristy-Eckerman, Stabin, and ICRP-133). RESULTS: Median MC-calculated dose/activity values are within 0.01-0.03 mGy/MBq for most organs, with higher doses delivered especially to the bladder wall, major vessels, and brain (medians of 0.058, 0.060, 0.066 mGy/MBq, respectively). Comparison with values obtained with nonindividualized phantom-based methods has shown important differences in many cases (ranging from -80% to + 260%). These differences are significant (p < 0.05) for several organs/tissues, namely, remaining tissues, adrenals, bladder wall, bones, upper large intestine, heart, pancreas, skin, and stomach wall. CONCLUSIONS: The methodology presented in this work is a viable and useful method to calculate internal dose distributions in patients undergoing medical procedures involving radiopharmaceuticals, individually, with higher accuracy than phantom-based methods, fulfilling the guidelines provided by the European Council directive 2013/59/Euratom.

eIMRT: a web platform for the verification and optimization of radiation treatment plans.

The eIMRT platform is a remote distributed computing tool that provides users with Internet access to three different services: Monte Carlo optimization of treatment plans, CRT & IMRT treatment optimization, and a database of relevant radiation treatments/clinical cases. These services are accessible through a user-friendly and platform independent web page. Its flexible and scalable design focuses on providing the final users with services rather than a collection of software pieces. All input and output data (CT, contours, treatment plans and dose distributions) are handled using the DICOM format. The design, implementation, and support of the verification and optimization algorithms are hidden to the user. This allows a unified, robust handling of the software and hardware that enables these computation-intensive services. The eIMRT platform is currently hosted by the Galician Supercomputing Center (CESGA) and may be accessible upon request (there is a demo version at http://eimrt.cesga.es:8080/eIMRT2/demo; request access in http://eimrt.cesga.es/signup.html). This paper describes all aspects of the eIMRT algorithms in depth, its user interface, and its services. Due to the flexible design of the platform, it has numerous applications including the intercenter comparison of treatment planning, the quality assurance of radiation treatments, the design and implementation of new approaches to certain types of treatments, and the sharing of information on radiation treatment techniques. In addition, the web platform and software tools developed for treatment verification and optimization have a modular design that allows the user to extend them with new algorithms. This software is not a commercial product. It is the result of the collaborative effort of different public research institutions and is planned to be distributed as an open source project. In this way, it will be available to any user; new releases will be generated with the new implemented codes or upgrades.

A Model of Indirect Cell Death Caused by Tumor Vascular Damage after High-Dose Radiotherapy.

There is increasing evidence that high doses of radiotherapy, like those delivered in stereotactic body radiotherapy (SBRT), trigger indirect mechanisms of cell death. Such effect seems to be two-fold. High doses may trigger an immune response and may cause vascular damage, leading to cell starvation and death. Development of mathematical response models, including indirect death, may help clinicians to design SBRT optimal schedules. Despite increasing experimental literature on indirect tumor cell death caused by vascular damage, efforts on modeling this effect have been limited. In this work, we present a biomathematical model of this effect. In our model, tumor oxygenation is obtained by solving the reaction-diffusion equation; radiotherapy kills tumor cells according to the linear-quadratic model, and also endothelial cells (EC), which can trigger loss of functionality of capillaries. Capillary death will affect tumor oxygenation, driving nearby tumor cells into severe hypoxia. Capillaries can recover functionality due to EC proliferation. Tumor cells entering a predetermined severe hypoxia status die according to a hypoxia-death model. This model fits recently published experimental data showing the effect of vascular damage on surviving fractions. It fits surviving fraction curves and qualitatively reproduces experimental values of percentages of functional capillaries 48 hours postirradiation, and hypoxic cells pre- and 48 hours postirradiation. This model is useful for exploring aspects of tumor and EC response to radiotherapy and constitutes a stepping stone toward modeling indirect tumor cell death caused by vascular damage and accounting for this effect during SBRT planning. SIGNIFICANCE: A novel biomathematical model of indirect tumor cell death caused by vascular radiation damage could potentially help clinicians interpret experimental data and design better radiotherapy schedules.

Impact of different biologically-adapted radiotherapy strategies on tumor control evaluated with a tumor response model.

Motivated by the capabilities of modern radiotherapy techniques and by the recent developments of functional imaging techniques, dose painting by numbers (DPBN) was proposed to treat tumors with heterogeneous biological characteristics. This work studies different DPBN optimization techniques for virtual head and neck tumors assessing tumor response in terms of cell survival and tumor control probability with a previously published tumor response model (TRM). Uniform doses of 2 Gy are redistributed according to the microscopic oxygen distribution and the density distribution of tumor cells in four virtual tumors with different biological characteristics. In addition, two different optimization objective functions are investigated, which: i) minimize tumor cell survival (OFsurv) or; ii) maximize the homogeneity of the density of surviving tumor cells (OFstd). Several adaptive schemes, ranging from single to daily dose optimization, are studied and the treatment response is compared to that of the uniform dose. The results show that the benefit of DPBN treatments depends on the tumor reoxygenation capability, which strongly differed among the set of virtual tumors investigated. The difference between daily (fraction by fraction) and three weekly optimizations (at the beginning of weeks 1, 3 and 4) was found to be small, and higher benefit was observed for the treatments optimized using OFsurv. This in silico study corroborates the hypothesis that DPBN may be beneficial for treatments of tumors which show reoxygenation during treatment, and that a few optimizations may be sufficient to achieve this therapeutic benefit.

Structural and functional identification of vasculogenic mimicry in vitro.

Vasculogenic mimicry (VM) describes a process by which cancer cells establish an alternative perfusion pathway in an endothelial cell-free manner. Despite its strong correlation with reduced patient survival, controversy still surrounds the existence of an in vitro model of VM. Furthermore, many studies that claim to demonstrate VM fail to provide solid evidence of true hollow channels, raising concerns as to whether actual VM is actually being examined. Herein, we provide a standardized in vitro assay that recreates the formation of functional hollow channels using ovarian cancer cell lines, cancer spheres and primary cultures derived from ovarian cancer ascites. X-ray microtomography 3D-reconstruction, fluorescence confocal microscopy and dye microinjection conclusively confirm the existence of functional glycoprotein-rich lined tubular structures in vitro and demonstrate that many of structures reported in the literature may not represent VM. This assay may be useful to design and test future VM-blocking anticancer therapies.

Modelling radiation-induced cell death and tumour re-oxygenation: local versus global and instant versus delayed cell death.

The resistance of hypoxic cells to radiation, due to the oxygen dependence of radiosensitivity, is well known and must be taken into account to accurately calculate the radiation induced cell death. A proper modelling of the response of tumours to radiation requires deriving the distribution of oxygen at a microscopic scale. This usually involves solving the reaction-diffusion equation in tumour voxels using a vascularization distribution model. Moreover, re-oxygenation arises during the course of radiotherapy, one reason being the increase of available oxygen caused by cell killing, which can turn hypoxic tumours into oxic. In this work we study the effect of cell death kinetics in tumour oxygenation modelling, analysing how it affects the timing of re-oxygenation, surviving fraction and tumour control. Two models of cell death are compared, an instantaneous cell killing, mimicking early apoptosis, and a delayed cell death scenario in which cells can die shortly after being damaged, as well as long after irradiation. For each of these scenarios, the decrease in oxygen consumption due to cell death can be computed globally (macroscopic voxel average) or locally (microscopic). A re-oxygenation model already used in the literature, the so called full re-oxygenation, is also considered. The impact of cell death kinetics and re-oxygenation on tumour responses is illustrated for two radiotherapy fractionation schemes: a conventional schedule, and a hypofractionated treatment. The results show large differences in the doses needed to achieve 50% tumour control for the investigated cell death models. Moreover, the models affect the tumour responses differently depending on the treatment schedule. This corroborates the complex nature of re-oxygenation, showing the need to take into account the kinetics of cell death in radiation response models.

Correction factors for ionization chamber dosimetry in CyberKnife: machine-specific, plan-class, and clinical fields.

PURPOSE: The aim of this work is the application of the formalism for ionization chamber reference dosimetry of small and nonstandard fields [R. Alfonso, P. Andreo, R. Capote, M. S. Huq, W. Kilby, P. Kjäll, T. R. Mackie, H. Palmans, K. Rosser, J. Seuntjens, W. Ullrich, and S. Vatnitsky, "A new formalism for reference dosimetry of small and nonstandard fields," Med. Phys. 35, 5179-5186 (2008)] to the CyberKnife robotic radiosurgery system. Correction factors for intermediate calibration fields, a machine-specific reference field (msr) and two plan-class specific reference fields (pcsr), have been studied. Furthermore, the applicability of the new formalism to clinical dosimetry has been analyzed through the investigation of two clinical treatments. METHODS: PTW31014 and Scanditronix-Wellhofer CC13 ionization chamber measurements were performed for the fields under investigation. Absorbed dose to water was determined using alanine reference dosimetry, and experimental correction factors were calculated from alanine to ionization chamber readings ratios. In addition, correction factors were calculated for the intermediate calibration fields and one of the clinical treatment fields using the Monte Carlo method and these were compared with the experimental values. RESULTS: Overall correction factors deviating from unity by approximately 2% were obtained from both measurements and simulations, with values below and above unity for the studied intermediate calibration fields and clinical fields for the ionization chambers under consideration. Monte Carlo simulations yielded correction factors comparable with those obtained from measurements for the machine-specific reference field, although differences from 1% to 3.3% were observed between measured and calculated correction factors for the composite intermediate calibration fields. Dose distribution inhomogeneities are thought to be responsible for such discrepancies. CONCLUSIONS: The differences found between overall correction factors associated with the proposed intermediate calibration fields and the clinical fields under investigation show that depending on the clinical field and the detector used, either a machine-specific reference field or a plan-class specific reference field is more representative for the clinical field. Given the experimental and numerical uncertainties and the small number of clinical fields considered in this study the significance of these observations is limited and it remains unclear for the CyberKnife if there would be a significant gain in using a pcsr field rather than a msr field as reference field for relative dosimetry.

A two-dimensional liquid-filled ionization chamber array prototype for small-field verification: characterization and first clinical tests.

In this work we present the design, characterization and first clinical tests of an in-house developed two-dimensional liquid-filled ionization chamber prototype for the verification of small radiotherapy fields and treatments containing such small fields as in radiosurgery, which consists of 2 mm × 2 mm pixels arranged on a 16×8 rectangular grid. The ionization medium is isooctane. The characterization of the device included the study of depth, field-size and dose-rate dependences, which are sufficiently moderate for a good operation at therapy radiation levels. However, the detector presents an important anisotropic response, up to ≃ 12% for front versus near-lateral incidence, which can impact the verification of full treatments with different incidences. In such a case, an anisotropy correction factor can be applied. Output factors of small square fields measured with the device show a small systematic over-response, less than 1%, when compared to unshielded diode measurements. An IMRT radiosurgery treatment has been acquired with the liquid-filled ionization chamber device and compared with film dosimetry by using the gamma method, showing good agreement: over 99% passing rates for 1.2% and 1.2 mm for an incidence-per-incidence analysis; 100% passing rates for tolerances 1.8% and 1.8 mm when the whole treatment is analysed and the anisotropy correction factor is applied. The point dose verification for each incidence of the treatment performed with the liquid-filled ionization chamber agrees within 1% with a CC01 ionization chamber. This prototype has shown the utility of this kind of technology for the verification of small fields/treatments. Currently, a larger device covering a 5 cm × 5 cm area is under development.