Monte-Carlo dosimetry on a realistic cell monolayer geometry exposed to alpha particles

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.

A GATE simulation study for dosimetry in cancer cell and micrometastasis from the 225Ac decay chain

BACKGROUND

Radiopharmaceutical therapy (RPT) with alpha-emitting radionuclides has shown great promise in treating metastatic cancers. The successive emission of four alpha particles in the ²²⁵Ac decay chain leads to highly targeted and effective cancer cell death. Quantifying cellular dosimetry for ²²⁵Ac RPT is essential for predicting cell survival and therapeutic success. However, the leading assumption that all ²²⁵Ac progeny remain localized at their target sites likely overestimates the absorbed dose to cancer cells. To address limitations in existing semi-analytic approaches, this work evaluates S-values for ²²⁵Ac's progeny radionuclides with GATE Monte Carlo simulations.

METHODS

The cellular geometries considered were an individual cell (10 µm diameter with a nucleus of 8 µm diameter) and a cluster of cells (micrometastasis) with radionuclides localized in four subcellular regions: cell membrane, cytoplasm, nucleus, or whole cell. The absorbed dose to the cell nucleus was scored, and self- and cross-dose S-values were derived. We also evaluated the total absorbed dose with various degrees of radiopharmaceutical internalization and retention of the progeny radionuclides ²²¹Fr (t_1/2 = 4.80 m) and ²¹³Bi (t_1/2 = 45.6 m).

RESULTS

For the cumulative ²²⁵Ac decay chain, our self- and cross-dose nuclear S-values were both in good agreement with S-values published by MIRDcell, with per cent differences ranging from - 2.7 to - 8.7% for the various radionuclide source locations. Source location had greater effects on self-dose S-values than the intercellular cross-dose S-values. Cumulative ²²⁵Ac decay chain self-dose S-values increased from 0.167 to 0.364 GyBq^-1 s^-1 with radionuclide internalization from the cell surface into the cell. When progeny migration from the target site was modelled, the cumulative self-dose S-values to the cell nucleus decreased by up to 71% and 21% for ²²¹Fr and ²¹³Bi retention, respectively.

CONCLUSIONS

Our GATE Monte Carlo simulations resulted in cellular S-values in agreement with existing MIRD S-values for the alpha-emitting radionuclides in the ²²⁵Ac decay chain. To obtain accurate absorbed dose estimates in ²²⁵Ac studies, accurate understanding of daughter migration is critical for optimized injected activities. Future work will investigate other novel preclinical alpha-emitting radionuclides to evaluate therapeutic potency and explore realistic cellular geometries corresponding to targeted cancer cell lines.

Recruitment Kinetics of XRCC1 and RNF8 Following MeV Proton and α-Particle Micro-Irradiation

Time-lapse fluorescence imaging coupled to micro-irradiation devices provides information on the kinetics of DNA repair protein accumulation, from a few seconds to several minutes after irradiation. Charged-particle microbeams are valuable tools for such studies since they provide a way to selectively irradiate micrometric areas within a cell nucleus, control the dose and the micro-dosimetric quantities by means of advanced detection systems and Monte Carlo simulations and monitor the early cell response by means of beamline microscopy. We used the charged-particle microbeam installed at the AIFIRA facility to perform micro-irradiation experiments and measure the recruitment kinetics of two proteins involved in DNA signaling and repair pathways following exposure to protons and α-particles. We developed and validated image acquisition and processing methods to enable a systematic study of the recruitment kinetics of GFP-XRCC1 and GFP-RNF8. We show that XRCC1 is recruited to DNA damage sites a few seconds after irradiation as a function of the total deposited energy and quite independently of the particle LET. RNF8 is recruited to DNA damage sites a few minutes after irradiation and its recruitment kinetics depends on the particle LET.

On the effectiveness of proton boron fusion therapy (PBFT) at cellular level

