Calculation of dose kernels for radiotherapy applications using the GEANT 4 Monte Carlo toolkit

Dose kernels may be used to calculate dose distributions in radiotherapy (as described by Ahnesjo et al., 1999). Their calculation requires use of Monte Carlo methods, usually by forcing interactions to occur at a point. The Geant4 Monte Carlo toolkit provides a capability to force interactions to occur in a particular volume. We have modified this capability and created a Geant4 application to calculate dose kernels in cartesian, cylindrical, and spherical scoring systems. The simulation considers monoenergetic photons incident at the origin of a 3 m x 3 x 9 3 m water volume. Photons interact via compton, photo-electric, pair production, and rayleigh scattering. By default, Geant4 models photon interactions by sampling a physical interaction length (PIL) for each process. The process returning the smallest PIL is then considered to occur. In order to force the interaction to occur within a given length, L_FIL, we scale each PIL according to the formula: PIL_forced = L_FIL 9 (1 - exp(-PIL/PILo)) where PILo is a constant. This ensures that the process occurs within L_FIL, whilst correctly modelling the relative probability of each process. Dose kernels were produced for an incident photon energy of 0.1, 1.0, and 10.0 MeV. In order to benchmark the code, dose kernels were also calculated using the EGSnrc Edknrc user code. Identical scoring systems were used; namely, the collapsed cone approach of the Edknrc code. Relative dose difference images were then produced. Preliminary results demonstrate the ability of the Geant4 application to reproduce the shape of the dose kernels; median relative dose differences of 12.6, 5.75, and 12.6 % were found for an incident photon energy of 0.1, 1.0, and 10.0 MeV respectively.

Hydrogen bonding in two ammonium salts of 5-sulfosalicylic acid : ammonium 3-carboxy-4-hydroxybenzenesulfonate monohydrate and triammonium 3-carboxy-4-hydroxybenzenesulfonate 3-carboxylato-4-hydroxybenzenesulfonate

The structures of two ammonium salts of 3-carboxy-4-hydroxybenzenesulfonic acid (5-sulfosalicylic acid, 5-SSA) have been determined at 200 K. In the 1:1 hydrated salt, ammonium 3-carboxy-4-hydroxybenzenesulfonate monohydrate, NH4+·C7H5O6S-·H2O, (I), the 5-SSA- monoanions give two types of head-to-tail laterally linked cyclic hydrogen-bonding associations, both with graph-set R44(20). The first involves both carboxylic acid O-HOwater and water O-HOsulfonate hydrogen bonds at one end, and ammonium N-HOsulfonate and N-HOcarboxy hydrogen bonds at the other. The second association is centrosymmetric, with end linkages through water O-HOsulfonate hydrogen bonds. These conjoined units form stacks down c and are extended into a three-dimensional framework structure through N-HO and water O-HO hydrogen bonds to sulfonate O-atom acceptors. Anhydrous triammonium 3-carboxy-4-hydroxybenzenesulfonate 3-carboxylato-4-hydroxybenzenesulfonate, 3NH4+·C7H4O6S2-·C7H5O6S-, (II), is unusual, having both dianionic 5-SSA2- and monoanionic 5-SSA- species. These are linked by a carboxylic acid O-HO hydrogen bond and, together with the three ammonium cations (two on general sites and the third comprising two independent half-cations lying on crystallographic twofold rotation axes), give a pseudo-centrosymmetric asymmetric unit. Cation-anion hydrogen bonding within this layered unit involves a cyclic R33(8) association which, together with extensive peripheral N-HO hydrogen bonding involving both sulfonate and carboxy/carboxylate acceptors, gives a three-dimensional framework structure. This work further demonstrates the utility of the 5-SSA- monoanion for the generation of stable hydrogen-bonded crystalline materials, and provides the structure of a dianionic 5-SSA2- species of which there are only a few examples in the crystallographic literature.

Dosimetric effects of a high-density spinal implant

