A feasibility study of X-ray phase-contrast mammographic tomography at the Imaging and Medical beamline of the Australian Synchrotron.

Results are presented of a recent experiment at the Imaging and Medical beamline of the Australian Synchrotron intended to contribute to the implementation of low-dose high-sensitivity three-dimensional mammographic phase-contrast imaging, initially at synchrotrons and subsequently in hospitals and medical imaging clinics. The effect of such imaging parameters as X-ray energy, source size, detector resolution, sample-to-detector distance, scanning and data processing strategies in the case of propagation-based phase-contrast computed tomography (CT) have been tested, quantified, evaluated and optimized using a plastic phantom simulating relevant breast-tissue characteristics. Analysis of the data collected using a Hamamatsu CMOS Flat Panel Sensor, with a pixel size of 100 µm, revealed the presence of propagation-based phase contrast and demonstrated significant improvement of the quality of phase-contrast CT imaging compared with conventional (absorption-based) CT, at medically acceptable radiation doses.

Review of chemical models and applications in Geant4-DNA: Report from the ESA BioRad III Project.

A chemistry module has been implemented in Geant4-DNA since Geant4 version 10.1 to simulate the radiolysis of water after irradiation. It has been used in a number of applications, including the calculation of G-values and early DNA damage, allowing the comparison with experimental data. Since the first version, numerous modifications have been made to the module to improve the computational efficiency and extend the simulation to homogeneous kinetics in bulk solution. With these new developments, new applications have been proposed and released as Geant4 examples, showing how to use chemical processes and models. This work reviews the models implemented and application developments for modeling water radiolysis in Geant4-DNA as reported in the ESA BioRad III Project.

A Geant4 based simulation platform of the HollandPTC R&D proton beamline for radiobiological studies.

A Geant4 based simulation platform of the Holland Proton Therapy Centre (HollandPTC, Netherlands) R&D beamline (G4HPTC-R&D) was developed to enable the planning, optimisation and advanced dosimetry for radiobiological studies. It implemented a six parameter non-symmetrical Gaussian pencil beam surrogate model to simulate the R&D beamline in both a pencil beam and passively scattered field configuration. Three different experimental proton datasets (70 MeV, 150 MeV, and 240 MeV) of the pencil beam envelope evolution in free air and depth-dose profiles in water were used to develop a set of individual parameter surrogate functions to enable the modelling of the non-symmetrical Gaussian pencil beam properties with only the ProBeam isochronous cyclotron mean extraction proton energy as input. This refined beam model was then benchmarked with respect to three independent experimental datasets of the R&D beamline operating in both a pencil beam configuration at 120 and 200 MeV, and passively scattered field configuration at 150 MeV. It was shown that the G4HPTC-R&D simulation platform can reproduce the pencil beam envelope evolution in free air and depth-dose profiles to within an accuracy on the order of ±5% for all tested energies, and that it was able to reproduce the 150 MeV passively scattered field to the specifications need for clinical and radiobiological applications.

Geant4-DNA simulation of human cancer cells irradiation with helium ion beams.

PURPOSE: This study aimed to develop a computational environment for the accurate simulation of human cancer cell irradiation using Geant4-DNA. New cell geometrical models were developed and irradiated by alpha particle beams to induce DNA damage. The proposed approach may help further investigation of the benefits of external alpha irradiation therapy. METHODS: The Geant4-DNA Monte Carlo (MC) toolkit allows the simulation of cancer cell geometries that can be combined with accurate modelling of physical, physicochemical and chemical stages of liquid water irradiation, including radiolytic processes. Geant4-DNA is used to calculate direct and non-direct DNA damage yields, such as single and double strand breaks, produced by the deposition of energy or by the interaction of DNA with free radicals. RESULTS: In this study, the "molecularDNA" example application of Geant4-DNA was used to quantify early DNA damage in human cancer cells upon irradiation with alpha particle beams, as a function of linear energy transfer (LET). The MC simulation results are compared to experimental data, as well as previously published simulation data. The simulation results agree well with the experimental data on DSB yields in the lower LET range, while the experimental data on DSB yields are lower than the results obtained with the "molecularDNA" example in the higher LET range. CONCLUSION: This study explored and demonstrated the possibilities of the Geant4-DNA toolkit together with the "molecularDNA" example to simulate the helium beam irradiation of cancer cell lines, to quantify the early DNA damage, or even the following DNA damage response.