The present work introduced a framework to investigate the effectiveness of proton boron fusion therapy (PBFT) at the cellular level. The framework consisted of a cell array generator program coupled with PHITS Monte Carlo package with a dedicated terminal-based code editor that was developed in this work. The framework enabled users to model large cell arrays with normal, all boron, and random boron filled cytoplasm, to investigate the underlying mechanism of PBFT. It was found that alpha particles and neutrons could be produced in absence of boron mainly because of nuclear reaction induced by proton interaction with ¹⁶O, ¹²C and ¹⁴N nuclei. The effectiveness of PBFT is highly dependent on the incident proton energy, source size, cell array size, buffer medium thickness layer, concentration and distribution of boron in the cell array. To quantitatively assess the effectiveness of PBFT, of the total energy deposition by alpha particle for different cases were determined. The number of alpha particle hits in cell cytoplasm and nucleus for normal and 100 ppm boron were determined. The obtained results and the developed tools would be useful for future development of PBFT to objectively determine the effectiveness of this treatment modality.

Small-scale dosimetry for alpha particle 241Am source cell irradiation and estimation of γ-H2AX foci distribution in prostate cancer cell line PC3

Background

The development of new targeted alpha therapies motivates improving alpha particle dosimetry. For alpha particles, microscopic targets must be considered to estimate dosimetric quantities that can predict the biological response. As double-strand breaks (DSB) on DNA are the main cause of cell death by ionizing radiation, cell nuclei are relevant volumes necessary to consider as targets. Since a large variance is expected of alpha particle hits in individual cell nuclei irradiated by an uncollimated alpha-emitting source, the damage induced should have a similar distribution. The induction of DSB can be measured by immunofluorescent γ-H2AX staining. The cell γ-H2AX foci distribution and alpha particle hits distribution should be comparable and thereby verify the necessity to consider the relevant dosimetric volumes.

Methods

A Monte Carlo simulation model of an ²⁴¹Am source alpha particle irradiation setup was combined with two versions of realistic cell nuclei phantoms. These were generated from DAPI-stained PC3 cells imaged with fluorescent microscopy, one consisting of elliptical cylinders and the other of segmented mesh volumes. PC3 cells were irradiated with the ²⁴¹Am source for 4, 8 and 12 min, and after 30 min fixated and stained with immunofluorescent γ-H2AX marker. The detected radiation-induced foci (RIF) were compared to simulated RIF.

Results

The mesh volume phantom detected a higher mean of alpha particle hits and energy imparted (MeV) per cell nuclei than the elliptical cylinder phantom, but the mean specific energy (Gy) was very similar. The mesh volume phantom detected a slightly larger variance between individual cells, stemming from the more extreme and less continuous distribution of cell nuclei sizes represented in this phantom. The simulated RIF distribution from both phantoms was in good agreement with the detected RIF, although the detected distribution had a zero-inflated shape not seen in the simulated distributions. An estimate of undetected foci was used to correct the detected RIF distribution and improved the agreement with the simulations.

Conclusion

Two methods to generate cell nuclei phantoms for Monte Carlo dosimetry simulations were tested and generated similar results. The simulated and detected RIF distributions from alpha particle-irradiated PC3 cells were in good agreement, proposing the necessity to consider microscopic targets in alpha particle dosimetry.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40658-022-00475-x.

CPOP: An open source C++ cell POPulation modeler for radiation biology applications

PURPOSE

Multicellular tumor spheroids are realistic in-vitro systems used in radiation biology research to study the effect of anticancer drugs or to evaluate the resistance of cancer cells under specific conditions. When combining the modeling of spheroids together with the simulation of radiation using Monte Carlo methods, one could estimate cell and DNA damage to be compared with experimental data. We developed a Cell Population (CPOP) modeler combined to Geant4 simulations in order to tackle how energy depositions are allocated to cells, especially when enhancing radiation outcomes using high-Z nanoparticles. CPOP manages to model large three-dimensional cell populations with independent deformable cells described with their nucleus, cytoplasm and membranes together with force law systems to manage cell-cell interactions.