In this study, a treatment plan for a spinal lesion, with all beams transmitted though a titanium vertebral reconstruction implant, was used to investigate the potential effect of a high-density implant on a three-dimensional dose distribution for a radiotherapy treatment. The BEAMnrc/DOSXYZnrc and MCDTK Monte Carlo codes were used to simulate the treatment using both a simplified, recltilinear model and a detailed model incorporating the full complexity of the patient anatomy and treatment plan. The resulting Monte Carlo dose distributions showed that the commercial treatment planning system failed to accurately predict both the depletion of dose downstream of the implant and the increase in scattered dose adjacent to the implant. Overall, the dosimetric effect of the implant was underestimated by the commercial treatment planning system and overestimated by the simplified Monte Carlo model. The value of performing detailed Monte Carlo calculations, using the full patient and treatment geometry, was demonstrated.

The use of electronic portal imaging to verify patient position during intensity-modulated radiotherapy delivered by the dynamic MLC technique

Purpose: The precise shape of the three-dimensional dose distributions created by intensity-modulated radiotherapy means that the verification of patient position and setup is crucial to the outcome of the treatment. In this paper, we investigate and compare the use of two different image calibration procedures that allow extraction of patient anatomy from measured electronic portal images of intensity-modulated treatment beams. Methods and Materials: Electronic portal images of the intensity-modulated treatment beam delivered using the dynamic multileaf collimator technique were acquired. The images were formed by measuring a series of frames or segments throughout the delivery of the beams. The frames were then summed to produce an integrated portal image of the delivered beam. Two different methods for calibrating the integrated image were investigated with the aim of removing the intensity modulations of the beam. The first involved a simple point-by-point division of the integrated image by a single calibration image of the intensity-modulated beam delivered to a homogeneous polymethyl methacrylate (PMMA) phantom. The second calibration method is known as the quadratic calibration method and required a series of calibration images of the intensity-modulated beam delivered to different thicknesses of homogeneous PMMA blocks. Measurements were made using two different detector systems: a Varian amorphous silicon flat-panel imager and a Theraview camera-based system. The methods were tested first using a contrast phantom before images were acquired of intensity-modulated radiotherapy treatment delivered to the prostate and pelvic nodes of cancer patients at the Royal Marsden Hospital. Results: The results indicate that the calibration methods can be used to remove the intensity modulations of the beam, making it possible to see the outlines of bony anatomy that could be used for patient position verification. This was shown for both posterior and lateral delivered fields. Conclusions: Very little difference between the two calibration methods was observed, so the simpler division method, requiring only the single extra calibration measurement and much simpler computation, was the favored method. This new method could provide a complementary tool to existing position verification methods, and it has the advantage that it is completely passive, requiring no further dose to the patient and using only the treatment fields.

The development of Monte Carlo techniques for the verification of radiotherapy treatments

Introduction: Recent advances in the planning and delivery of radiotherapy treatments have resulted in improvements in the accuracy and precision with which therapeutic radiation can be administered. As the complexity of the treatments increases it becomes more difficult to predict the dose distribution in the patient accurately. Monte Carlo (MC) methods have the potential to improve the accuracy of the dose calculations and are increasingly being recognised as the ‘gold standard’ for predicting dose deposition in the patient [1]. This project has three main aims: 1. To develop tools that enable the transfer of treatment plan information from the treatment planning system (TPS) to a MC dose calculation engine. 2. To develop tools for comparing the 3D dose distributions calculated by the TPS and the MC dose engine. 3. To investigate the radiobiological significance of any errors between the TPS patient dose distribution and the MC dose distribution in terms of Tumour Control Probability (TCP) and Normal Tissue Complication Probabilities (NTCP). The work presented here addresses the first two aims. Methods: (1a) Plan Importing: A database of commissioned accelerator models (Elekta Precise and Varian 2100CD) has been developed for treatment simulations in the MC system (EGSnrc/BEAMnrc). Beam descriptions can be exported from the TPS using the widespread DICOM framework, and the resultant files are parsed with the assistance of a software library (PixelMed Java DICOM Toolkit). The information in these files (such as the monitor units, the jaw positions and gantry orientation) is used to construct a plan-specific accelerator model which allows an accurate simulation of the patient treatment field. (1b) Dose Simulation: The calculation of a dose distribution requires patient CT images which are prepared for the MC simulation using a tool (CTCREATE) packaged with the system. Beam simulation results are converted to absolute dose per- MU using calibration factors recorded during the commissioning process and treatment simulation. These distributions are combined according to the MU meter settings stored in the exported plan to produce an accurate description of the prescribed dose to the patient. (2) Dose Comparison: TPS dose calculations can be obtained using either a DICOM export or by direct retrieval of binary dose files from the file system. Dose difference, gamma evaluation and normalised dose difference algorithms [2] were employed for the comparison of the TPS dose distribution and the MC dose distribution. These implementations are spatial resolution independent and able to interpolate for comparisons. Results and Discussion: The tools successfully produced Monte Carlo input files for a variety of plans exported from the Eclipse (Varian Medical Systems) and Pinnacle (Philips Medical Systems) planning systems: ranging in complexity from a single uniform square field to a five-field step and shoot IMRT treatment. The simulation of collimated beams has been verified geometrically, and validation of dose distributions in a simple body phantom (QUASAR) will follow. The developed dose comparison algorithms have also been tested with controlled dose distribution changes. Conclusion: The capability of the developed code to independently process treatment plans has been demonstrated. A number of limitations exist: only static fields are currently supported (dynamic wedges and dynamic IMRT will require further development), and the process has not been tested for planning systems other than Eclipse and Pinnacle. The tools will be used to independently assess the accuracy of the current treatment planning system dose calculation algorithms for complex treatment deliveries such as IMRT in treatment sites where patient inhomogeneities are expected to be significant. Acknowledgements: Computational resources and services used in this work were provided by the HPC and Research Support Group, Queensland University of Technology, Brisbane, Australia. Pinnacle dose parsing made possible with the help of Paul Reich, North Coast Cancer Institute, North Coast, New South Wales.