Prediction of DNA rejoining kinetics and cell survival after proton irradiation for V79 cells using Geant4-DNA.

PURPOSE: Track structure Monte Carlo (MC) codes have achieved successful outcomes in the quantitative investigation of radiation-induced initial DNA damage. The aim of the present study is to extend a Geant4-DNA radiobiological application by incorporating a feature allowing for the prediction of DNA rejoining kinetics and corresponding cell surviving fraction along time after irradiation, for a Chinese hamster V79 cell line, which is one of the most popular and widely investigated cell lines in radiobiology. METHODS: We implemented the Two-Lesion Kinetics (TLK) model, originally proposed by Stewart, which allows for simulations to calculate residual DNA damage and surviving fraction along time via the number of initial DNA damage and its complexity as inputs. RESULTS: By optimizing the model parameters of the TLK model in accordance to the experimental data on V79, we were able to predict both DNA rejoining kinetics at low linear energy transfers (LET) and cell surviving fraction. CONCLUSION: This is the first study to demonstrate the implementation of both the cell surviving fraction and the DNA rejoining kinetics with the estimated initial DNA damage, in a realistic cell geometrical model simulated by full track structure MC simulations at DNA level and for various LET. These simulation and model make the link between mechanistic physical/chemical damage processes and these two specific biological endpoints.

Feasibility Study on the Radiation Dose by Radioactive Magnetic Core-Shell Nanoparticles for Open-Source Brachytherapy.

BACKGROUND: Treatment of early-stage breast cancer currently includes surgical removal of the tumor and (partial) breast irradiation of the tumor site performed at fractionated dose. Although highly effective, this treatment is exhaustive for both patient and clinic. In this study, the theoretical potential of an alternative treatment combining thermal ablation with low dose rate (LDR) brachytherapy using radioactive magnetic nanoparticles (RMNPs) containing 103-palladium was researched. METHODS: The radiation dose characteristics and emission spectra of a single RMNP were calculated, and dose distributions of a commercial brachytherapy seed and an RMNP brachytherapy seed were simulated using Geant4 Monte Carlo toolkit. RESULTS: It was found that the RMNP seeds deliver a therapeutic dose similar to currently used commercial seed, while the dose distribution shows a spherical fall off compared to the more inhomogeneous dose distribution of the commercial seed. Changes in shell thickness only changed the dose profile between 2 × 10-4 mm and 3 × 10-4 mm radial distance to the RMNP, not effecting long-range dose. CONCLUSION: The dose distribution of the RMNP seed is comparable with current commercial brachytherapy seeds, while anisotropy of the dose distribution is reduced. Because this reduces the dependency of the dose distribution on the orientation of the seed, their surgical placement is easier. This supports the feasibility of the clinical application of the proposed novel treatment modality.

Full Field X-Ray Scatter Tomography.

In X-ray imaging, photons are transmitted through and absorbed by the target object, but are also scattered in significant quantities. Previous attempts to use scattered X-ray photons for imaging applications used pencil or fan beam illumination. Here we present 3D X-ray Scatter Tomography using full-field illumination for small-animal imaging. Synchrotron imaging experiments were performed on a phantom and the chest of a juvenile rat. Transmitted and scattered photons were simultaneously imaged with separate cameras; a scientific camera directly downstream of the sample stage, and a pixelated detector with a pinhole imaging system placed at 45° to the beam axis. We obtained scatter tomogram feature fidelity sufficient for segmentation of the lungs and major airways in the rat. The image contrast in the scatter tomogram slices approached that of transmission imaging, indicating robustness to the amount of multiple scattering present in our case. This opens the possibility of augmenting full-field 2D imaging systems with additional scatter detectors to obtain complementary modes or to improve the fidelity of existing images without additional dose, potentially leading to single-shot or reduced-angle tomography or overall dose reduction for live animal studies.