METHODS

CPOP is an opensource platform written in C++. It is divided into two main libraries: a "Modeler" library, for cell geometry modeling using meshes, and a Multi Agent System (MAS) library, simulating all agent (cell) interactions among the population. CPOP is fully interfaced with the Geant4 Monte Carlo toolkit and is able to directly launch Geant4 simulations after compilation. We modeled a full and realistic 3D cell population from SK-MEL28 melanoma cell population cultured experimentally. The spheroid diameter of 550 ± 40 µm corresponds to a population of approximately 1000 cells having a diameter of 17.2 ± 2.5 µm and a nucleus diameter of 11.2 ± 2.0 µm. We decided to reproduce cell irradiations performed with a X-RAD 320 Biological Irradiator (Precision XRay Inc., North Branford, CT).

RESULTS

We simulated the energy spectrum of secondary particles generated in the vicinity of the spheroid and plotted the different energy spectra recovered internally to the spheroid. We evaluated also the impact of AGuIX (Gadolinium) nanoparticles modeled into the spheroid with their corresponding secondary energy spectra.

CONCLUSIONS

We succeeded into modeling cell populations and combined them with Geant4 simulations. The next step will be to integrate DNA geometrical models into cell nuclei and to use the Geant4-DNA physics and radiolysis modeling capabilities in order to evaluate early strand breaks induced on DNA.

Study on Caenorhabditis Elegans as a Combined Model of Microdosimetry and Biology

Microdosimetry is a tool for the investigation of microscopic energy deposition of ionizing radiation. This work used Caenorhabditis elegans as a model to estimate the microdosimetric deposition level at the ⁶⁰Co gamma radiation. Monte Carlo software PHITS was employed to establish irradiated nematodes model. The dose deposition of the entire body and gonad irradiated to 100 Gy was calculated. The injury levels of radiation were evaluated by the detection of biological indicators. The result of microdosimetric experiment suggested that the dose of whole body of nematodes was estimated to be 99.9 ± 57.8 Gy, ranging from 19.6 to 332.2 Gy. The dose of gonad was predicted to be 129.4 ± 558.8 Gy (9.5-6597 Gy). The result of biological experiment suggested that there were little changes in the length of nematodes after irradiation. However, times of head thrash per minute and the spawning yield in 3 consecutive days decreased 27.1% and 94.7%, respectively. Nematodes in the irradiated group displayed heterogeneity. Through contour analysis, trends of behavior kinematics and reproductive capacity of irradiated nematodes proved to be consistent with the dose distribution levels estimated by microdosimetric model. Finally, C. elegans presented a suitable combined model of microdosimetry and biology for studying radiation.

Dosimetric uncertainties impact on cell survival curve with low energy proton

There is an increasing number of radiobiological experiments being conducted with low energy protons (less than 5 MeV) for radiobiological studies due to availability of sub-millimetre focused beam. However, low energy proton has broad microdosimetric spectra which can introduce dosimetric uncertainty. In this work, we quantify the impact of this dosimetric uncertainties on the cell survival curve and how it affects the estimation of the alpha and beta parameters in the LQ formalism. Monte Carlo simulation is used to generate the microdosimetric spectra in a micrometer-sized water sphere under proton irradiation. This is modelled using radiobiological experiment set-up at the Centre of Ion Beam Application (CIBA) in National University of Singapore. Our results show that the microdosimetric spectra can introduce both systematic and random shifts in dose and cell survival; this effect is most pronounced with low energy protons. The alpha and beta uncertainties can be up to 10% and above 30%, respectively for low energy protons passing through thin cell target (about 10 microns). These uncertainties are non-negligible and show that care must be taken in using the cell survival curve and its derived parameters for radiobiological models.

