Desenvolvimento e avaliação de um sistema de cálculo de dose independente para controle de qualidade de IMRT do tipo jaws-only

Dissertação (mestrado)—Universidade de Brasília, Faculdade Gama, Programa de Pós-Graduação em Engenharia Biomédica, 2015.

On The Consistency Of Monte Carlo Track Structure Dna Damage Simulations.

Monte Carlo track structures (MCTS) simulations have been recognized as useful tools for radiobiological modeling. However, the authors noticed several issues regarding the consistency of reported data. Therefore, in this work, they analyze the impact of various user defined parameters on simulated direct DNA damage yields. In addition, they draw attention to discrepancies in published literature in DNA strand break (SB) yields and selected methodologies. The MCTS code Geant4-DNA was used to compare radial dose profiles in a nanometer-scale region of interest (ROI) for photon sources of varying sizes and energies. Then, electron tracks of 0.28 keV-220 keV were superimposed on a geometric DNA model composed of 2.7 × 10(6) nucleosomes, and SBs were simulated according to four definitions based on energy deposits or energy transfers in DNA strand targets compared to a threshold energy ETH. The SB frequencies and complexities in nucleosomes as a function of incident electron energies were obtained. SBs were classified into higher order clusters such as single and double strand breaks (SSBs and DSBs) based on inter-SB distances and on the number of affected strands. Comparisons of different nonuniform dose distributions lacking charged particle equilibrium may lead to erroneous conclusions regarding the effect of energy on relative biological effectiveness. The energy transfer-based SB definitions give similar SB yields as the one based on energy deposit when ETH ≈ 10.79 eV, but deviate significantly for higher ETH values. Between 30 and 40 nucleosomes/Gy show at least one SB in the ROI. The number of nucleosomes that present a complex damage pattern of more than 2 SBs and the degree of complexity of the damage in these nucleosomes diminish as the incident electron energy increases. DNA damage classification into SSB and DSB is highly dependent on the definitions of these higher order structures and their implementations. The authors' show that, for the four studied models, different yields are expected by up to 54% for SSBs and by up to 32% for DSBs, as a function of the incident electrons energy and of the models being compared. MCTS simulations allow to compare direct DNA damage types and complexities induced by ionizing radiation. However, simulation results depend to a large degree on user-defined parameters, definitions, and algorithms such as: DNA model, dose distribution, SB definition, and the DNA damage clustering algorithm. These interdependencies should be well controlled during the simulations and explicitly reported when comparing results to experiments or calculations.

Evaluation of IMRT and VMAT techniques with and without flattening filter using pareto optimal fronts

Mestrado em Radiações Aplicadas às Tecnologias da Saúde.

Comparison of different breast planning techniques and algorithms for radiation therapy treatment

This work aims at investigating the impact of treating breast cancer using different radiation therapy (RT) techniques – forwardly-planned intensity-modulated, f-IMRT, inversely-planned IMRT and dynamic conformal arc (DCART) RT – and their effects on the whole-breast irradiation and in the undesirable irradiation of the surrounding healthy tissues. Two algorithms of iPlan BrainLAB treatment planning system were compared: Pencil Beam Convolution (PBC) and commercial Monte Carlo (iMC). Seven left-sided breast patients submitted to breast-conserving surgery were enrolled in the study. For each patient, four RT techniques – f-IMRT, IMRT using 2-fields and 5-fields (IMRT2 and IMRT5, respectively) and DCART – were applied. The dose distributions in the planned target volume (PTV) and the dose to the organs at risk (OAR) were compared analyzing dose–volume histograms; further statistical analysis was performed using IBM SPSS v20 software. For PBC, all techniques provided adequate coverage of the PTV. However, statistically significant dose differences were observed between the techniques, in the PTV, OAR and also in the pattern of dose distribution spreading into normal tissues. IMRT5 and DCART spread low doses into greater volumes of normal tissue, right breast, right lung and heart than tangential techniques. However, IMRT5 plans improved distributions for the PTV, exhibiting better conformity and homogeneity in target and reduced high dose percentages in ipsilateral OAR. DCART did not present advantages over any of the techniques investigated. Differences were also found comparing the calculation algorithms: PBC estimated higher doses for the PTV, ipsilateral lung and heart than the iMC algorithm predicted.

