Performance of dose calculation algorithms from three generations in lung SBRT: comparison with full Monte Carlo-based dose distributions.

The accuracy of dose calculation is a key challenge in stereotactic body radiotherapy (SBRT) of the lung. We have benchmarked three photon beam dose calculation algorithms--pencil beam convolution (PBC), anisotropic analytical algorithm (AAA), and Acuros XB (AXB)--implemented in a commercial treatment planning system (TPS), Varian Eclipse. Dose distributions from full Monte Carlo (MC) simulations were regarded as a reference. In the first stage, for four patients with central lung tumors, treatment plans using 3D conformal radiotherapy (CRT) technique applying 6 MV photon beams were made using the AXB algorithm, with planning criteria according to the Nordic SBRT study group. The plans were recalculated (with same number of monitor units (MUs) and identical field settings) using BEAMnrc and DOSXYZnrc MC codes. The MC-calculated dose distributions were compared to corresponding AXB-calculated dose distributions to assess the accuracy of the AXB algorithm, to which then other TPS algorithms were compared. In the second stage, treatment plans were made for ten patients with 3D CRT technique using both the PBC algorithm and the AAA. The plans were recalculated (with same number of MUs and identical field settings) with the AXB algorithm, then compared to original plans. Throughout the study, the comparisons were made as a function of the size of the planning target volume (PTV), using various dose-volume histogram (DVH) and other parameters to quantitatively assess the plan quality. In the first stage also, 3D gamma analyses with threshold criteria 3%/3mm and 2%/2 mm were applied. The AXB-calculated dose distributions showed relatively high level of agreement in the light of 3D gamma analysis and DVH comparison against the full MC simulation, especially with large PTVs, but, with smaller PTVs, larger discrepancies were found. Gamma agreement index (GAI) values between 95.5% and 99.6% for all the plans with the threshold criteria 3%/3 mm were achieved, but 2%/2 mm threshold criteria showed larger discrepancies. The TPS algorithm comparison results showed large dose discrepancies in the PTV mean dose (D50%), nearly 60%, for the PBC algorithm, and differences of nearly 20% for the AAA, occurring also in the small PTV size range. This work suggests the application of independent plan verification, when the AAA or the AXB algorithm are utilized in lung SBRT having PTVs smaller than 20-25 cc. The calculated data from this study can be used in converting the SBRT protocols based on type 'a' and/or type 'b' algorithms for the most recent generation type 'c' algorithms, such as the AXB algorithm.

Evaluation of dose-response models and parameters predicting radiation induced pneumonitis using clinical data from breast cancer radiotherapy.

The purpose of this work is to evaluate the predictive strength of the relative seriality, parallel and LKB normal tissue complication probability (NTCP) models regarding the incidence of radiation pneumonitis, in a large group of patients following breast cancer radiotherapy, and furthermore, to illustrate statistical methods for examining whether certain published radiobiological parameters are compatible with a clinical treatment methodology and patient group characteristics. The study is based on 150 consecutive patients who received radiation therapy for breast cancer. For each patient, the 3D dose distribution delivered to lung and the clinical treatment outcome were available. Clinical symptoms and radiological findings, along with a patient questionnaire, were used to assess the manifestation of radiation-induced complications. Using this material, different methods of estimating the likelihood of radiation effects were evaluated. This was attempted by analysing patient data based on their full dose distributions and associating the calculated complication rates with the clinical follow-up records. Additionally, the need for an update of the criteria that are being used in the current clinical practice was also examined. The patient material was selected without any conscious bias regarding the radiotherapy treatment technique used. The treatment data of each patient were applied to the relative seriality, LKB and parallel NTCP models, using published parameter sets. Of the 150 patients, 15 experienced radiation-induced pneumonitis (grade 2) according to the radiation pneumonitis scoring criteria used. Of the NTCP models examined, the relative seriality model was able to predict the incidence of radiation pneumonitis with acceptable accuracy, although radiation pneumonitis was developed by only a few patients. In the case of modern breast radiotherapy, radiobiological modelling appears to be very sensitive to model and parameter selection giving clinically acceptable results in certain cases selectively (relative seriality model with Seppenwoolde et al and Gagliardi et al parameter sets). The use of published parameters should be considered as safe only after their examination using local clinical data. The variation of inter-patient radiosensitivity seems to play a significant role in the prediction of such low incidence rate complications. Scoring grades were combined to give stronger evidence of radiation pneumonitis since their differences could not be strictly associated with dose. This obviously reveals a weakness of the scoring related to this endpoint, and implies that the probability of radiation pneumonitis induction may be too low to be statistically analysed with high accuracy, at least with the latest advances of dose delivery in breast radiotherapy.

Effects of positioning uncertainty and breathing on dose delivery and radiation pneumonitis prediction in breast cancer.

The quality of the radiation therapy delivered in the treatment of breast cancer is susceptible to setup errors and organ motion uncertainties. For 60 breast cancer patients (24 resected with negative node involvement, 13 resected with positive node involvement and 23 ablated) who were treated with three different irradiation techniques. these uncertainties are simulated. The delivered dose distributions in the lung were recalculated taking positioning uncertainty and breathing effects into account. In this way the real dose distributions delivered to the patients are more closely determined. The positioning uncertainties in the anteroposterior (AP) and the craniocaudal (CC) directions are approximated by Gaussian distributions based on the fact that setup errors are random. Breathing is assumed to have a linear behavior because of the chest wall movement during expiration and inspiration. The combined frequency distribution of the positioning and breathing distributions is obtained by convolution. By integrating the convolved distribution over a number of intervals, the positions and the weights of the fields that simulate the original 'effective fields' are calculated. Opposed tangential fields are simulated by a set of 5 pairs of fields in the AP direction and 3 such sets in the CC direction. Opposed AP + PA fields are simulated by a set of 3 pairs of fields in the AP direction and 3 such sets in the CC direction. Single frontal fields are simulated by a set of 5 fields. In radiotherapy for breast cancer, the lung is often partly within the irradiated volume even though it is a sensitive organ at risk. The influence of the deviation in the dose delivered by the original and the adjusted treatment plans on the clinical outcome is estimated by using the relative seriality model and the biologically effective uniform dose concept. Radiation pneumonitis is used as the clinical endpoint for lung complications. The adjusted treatment plans show larger lung complication probabilities than the original plans. This means that the true expected complications are often underestimated in clinical practice. The lung density variation during breathing is calculated from the maximal change in average density during tidal breathing. The change in density in the lung due to breathing is shown to have almost no influence on the dose distribution in the lung. The proposed treatment-plan adjustments taking positioning uncertainty and breathing effects into account indicate significant deviations in the dose delivery and the predicted lung complications.