Monte Carlo dosimetry of a realistic multicellular model of follicular lymphoma in a context of radioimmunotherapy

PURPOSE

Small-scale dosimetry studies generally consider an artificial environment where the tumors are spherical and the radionuclides are homogeneously biodistributed. However, tumor shapes are irregular and radiopharmaceutical biodistributions are heterogeneous, impacting the energy deposition in targeted radionuclide therapy. To bring realism, we developed a dosimetric methodology based on a three-dimensional in vitro model of follicular lymphoma incubated with rituximab, an anti-CD20 monoclonal antibody used in the treatment of non-Hodgkin lymphomas, which might be combined with a radionuclide. The effects of the realistic geometry and biodistribution on the absorbed dose were highlighted by comparison with literature data. Additionally, to illustrate the possibilities of this methodology, the effect of different radionuclides on the absorbed dose distribution delivered to the in vitro tumor were compared.

METHODS

The starting point was a model named multicellular aggregates of lymphoma cells (MALC). Three MALCs of different dimensions and their rituximab biodistribution were considered. Geometry, antibody location and concentration were extracted from selective plane illumination microscopy. Assuming antibody radiolabeling with Auger electron (¹²⁵ I and ¹¹¹ In) and β^- particle emitters (¹⁷⁷ Lu, ¹³¹ I and ⁹⁰ Y), we simulated energy deposition in MALCs using two Monte Carlo codes: Geant4-DNA with "CPA100" physics models for Auger electron emitters and Geant4 with "Livermore" physics models for β^- particle emitters.

RESULTS

MALCs had ellipsoid-like shapes with major radii, r, of ~0.25, ~0.5 and ~1.3 mm. Rituximab was concentrated in the periphery of the MALCs. The absorbed doses delivered by ¹⁷⁷ Lu, ¹³¹ I and ⁹⁰ Y in MALCs were compared with literature data for spheres with two types of homogeneous biodistributions (on the surface or throughout the volume). Compared to the MALCs, the mean absorbed doses delivered in spheres with surface biodistributions were between 18% and 38% lower, while with volume biodistribution they were between 15% and 29% higher. Regarding the radionuclides comparison, the relationship between MALC dimensions, rituximab biodistribution and energy released per decay impacted the absorbed doses. Despite releasing less energy, ¹²⁵ I delivered a greater absorbed dose per decay than ¹¹¹ In in the r ~ 0.25 mm MALC (6.78·10^-2 vs 6.26·10^-2  µGy·Bq^-1 ·s^-1 ). Similarly, the absorbed doses per decay in the r ~ 0.5 mm MALC for ¹⁷⁷ Lu (2.41·10^-2  µGy·Bq^-1 ·s^-1 ) and ¹³¹ I (2.46·10^-2  µGy·Bq^-1 ·s^-1 ) are higher than for ⁹⁰ Y (1.98·10^-2  µGy·Bq^-1 ·s^-1 ). Furthermore, radionuclides releasing more energy per decay delivered absorbed dose more uniformly through the MALCs. Finally, when considering the radiopharmaceutical effective half-life, due to the biological half-life of rituximab being best matched by the physical half-life of ¹⁷⁷ Lu and ¹³¹ I compared to ⁹⁰ Y, the first two radionuclides delivered higher absorbed doses.

CONCLUSION

In the simulated configurations, β^- emitters delivered higher and more uniform absorbed dose than Auger electron emitters. When considering radiopharmaceutical half-lives, ¹⁷⁷ Lu and ¹³¹ I delivered absorbed doses higher than ⁹⁰ Y. In view of real irradiation of MALCs, such a work may be useful to select suited radionuclides and to help explain the biological effects.

Monte-Carlo dosimetry and real-time imaging of targeted irradiation consequences in 2-cell stage Caenorhabditis elegans embryo