Comparison of the pencil beam convolution algorithm and the analytical anisotropic algorithm in breast tumors

The calculation of the dose is one of the key steps in radiotherapy planning1-5. This calculation should be as accurate as possible, and over the years it became feasible through the implementation of new algorithms to calculate the dose on the treatment planning systems applied in radiotherapy. When a breast tumour is irradiated, it is fundamental a precise dose distribution to ensure the planning target volume (PTV) coverage and prevent skin complications. Some investigations, using breast cases, showed that the pencil beam convolution algorithm (PBC) overestimates the dose in the PTV and in the proximal region of the ipsilateral lung. However, underestimates the dose in the distal region of the ipsilateral lung, when compared with analytical anisotropic algorithm (AAA). With this study we aim to compare the performance in breast tumors of the PBC and AAA algorithms.

Comparison of inverse planning systems based on biological or physical factors: Pinnacle®, Corvus® and Monaco®

Radiotherapy (RT) is one of the most important approaches in the treatment of cancer and its performance can be improved in three different ways: through the optimization of the dose distribution, by the use of different irradiation techniques or through the study of radiobiological initiatives. The first is purely physical because is related to the physical dose distributiuon. The others are purely radiobiological because they increase the differential effect between the tumour and the health tissues. The Treatment Planning Systems (TPS) are used in RT to create dose distributions with the purpose to maximize the tumoral control and minimize the complications in the healthy tissues. The inverse planning uses dose optimization techniques that satisfy the criteria specified by the user, regarding the target and the organs at risk (OAR’s). The dose optimization is possible through the analysis of dose-volume histograms (DVH) and with the use of computed tomography, magnetic resonance and other digital image techniques.

Dosimetric effect of tissue heterogeneity for 125I prostate implants

Aim - To use Monte Carlo (MC) together with voxel phantoms to analyze the tissue heterogeneity effect in the dose distributions and equivalent uniform dose (EUD) for (125)I prostate implants. Background - Dose distribution calculations in low dose-rate brachytherapy are based on the dose deposition around a single source in a water phantom. This formalism does not take into account tissue heterogeneities, interseed attenuation, or finite patient dimensions effects. Tissue composition is especially important due to the photoelectric effect. Materials and Methods - The computed tomographies (CT) of two patients with prostate cancer were used to create voxel phantoms for the MC simulations. An elemental composition and density were assigned to each structure. Densities of the prostate, vesicles, rectum and bladder were determined through the CT electronic densities of 100 patients. The same simulations were performed considering the same phantom as pure water. Results were compared via dose-volume histograms and EUD for the prostate and rectum. Results - The mean absorbed doses presented deviations of 3.3-4.0% for the prostate and of 2.3-4.9% for the rectum, when comparing calculations in water with calculations in the heterogeneous phantom. In the calculations in water, the prostate D 90 was overestimated by 2.8-3.9% and the rectum D 0.1cc resulted in dose differences of 6-8%. The EUD resulted in an overestimation of 3.5-3.7% for the prostate and of 7.7-8.3% for the rectum. Conclusions - The deposited dose was consistently overestimated for the simulation in water. In order to increase the accuracy in the determination of dose distributions, especially around the rectum, the introduction of the model-based algorithms is recommended.

Hand exposure in diagnostic nuclear medicine with 18F- and 99mTc-labelled radiopharmaceuticals - Results of the ORAMED project