Developing an adaptive radiotherapy technique for virally mediated head and neck cancer

Purpose Virally mediated head and neck cancers (VMHNC) often present with nodal involvement and are highly radioresponsive, meaning that treatment plan adaptation during radiotherapy (RT) in a subset of patients is required. We sought to determine potential risk profiles and a corresponding adaptive treatment strategy for these patients. Methodology 121 patients with virally mediated, node positive nasopharyngeal (Epstein Barr Virus positive) or oropharyngeal (Human Papillomavirus positive) cancers, receiving curative intent RT were reviewed. The type, frequency and timing of adaptive interventions, including source-to-skin distance (SSD) corrections, re-scanning and re-planning, were evaluated. Patients were reviewed based on the maximum size of the dominant node to assess the need for plan adaptation. Results Forty-six patients (38%) required plan adaptation during treatment. The median fraction at which the adaptive intervention occurred was 26 for SSD corrections and 22 for re-planning CTs. A trend toward 3 risk profile groupings was discovered: 1) Low risk with minimal need (< 10%) for adaptive intervention (dominant pre-treatment nodal size of ≤ 35 mm), 2) Intermediate risk with possible need (< 20%) for adaptive intervention (dominant pre-treatment nodal size of 36 mm – 45 mm) and 3) High-risk with increased likelihood (> 50%) for adaptive intervention (dominant pre-treatment nodal size of ≥ 46 mm). Conclusion In this study, patients with VMHNC and a maximum dominant nodal size of > 46 mm were identified at a higher risk of requiring re-planning during a course of definitive RT. Findings will be tested in a future prospective adaptive RT study.

Molecular cocrystals of 5-(4-bromophenyl)-1,3,4-thiadiazol-2-amine: hydrogen bonding in the structures of the 1:1 adduct with 2-(naphthalen-2-yloxy)acetic acid and the salt with 3,5-dinitrobenzoic acid