In-silico optimisation of tileable philips digital SiPM based thin monolithic scintillator detectors for SPECT applications.

Over the last decade one of the most significant technological advances made in the field of radiation detectors for nuclear medicine was the development of Silicon Photomultipler (SiPM) sensors. At present only a small number of SiPM based radiation detectors for Single Photon Emission Computed Tomography (SPECT) applications have been explored, and even fewer experimental prototypes developed. An in-silico investigation into the optimal design of a Philips DPC3200 SiPM photosensor-based thin monolithic scintillator detector for SPECT applications was undertaken using the Monte Carlo radiation transport modelling toolkit Geant4 version 10.5. The performance of the 20 different SPECT radiation detector configurations, 4 scintillator materials (NaI(Tl), GAGG(Ce), CsI(Tl) and LaBr3(Ce)) and 5 thicknesses (1-5 mm), were determined through the use of seven figures of merit. It was found that a crystal thickness range of 4-5 mm was required for all four materials to ensure acceptable energy resolution, sensitivity and spatial resolution performance with the Philips DPC3200 SiPM. Any thinner than this and the performance of all four materials was found to degrade rapidly due to a high probability of material specific fluorescence x-ray escape after incident gamma/x-ray photoelectric absorption. When factoring in each material's magnetic resonance imaging compatibility, hygroscopy, and cost, it was found that CsI(Tl) represents the most promising material to construct tileable Philips digital SiPM based thin monolithic scintillator detectors for SPECT applications.

Evaluation of early radiation DNA damage in a fractal cell nucleus model using Geant4-DNA.

The advancement of multidisciplinary research fields dealing with ionising radiation induced biological damage - radiobiology, radiation physics, radiation protection and, in particular, medical physics - requires a clear mechanistic understanding of how cellular damage is induced by ionising radiation. Monte Carlo (MC) simulations provide a promising approach for the mechanistic simulation of radiation transport and radiation chemistry, towards the in silico simulation of early biological damage. We have recently developed a fully integrated MC simulation that calculates early single strand breaks (SSBs) and double strand breaks (DSBs) in a fractal chromatin based human cell nucleus model. The results of this simulation are almost equivalent to past MC simulations when considering direct/indirect strand break fraction, DSB yields and fragment distribution. The simulation results agree with experimental data on DSB yields within 13.6% on average and fragment distributions agree within an average of 34.8%.

Mechanistic DNA damage simulations in Geant4-DNA Part 2: Electron and proton damage in a bacterial cell.

We extended a generic Geant4 application for mechanistic DNA damage simulations to an Escherichia coli cell geometry, finding electron damage yields and proton damage yields largely in line with experimental results. Depending on the simulation of radical scavenging, electrons double strand breaks (DSBs) yields range from 0.004 to 0.010 DSB Gy-1 Mbp-1, while protons have yields ranging from 0.004 DSB Gy-1 Mbp-1 at low LETs and with strict assumptions concerning scavenging, up to 0.020 DSB Gy-1 Mbp-1 at high LETs and when scavenging is weakest. Mechanistic DNA damage simulations can provide important limits on the extent to which physical processes can impact biology in low background experiments. We demonstrate the utility of these studies for low dose radiation biology calculating that in E. coli, the median rate at which the radiation background induces double strand breaks is 2.8 × 10-8 DSB day-1, significantly less than the mutation rate per generation measured in E. coli, which is on the order of 10-3.

Mechanistic DNA damage simulations in Geant4-DNA part 1: A parameter study in a simplified geometry.