Workers performing preparation and administration of radiopharmaceuticals in NM departments are likely to receive high local skin doses to the hands which may even surpass the dose limit of 500 mSv whenever radiation protection standards are insufficient. A large measurement campaign was organised within the framework of the ORAMED project to determine the dose distribution across the hands received during preparation and administration of 18F- and 99mTc-labelled radiopharmaceuticals. The final data, collected over almost 3 years, include 641 measurements from 96 workers in 30 NM departments from 6 European countries. Results have provided levels of reference doses for the considered standard NM diagnostic procedures (mean maximum normalised skin dose of 230 μSv/GBq, 430 μSv/GBq, 930 μSv/GBq and 1200 μSv/GBq for the administration of 99mTc, preparation of 99mTc, administration of 18F and preparation of 18F, respectively). Finger dose was analysed as a function of the potential parameters of influence showing that shielding is the most efficient means of radiation protection to reduce skin dose. An appropriate method for routine monitoring of the extremities is also proposed: the base of the index finger of the non-dominant hand is a suitable position to place the ring dosemeter, with its sensitive part oriented towards the palm side; its reading may be multiplied by a factor of 6 to estimate the maximum local skin dose. Finally, results were compared to earlier published data, which correspond mostly to individual works with a reduced number of workers and measurements.

Recurrence pattern after [(18)F]fluoroethyltyrosine-positron emission tomography-guided radiotherapy for high-grade glioma: a prospective study.

PURPOSE: To assess the failure pattern observed after (18)F fluoroethyltyrosine (FET) planning after chemo- and radiotherapy (RT) for high-grade glioma. METHODS: All patients underwent prospectively RT planning using morphological gross tumour volumes (GTVs) and biological tumour volumes (BTVs). The post-treatment recurrence tumour volumes (RTVs) of 10 patients were transferred on their CT planning. First, failure patterns were defined in terms of percentage of RTV located outside the GTV and BTV. Second, the location of the RTV with respect to the delivered dose distribution was assessed using the RTV's DVHs. Recurrences with >95% of their volume within 95% isodose line were considered as central recurrences. Finally, the relationship between survival and GTV/BTV mismatches was assessed. RESULTS: The median percentages of RTV outside the GTV and BTV were 41.8% (range, 10.5-92.4) and 62.8% (range, 34.2-81.1), respectively. The majority of recurrences (90%) were centrally located. Using a composite target volume planning formalism, the degree of GTV and BTV mismatch did not correlate with survivorship. CONCLUSIONS: The observed failure pattern after FET-PET planning and chemo-RT is primarily central. The target mismatch-survival data suggest that using FET-PET planning may counteract the possibility of BTV-related progression, which may have a detrimental effect on survival.

Symmetry-breaking in cumulative measures of shapes of polymer models.

Using numerical simulations we investigate shapes of random equilateral open and closed chains, one of the simplest models of freely fluctuating polymers in a solution. We are interested in the 3D density distribution of the modeled polymers where the polymers have been aligned with respect to their three principal axes of inertia. This type of approach was pioneered by Theodorou and Suter in 1985. While individual configurations of the modeled polymers are almost always nonsymmetric, the approach of Theodorou and Suter results in cumulative shapes that are highly symmetric. By taking advantage of asymmetries within the individual configurations, we modify the procedure of aligning independent configurations in a way that shows their asymmetry. This approach reveals, for example, that the 3D density distribution for linear polymers has a bean shape predicted theoretically by Kuhn. The symmetry-breaking approach reveals complementary information to the traditional, symmetrical, 3D density distributions originally introduced by Theodorou and Suter.

Dummy run and conformity indices in the ongoing EORTC low-grade glioma trial 22033-26033: First evaluation of quality of radiotherapy planning.