The structures of the anhydrous products from the interaction of 2-amino-5-(4-bromophenyl)-1,3,4-thiadiazole with (2-naphthoxy)acetic acid, the 1:1 adduct C8H6BrN3S . C12H10O3 (I) and 3,5-dinitrobenzoic acid, the salt C8H7BrN3S+ C7H3N2O6- (II) have been determined. In the adduct (I), a heterodimer is formed through a cyclic hydrogen-bonding motif [graph set R2/2(8)], involving carboxylic acid O-H...N(hetero)and amine N-H...O(carboxyl) interactions. The heterodimers are essentially planar with a thiadiazole to naphthyl ring dihedral angle of 15.9(2)deg. and the intramolecular thiadiazole to phenyl ring angle of 4.7(2)deg. An amine N-H...N(hetero) hydrogen bond between the heterodimers generates a one-dimensional chain structure extending down [001]. Also present are weak benzene-benzene and naphthalene-naphthalene pi-pi stacking interactions down the b axis [minimum ring centroid separation, 3.936(3) Ang.]. With the salt (II), the cation-anion association is also through a cyclic R2/2(8) motif but involving duplex N-H...O(carboxyl) hydrogen bonds, giving a heterodimer which is close to planar [dihedral angles between the thiadiazole ring and the two benzene rings, 5.00(16)deg. (intra) and 7.23(15)deg. (inter)]. A secondary centrosymmetric cyclic N-H...O(carboxyl) hydrogen-bonding association involving the second amino H-atom generates a heterotetramer. Also present in the crystal are weak pi-pi i-\p interactions between thiadiazolium rings [minimum ring centroid separation, 3.936(3)Ang.], as well as a short Br...O(nitro) interaction [3.314(4)Ang.]. The two structures reported here now provide a total of three crystallographically characterized examples of co-crystalline products from the interaction of 2-amino-5-(4-bromophenyl)-1,3,4-thiadiazole with carboxylic acids, of which only one involves proton-transfer.

A simple method for fabricating 3-D multilayered composite scaffolds

One limitation of electrospinning stems from the charge build-up that occurs during processing, preventing further fibre deposition and limiting the scaffold overall thickness and hence their end-use in tissue engineering applications targeting large tissue defect repair. To overcome this, we have developed a technique in which thermally induced phase separation (TIPS) and electrospinning are combined. Thick three-dimensional, multilayered composite scaffolds were produced by simply stacking individual polycaprolactone (PCL) microfibrous electrospun discs into a cylindrical holder that was filled with a 3% poly(lactic-co-glycolic acid) (PLGA) solution in dimethylsulfoxide (a good solvent for PLGA but a poor one for PCL). The construct was quenched in liquid nitrogen and the solvent removed by leaching out in cold water. This technique enables the fabrication of scaffolds composed principally of electrospun membranes that have no limit to their thickness. The mechanical properties of these scaffolds were assessed under both quasi-static and dynamic conditions. The multilayered composite scaffolds had similar compressive properties to 5% PCL scaffolds fabricated solely by the TIPS methodology. However, tensile tests demonstrated that the multilayered construct outperformed a scaffold made purely by TIPS, highlighting the contribution of the electrospun component of the composite scaffold to enhancing the overall mechanical property slate. Cell studies revealed cell infiltration principally from the scaffold edges under static seeding conditions. This fabrication methodology permits the rapid construction of thick, strong scaffolds from a range of biodegradable polymers often used in tissue engineering, and will be particularly useful when large dimension electrospun scaffolds are required.

Anodic formation of a thick three-dimensional nanoporous WO3 film and its photocatalytic property

We report an efficient solar-light-driven photocatalyst based on three-dimensional nanoporous tungsten trioxide (WO3) films. These films are obtained by anodizing W foils in fluoride-containing electrolytes at room temperature and under low applied voltages with an efficient growth rate of 2 μm h− 1. The maximum thickness of the films is ~ 3 μm that exceeds those of previously reported anodized WO3 films in fluoride-containing electrolytes. By investigating the photocatalytic properties of the films with thicknesses ranging from ~ 0.5 to ~ 3 μm, the optimum thickness of the nanoporous film is found to be ~ 1 μm, which demonstrates an impressive 120% improvement in the photocatalytic performance compared to that of a RF-sputtered nanotextured film with similar weights. We mainly ascribe this to large surface area and smaller bandgap.

Experiences of regional and rural people with cancer being treated with radiotherapy in a metropolitan centre

This paper explores the issues related to rural people with cancer whose choice of radiotherapy treatment necessitated travel and accommodation in a metropolitan centre. Semi-structured interviews with 46 participants, from the Toowoomba and Darling Downs region of Queensland, Australia, were conducted and the data thematically analysed. The specific themes identified were: being away from loved ones, maintaining responsibilities whilst undergoing treatment, emotional stress, burden on significant others, choice about radiotherapy as a treatment, travel and accommodation, and financial burden. This study supports the need for a radiotherapy centre in the location of Toowoomba as a way of providing some equity and access to such treatment for the rural people of Queensland.