Charged-particle microbeams (CPMs) provide a unique opportunity to investigate the effects of ionizing radiation on living biological specimens with a precise control of the delivered dose, i.e. the number of particles per cell. We describe a methodology to manipulate and micro-irradiate early stage C. elegans embryos at a specific phase of the cell division and with a controlled dose using a CPM. To validate this approach, we observe the radiation-induced damage, such as reduced cell mobility, incomplete cell division and the appearance of chromatin bridges during embryo development, in different strains expressing GFP-tagged proteins in situ after irradiation. In addition, as the dosimetry of such experiments cannot be extrapolated from random irradiations of cell populations, realistic three-dimensional models of 2 cell-stage embryo were imported into the Geant4 Monte-Carlo simulation toolkit. Using this method, we investigate the energy deposit in various chromatin condensation states during the cell division phases. The experimental approach coupled to Monte-Carlo simulations provides a way to selectively irradiate a single cell in a rapidly dividing multicellular model with a reproducible dose. This method opens the way to dose-effect investigations following targeted irradiation.

A Monte Carlo study of macroscopic and microscopic dose descriptors for kilovoltage cellular dosimetry

This work investigates how doses to cellular targets depend on cell morphology, as well as relations between cellular doses and doses to bulk tissues and water. Multicellular models of five healthy and cancerous soft tissues are developed based on typical values of cell compartment sizes, elemental compositions and number densities found in the literature. Cells are modelled as two concentric spheres with nucleus and cytoplasm compartments. Monte Carlo simulations are used to calculate the absorbed dose to the nucleus and cytoplasm for incident photon energies of 20-370 keV, relevant for brachytherapy, diagnostic radiology, and out-of-field radiation in higher-energy external beam radiotherapy. Simulations involving cell clusters, single cells and single nuclear cavities are carried out for cell radii between 5 and [Formula: see text]m, and nuclear radii between 2 and [Formula: see text]m. Seven nucleus and cytoplasm elemental compositions representative of animal cells are considered. The presence of a cytoplasm, extracellular matrix and surrounding cells can affect the nuclear dose by up to [Formula: see text]. Differences in cell and nucleus size can affect dose to the nucleus (cytoplasm) of the central cell in a cluster of 13 cells by up to [Formula: see text] ([Formula: see text]). Furthermore, the results of this study demonstrate that neither water nor bulk tissue are reliable substitutes for subcellular targets for incident photon energies  <50 keV: nuclear (cytoplasm) doses differ from dose-to-medium by up to [Formula: see text] ([Formula: see text]), and from dose-to-water by up to [Formula: see text] ([Formula: see text]). The largest differences between dose descriptors are seen for the lowest incident photon energies; differences are less than [Formula: see text] for energies [Formula: see text]90 keV. The sensitivity of results with regard to the parameters of the microscopic tissue structure model and cell model geometry, and the importance of the nucleus and cytoplasm as targets for radiation-induced cell death emphasize the importance of accurate models for cellular dosimetry studies.

Single α-particle irradiation permits real-time visualization of RNF8 accumulation at DNA damaged sites

As well as being a significant source of environmental radiation exposure, α-particles are increasingly considered for use in targeted radiation therapy. A better understanding of α-particle induced damage at the DNA scale can be achieved by following their tracks in real-time in targeted living cells. Focused α-particle microbeams can facilitate this but, due to their low energy (up to a few MeV) and limited range, α-particles detection, delivery, and follow-up observations of radiation-induced damage remain difficult. In this study, we developed a thin Boron-doped Nano-Crystalline Diamond membrane that allows reliable single α-particles detection and single cell irradiation with negligible beam scattering. The radiation-induced responses of single 3 MeV α-particles delivered with focused microbeam are visualized in situ over thirty minutes after irradiation by the accumulation of the GFP-tagged RNF8 protein at DNA damaged sites.

Review of Geant4-DNA applications for micro and nanoscale simulations