PURPOSE: Early assessment of radiotherapy (RT) quality in the ongoing EORTC trial comparing primary temozolomide versus RT in low-grade gliomas. MATERIALS AND METHODS: RT plans provided for dummy cases were evaluated and compared against expert plans. We analysed: (1) tumour and organs-at-risk delineation, (2) geometric and dosimetric characteristics, (3) planning parameters, compliance with dose prescription and Dmax for OAR (4) indices: RTOG conformity index (CI), coverage factor (CF), tissue protection factor (PF); conformity number (CN = PF x CF); dose homogeneity in PTV (U). RESULTS: Forty-one RT plans were evaluated. Only two (5%) centres were requested to repeat CTV-PTV delineations. Three (7%) plans had a significant under-dosage and dose homogeneity in one deviated > 10%. Dose distribution was good with mean values of 1.5, 1, 0.68, and 0.68 (ideal values = 1) for CI, CF, PF, and CN, respectively. CI and CN strongly correlated with PF and they correlated with PTV. Planning with more beams seems to increase PTV(Dmin), improving CF. U correlated with PTV(Dmax). CONCLUSION: Preliminary results of the dummy run procedure indicate that most centres conformed to protocol requirements. To quantify plan quality we recommend systematic calculation of U and either CI or CN, both of which measure the amount of irradiated normal brain tissue.

Dosimetric comparison of different treatment modalities for stereotactic radiosurgery of arteriovenous malformations and acoustic neuromas.

PURPOSE: We investigated the influence of beam modulation on treatment planning by comparing four available stereotactic radiosurgery (SRS) modalities: Gamma-Knife-Perfexion, Novalis-Tx Dynamic-Conformal-Arc (DCA) and Dynamic-Multileaf-Collimation-Intensity-Modulated-radiotherapy (DMLC-IMRT), and Cyberknife. MATERIAL AND METHODS: Patients with arteriovenous malformation (n = 10) or acoustic neuromas (n = 5) were planned with different treatment modalities. Paddick conformity index (CI), dose heterogeneity (DH), gradient index (GI) and beam-on time were used as dosimetric indices. RESULTS: Gamma-Knife-Perfexion can achieve high degree of conformity (CI = 0.77 ± 0.04) with limited low-doses (GI = 2.59 ± 0.10) surrounding the inhomogeneous dose distribution (D(H) = 0.84 ± 0.05) at the cost of treatment time (68.1 min ± 27.5). Novalis-Tx-DCA improved this inhomogeneity (D(H) = 0.30 ± 0.03) and treatment time (16.8 min ± 2.2) at the cost of conformity (CI = 0.66 ± 0.04) and Novalis-TX-DMLC-IMRT improved the DCA CI (CI = 0.68 ± 0.04) and inhomogeneity (D(H) = 0.18 ± 0.05) at the cost of low-doses (GI = 3.94 ± 0.92) and treatment time (21.7 min ± 3.4) (p<0.01). Cyberknife achieved comparable conformity (CI = 0.77 ± 0.06) at the cost of low-doses (GI = 3.48 ± 0.47) surrounding the homogeneous (D(H) = 0.22 ± 0.02) dose distribution and treatment time (28.4min±8.1) (p<0.01). CONCLUSIONS: Gamma-Knife-Perfexion will comply with all SRS constraints (high conformity while minimizing low-dose spread). Multiple focal entries (Gamma-Knife-Perfexion and Cyberknife) will achieve better conformity than High-Definition-MLC of Novalis-Tx at the cost of treatment time. Non-isocentric beams (Cyberknife) or IMRT-beams (Novalis-Tx-DMLC-IMRT) will spread more low-dose than multiple isocenters (Gamma-Knife-Perfexion) or dynamic arcs (Novalis-Tx-DCA). Inverse planning and modulated fluences (Novalis-Tx-DMLC-IMRT and CyberKnife) will deliver the most homogeneous treatment. Furthermore, Linac-based systems (Novalis and Cyberknife) can perform image verification at the time of treatment delivery.

A Monte Carlo-based procedure for independent monitor unit calculation in IMRT treatment plans.