Carbonic anhydrase IX expression and outcome after radiotherapy for muscle-invasive bladder cancer

Aims: Carbonic anhydrase IX (CA IX) expression has been described as an endogenous marker of hypoxia in solid neoplasms. Furthermore, CA IX expression has been associated with an aggressive phenotype and resistance to radiotherapy. We assessed the prognostic significance of CA IX expression in patients with muscle-invasive bladder cancer treated with radiotherapy. Materials and methods: A standard immunohistochemistry technique was used to show CA IX expression in 110 muscle-invasive bladder tumours treated with radiotherapy. Clinicopathological data were obtained from medical case notes. Results: CA IX immunostaining was detected in 89 (∼81%) patients. Staining was predominantly membranous, with areas of concurrent cytoplasmic and nuclear staining and was abundant in luminal and perinecrotic areas. No significant correlation was shown between the overall CA IX status and the initial response to radiotherapy, 5-year bladder cancer-specific survival or the time to local recurrence. Conclusions: The distribution of CA IX expression in paraffin-embedded tissue sections seen in this series is consistent with previous studies in bladder cancer, but does not provide significant prognostic information with respect to the response to radiotherapy at 3 months and disease-specific survival after radical radiotherapy. © 2007 The Royal College of Radiologists.

Weekly cisplatin concurrently with radiotherapy in head and neck squamous cell cancer : a retrospective analysis of a tertiary institute experience

Radiotherapy combined with three weekly 100 mg/m2 of cisplatin is the accepted standard of care in head and neck squamous cell carcinoma. However, this regimen is associated with severe toxicities with devastating effects on patients. Alternative protocols like weekly 40 mg/m2 have been used in an attempt to reduce toxicities. The main objective of the present study is to identify the dose intensities and toxicities of weekly cisplatin in patients treated in a tertiary centre over a 12 month period. Included patients had squamous cell carcinoma arising in the oral cavity, oropharynx, larynx, or hypopharynx. Patients were excluded if they had nasopharyngeal squamous cell carcinoma, distant metastasis or if they had prior treatment for head and neck cancer excluding neck dissection. During the study period, 52 patients met the inclusion criteria and their data were retrospectively obtained from the patients' database of St James hospital, Dublin. The median age of the study cohort was 54 years (range 33-73). Of the patients, 40 (76.9 %) were male and 12 (20.1 %) were female. The primary tumour sites were as follows: oral cavity and oropharynx in 38 (73 %), larynx in 10 (19 %), and hypopharynx in 4 (8 %). In total, 33 (63.5 %) patients had stage IV disease, while 19 (36.5 %) had stage III disease. Treatment was definitive in 35 (67 %) patients and adjuvant in 17 (35 %). Full-dose radiotherapy was achieved in 50 (96 %) patients. Only 22 (42.3 %) patients completed the intended six cycles of chemotherapy. Cumulative dose of 200 mg/m2 or more was reached in 37 (71 %) patients. The acute adverse effects included grades 3 and 4 mucositis, which occurred in 22 (43.3 %) and 6 patients (12 %), respectively. Grade 3 and 4 neutropenia occurred in six (11.5 %) and three (5.7 %) patients, respectively. The only other haematological toxicity was grade 3 anaemia in 20 (38.4 %) patients. There was no grade 3 or 4 renal toxicity among the study cohort, although grade 2 was observed in six (11.5 %) patients. Death occurred in one patient due to neutropenic septicaemia. In conclusion, weekly cisplatin is associated with moderate to severe toxicities and might lead to suboptimal chemotherapy delivery. More prospective clinical studies are required to determine the optimal chemoradiation regimen in head and neck squamous cell carcinoma.

Constructing a clinical decision-making framework for image-guided radiotherapy using a Bayesian Network