Mechanistic modelling of DNA damage in Monte Carlo simulations is highly sensitive to the parameters that define DNA damage. In this work, we use a simple testing geometry to investigate how different choices of physics models and damage model parameters can change the estimation of DNA damage in a mechanistic DNA damage simulation built in Geant4-DNA. The choice of physics model can lead to variations by up to a factor of two in the yield of physically induced strand breaks, and the parameters that determine scavenging, and physical and chemical single strand break induction can have even larger consequences. Using low energy electrons as primary particles, a variety of parameters are tested in this geometry in order to arrive at a parameter set consistent with past simulation studies. We find that the modelling of scavenging can play an important role in determining results, and speculate that high-scavenging regimes, where only chemical radicals within 1 nm of DNA are simulated, could provide a good means of testing mechanistic DNA simulations.

Towards photon radiotherapy treatment planning with high Z nanoparticle radiosensitisation agents: the Relative Biological Effective Dose (RBED) framework.

A novel treatment planning framework, the Relative Biological Effective Dose (RBED), for high Z nanoparticle (NP)-enhanced photon radiotherapy is developed and tested in silico for the medical exemplar of neoadjuvant (preoperative) breast cancer MV photon radiotherapy. Two different treatment scenarios, conventional and high Z NP enhanced, were explored with a custom Geant4 application that was developed to emulate the administration of a single 2 Gy fraction as part of a 50 Gy radiotherapy treatment plan. It was illustrated that there was less than a 1% difference in the dose deposition throughout the standard and high Z NP-doped adult female phantom. Application of the RBED framework found that the extent of possible biological response with high Z NP doping was great than expected via the dose deposition alone. It is anticipated that this framework will assist the scientific community in future high Z NP-enhanced in-silico, pre-clinical and clinical trials.

Observation of dose-rate dependence in a Fricke dosimeter irradiated at low dose rates with monoenergetic X-rays.

Absolute measurements of the radiolytic yield of Fe3+ in a ferrous sulphate dosimeter formulation (6 mM Fe2+), with a 20 keV x-ray monoenergetic beam, are reported. Dose-rate suppression of the radiolytic yield was observed at dose rates lower than and different in nature to those previously reported with x-rays. We present evidence that this effect is most likely to be due to recombination of free radicals radiolytically produced from water. The method used to make these measurements is also new and it provides radiolytic yields which are directly traceable to the SI standards system. The data presented provides new and exacting tests of radiation chemistry codes.

A local effect model-based interpolation framework for experimental nanoparticle radiosensitisation data.

A local effect model (LEM)-based framework capable of interpolating nanoparticle-enhanced photon-irradiated clonogenic cell survival fraction measurements as a function of nanoparticle concentration was developed and experimentally benchmarked for gold nanoparticle (AuNP)-doped bovine aortic endothelial cells (BAECs) under superficial kilovoltage X-ray irradiation. For three different superficial kilovoltage X-ray spectra, the BAEC survival fraction response was predicted for two different AuNP concentrations and compared to experimental data. The ability of the developed framework to predict the cell survival fraction trends is analysed and discussed. This developed framework is intended to fill in the existing gaps of individual cell line response as a function of NP concentration under photon irradiation and assist the scientific community in planning future pre-clinical trials of high Z nanoparticle-enhanced photon radiotherapy.

Simulating the Impact of the Natural Radiation Background on Bacterial Systems: Implications for Very Low Radiation Biological Experiments.

At very low radiation dose rates, the effects of energy depositions in cells by ionizing radiation is best understood stochastically, as ionizing particles deposit energy along tracks separated by distances often much larger than the size of cells. We present a thorough analysis of the stochastic impact of the natural radiative background on cells, focusing our attention on E. coli grown as part of a long term evolution experiment in both underground and surface laboratories. The chance per day that a particle track interacts with a cell in the surface laboratory was found to be 6 × 10-5 day-1, 100 times less than the expected daily mutation rate for E. coli under our experimental conditions. In order for the chance cells are hit to approach the mutation rate, a gamma background dose rate of 20 µGy hr-1 is predicted to be required.