Intensity-modulated radiotherapy (IMRT) treatment plan verification by comparison with measured data requires having access to the linear accelerator and is time consuming. In this paper, we propose a method for monitor unit (MU) calculation and plan comparison for step and shoot IMRT based on the Monte Carlo code EGSnrc/BEAMnrc. The beamlets of an IMRT treatment plan are individually simulated using Monte Carlo and converted into absorbed dose to water per MU. The dose of the whole treatment can be expressed through a linear matrix equation of the MU and dose per MU of every beamlet. Due to the positivity of the absorbed dose and MU values, this equation is solved for the MU values using a non-negative least-squares fit optimization algorithm (NNLS). The Monte Carlo plan is formed by multiplying the Monte Carlo absorbed dose to water per MU with the Monte Carlo/NNLS MU. Several treatment plan localizations calculated with a commercial treatment planning system (TPS) are compared with the proposed method for validation. The Monte Carlo/NNLS MUs are close to the ones calculated by the TPS and lead to a treatment dose distribution which is clinically equivalent to the one calculated by the TPS. This procedure can be used as an IMRT QA and further development could allow this technique to be used for other radiotherapy techniques like tomotherapy or volumetric modulated arc therapy.

The use of different types of thermoluminescent dosimeters to measure extremity doses in nuclear medicine

Depth-dose curves in LiF detectors of different effective thicknesses, together with their responses, were calculated for typical nuclear medicine radiation fields with 99mTc, 18F and 90Y sources. Responses were analysed in function of the radionuclide, detector effective thickness and irradiation geometry. On the other hand the results of the nuclear medicine measurement campaign of the ORAMED project were presented focussing on the dose distribution across the hand and on the appropriate position to wear the dosimeter.According to the results, thin LiF detectors provide better responses in all cases. Its use is essential for 18F, since thick dosimeters can underestimate Hp(0.07) up to a 50% because of the very inhomogeneous dose deposition on the active layer. The preliminary results of the measurement campaign showed that the index tip of the non-dominant hand is usually the most exposed position among the 22 monitored positions. It was also found that, in average, wrist dosimeters are likely to underestimate the maximum skin dose by a factor of the order of 20. This factor is reduced to around 6 for a ring dosimeter worn on the base of the index of the non-dominant hand. Thus, for typical nuclear medicine procedures, the base of the index of the non-dominant hand is recommended as the best monitoring option.

Modélisation radiobiologique pour la planification des traitements en radiothérapie à partir de données d’imagerie spécifiques aux patients

Un modèle de croissance et de réponse à la radiothérapie pour le glioblastome multiforme (GBM) basé le formalisme du modèle de prolifération-invasion (PI) et du modèle linéaire-quadratique a été développé et implémenté. La géométrie spécifique au patient est considérée en modélisant, d'une part, les voies d'invasion possibles des GBM avec l'imagerie du tenseur de diffusion (DTI) et, d'autre part, les barrières à la propagation à partir des images anatomiques disponibles. La distribution de dose réelle reçue par un patient donné est appliquée telle quelle dans les simulations, en respectant l'horaire de traitement. Les paramètres libres du modèle (taux de prolifération, coefficient de diffusion, paramètres radiobiologiques) sont choisis aléatoirement à partir de distributions de valeurs plausibles. Un total de 400 ensembles de valeurs pour les paramètres libres sont ainsi choisis pour tous les patients, et une simulation de la croissance et de la réponse au traitement est effectuée pour chaque patient et chaque ensemble de paramètres. Un critère de récidive est appliqué sur les résultats de chaque simulation pour identifier un lieu probable de récidive (SPR). La superposition de tous les SPR obtenus pour un patient donné permet de définir la probabilité d'occurrence (OP). Il est démontré qu'il existe des valeurs de OP élevées pour tous les patients, impliquant que les résultats du modèle PI ne sont pas très sensibles aux valeurs des paramètres utilisés. Il est également démontré comment le formalisme développé dans cet ouvrage pourrait permettre de définir un volume cible personnalisé pour les traitements de radiothérapie du GBM.