Emerging radiotherapy treatments including targeted particle therapy, hadron therapy or radiosensitisation of cells by high-Z nanoparticles demand the theoretical determination of radiation track structure at the nanoscale. This is essential in order to evaluate radiation damage at the cellular and DNA level. Since 2007, Geant4 offers physics models to describe particle interactions in liquid water at the nanometre level through the Geant4-DNA Package. This package currently provides a complete set of models describing the event-by-event electromagnetic interactions of particles with liquid water, as well as developments for the modelling of water radiolysis. Since its release, Geant4-DNA has been adopted as an investigational tool in kV and MV external beam radiotherapy, hadron therapies using protons and heavy ions, targeted therapies and radiobiology studies. It has been benchmarked with respect to other track structure Monte Carlo codes and, where available, against reference experimental measurements. While Geant4-DNA physics models and radiolysis modelling functionalities have already been described in detail in the literature, this review paper summarises and discusses a selection of representative papers with the aim of providing an overview of a) geometrical descriptions of biological targets down to the DNA size, and b) the full spectrum of current micro- and nano-scale applications of Geant4-DNA.

Applications of High-Throughput Clonogenic Survival Assays in High-LET Particle Microbeams

Charged particle therapy is increasingly becoming a valuable tool in cancer treatment, mainly due to the favorable interaction of particle radiation with matter. Its application is still limited due, in part, to lack of data regarding the radiosensitivity of certain cell lines to this radiation type, especially to high-linear energy transfer (LET) particles. From the earliest days of radiation biology, the clonogenic survival assay has been used to provide radiation response data. This method produces reliable data but it is not optimized for high-throughput microbeam studies with high-LET radiation where high levels of cell killing lead to a very low probability of maintaining cells' clonogenic potential. A new method, therefore, is proposed in this paper, which could potentially allow these experiments to be conducted in a high-throughput fashion. Cells are seeded in special polypropylene dishes and bright-field illumination provides cell visualization. Digital images are obtained and cell detection is applied based on corner detection, generating individual cell targets as x-y points. These points in the dish are then irradiated individually by a micron field size high-LET microbeam. Post-irradiation, time-lapse imaging follows cells' response. All irradiated cells are tracked by linking trajectories in all time-frames, based on finding their nearest position. Cell divisions are detected based on cell appearance and individual cell temporary corner density. The number of divisions anticipated is low due to the high probability of cell killing from high-LET irradiation. Survival curves are produced based on cell's capacity to divide at least four to five times. The process is repeated for a range of doses of radiation. Validation shows the efficiency of the proposed cell detection and tracking method in finding cell divisions.

Effects of radiation type and delivery mode on a radioresistant eukaryote Cryptococcus neoformans

INTRODUCTION

Most research on radioresistant fungi, particularly on human pathogens such as Cryptococcus neoformans, involves sparsely-ionizing radiation. Consequently, fungal responses to densely-ionizing radiation, which can be harnessed to treat life-threatening fungal infections, remain incompletely understood.

METHODS

We addressed this issue by quantifying and comparing the effects of densely-ionizing α-particles (delivered either by external beam or by (213)Bi-labeled monoclonal antibodies), and sparsely-ionizing (137)Cs γ-rays, on Cryptococcus neoformans.

RESULTS

The best-fit linear-quadratic parameters for clonogenic survival were the following: α = 0.24 × 10(-2) Gy(-1) for γ-rays and 1.07 × 10(-2) Gy(-1) for external-beam α-particles, and β = 1.44 × 10(-5) Gy(-2) for both radiation types. Fungal cell killing by radiolabeled antibodies was consistent with predictions based on the α-particle dose to the cell nucleus and the linear-quadratic parameters for external-beam α-particles. The estimated RBE (for α-particles vs. γ-rays) at low doses was 4.47 for the initial portion of the α-particle track, and 7.66 for the Bragg peak. Non-radiological antibody effects accounted for up to 23% of cell death.

CONCLUSIONS