A decision-making framework for image-guided radiotherapy (IGRT) is being developed using a Bayesian Network (BN) to graphically describe, and probabilistically quantify, the many interacting factors that are involved in this complex clinical process. Outputs of the BN will provide decision-support for radiation therapists to assist them to make correct inferences relating to the likelihood of treatment delivery accuracy for a given image-guided set-up correction. The framework is being developed as a dynamic object-oriented BN, allowing for complex modelling with specific sub-regions, as well as representation of the sequential decision-making and belief updating associated with IGRT. A prototype graphic structure for the BN was developed by analysing IGRT practices at a local radiotherapy department and incorporating results obtained from a literature review. Clinical stakeholders reviewed the BN to validate its structure. The BN consists of a sub-network for evaluating the accuracy of IGRT practices and technology. The directed acyclic graph (DAG) contains nodes and directional arcs representing the causal relationship between the many interacting factors such as tumour site and its associated critical organs, technology and technique, and inter-user variability. The BN was extended to support on-line and off-line decision-making with respect to treatment plan compliance. Following conceptualisation of the framework, the BN will be quantified. It is anticipated that the finalised decision-making framework will provide a foundation to develop better decision-support strategies and automated correction algorithms for IGRT.

Gelatine methacrylamide-based hydrogels : an alternative three-dimensional cancer cell culture system

Modern cancer research requires physiological, three-dimensional (3-D) cell culture platforms, wherein the physical and chemical characteristics of the extracellular matrix (ECM) can be modified. In this study, gelatine methacrylamide (GelMA)-based hydrogels were characterized and established as in vitro and in vivo spheroid-based models for ovarian cancer, reflecting the advanced disease stage of patients, with accumulation of multicellular spheroids in the tumour fluid (ascites). Polymer concentration (2.5-7% w/v) strongly influenced hydrogel stiffness (0.5±0.2kPa to 9.0±1.8kPa) but had little effect on solute diffusion. The diffusion coefficient of 70kDa fluorescein isothiocyanate (FITC)-labelled dextran in 7% GelMA-based hydrogels was only 2.3 times slower compared to water. Hydrogels of medium concentration (5% w/v GelMA) and stiffness (3.4kPa) allowed spheroid formation and high proliferation and metabolic rates. The inhibition of matrix metalloproteinases and consequently ECM degradability reduced spheroid formation and proliferation rates. The incorporation of the ECM components laminin-411 and hyaluronic acid further stimulated spheroid growth within GelMA-based hydrogels. The feasibility of pre-cultured GelMA-based hydrogels as spheroid carriers within an ovarian cancer animal model was proven and led to tumour development and metastasis. These tumours were sensitive to treatment with the anti-cancer drug paclitaxel, but not the integrin antagonist ATN-161. While paclitaxel and its combination with ATN-161 resulted in a treatment response of 33-37.8%, ATN-161 alone had no effect on tumour growth and peritoneal spread. The semi-synthetic biomaterial GelMA combines relevant natural cues with tunable properties, providing an alternative, bioengineered 3-D cancer cell culture in in vitro and in vivo model systems.

A simple method for EPID-based in vivo dosimetry for radiotherapy treatments of head-and-neck cancers

This study used a homogeneous water-equivalent model of an electronic portal imaging device (EPID), contoured as a structure in a radiotherapy treatment plan, to produce reference dose images for comparison with in vivo EPID dosimetry images. Head and neck treatments were chosen as the focus of this study, due to the heterogeneous anatomies involved and the consequent difficulty of rapidly obtaining reliable reference dose images by other means. A phantom approximating the size and heterogeneity of a typical neck, with a maximum radiological thickness of 8.5 cm, was constructed for use in this study. This phantom was CT scanned and a simple treatment including five square test fields and one off-axis IMRT field was planned. In order to allow the treatment planning system to calculate dose in a model EPID positioned a distance downstream from the phantom to achieve a source-to-detector distance (SDD) of 150 cm, the CT images were padded with air and the phantom’s “body” contour was extended to encompass the EPID contour. Comparison of dose images obtained from treatment planning calculations and experimental irradiations showed good agreement, with more than 90% of points in all fields passing a gamma evaluation, at γ (3%, 3mm )Similar agreement was achieved when the phantom was over-written with air in the treatment plan and removed from the experimental beam, suggesting that water EPID model at 150 cm SDD is capable of providing accurate reference images for comparison with clinical IMRT treatment images, for patient anatomies with radiological thicknesses ranging from 0 up to approximately 9 cm. This methodology therefore has the potential to be used for in vivo dosimetry during treatments to tissues in the neck as well as the oral and nasal cavities, in the head-and-neck region.