These results quantify the degree of C. neoformans resistance to densely-ionizing radiations, and show how this resistance can be overcome with fungus-specific radiolabeled antibodies.

Voxel model of individual cells and its implementation in microdosimetric calculations using GEANT4

Accurate dosimetric calculations at cellular and sub-cellular levels are crucial to obtain an increased understanding of the interactions of ionizing radiation with a cell and its nucleus and cytoplasm. Ion microbeams provide a superior opportunity to irradiate small biological samples, e.g., DNA, cells, and to compare their response to computer simulations. However, the phantoms used to simulate small biological samples at cellular levels are often simplified as simple volumes filled with water. As a first step to improve the situation in comparing measurements of cell response to ionizing radiation with model calculations, a realistic voxel model of a KB cell was constructed and used together with an already constructed geometry and tracking 4 (GEANT4) model of the horizontal microbeam line of the Centre d'Etudes Nucléaires de Bordeaux-Gradignan (CENBG) 3.5 MV Van de Graaf accelerator at the CENBG, France. The microbeam model was then implemented into GEANT4 for simulations of the average number of particles hitting an irradiated cell when a specified number of particles are produced in the beam line. The result shows that when irradiating the developed voxel model of a KB cell with 200 α particles, with a nominal energy of 3 MeV in the beam line and 2.34 MeV at the cell entrance, 100 particles hit the cell on average. The mean specific energy is 0.209 ± 0.019 Gy in the nucleus and 0.044 ± 0.001 Gy in the cytoplasm. These results are in agreement with previously published data, which indicates that this model could act as a reference model for dosimetric calculations of radiobiological experiments, and that the proposed method could be applied to build a cell model database.

Response functions of imaging plates to photons, electrons and 4He particles

Imaging plates from Fuji (BAS-SR, MS, and TR types) are phosphor films routinely used in ultra high intensity laser experiments. However, few data are available on the absolute IP response functions to ionizing particles. We have previously measured and modeled the IP response functions to protons. We focus here on the determination of the responses to photons, electrons, and (4)He particles. The response functions are obtained on an energy range going from a few tens of keV to a few tens of MeV and are compared to available data. The IP sensitivities to the different ionizing particles demonstrate a quenching effect depending on the particle stopping power.

Absorbed fractions for alpha particles in ellipsoidal volumes

Internal dosimetry of alpha particles is gaining attention due to the increasing applications in cancer treatment and also for the assessment of environmental contamination from radionuclides. We developed a Monte Carlo simulation in GEANT4 in order to calculate the absorbed fractions for monoenergetic alpha particles in the energy interval between 0.1 and 10 MeV, uniformly distributed in ellipsoids made of soft tissue. For each volume, we simulated a spherical shape, three oblate and three prolate ellipsoids, and one scalene shape. For each energy and for every geometrical configuration, an analytical relationship between the absorbed fraction and a 'generalized radius' was found; and the dependence of the fit parameters on the alpha energy is discussed and fitted by parametric functions. With the proposed formulation, the absorbed fraction for alpha particles in the energy range explored can be calculated for volumes and for ellipsoidal shapes of practical interest. This method can be applied to the evaluation of absorbed fraction from alpha-emitting radionuclides. The contribution to the deposited energy coming from electron and photon emissions can be accounted for exploiting the specific formulations previously introduced. As an example of application, the dosimetry of (213)Bi and its decay chain in ellipsoids is reported.

In silico modelling of treatment-induced tumour cell kill: developments and advances

Mathematical and stochastic computer (in silico) models of tumour growth and treatment response of the past and current eras are presented, outlining the aims of the models, model methodology, the key parameters used to describe the tumour system, and treatment modality applied, as well as reported outcomes from simulations. Fractionated radiotherapy, chemotherapy, and combined therapies are reviewed, providing a comprehensive overview of the modelling literature for current modellers and radiobiologists to ignite the interest of other computational scientists and health professionals of the ever evolving and clinically relevant field of tumour modelling.