A microdosimetry analysis of reversible electroporation in scattered, overlapping, and cancerous cervical cells

We present a numerical method for studying reversible electroporation on normal and cancerous cervical cells. This microdosimetry analysis builds on a unique approach for extracting contours of free and overlapping cervical cells in the cluster from the Extended Depth of Field (EDF) images. The algorithm used for extracting the contours is a joint optimization of multiple-level set function along with the Gaussian mixture model and Maximally Stable Extremal Regions. These contours are then exported to a multi-physics domain solver, where a variable frequency pulsed electric field is applied. The trans-Membrane voltage (TMV) developed across the cell membrane is computed using the Maxwell equation coupled with a statistical approach, employing the asymptotic Smoluchowski equation. The numerical model was validated by successful replication of existing experimental configurations that employed low-frequency uni-polar pulses on the overlapping cells to obtain reversible electroporation, wherein, several overlapping clumps of cervical cells were targeted. For high-frequency calculation, a combination of normal and cancerous cells is introduced to the computational domain. The cells are assumed to be dispersive and the Debye dispersion equation is used for further calculations. We also present the resulting strength-duration relationship for achieving the threshold value of electroporation between the normal and cancerous cervical cells due to their size and conductivity differences. The dye uptake modulation during the high-frequency electric field electroporation is further advocated by a mathematical model.

Calculation of biological effectiveness of SOBP proton beams: a TOPAS Monte Carlo study

Objective.This study aims to investigate the biological effectiveness of Spread-Out Bragg-Peak (SOBP) proton beams with initial kinetic energies 50-250 MeV at different depths in water using TOPAS Monte Carlo code.Approach.The study modelled SOBP proton beams using TOPAS time feature. Various LET-based models and Repair-Misrepair-Fixation model were employed to calculate Relative Biological Effectiveness (RBE) for V79 cell lines at different on-axis depths based on TOPAS. Microdosimetric Kinetic Model and biological weighting function-based models, which utilize microdosimetric distributions, were also used to estimate the RBE. A phase-space-based method was adopted for calculating microdosimetric distributions.Main results.The trend of variation of RBE with depth is similar in all the RBE models, but the absolute RBE values vary based on the calculation models. RBE sharply increases at the distal edge of SOBP proton beams. In the entrance region of all the proton beams, RBE values at 4 Gy i.e. RBE(4 Gy) resulting from different models are in the range of 1.04-1.07, comparable to clinically used generic RBE of 1.1. Moving from the proximal to distal end of the SOBP, RBE(4 Gy) is in the range of 1.15-1.33, 1.13-1.21, 1.11-1.17, 1.13-1.18 and 1.17-1.21, respectively for 50, 100, 150, 200 and 250 MeV SOBP beams, whereas at the distal dose fall-off region, these values are 1.68, 1.53, 1.44, 1.42 and 1.40, respectively.Significance.The study emphasises application of depth-, dose- and energy- dependent RBE values in clinical application of proton beams.

Integrating microdosimetricin vitroRBE models for particle therapy into TOPAS MC using the MicrOdosimetry-based modeliNg for RBE ASsessment (MONAS) tool

Objective.In this paper, we present MONAS (MicrOdosimetry-based modelliNg for relative biological effectiveness (RBE) ASsessment) toolkit. MONAS is a TOPAS Monte Carlo extension, that combines simulations of microdosimetric distributions with radiobiological microdosimetry-based models for predicting cell survival curves and dose-dependent RBE.Approach.MONAS expands TOPAS microdosimetric extension, by including novel specific energy scorers to calculate the single- and multi-event specific energy microdosimetric distributions at different micrometer scales. These spectra are used as physical input to three different formulations of themicrodosimetric kinetic model, and to thegeneralized stochastic microdosimetric model(GSM2), to predict dose-dependent cell survival fraction and RBE. MONAS predictions are then validated against experimental microdosimetric spectra andin vitrosurvival fraction data. To show the MONAS features, we present two different applications of the code: (i) the depth-RBE curve calculation from a passively scattered proton SOBP and monoenergetic12C-ion beam by using experimentally validated spectra as physical input, and (ii) the calculation of the 3D RBE distribution on a real head and neck patient geometry treated with protons.Main results.MONAS can estimate dose-dependent RBE and cell survival curves from experimentally validated microdosimetric spectra with four clinically relevant radiobiological models. From the radiobiological characterization of a proton SOBP and12C fields, we observe the well-known trend of increasing RBE values at the distal edge of the radiation field. The 3D RBE map calculated confirmed the trend observed in the analysis of the SOBP, with the highest RBE values found in the distal edge of the target.Significance.MONAS extension offers a comprehensive microdosimetry-based framework for assessing the biological effects of particle radiation in both research and clinical environments, pushing closer the experimental physics-based description to the biological damage assessment, contributing to bridging the gap between a microdosimetric description of the radiation field and its application in proton therapy treatment with variable RBE.

Multiscale Monte Carlo simulations for dosimetry in x-ray breast imaging: Part II - Microscopic scales

BACKGROUND

Although the benefits of breast screening and early diagnosis are known for reducing breast cancer mortality rates, the effects and risks of low radiation doses to the cells in the breast are still ongoing topics of study.

PURPOSE

To study specific energy distributions (f ( z , D g ) $f(z,D_{g})$ ) in cytoplasm and nuclei of cells corresponding to glandular tissue for different x-ray breast imaging modalities.

METHODS

A cubic lattice (500 μm length side) containing 4064 spherical cells was irradiated with photons loaded from phase space files with varying glandular voxel doses (D g $D_{g}$ ). Specific energy distributions were scored for nucleus and cytoplasm compartments using the PENELOPE (v. 2018) + penEasy (v. 2020) Monte Carlo (MC) code. The phase space files, generated in part I of this work, were obtained from MC simulations in a voxelized anthropomorphic phantom corresponding to glandular voxels for different breast imaging modalities, including digital mammography (DM), digital breast tomosynthesis (DBT), contrast enhanced digital mammography (CEDM) and breast CT (BCT).

RESULTS

In general, the average specific energy in nuclei is higher than the respective glandular dose scored in the same region, by up to 10%. The specific energy distributions for nucleus and cytoplasm are directly related to the magnitude of the glandular dose in the voxel (D g $D_{g}$ ), with little dependence on the spatial location. For similarD g $D_{g}$ values,f ( z , D g ) $f(z,D_{g})$ for nuclei is different between DM/DBT and CEDM/BCT, indicating that distinct x-ray spectra play significant roles inf ( z , D g ) $f(z,D_{g})$ . In addition, this behavior is also present when the specific energy distribution (F g ( z ) $F_{g}(z)$ ) is considered taking into account the GDD in the breast.

CONCLUSIONS

Microdosimetry studies are complementary to the traditional macroscopic breast dosimetry based on the mean glandular dose (MGD). For the same MGD, the specific energy distribution in glandular tissue varies between breast imaging modalities, indicating that this effect could be considered for studying the risks of exposing the breast to ionizing radiation.

Boron-Containing MOF Nanoparticles with Stable Metabolism in U87-MG Cells Combining Microdosimetry To Evaluate Relative Biological Effectiveness of Boron Neutron Capture Therapy

Accurate prediction of the relative biological effectiveness (RBE) of boron neutron capture therapy (BNCT) is challenging. The therapy is different from other radiotherapy; the dynamic distribution of boron-containing compounds in tumor cells affects the therapeutic outcome considerably and hampers accurate measurement of the neutron-absorbed dose. Herein, we used boron-containing metal-organic framework nanoparticles (BMOFs) with high boron content to target U87-MG cells and maintain the concentration of the 10B isotope in cells. The content of boron in the cells could maintain 90% (60 ppm) within 20 min compared with that at the beginning; therefore, the accurate RBE of BNCT can be acquired. The effects of BNCT upon cells after neutron irradiation were observed, and the neutron-absorbed dose was obtained by Monte Carlo simulations. The RBE of BMOFs was 6.78, which was 4.1-fold higher than that of a small-molecule boron-containing agent (boric acid). The energy spectrum of various particles was analyzed by Monte Carlo simulations, and the RBE was verified theoretically. Our results suggested that the use of nanoparticle-based boron carriers in BNCT may have many advantages and that maintaining a stable boron distribution within cells may significantly improve the efficiency of BNCT.

Diamond based integrated detection system for dosimetric and microdosimetric characterization of radiotherapy ion beams

BACKGROUND

Ion beam therapy allows for a substantial sparing of normal tissues and higher biological efficacy. Synthetic single crystal diamond is a very good material to produce high-spatial-resolution and highly radiation hard detectors for both dosimetry and microdosimetry in ion beam therapy.

PURPOSE

The aim of this work is the design, fabrication and test of an integrated waterproof detector based on synthetic single crystal diamond able to simultaneously perform dosimetric and microdosimetric characterization of clinical ion beams.

METHODS

The active elements of the integrated diamond device, that is, dosimeter and microdosimeter, were both realized in a Schottky diode configuration featured by different area, thickness, and shape by means of photolithography technologies for the selective growth of intrinsic and boron-doped CVD diamond. The cross-section of the sensitive volume of the dosimetric element is 4 mm2 and 1 μm-thick, while the microdosimetric one has an active cross-sectional area of 100 × 100 μm2 and a thickness of about 6.2 μm. The dosimetric and microdosimetric performance of the developed device was assessed at different depths in a water phantom at the MedAustron ion beam therapy facility using a monoenergetic uniformly scanned carbon ion beam of 284.7 MeV/u and proton beam of 148.7 MeV. The particle flux in the region of the microdosimeter was 6·107  cm2 /s for both irradiation fields. At each depth, dose and dose distributions in lineal energy were measured simultaneously and the dose mean lineal energy values were then calculated. Monte Carlo simulations were also carried out by using the GATE-Geant4 code to evaluate the relative dose, dose averaged linear energy transfer (LETd ), and microdosimetric spectra at various depths in water for the radiation fields used, by considering the contribution from the secondary particles generated in the ion interaction processes as well.

RESULTS

Dosimetric and microdosimetric quantities were measured by the developed prototype with relatively low noise (∼2 keV/μm). A good agreement between the measured and simulated dose profiles was found, with discrepancies in the peak to plateau ratio of about 3% and 4% for proton and carbon ion beams respectively, showing a negligible LET dependence of the dosimetric element of the device. The microdosimetric spectra were validated with Monte Carlo simulations and a good agreement between the spectra shapes and positions was found. Dose mean lineal energy values were found to be in close agreement with those reported in the literature for clinical ion beams, showing a sharp increase along the Bragg curve, being also consistent with the calculated LETd for all depths within the experimental error of 10%.

CONCLUSIONS

The experimental indicate that the proposed device can allow enhanced dosimetry in particle therapy centers, where the absorbed dose measurement is implemented by the microdosimetric characterization of the radiation field, thus providing complementary results. In addition, the proposed device allows for the reduction of the experimental uncertainties associated with detector positioning and could facilitate the partial overcoming of some drawbacks related to the low sensitivity of diamond microdosimeters to low LET radiation.

Variable RBE in proton radiotherapy: a comparative study with the predictive Mayo Clinic Florida microdosimetric kinetic model and phenomenological models of cell survival

Objectives. (1) To examine to what extent the cell- and exposure- specific information neglected in the phenomenological proton relative biological effectiveness (RBE) models could influence the computed RBE in proton therapy. (2) To explore similarities and differences in the formalism and the results between the linear energy transfer (LET)-based phenomenological proton RBE models and the microdosimetry-based Mayo Clinic Florida microdosimetric kinetic model (MCF MKM). (3) To investigate how the relationship between the RBE and the dose-mean proton LET is affected by the proton energy spectrum and the secondary fragments.Approach. We systematically compared six selected phenomenological proton RBE models with the MCF MKM in track-segment simulations, monoenergetic proton beams in a water phantom, and two spread-out Bragg peaks. A representative comparison within vitrodata for human glioblastoma cells (U87 cell line) is also included.Main results. Marked differences were observed between the results of the phenomenological proton RBE models, as reported in previous studies. The dispersion of these models' results was found to be comparable to the spread in the MCF MKM results obtained by varying the cell-specific parameters neglected in the phenomenological models. Furthermore, while single cell-specific correlation between RBE and the dose-mean proton LET seems reasonable above 2 keVμm-1, caution is necessary at lower LET values due to the relevant contribution of secondary fragments. The comparison within vitrodata demonstrates comparable agreement between the MCF MKM predictions and the results of the phenomenological models.Significance. The study highlights the importance of considering cell-specific characteristics and detailed radiation quality information for accurate RBE calculations in proton therapy. Furthermore, these results provide confidence in the use of the MCF MKM for clonogenic survival RBE calculations in proton therapy, offering a more mechanistic approach compared to phenomenological models.

Microdosimetry-based investigation of biological effectiveness of252Cf brachytherapy source: TOPAS Monte Carlo study

Objective.To investigate biological effectiveness of252Cf brachytherapy source using Monte Carlo-calculated microdosimetric distributions.Approach.252Cf source capsule was placed at the center of the spherical water phantom and phase-space data were scored as a function of radial distance in water (R= 1-5 cm) using TOPAS Monte Carlo code. The phase-space data were used to calculate microdosimetric distributions at 1μm site size. Using these distributions, Relative Biological Effectiveness (RBE), mean quality factor (Q̅) and Oxygen Enhancement Ratio (OER) were calculated as a function ofR.Main results.The overall shapes of the microdosimetric distributions are comparable at all the radial distances in water. However, slight variation in the bin-wise yield is observed withR. RBE,Q̅and OER are insensitive to R over the range 1-5 cm. Microdosimetric kinetic model based RBE values are about 2.3 and 2.8 for HSG tumour cells and V79 cells, respectively, whereas biological weighting function-based RBE is about 2.8. ICRP 60 and ICRU 40 recommendation-basedQ̅values are about 14.5 and 16, respectively. Theory of dual radiation action based RBE is 11.4. The calculated value of OER is 1.6.Significance.This study demonstrates the relative insensitivity of RBE,Q̅and OER radially away from the252Cf source along the distances of 1-5 cm in water.

Haralick texture analysis for microdosimetry: characterization of Monte Carlo generated 3D specific energy distributions

Objective.Explore the application of Haralick textural analysis to 3D distributions of specific energy (energy imparted per unit mass) scored in cell-scale targets considering varying mean specific energy (absorbed dose), target volume, and incident spectrum.Approach.Monte Carlo simulations are used to generate specific energy distributions in cell-scale water voxels ((1μm)3-(15μm)3) irradiated by photon sources (mean energies: 0.02-2 MeV) to varying mean specific energies (10-400 mGy). Five Haralick features (homogeneity, contrast, entropy, correlation, local homogeneity) are calculated using an implementation of Haralick analysis designed to reduce sensitivity to grey level quantization and are interpreted using fundamental radiation physics.Main results.Haralick measures quantify differences in 3D specific energy distributions observed with varying voxel volume, absorbed dose magnitude, and source spectrum. For example, specific energy distributions in small (1-3μm) voxels with low magnitudes of absorbed dose (10 mGy) have relatively high measures of homogeneity and local homogeneity and relatively low measures of contrast and entropy (all relative to measures for larger voxels), reflecting the many voxels with zero specific energy in an otherwise sporadic distribution. With increasing target size, energy is shared across more target voxels, and trends in Haralick measures, such as decreasing homogeneity and increasing contrast and entropy, reflect characteristics of each 3D specific energy distribution. Specific energy distributions for sources of differing mean energy are characterized by Haralick measures, e.g. contrast generally decreases with increasing source energy, correlation and homogeneity are often (not always) higher for higher energy sources.Significance.Haralick texture analysis successfully quantifies spatial trends in 3D specific energy distributions characteristic of radiation source, target size, and absorbed dose magnitude, thus offering new avenues to quantify microdosimetric data beyond first order histogram features. Promising future directions include investigations of multiscale tissue models, targeted radiation therapy techniques, and biological response to radiation.

Multiscale Monte Carlo simulations of gold nanoparticle dose-enhanced radiotherapy I: Cellular dose enhancement in microscopic models

BACKGROUND

The introduction of Gold NanoParticles (GNPs) in radiotherapy treatments necessitates considerations such as GNP size, location, and quantity, as well as patient geometry and beam quality. Physics considerations span length scales across many orders of magnitude (nanometer-to-centimeter), presenting challenges that often limit the scope of dosimetric studies to either micro- or macroscopic scales.

PURPOSE

To investigate GNP dose-enhanced radiation Therapy (GNPT) through Monte Carlo (MC) simulations that bridge micro-to-macroscopic scales. The work is presented in two parts, with Part I (this work) investigating accurate and efficient MC modeling at the single cell level to calculate nucleus and cytoplasm Dose Enhancement Factors (n,cDEFs), considering a broad parameter space including GNP concentration, GNP intracellular distribution, cell size, and incident photon energy. Part II then evaluates cell dose enhancement factors across macroscopic (tumor) length scales.

METHODS

Different methods of modeling gold within cells are compared, from a contiguous volume of either pure gold or gold-tissue mixture to discrete GNPs in a hexagonal close-packed lattice. MC simulations with EGSnrc are performed to calculate n,cDEF for a cell with radiusr cell = 7.35 $r_{\rm cell}=7.35$  µm and nucleusr nuc = 5 $r_{\rm nuc} = 5$  µm considering 10 to 370 keV incident photons, gold concentrations from 4 to 24 mgAu /gtissue , and three different GNP configurations within the cell: GNPs distributed around the surface of the nucleus (perinuclear) or GNPs packed into one (or four) endosome(s). Select simulations are extended to cells with different cell (and nucleus) sizes: 5 µm (2, 3, and 4 µm), 7.35 µm (4 and 6 µm), and 10 µm (7, 8, and 9 µm).

RESULTS

n,cDEFs are sensitive to the method of modeling gold in the cell, with differences of up to 17% observed; the hexagonal lattice of GNPs is chosen (as the most realistic model) for all subsequent simulations. Across cell/nucleus radii, source energies, and gold concentrations, both nDEF and cDEF are highest for GNPs in the perinuclear configuration, compared with GNPs in one (or four) endosome(s). Across all simulations of the (rcell , rnuc ) = (7.35, 5) µm cell, nDEFs and cDEFs range from unity to 6.83 and 3.87, respectively. Including different cell sizes, nDEFs and cDEFs as high as 21.5 and 5.5, respectively, are observed. Both nDEF and cDEF are maximized at photon energies above the K- or L-edges of gold by 10 to 20 keV.

CONCLUSIONS

Considering 5000 unique simulation scenarios, this work comprehensively investigates many physics trends on DEFs at the cellular level, including demonstrating that cellular DEFs are sensitive to gold modeling approach, intracellular GNP configuration, cell/nucleus size, gold concentration, and incident source energy. These data should prove especially useful in research as well as treatment planning, allowing one to optimize or estimate DEF using not only GNP uptake, but also account for average tumor cell size, incident photon energy, and intracellular configuration of GNPs. Part II will expand the investigation, taking the Part I cell model and applying it in cm-scale phantoms.

GPU-accelerated calculation of proton microdosimetric spectra as a function of target size, proton energy, and bounding volume size

Objective.Shortcomings of dose-averaged linear energy transfer (LETD), the quantity which is most commonly used to quantify proton relative biological effectiveness, have long been recognized. Microdosimetric spectra may overcome the limitations of LETDbut are extremely computationally demanding to calculate. A systematic library of lineal energy spectra for monoenergetic protons could enable rapid determination of microdosimetric spectra in a clinical environment. The objective of this work was to calculate and validate such a library of lineal energy spectra.Approach. SuperTrack, a GPU-accelerated CUDA/C++ based application, was developed to superimpose tracks calculated using Geant4 onto targets of interest and to compute microdosimetric spectra. Lineal energy spectra of protons with energies from 0.1 to 100 MeV were determined in spherical targets of diameters from 1 nm to 10μm and in bounding voxels with side lengths of 5μm and 3 mm.Main results.Compared to an analogous Geant4-based application, SuperTrack is up to 3500 times more computationally efficient if each track is resampled 1000 times. Dose spectra of lineal energy and dose-mean lineal energy calculated with SuperTrack were consistent with values published in the literature and with comparison to a Geant4 simulation. Using SuperTrack, we developed the largest known library of proton microdosimetric spectra as a function of primary proton energy, target size, and bounding volume size.Significance. SuperTrack greatly increases the computational efficiency of the calculation of microdosimetric spectra. The elevated lineal energy observed in a 3 mm side length bounding volume suggests that lineal energy spectra determined experimentally or computed in small bounding volumes may not be representative of the lineal energy spectra in voxels of a dose calculation grid. The library of lineal energy spectra calculated in this work could be integrated with a treatment planning system for rapid determination of lineal energy spectra in patient geometries.

The complexity of DNA damage by radiation follows a Gamma distribution: insights from the Microdosimetric Gamma Model

Introduction

DNA damage is the main predictor of response to radiation therapy for cancer. Its Q8 quantification and characterization are paramount for treatment optimization, particularly in advanced modalities such as proton and alpha-targeted therapy.

Methods

We present a novel approach called the Microdosimetric Gamma Model (MGM) to address this important issue. The MGM uses the theory of microdosimetry, specifically the mean energy imparted to small sites, as a predictor of DNA damage properties. MGM provides the number of DNA damage sites and their complexity, which were determined using Monte Carlo simulations with the TOPAS-nBio toolkit for monoenergetic protons and alpha particles. Complexity was used together with a illustrative and simplistic repair model to depict the differences between high and low LET radiations.

Results

DNA damage complexity distributions were were found to follow a Gamma distribution for all monoenergetic particles studied. The MGM functions allowed to predict number of DNA damage sites and their complexity for particles not simulated with microdosimetric measurements (yF) in the range of those studied.

Discussion

Compared to current methods, MGM allows for the characterization of DNA damage induced by beams composed of multi-energy components distributed over any time configuration and spatial distribution. The output can be plugged into ad hoc repair models that can predict cell killing, protein recruitment at repair sites, chromosome aberrations, and other biological effects, as opposed to current models solely focusing on cell survival. These features are particularly important in targeted alpha-therapy, for which biological effects remain largely uncertain. The MGM provides a flexible framework to study the energy, time, and spatial aspects of ionizing radiation and offers an excellent tool for studying and optimizing the biological effects of these radiotherapy modalities.

Haralick texture feature analysis for characterization of specific energy and absorbed dose distributions across cellular to patient length scales

Objective.To investigate an approach for quantitative characterization of the spatial distribution of dosimetric data by introducing Haralick texture feature analysis in this context.Approach.Monte Carlo simulations are used to generate 3D arrays of dosimetric data for 2 scenarios: (1) cell-scale microdosimetry: specific energy (energy imparted per unit mass) in cell-scale targets irradiated by photon spectra (¹²⁵I,¹⁹²Ir, 6 MV); (2) tumour-scale dosimetry: absorbed dose in voxels for idealized models of¹²⁵I permanent implant prostate brachytherapy, considering 'TG186' (realistic tissues including 0% to 5% intraprostatic calcifications; interseed attenuation) and 'TG43' (water model, no interseed attenuation) conditions. Five prominent Haralick features (homogeneity, contrast, correlation, local homogeneity, entropy) are computed and trends are interpreted using fundamental radiation physics.Main results.In the cell-scale scenario, the Haralick measures quantify differences in 3D specific energy distributions due to source spectra. For example, contrast and entropy are highest for¹²⁵I reflecting the large variations in specific energy in adjacent voxels (photoelectric interactions; relatively short range of electrons), while 6 MV has the highest homogeneity with smaller variations in specific energy between voxels (Compton scattering dominates; longer range of electrons). For the tumour-scale scenario, the Haralick measures quantify differences due to TG186/TG43 simulation conditions and the presence of calcifications. For example, as calcifications increase from 0% to 5%, contrast increases while correlation decreases, reflecting the large differences in absorbed dose in adjacent voxels (higher absorbed dose in voxels with calcification due to photoelectric interactions).Significance.Haralick texture analysis provides a quantitative method for the characterization of 3D dosimetric distributions across cellular to tumour length scales, with promising future applications including analyses of multiscale tissue models, patient-specific data, and comparison of treatment approaches.

The Importance of Subcellular Structures to the Modeling of Biological Cells in the Context of Computational Bioelectromagnetics Simulations

Numerical investigation of the interaction of electromagnetic fields with eukaryotic cells requires specifically adapted computer models. Virtual microdosimetry, used to investigate exposure, requires volumetric cell models, which are numerically challenging. For this reason, a method is presented here to determine the current and volumetric loss densities occurring in single cells and their distinct compartments in a spatially accurate manner as a first step toward multicellular models within the microstructure of tissue layers. To achieve this, 3D models of the electromagnetic exposure of generic eukaryotic cells of different shape (i.e. spherical and ellipsoidal) and internal complexity (i.e. different organelles) are performed in a virtual, finite element method-based capacitor experiment in the frequency range from 10 Hz to 100 GHz. In this context, the spectral response of the current and loss distribution within the cell compartments is investigated and any effects that occur are attributed either to the dispersive material properties of these compartments or to the geometric characteristics of the cell model investigated in each case. In these investigations, the cell is represented as an anisotropic body with an internal distributed membrane system of low conductivity that mimics the endoplasmic reticulum in a simplified manner. This will be used to determine which details of the cell interior need to be modeled, how the electric field and the current density will be distributed in this region, and where the electromagnetic energy is absorbed in the microstructure regarding electromagnetic microdosimetry. Results show that for 5 G frequencies, membranes make a significant contribution to the absorption losses. © 2023 The Authors. Bioelectromagnetics published by Wiley Periodicals LLC on behalf of Bioelectromagnetics Society.

MINAS TIRITH: a new tool for simulating radiation-induced DNA damage at the cell population level

Objective. The mechanisms of radiation-induced DNA damage can be understood via the fundamental acquisition of knowledge through a combination of experiments and modeling. Currently, most biological experiments are performed by irradiating an entire cell population, whereas modeling of radiation-induced effects is usually performed via Monte Carlo simulations with track structure codes coupled to realistic DNA geometries of a single-cell nucleus. However, the difference in scale between the two methods hinders a direct comparison because the dose distribution in the cell population is not necessarily uniform owing to the stochastic nature of the energy deposition. Thus, this study proposed the MINAS TIRITH tool to model the distribution of radiation-induced DNA damage in a cell population.Approach. The proposed method is based on precomputed databases of microdosimetric parameters and DNA damage distributions generated using the Geant4-DNA Monte Carlo Toolkit. First, a specific energyzwas assigned to each cell of an irradiated population for a particular absorbed doseDabs,following microdosimetric formalism. Then, each cell was assigned a realistic number of DNA damage events according to the specific energyz,respecting the stochastic character of its occurrence.Main results. This study validated the MINAS TIRITH tool by comparing its results with those obtained using the Geant4-DNA track structure code and a Geant4-DNA based simulation chain for DNA damage calculation. The different elements of comparison indicated consistency between MINAS TIRITH and the Monte Carlo simulation in case of the dose distribution in the population and the calculation of the amount of DNA damage.Significance. MINAS TIRITH is a new approach for the calculation of radiation-induced DNA damage at the cell population level that facilitates reasonable simulation times compared to those obtained with track structure codes. Moreover, this tool enables a more direct comparison between modeling and biological experimentation.

TOPAS simulations of the response of a mini-TEPC: benchmark with experimental data

Objective. Microdosimetry offers a fast tool for radiation quality (RQ) verification to be implemented in treatment planning systems in proton therapy based on variable LET or RBE to move forward from the use of a fixed RBE of 1.1. It is known that the RBE of protons can increase up to 50% higher than that value in the last few millimetres of their range. Microdosimetry can be performed both experimentally and by means of Monte Carlo (MC) simulations. This paper has the aim of comparing the two approaches.Approach. Experimental measurements have been performed using a miniaturized Tissue equivalent proportional counter developed at the Legnaro National Laboratories of the Italian National Institute for Nuclear Physics with the aim of being used as RQ monitors for high intensity beams. MC simulations have been performed using the microdosimetric extension of TOPAS which provides optimized parameters and scorers for this application.Main results. Simulations were compared with experimental microdosimetric spectra in terms of shape of the spectra and their average values. Moreover, the latter have been investigated as possible estimators of LET obtained with the same MC code. The shape of the spectra is in general consistent with the experimental distributions and the average values of the distributions in both cases can predict the RQ increase with depth.Significance. This study aims at the comparison of microdosimetric spectra obtained from both experimental measurements and the microdosimetric extension of TOPAS in the same radiation field.

First experimental measurements of 2D microdosimetry maps in proton therapy

BACKGROUND

Empirical data in proton therapy indicate that relative biological effectiveness (RBE) is not constant, and it is directly related to the linear energy transfer (LET). The experimental assessment of LET with high resolution would be a powerful tool for minimizing the LET hot spots in intensity-modulated proton therapy, RBE- or LET-guided evaluation and optimization to achieve biologically optimized proton plans, verifying the theoretical predictions of variable proton RBE models, and so on. This could impact clinical outcomes by reducing toxicities in organs at risk.

PURPOSE

The present work shows the first 2D LET maps obtained at a proton therapy facility using the double scattering delivery mode in clinical conditions by means of new silicon 3D-cylindrical microdetectors.

METHODS

The device consists of a matrix of 121 independent silicon-based detectors that have 3D-cylindrical electrodes of 25-µm diameter and 20-µm depth, resulting each one of them in a well-defined micrometric radiation sensitive volume etched inside the silicon. They have been specifically designed for a hadron therapy, improving the performance of current silicon-based microdosimeters. Microdosimetry spectra were obtained at different positions of the Bragg curve by using a water-equivalent phantom along an 89-MeV pristine proton beam generated in the Y1 proton passive scattering beamline of the Orsay Proton Therapy Centre (Institut Curie, France).

RESULTS

Microdosimetry 2D-maps showing the variation of the lineal energy with depth in the three dimensions were obtained in situ during irradiation at clinical fluence rates (∼10⁸  s^-1  cm^-2 ) for the first time with a spatial resolution of 200 µm, the highest achieved in the transverse plane so far. The experimental results were cross-checked with Monte Carlo simulations and a good agreement between the spectra shapes was found. The experimental frequency-mean lineal energy values in silicon were 1.858 ± 0.019 keV µm^-1 at the entrance, 2.61 ± 0.03 keV µm^-1 at the proximal distance, 4.97 ± 0.05 keV µm^-1 close to the Bragg peak, and 8.6 ± 0.1 keV µm^-1 at the distal edge. They are in good agreement with the expected trends in the literature in clinical proton beams.

CONCLUSIONS

We present the first 2D microdosimetry maps obtained in situ during irradiation at clinical fluence rates in proton therapy. Our results show that the arrays of 3D-cylindrical microdetectors are a reliable microdosimeter to evaluate LET maps not only in the longitudinal axis of the beam, but also in the transverse plane allowing for LET characterization in three dimensions. This work is a proof of principle showing the capacity of our system to deliver LET 2D maps. This kind of experimental data is needed to validate variable proton RBE models and to optimize LET-guided plans.

On the measurement uncertainty of microdosimetric quantities using diamond and silicon microdosimeters in carbon-ion beams

PURPOSE

The purpose of this paper is to compare the response of two different types of solid-state microdosimeters, that is, silicon and diamond, and their uncertainties. A study of the conversion of silicon microdosimetric spectra to the diamond equivalent for microdosimeters with different geometry of the sensitive volumes is performed, including the use of different stopping power databases.

METHOD

Diamond and silicon microdosimeters were irradiated under the same conditions, aligned at the same depth in a carbon-ion beam at the MedAustron ion therapy center. In order to estimate the microdosimetric quantities, the readout electronic linearity was investigated with three different methods, that is, the first being a single linear regression, the second consisting of a double linear regression with a channel transition and last a multiple linear regression by splitting the data into odd and even groups. The uncertainty related to each of these methods was estimated as well. The edge calibration was performed using the intercept with the horizontal axis of the tangent through the inflection point of the Fermi function approximation multi-channel analyzer spectrum. It was assumed that this point corresponds to the maximum energy difference of particle traversing the sensitive volume (SV) for which the residual range difference in the continuous slowing down approximation is equal to the thickness of the SV of the microdosimeter. Four material conversion methods were explored, the edge method, the density method, the maximum-deposition energy method and the bin-by-bin transformation method. The uncertainties of the microdosimetric quantities resulting from the linearization, the edge calibration and the detectors thickness were also estimated.

RESULTS

It was found that the double linear regression had the lowest uncertainty for both microdosimeters. The propagated standard (k = 1) uncertainties on the frequency-mean lineal energy y ¯ F ${\bar{y}}_{\rm{F}}$ and the dose-mean lineal energy y ¯ D ${\bar{y}}_{\rm{D}}$ values from the marker point, in the spectra, in the plateau were 0.1% and 0.2%, respectively, for the diamond microdosimeter, whilst for the silicon microdosimeter data converted to diamond, the uncertainty was estimated to be 0.1%. In the range corresponding to the 90% of the amplitude of the Bragg Peak at the distal part of the Bragg curve (R₉₀ ) the uncertainty was found to be 0.1%. The uncertainty propagation from the stopping power tables was estimated to be between 5% and 7% depending on the method. The uncertainty on the y ¯ F ${\bar{y}}_{\rm{F}}$ and y ¯ D ${\bar{y}}_{\rm{D}}$ coming from the thickness of the detectors varied between 0.3% and 0.5%.

CONCLUSION

This article demonstrate that the linearity of the readout electronics affects the microdosimetric spectra with a difference in y ¯ F ${\bar{y}}_{\rm{F}}$ values between the different linearization methods of up to 17.5%. The combined uncertainty was dominated by the uncertainty of stopping power on the edge.

Microdosimetric and radiobiological effects of gold nanoparticles at therapeutic radiation energies

PURPOSE

The purpose of this study was to quantify the microscopic dose distribution surrounding gold nanoparticles (GNPs) irradiated at therapeutic energies and to measure the changes in cell survival in vitro caused by this dose enhancement.

METHODS

The dose distributions from secondary electrons surrounding a single gold nanosphere and single gold nanocube of equal volume were both simulated using MCNP6. Dose enhancement factors (DEFs) in the 1 μm3 volume surrounding a GNP were calculated and compared between a nanosphere and nanocube and between 6 and 18 MV energies. This microscopic effect was explored further by experimentally measuring the cell survival of C-33a cervical cancer cells irradiated at 18 MV with varying doses of energy and concentrations of GNPs. Survival of cells receiving no irradiation, a 3 Gy dose, and a 6 Gy dose of 18 MV energy were determined for each concentration of GNPs.

RESULTS

It was observed that the dose from electrons surrounding the gold nanocube surpasses that of a gold nanosphere up to a distance of 1.1 μm by 18.5% for the 18 MV energy spectrum and by 23.1% for the 6 MV spectrum. DEFs ranging from ∼2 to 8 were found, with the maximum DEF resulting from the case of the gold nanocube irradiated at 6 MV energy. Experimentally, for irradiation at 18 MV, incubating cells with 6 nM (0.10% gold by mass) GNPs produces an average 6.7% decrease in cell survival, and incubating cells with 9 nM (0.15% gold by mass) GNPs produces an average 14.6% decrease in cell survival, as compared to cells incubated and irradiated without GNPs.

CONCLUSION

We have successfully demonstrated the potential radiation dose enhancing effects in vitro and microdosimetrically from gold nanoparticles.

Impact of Radiation Quality on Microdosimetry and Chromosome Aberrations for High-Energy (>250 MeV/n) Ions

Studying energy deposition by space radiation at the cellular scale provides insights on health risks to astronauts. Using the Monte Carlo track structure code RITRACKS, and the chromosome aberrations code RITCARD, we performed a modeling study of single-ion energy deposition spectra and chromosome aberrations for high-energy (>250 MeV/n) ion beams with linear energy transfer (LET) varying from 0.22 to 149.2 keV/µm. The calculations were performed using cells irradiated directly by mono-energetic ion beams, and by poly-energetic beams after particle transport in a digital mouse model, representing the radiation exposure of a cell in a tissue. To discriminate events from ion tracks directly traversing the nucleus, to events from δ-electrons emitted by distant ion tracks, we categorized ion contributions to microdosimetry or chromosome aberrations into direct and indirect contributions, respectively. The ions were either ions of the mono-energetic beam or secondary ions created in the digital mouse due to interaction of the beam with tissues. For microdosimetry, the indirect contribution is largely independent of the beam LET and minimally impacted by the beam interactions in mice. In contrast, the direct contribution is strongly dependent on the beam LET and shows increased probabilities of having low and high-energy deposition events when considering beam transport. Regarding chromosome aberrations, the indirect contribution induces a small number of simple exchanges, and a negligible number of complex exchanges. The direct contribution is responsible for most simple and complex exchanges. The complex exchanges are significantly increased for some low-LET ion beams when considering beam transport.

Microdosimetric characterization of a clinical proton therapy beam: comparison between simulated lineal energy distributions in spherical water targets and experimental measurements with a silicon detector

Objective Treatment planning based on computer simulations were proposed to account for the increase in the relative biological effectiveness (RBE) of proton radiotherapy beams near to the edges of the irradiated volume. Since silicon detectors could be used to validate the results of these simulations, it is important to explore the limitations of this comparison. Approach Microdosimetric measurements with a MicroPlus Bridge V2 silicon detector (thickness = 10 µm) were performed along the Bragg peak of a clinical proton beam. The lineal energy distributions, the dose mean values, and the RBE calculated with a biological weighting function were compared with simulations with PHITS (microdosimetric target = 1 µm water sphere), and published clonogenic survival in vitro RBE data for the V79 cell line. The effect of the silicon-to-water conversion was also investigated by comparing three different methodologies (conversion based on a single value, novel bin-to-bin conversions based on SRIM and PSTAR). Main results Mainly due to differences in the microdosimetric targets, the experimental dose-mean lineal energy and RBE values at the distal edge were respectively up to 53% and 28% lower than the simulated ones. Furthermore, the methodology chosen for the silicon-to-water conversion was proven to affect the dose mean lineal energy and the RBE10 up to 32% and 11% respectively. The best methodology to compensate for this underestimation was the bin-to-bin silicon-to-water conversion based on PSTAR. Significance This work represents the first comparison between PHITS-simulated lineal energy distributions in water targets and corresponding experimental spectra measured with silicon detectors. Furthermore, the effect of the silicon-to-water conversion on the RBE was explored for the first time. The proposed methodology based on the PSTAR bin-to-bin conversion appears to provide superior results with respect to commonly used single scaling factors and is recommended for future studies.

Track-structure modes in particle and heavy ion transport code system (PHITS): application to radiobiological research

PURPOSE

In radiation physics, Monte Carlo radiation transport simulations are powerful tools to evaluate the cellular responses after irradiation. When investigating such radiation-induced biological effects, it is essential to perform track structure simulations by explicitly considering each atomic interaction in liquid water at the sub-cellular and DNA scales. The Particle and Heavy-Ion Transport code System (PHITS) is a Monte Carlo code which enables to calculate track structure at DNA scale by employing the track-structure modes for electrons, protons and carbon ions. In this paper, we review the recent development status and future prospects of the track-structure modes in the PHITS code.

CONCLUSIONS

To date, the physical features of these modes have been verified using the available experimental data and Monte Carlo simulation results reported in literature. These track-structure modes can be used for calculating microdosimetric distributions to estimate cell survival and for estimating initial DNA damage yields. The use of PHITS track-structure mode is expected not only to clarify the underlying mechanisms of radiation effects but also to predict curative effects in radiation therapy. The results of PHITS simulations coupled with biophysical models will contribute to the radiobiological studies by precisely predicting radiation-induced biological effects based on the Monte Carlo approach.

Microdosimetric Investigation and a Novel Model of Radiosensitization in the Presence of Metallic Nanoparticles

Auger cascades generated in high atomic number nanoparticles (NPs) following ionization were considered a potential mechanism for NP radiosensitization. In this work, we investigated the microdosimetric consequences of the Auger cascades using the theory of dual radiation action (TDRA), and we propose the novel Bomb model as a general framework for describing NP-related radiosensitization. When triggered by an ionization event, the Bomb model considers the NPs that are close to a radiation sensitive cellular target, generates dense secondary electrons and kills the cells according to a probability distribution, acting like a “bomb.” TDRA plus a distance model were used as the theoretical basis for calculating the change in α of the linear-quadratic survival model and the relative biological effectiveness (RBE). We calculated these quantities for SQ20B and Hela human cancer cells under 250 kVp X-ray irradiation with the presence of gadolinium-based NPs (AGuIX^TM), and 220 kVp X-ray irradiation with the presence of 50 nm gold NPs (AuNPs), respectively, and compared with existing experimental data. Geant4-based Monte Carlo (MC) simulations were used to (1) generate the electron spectrum and the phase space data of photons entering the NPs and (2) calculate the proximity functions and other related parameters for the TDRA and the Bomb model. The Auger cascade electrons had a greater proximity function than photoelectric and Compton electrons in water by up to 30%, but the resulting increases in α were smaller than those derived from experimental data. The calculated RBEs cannot explain the experimental findings. The relative increase in α predicted by TDRA was lower than the experimental result by a factor of at least 45 for SQ20B cells with AGuIX under 250 kVp X-ray irradiation, and at least four for Hela cells with AuNPs under 220 kVp X-ray irradiation. The application of the Bomb model to Hela cells with AuNPs under 220 kVp X-ray irradiation indicated that a single ionization event for NPs caused by higher energy photons has a higher probability of killing a cell. NPs that are closer to the cell nucleus are more effective for radiosensitization. Microdosimetric calculations of the RBE for cell death of the Auger electron cascade cannot explain the experimentally observed radiosensitization by AGuIX or AuNP, while the proposed Bomb model is a potential candidate for describing NP-related radiosensitization at low NP concentrations.

Dose-mean lineal energy values for electrons by different Monte Carlo codes: Consequences for estimates of radiation quality in photon beams

BACKGROUND

The microdosimetric quantity lineal energy and its mean values have proven useful for quantifying radiation quality in many situations. The ratio of dose-mean lineal energies is perhaps the simplest quantity for quantifying differences between two radiation qualities. However, published dose-mean lineal energy values from different codes may differ significantly with potential influence on radiation quality estimates.

PURPOSE

The purpose was to compare dose-mean lineal energy values from different track-structure data sets for condensed water vapor and liquid water, and to evaluate the influence on radiation quality estimations for some photon sources.

METHODS

Published dose-mean lineal energy values for 0.1 keV to 1 MeV electrons in spheres with diameters 2 nm to 1 μm, calculated with water vapor and liquid water track structure codes and proximity functions, were collected, analyzed, and compared. Data for cylinders were converted to spheres using a theoretical transformation published by Kellerer. A new set of dose-mean lineal energy values was calculated to cover the whole range of volumes of interest here using the GEANT4-DNA code. The influence from the differences between codes on radiation quality calculations was estimated using dose-mean lineal energy ratios for the photon sources ¹²⁵ I, ¹⁶⁹ Yb, and ¹⁹² Ir relative to ⁶⁰ Co.

RESULTS

The theoretical relation for converting the dose-mean lineal energy between different geometrical volumes, results in differences up to 10% between cylinders and spheres depending on electron energy and target size, in agreement with published simulated results. For spheres with diameter above 100 nm, dose-mean lineal energy values for condensed water vapor and liquid water are with few exceptions within ±10%. Below 100 nm, the difference increases with decreasing diameter reaching a factor of two at 2 nm. The values from water vapor codes are in general larger than from liquid water codes. If the dose-mean lineal energy ratio is based on condensed water vapor instead of liquid water, the ratio differs less than 9% for the nuclides ¹²⁵ I, ¹⁶⁹ Yb, and ¹⁹² Ir relative to ⁶⁰ Co independent of the volume simulated. However, a specific value of the dose-mean lineal energy ratio, is found at a larger target diameter in liquid water than in condensed water vapor.

CONCLUSIONS

When ratios of the dose-mean lineal energy are used as a measure of the radiation quality it is important to compare values for geometrically equal target shapes. A practical method of converting values for cylinders of equal diameter and height to spheres was demonstrated. Although dose-mean lineal energy values calculated with water vapor and liquid water codes may differ significantly, the radiation quality, in terms of ratios of dose-mean lineal energy, for the three photon sources ¹⁹² Ir, ¹⁶⁹ Yb, and ¹²⁵ I relative to ⁶⁰ Co, agree within 9%. The same ratio appears at a larger diameter when a liquid water code is used. It is therefore important to use the same code in radiation quality investigations. The present findings may be of special interest in studies related to the relative biological effectiveness (RBE).

Comparison between the results of a recently-developed biological weighting function (V79-RBE10BWF) and thein vitroclonogenic survival RBE10of other repair-competent asynchronized normoxic mammalian cell lines and ions not used for the development of the model

728 simulated microdosimetric lineal energy spectra (26 different ions between 1H and 238U, 28 energy points from 1 to 1000 MeV/n) were used in combination with a recently-developed biological weighting function (Parisi et al., 2020) and 571 published in vitro clonogenic survival curves in order to: 1) assess prediction intervals for the in silico results by deriving an empirical indication of the experimental uncertainty from the dispersion in the in vitro hamster lung fibroblast (V79) data used for the development of the biophysical model; 2) explore the possibility of modeling the relative biological effectiveness (RBE) of the 10% clonogenic survival of asynchronized normoxic repair-competent mammalian cell lines other than the one used for the development of the model (V79); 3) investigate the predictive power of the model through a comparison between in silico results and in vitro data for 10 ions not used for the development of the model. At first, different strategies for the assessment of the in silico prediction intervals were compared. The possible sources of uncertainty responsible for the dispersion in the in vitro data were also shortly reviewed. Secondly, also because of the relevant scatter in the in vitro data, no statistically-relevant differences were found between the RBE10 of the investigated different asynchronized normoxic repair-competent mammalian cell lines. The only exception (Chinese Hamster peritoneal fibroblasts, B14FAF28), is likely due to the limited dataset (all in vitro ion data were extracted from a single publication), systematic differences in the linear energy transfer (LET) calculations for the employed very-heavy ions, and the use of reference photon survival curves extracted from a different publication. Finally, the in silico predictions for the 10 ions not used for the model development were in good agreement with the corresponding in vitro data.

Response of SOI microdosimeter in fast neutron beams: experiment and Monte Carlo simulations

In this study, Monte Carlo codes, Geant4 and MCNP6, were used to characterize the fast neutron therapeutic beam produced at iThemba LABS in South Africa. Experimental and simulation results were compared using the latest generation of Silicon on Insulator (SOI) microdosimeters from the Centre for Medical Radiation Physics (CMRP). Geant4 and MCNP6 were able to successfully model the neutron gantry and simulate the expected neutron energy spectrum produced from the reaction by protons bombarding a ⁹Be target. The neutron beam was simulated in a water phantom and its characteristics recorded by the silicon microdosimeters; bare and covered by a ¹⁰B enriched boron carbide converter, at different positions. The microdosimetric quantities calculated using Geant4 and MCNP6 are in agreement with experimental measurements. The thermal neutron sensitivity and production of ¹⁰B capture products in the p+ boron-implanted dopant regions of the Bridge microdosimeter is investigated. The obtained results are useful for the future development of dedicated SOI microdosimeters for Boron Neutron Capture Therapy (BNCT). This paper provides a benchmark comparison of Geant4 and MCNP6 capabilities in the context of further applications of these codes for neutron microdosimetry.

Uncertainty in tissue equivalent proportional counter assessments of microdosimetry and RBE estimates in carbon radiotherapy

Microdosimetry is an important tool for assessing energy deposition distributions from ionizing radiation at cellular and cellular nucleus scales. It has served as an input parameter for multiple common mathematical models, including evaluation of relative biological effectiveness (RBE) of carbon ion therapy. The most common detector used for microdosimetry is the tissue-equivalent proportional counter (TEPC). Although it is widely applied, TEPC has various inherent uncertainties. Therefore, this work quantified the magnitude of TEPC measurement uncertainties and their impact on RBE estimates for therapeutic carbon beams. Microdosimetric spectra and frequency-, dose-, and saturation-corrected dose-mean lineal energy (****) were calculated using the Monte Carlo toolkit Geant4 for five monoenergetic and three spread-out Bragg peak carbon beams in water at every millimeter along the central beam axis. We simulated the following influences on these spectra from eight sources of uncertainty: wall effects, pulse pile-up, electronics, gas pressure, W-value, gain instability, low energy cut-off, and counting statistics. Statistic uncertainty was quantified as the standard deviation of perturbed values for each source. Bias was quantified as the difference between default lineal energy values and the mean of perturbed values for each systematic source. Uncertainties were propagated to RBE using the modified microdosimetric kinetic model (MKM). Variance introduced by statistic sources iny¯Fandy¯Daveraged 3.8% and 3.4%, respectively, and 1.5% iny*across beam depths and energies. Bias averaged 6.2% and 7.3% iny¯Fandy¯D,and 4.8% iny*.These uncertainties corresponded to 1.2 ± 0.9% on average in RBE_MKM. The largest contributors to variance and bias were pulse pile-up and wall effects. This study established an error budget for microdosimetric carbon measurements by quantifying uncertainty inherent to TEPC measurements. It is necessary to understand how robust the measurement of RBE model input parameters are against this uncertainty in order to verify clinical model implementation.

The relation between microdosimetry and induction of direct damage to DNA by alpha particles

In radiopharmaceutical treatmentsα-particles are employed to treat tumor cells. However, the mechanism that drives the biological effect induced is not well known. Being ionizing radiation,α-particles can affect biological organisms by producing damage to the DNA, either directly or indirectly. Following the principle that microdosimetry theory accounts for the stochastic way in which radiation deposits energy in sub-cellular sized volumes via physical collisions, we postulate that microdosimetry represents a reasonable framework to characterize the statistical nature of direct damage induction byα-particles to DNA. We used the TOPAS-nBio Monte Carlo package to simulate direct damage produced by monoenergetic alpha particles to different DNA structures. In separate simulations, we obtained the frequency-mean lineal energy (yF) and dose-mean lineal energy (yD) of microdosimetric distributions sampled with spherical sites of different sizes. The total number of DNA strand breaks, double strand breaks (DSBs) and complex strand breaks per track were quantified and presented as a function of eitheryForyD.The probability of interaction between a track and the DNA depends on how the base pairs are compacted. To characterize this variability on compactness, spherical sites of different size were used to match these probabilities of interaction, correlating the size-dependent specific energy (z) with the damage induced. The total number of DNA strand breaks per track was found to linearly correlate withyFandzFwhen using what we defined an effective volume as microdosimetric site, while the yield of DSB per unit dose linearly correlated withyDorzD,being larger for compacted than for unfolded DNA structures. The yield of complex breaks per unit dose exhibited a quadratic behavior with respect toyDand a greater difference among DNA compactness levels. Microdosimetric quantities correlate with the direct damage imparted on DNA.

Elke Bräuer-Krisch: dedication, creativity and generosity: May 17, 1961-September 10, 2018

PURPOSE

This extraordinary woman worked her professional way from a radiation protection engineer to become the successful principal investigator of a prestigious international European project for a new radiation therapy (ERC Synergy grant, HORIZON 2020). The evaluation of the submitted proposal was very positive. The panel proposed that it be funded. Elke tragically passed away a few days before this conclusion of the panel. The present account describes her gradual career development; it includes many episodes that Elke personally chronicled in her curriculum of 2017.

METHODS

An internet literature search was performed using Google Scholar and other sources to assist in the writing of this narrative review and account.

CONCLUSIONS

In parallel to the development of the new Biomedical Beamline ID17 at the European Synchrotron Radiation Facility in Grenoble in the late nineties, Elke focused her interest and her personal and professional priorities on MRT, particularly on its clinical goals. She outlined her main objectives in several documents: (1) develop a new paradigm of cancer care by broadening the foundation for MRT. (2) Filling the gaps in basic biological knowledge about the mechanisms of MRT effects on normal and neoplastic tissues. (3) Broaden the preclinical level of evidence for the low normal organ toxicity of MRT versus standard X-ray irradiations; preclinical experiments involved the application of MRT to animal tumor patients, to animals of larger size than laboratory rodents, using larger radiation field sizes, and irradiating in a real-time scenario comparable to the one planned for human patients. (4) To foster the specific purpose of radiosurgical MRT of tumor patients at the ESRF that required development of new, specific state of the art modalities and tools for treatment planning, dosimetry, dose calculation, patient positioning and, of particular importance, redundant levels of patient safety. Just as she was about to take responsibility as principal investigator for a prestigious international European project on a new radiation therapy, death called Elke in.

In-field and out-of-field microdosimetric characterisation of a 62 MeV proton beam at CATANA

PURPOSE

A 5 and 10 μm thin silicon on insulator (SOI) 3D mushroom microdosimeter was used to characterize both the in-field and out-of-field of a 62 MeV proton beam.

METHODS

The SOI mushroom microdosimeter consisted of an array of cylindrical sensitive volumes (SVs), developed by the Centre for Medical Radiation Physics, University of Wollongong, was irradiated with 62 MeV protons at the CATANA (Centro di AdroTerapia Applicazioni Nucleari Avanzate) facility in Catania, Italy, a facility dedicated to the radiation treatment of ocular melanomas. Dose mean lineal energy, ( y D ¯ ), values were obtained at various depths in PMMA along a pristine and spread out Bragg peak (SOBP). The measured microdosimetric spectra at each position were then used as inputs into the modified Microdosimetric Kinetic Model (MKM) to derive the RBE for absorbed dose in a middle of the SOBP 2Gy (RBE_D ). Microdosimetric spectra were obtained with both the 5 and 10 μm 3D SOI microdosimeters, with a focus on the distal part of the BP. The in-field and out-of-field measurement configurations along the Bragg curve were modeled in Geant4 for comparison with experimental results. Lateral out-of-field measurements were performed to study secondary particles' contribution to normal tissue's dose, up to 12 mm from the edge of the beam field, and quality factor and dose equivalent results were obtained.

RESULTS

Comparison between experimental and simulation results showed good agreement between one another for both the pristine and SOBP beams in terms of y D ¯ and RBED. Though a small discrepancy between experiment and simulation was seen at the entrance of the Bragg curve, where experimental results were slightly lower than Geant4. The dose equivalent value measured 12 mm from the edge of the target volume was 1.27 ± 0.15 mSv/Gy with a Q ¯ value of 2.52 ± 0.30, both of which agree within uncertainty with Geant4 simulation.

CONCLUSIONS

These results demonstrate that SOI microdosimeters are an effective tool to predict RBED in-field as well as dose equivalent monitoring out-of-field to provide insight to probability of second cancer generation.

Microdosimetry in low energy proton beam at therapeutic-equivalent fluence rate with silicon 3D-cylindrical microdetectors

In this work we show the first microdosimetry measurements on a low energy proton beam with therapeutic-equivalent fluence rates by using the second generation of 3D-cylindrical microdetectors. The sensors belong to an improved version of a novel silicon-based 3D-microdetector design with electrodes etched inside silicon, which were manufactured at the National Microelectronics Centre (IMB-CNM, CSIC) in Spain. A new microtechnology has been employed using quasi-toroid electrodes of 25μm diameter and a depth of 20μm within the silicon bulk, resulting in a well-defined cylindrical radiation sensitive volume. These detectors were tested at the 18 MeV proton beamline of the cyclotron at the National Accelerator Centre (CNA, Spain). They were assembled into an in-house low-noise readout electronics to assess their performance at a therapeutic-equivalent fluence rate. Microdosimetry spectra of lineal energy were recorded at several proton energies starting from 18 MeV by adding 50μm thick tungsten foils gradually at the exit-window of the cyclotron external beamline, which corresponds to different depths along the Bragg curve. The experimentalyF¯values in silicon cover from (5.7 ± 0.9) to (8.5 ± 0.4) keV μm^-1in the entrance to (27.4 ± 2.3) keV μm^-1in the distal edge. Pulse height energy spectra were crosschecked with Monte Carlo simulations and an excellent agreement was obtained. This work demonstrates the capability of the second generation 3D-microdetectors to assess accurate microdosimetric distributions at fluence rates as high as those used in clinical centers in proton therapy.

IDDRRA: A novel platform, based on Geant4-DNA to quantify DNA damage by ionizing radiation

PURPOSE

This study proposes a novel computational platform that we refer to as IDDRRA (DNA Damage Response to Ionizing RAdiation), which uses Monte Carlo (MC) simulations to score radiation induced DNA damage. MC simulations provide results of high accuracy on the interaction of radiation with matter while scoring the energy deposition based on state-of-the-art physics and chemistry models and probabilistic methods.

METHODS

The IDDRRA software is based on the Geant4-DNA toolkit together with new tools that were developed for the purpose of this study, including a new algorithm that was developed in Python for the design of the DNA molecules. New classes were developed in C++ to integrate the GUI and produce the simulation's output in text format. An algorithm was also developed to analyze the simulation's output in terms of energy deposition, Single Strand Breaks (SSB), Double Strand Breaks (DSB) and Cluster Damage Sites (CDS). Finally, a new tool was developed to implement probabilistic SSB and DSB repair models using MC techniques.

RESULTS

This article provides the first benchmarks that the user of the IDDRRA tool can use to validate the functionality of the software as well as to provide a starting point to produce different types of DNA simulations. These benchmarks incorporate different kind of particles (e-, e+, protons, electron spectrum) and DNA molecules.

CONCLUSION

We have developed the IDDRRA tool and demonstrated its use to study various aspects of the modeling and simulation of a DNA irradiation experiment. The tool is expandable and can be expanded by other users with new benchmarks and applications based on the user's needs and experience. New functionality will be added over time, including the quantification of the indirect damage.

Simultaneous Measurements of Dose and Microdosimetric Spectra in a Clinical Proton Beam Using a scCVD Diamond Membrane Microdosimeter

A single crystal chemical vapor deposition (scCVD) diamond membrane-based microdosimetric system was used to perform simultaneous measurements of dose profile and microdosimetric spectra with the Y1 proton passive scattering beamline of the Center of Proton Therapy, Institute Curie in Orsay, France. To qualify the performance of the set-up in clinical conditions of hadrontherapy, the dose, dose rate and energy loss pulse-height spectra in a diamond microdosimeter were recorded at multiple points along depth of a water-equivalent plastic phantom. The dose-mean lineal energy (y¯D) values were computed from experimental data and compared to silicon on insulator (SOI) microdosimeter literature results. In addition, the measured dose profile, pulse height spectra and y¯D values were benchmarked with a numerical simulation using TOPAS and Geant4 toolkits. These first clinical tests of a novel system confirm that diamond is a promising candidate for a tissue equivalent, radiation hard, high spatial resolution microdosimeter in beam quality assurance of proton therapy.

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.

Silicon 3D Microdetectors for Microdosimetry in Hadron Therapy

The present overview describes the evolution of new microdosimeters developed in the National Microelectronics Center in Spain (IMB-CNM, CSIC), ranging from the first ultra-thin 3D diodes (U3DTHINs) to the advanced 3D-cylindrical microdetectors, which have been developed over the last 10 years. In this work, we summarize the design, main manufacture processes, and electrical characterization of these devices. These sensors were specifically customized for use in particle therapy and overcame some of the technological challenges in this domain, namely the low noise capability, well-defined sensitive volume, high spatial resolution, and pile-up robustness. Likewise, both architectures reduce the loss of charge carriers due to trapping effects, the charge collection time, and the voltage required for full depletion compared to planar silicon detectors. In particular, a 3D‒cylindrical architecture with electrodes inserted into the silicon bulk and with a very well‒delimited sensitive volume (SV) mimicked a cell array with shapes and sizes similar to those of mammalian cells for the first time. Experimental tests of the carbon beamlines at the Grand Accélérateur National d’Lourds (GANIL, France) and Centro Nazionale Adroterapia Oncologica (CNAO, Italy) showed the feasibility of the U3DTHINs in hadron therapy beams and the good performance of the 3D‒cylindrical microdetectors for assessing linear energy distributions of clinical beams, with clinical fluence rates of 5 × 10⁷ s⁻¹cm⁻² without saturation. The dose-averaged lineal energies showed a generally good agreement with Monte Carlo simulations. The results indicated that these devices can be used to characterize the microdosimetric properties in hadron therapy, even though the charge collection efficiency (CCE) and electronic noise may pose limitations on their performance, which is studied and discussed herein. In the last 3D‒cylindrical microdetector generation, we considerably improved the CCE due to the microfabrication enhancements, which have led to shallower and steeper dopant profiles. We also summarize the successive microdosimetric characterizations performed with both devices in proton and carbon beamlines.

Experimental validation of an analytical microdosimetric model based on Geant4-DNA simulations by using a silicon-based microdosimeter

PURPOSE

To study the agreement between proton microdosimetric distributions measured with a silicon-based cylindrical microdosimeter and a previously published analytical microdosimetric model based on Geant4-DNA in-water Monte Carlo simulations for low energy proton beams.

METHODS AND MATERIAL

Distributions for lineal energy (y) are measured for four proton monoenergetic beams with nominal energies from 2.0 MeV to 4.5 MeV, with a tissue equivalent proportional counter (TEPC) and a silicon-based microdosimeter. The actual energy for protons traversing the silicon-based microdosimeter is simulated with SRIM. Monoenergetic beams with these energies are simulated with Geant4-DNA code by simulating a water cylinder site of dimensions equal to those of the microdosimeter. The microdosimeter response is calibrated by using the distribution peaks obtained from the TEPC. Analytical calculations fory ¯ F andy ¯ D using our methodology based on spherical sites are also performed choosing the equivalent sphere to be checked against experimental results.

RESULTS

Distributions for y at silicon are converted into tissue equivalent and compared to the Geant4-DNA simulated, yielding maximum deviations of 1.03% fory ¯ F and 1.17% fory ¯ D . Our analytical method generates maximum deviations of 1.29% and 3.33%, respectively, with respect to experimental results.

CONCLUSION

Simulations in Geant4-DNA with ideal cylindrical sites in liquid water produce similar results to the measurements in an actual silicon-based cylindrical microdosimeter properly calibrated. The found agreement suggests the possibility to experimentally verify the calculated clinicaly ¯ D with our analytical method.

Confocal Microscopy Improves 3D Microdosimetry Applied to Nanoporation Experiments Targeting Endoplasmic Reticulum

In the last years, microdosimetric numerical models of cells including intracellular compartments have been proposed, aiming to investigate the poration induced by the application of nanosecond pulsed electric fields (nsPEFs). A limitation of such models was the extremely approximate cell and organelle shapes, leading to an incorrect estimation of the electric field or transmembrane potential distribution in the studied domain. In order to obtain a reliable model of in vitro experiments and a one-to-one comparison between experimental and simulated results, here, a realistic model of 12 human mesenchymal stem cells was built starting from their optical microscopy images where different cell compartments were highlighted. The microdosimetric analysis of the cells group was quantified in terms of electric field and transmembrane potentials (TMPs) induced by an externally applied 10-ns trapezoidal pulse with rise and fall times of 2 ns, with amplitudes ranging from 2 to 30 MV/m. The obtained results showed that the plasma and endoplasmic reticulum (ER) membrane of each cell respond in a different way to the same electric field amplitude, depending on differences in shape, size, and position of the single cell with respect to the applied electric field direction. Therefore, also the threshold for an efficient electroporation is highly different from cell to cell. This difference was quantitatively estimated through the cumulative distribution function of the pore density for the plasma and ER membrane of each cell, representing the probability that a certain percentage of membrane has reached a specific value of pore density. By comparing the dose-response curves resulted from the simulations and those from the experimental study of De Menorval et al. (2016), we found a very good matching of results for plasma and ER membrane when 2% of the porated area is considered sufficient for permeabilizing the membrane. This result is worth of noting as it highlights the possibility to effectively predict the behavior of a cell (or of a population of cells) exposed to nsPEFs. Therefore, the microdosimetric realistic model described here could represent a valid tool in setting up more efficient and controlled electroporation protocols.

A new approach to modeling the microdosimetry of proton therapy beams

INTRODUCTION

To revisit the formulation of the mean chord length in microdosimetry and replace it by the particle mean free path appropriate for modelings in radiobiology.

METHODS

We perform a collision-by-collision followed by event-by-event Geant4 Monte Carlo simulation and calculate double-averaged stepping length, 〈〈l〉〉, for a range of target sizes from mm down to μm and depth in water. We consider 〈〈l〉〉 to represent the particle mean free path.

RESULTS

We show that 〈〈l〉〉 continuously drops as a function of depth and asymptotically saturates to a minimum value in low energies, where it exhibits a universal scaling behavior, independent of particle nominal beam energy. We correlate 〈〈l〉〉 to linear density of DNA damage, complexities of initial lethal lesions and illustrate a relative difference between predictive RBEs in model calculations using mean chord length vs the proposed mean free path. We demonstrate consistency between rapid increase in RBE within and beyond the Bragg peak and 〈〈l〉〉, a decreasing function of depth.

DISCUSSION AND CONCLUSION

An interplay between localities in imparted energy at nanometer scale and subsequent physio-chemical processes, causalities and pathways in DNA damage requires substitution of geometrical chord length of cell nuclei by mean-free path of proton and charged particles to account for a mean distance among sequential collisions in DNA materials. To this end, the event averaging over cell volume in the current microdosimetry formalism must be superseded by the collision averaging scored within the volume. The former, is fundamentally a global attribute of the cell nuclei surfaces and boundaries and is characterized by their membrane diameters, hence such global indices are not appropriate to quantitatively represent the radiobiological strength of the particles and their RBE variabilities that is associated with the sensitivities to local structure of the collisions and their spatio-temporal collective patterns in DNA materials.

A Model for Estimating Dose-Rate Effects on Cell-Killing of Human Melanoma after Boron Neutron Capture Therapy

Boron neutron capture therapy (BNCT) is a type of radiation therapy for eradicating tumor cells through a ¹⁰B(n,α)⁷Li reaction in the presence of ¹⁰B in cancer cells. When delivering a high absorbed dose to cancer cells using BNCT, both the timeline of ¹⁰B concentrations and the relative long dose-delivery time compared to photon therapy must be considered. Changes in radiosensitivity during such a long dose-delivery time can reduce the probability of tumor control; however, such changes have not yet been evaluated. Here, we propose an improved integrated microdosimetric-kinetic model that accounts for changes in microdosimetric quantities and dose rates depending on the ¹⁰B concentration and investigate the cell recovery (dose-rate effects) of melanoma during BNCT irradiation. The integrated microdosimetric–kinetic model used in this study considers both sub-lethal damage repair and changes in microdosimetric quantities during irradiation. The model, coupled with the Monte Carlo track structure simulation code of the Particle and Heavy Ion Transport code System, shows good agreement with in vitro experimental data for acute exposure to ⁶⁰Co γ-rays, thermal neutrons, and BNCT with ¹⁰B concentrations of 10 ppm. This indicates that microdosimetric quantities are important parameters for predicting dose-response curves for cell survival under BNCT irradiations. Furthermore, the model estimation at the endpoint of the mean activation dose exhibits a reduced impact of cell recovery during BNCT irradiations with high linear energy transfer (LET) compared to ⁶⁰Co γ-rays irradiation with low LET. Throughout this study, we discuss the advantages of BNCT for enhancing the killing of cancer cells with a reduced dose-rate dependency. If the neutron spectrum and the timelines for drug and dose delivery are provided, the present model will make it possible to predict radiosensitivity for more realistic dose-delivery schemes in BNCT irradiations.

Microdosimetric quantities of an accelerator-based neutron source used for boron neutron capture therapy measured using a gas-filled proportional counter

Boron neutron capture therapy (BNCT) is an emerging radiation treatment modality, exhibiting the potential to selectively destroy cancer cells. Currently, BNCT is conducted using a nuclear reactor. However, the future trend is to move toward an accelerator-based system for use in hospital environments. A typical BNCT radiation field has several different types of radiation. The beam quality should be quantified to accurately determine the dose to be delivered to the target. This study utilized a tissue equivalent proportional counter (TEPC) to measure microdosimetric and macrodosimetric quantities of an accelerator-based neutron source. The micro- and macro-dosimetric quantities measured with the TEPC were compared with those obtained via the the particle and heavy ion transport code system (PHITS) Monte Carlo simulation. The absorbed dose from events >20 keV/μm measured free in air for a 1-h irradiation was calculated as 1.31 ± 0.02 Gy. The simulated result was 1.41 ± 0.07 Gy. The measured and calculated values exhibit good agreement. The relative biological effectiveness (RBE) that was evaluated from the measured microdosimetric spectrum was calculated as 3.7 ± 0.02, similar to the simulated value of 3.8 ± 0.1. These results showed the PHITS Monte Carlo simulation can simulate both micro- and macro-dosimetric quantities accurately. The RBE was calculated using a single-response function, and the results were compared with those of several other institutes that used a similar method. However, care must be taken when using such a single-response function for clinical application, as it is only valid for low doses. For clinical dose ranges (i.e., high doses), multievent distribution inside the target needs to be considered.

Radiation Track Structure: How the Spatial Distribution of Energy Deposition Drives Biological Response

Ionising radiation is incredibly effective at causing biological effects. This is due to the unique way energy is deposited along highly structured tracks of ionisation and excitation events, which results in correlation with sites of DNA damage from the nanometre to the micrometre scale. Correlation of these events along the track on the nanometre scale results in clustered damage, which not only results in the formation of DNA double-strand breaks (DSB), but also more difficult to repair complex DSB, which include additional damage within a few base pairs. The track structure varies significantly with radiation quality and the increase in relative biological effectiveness observed with increasing linear energy transfer in part corresponds to an increase in the probability and complexity of clustered DNA damage produced. Likewise, correlation over larger scales, associated with packing of DNA and associated chromosomes within the cell nucleus, can also have a major impact on the biological response. The proximity of the correlated damage along the track increases the probability of miss-repair through pairwise interactions resulting in an increase in probability and complexity of DNA fragments/deletions, mutations and chromosomal rearrangements. Understanding the mechanisms underlying the biological effectiveness of ionising radiation can provide an important insight into ways of increasing the efficacy of radiotherapy, as well as the risks associated with exposure. This requires a multi-scale approach for modelling, not only considering the physics of the track structure from the millimetre scale down to the nanometre scale, but also the structural packing of the DNA within the nucleus, the resulting chemistry in the context of the highly reactive environment of the nucleus, together with the subsequent biological response.

The influence of Geant4-DNA toolkit parameters on electron microdosimetric track structure

The influence of different physical process factors on tracks of low-energy electrons in liquid water was analyzed and evaluated based on the Geant4-DNA toolkit of Geant4 version 10.4, and it provides theoretical support for obtaining the basic parameters of microdosimetry concerned with radiotherapy and radiation protection. According to the characteristics of different models, five physics constructors of Geant4-DNA toolkit were selected to simulate monoenergetic electrons in microscopic scale. Details of track structure of different Geant4-DNA physics constructors were compared, including total number of interaction processes, number and energy percentage of excitation and ionization; analyzing the impacts of mean lineal energy of several factors, including Geant4-DNA physics constructors, initial energy, radius of scoring spheres, interaction processes and cut-off energy. Firstly, 'G4EmDNAPhysics' (hereinafter referred to as 'dna') is well consistent with 'G4EmDNAPhysics_option 2' (hereinafter referred to as 'option 2'), and 'G4EmDNAPhysics_option 4' (hereinafter referred to as 'option 4') is well consistent with 'G4EmDNAPhysics_option 5' (hereinafter referred to as 'option 5'); secondly, there are differences for the information of track structure and mean lineal energy between 'option 2' 'option 4' and 'G4EmDNAPhysics_option 6' (hereinafter referred to as 'option 6'); thirdly, the influence of the model on the mean lineal energy decreases with the increase of the radius of the scoring spheres, whereas mean lineal energy increases as the tracking cut increases. Several alternative discrete physics constructors of Geant4-DNA are comprehensively discussed overlaying multiple perspectives under different conditions in this work.

Assessment of Radio-Induced Damage in Endothelial Cells Irradiated with 40 kVp, 220 kVp, and 4 MV X-rays by Means of Micro and Nanodosimetric Calculations

The objective of this work was to study the differences in terms of early biological effects that might exist between different X-rays energies by using a mechanistic approach. To this end, radiobiological experiments exposing cell monolayers to three X-ray energies were performed in order to assess the yields of early DNA damage, in particular of double-strand breaks (DSBs). The simulation of these irradiations was set in order to understand the differences in the obtained experimental results. Hence, simulated results in terms of microdosimetric spectra and early DSB induction were analyzed and compared to the experimental data. Human umbilical vein endothelial cells (HUVECs) were irradiated with 40, 220 kVp, and 4 MV X-rays. The Geant4 Monte Carlo simulation toolkit and its extension Geant4-DNA were used for the simulations. Microdosimetric calculations aiming to determine possible differences in the variability of the energy absorbed by the irradiated cell population for those photon spectra were performed on 10,000 endothelial cell nuclei representing a cell monolayer. Nanodosimetric simulations were also carried out using a computation chain that allowed the simulation of physical, physico-chemical, and chemical stages on a single realistic endothelial cell nucleus model including both heterochromatin and euchromatin. DNA damage was scored in terms of yields of prompt DSBs per Gray (Gy) and per giga (10⁹) base pair (Gbp) and DSB complexity was derived in order to be compared to experimental data expressed as numbers of histone variant H2AX (γ-H2AX) foci per cell. The calculated microdosimetric spread in the irradiated cell population was similar when comparing between 40 and 220 kVp X-rays and higher when comparing with 4 MV X-rays. Simulated yields of induced DSB/Gy/Gbp were found to be equivalent to those for 40 and 220 kVp but larger than those for 4 MV, resulting in a relative biological effectiveness (RBE) of 1.3. Additionally, DSB complexity was similar between the considered photon spectra. Simulated results were in good agreement with experimental data obtained by IRSN (Institut de radioprotection et de sûreté nucléaire) radiobiologists. Despite differences in photon energy, few differences were observed when comparing between 40 and 220 kVp X-rays in microdosimetric and nanodosimetric calculations. Nevertheless, variations were observed when comparing between 40/220 kVp and 4 MV X-rays. Thanks to the simulation results, these variations were able to be explained by the differences in the production of secondary electrons with energies below 10 keV.

Characterization of radiochromic films as a micrometer-resolution dosimeter by confocal Raman spectroscopy

PURPOSE

Micrometer spatial resolution dosimetry has become inevitable for advanced radiotherapy techniques. A new approach using radiochromic films was developed to measure a radiation dose at a micrometer spatial resolution by confocal Raman spectroscopy.

METHODS

The commercial radiochromic films (RCF), EBT3 and EBT-XD, were irradiated with known doses using 50, 100, 200, and 300 kVp, and 6-MV x rays. The dose levels ranged from 0.3 to 50 Gy. The Raman mapping technique developed in our early study was used to readout an area of 100 × 100 µm² on RCF with improved lateral and depth resolutions with confocal Raman spectrometry. The variation in Raman spectra of C-C-C deformation and C≡C stretching modes of diacetylene polymers around 676 and 2060 cm^-1 , respectively, as a function of therapeutic x-ray doses, was measured. The single peak (SP) of C≡C and the peak ratio (PR) of C≡C band height to C-C-C band height with a spatial resolution of 10 µm on both types of RCF were evaluated, averaged, and plotted as a function of dose. An achievable spatial resolution, clinically useful dose range, dosimetric sensitivity, dose uniformity, and postirradiation stability as well as the orientation, energy, and dose rate dependence, of both types of RCFs, were characterized by the technique developed in this study.

RESULTS

A spatial resolution on RCF achieved by SP and PR methods was ~4.5 and ~2.9 µm, respectively. Raman spectroscopy data showed dose nonuniformity of ~11% in SP method and <3% in PR method. The SP method provided dose ranges of up to ~10 and ~20 Gy for EBT3 and EBT-XD films, respectively while the PR method up to ~30 and ~50 Gy. The PR method diminished the orientation effect. The percent difference between landscape and portrait orientations for the EBT3 and the EBT-XD films at 4 Gy had an acceptable level of 1.2% and 2.4%, respectively. With both SP and PR methods, the EBT3 and the EBT-XD films showed weak energy (within ~10% and ~3% for SP and PR methods, respectively) and dose rate dependence (within ~5% and ~3% for SP and PR methods, respectively) and had a stable response after 24-h postirradiation.

CONCLUSIONS

A technique for micrometer-resolution dosimetry was successfully developed by detecting radiation-induced Raman shift on EBT3 and EBT-XD. Both types of RCFs were suitable for micrometer-resolution dosimetry using CRS. With CRS both lateral and depth resolutions on RCF were improved. The PR method provided superior characteristics in dose uniformity, dose ranges, orientation dependence, and laser effect for both types of RCFs. The overall dosimetric characteristics of the RCFs determined by this technique were similar to those known by optical density scanning. The CRS with the PR method is advantageous over other the traditional scanning systems as a spatial resolution of <10 µm on RCF can be achieved with less deviations.

Technical Note: Effect of magnetic fields on the microdosimetry of carbon-ion beams

PURPOSE

To investigate the influence of magnetic fields on the microdosimetry of carbon-ion beams and the scaling effect of tissue equivalent proportional counter (TEPC) defined as the change of energy deposition in a TEPC from that in a microscopic scale region of interest due to the presence of a magnetic field in combination with the TEPC larger physical dimensions.

METHODS

Geant4-based Monte Carlo simulations were conducted to calculate the microdosimetric quantities for carbon-ion beams with different initial energies (10-290 MeV/u) under magnetic fields of various strengths (0.5-3 T). The calculations were performed for a 1 μm spherical volume made of tissue, and for spherical TEPCs of 1 and 10 mm in diameter. Then, values of dose-averaged lineal energy (y_D ) were acquired for the different scenarios to analyze the effect of magnetic fields on the microdosimetry of carbon-ion beams and the scaling effect of TEPC.

RESULTS

The y_D values and lineal energy spectra in the 1 μm spherical tissue volume for the scenarios without magnetic field and with magnetic fields of different strengths and directions remained nearly the same for the various energy carbon-ion beams. However, compared with those of the 1 μm spherical tissue volume, an increase of of y_D values and an obvious shift of the lineal energy spectra for the TEPCs of 1 and 10 mm in diameter under magnetic fields were found.

CONCLUSIONS

The application of magnetic fields under 3 T has no significant influence on the microdosimetric results of carbon-ion beams. However, there is definitely a scaling effect when using TEPC for microdosimetric study, which alters the reading of TEPC in the presence of magnetic fields. Novel methods to correct the reading of TEPC or scaling effect-resistant microdosimetric measurement detectors are urgently needed to perform experimental microdosimetric studies under magnetic fields.

MPEXS-DNA, a new GPU-based Monte Carlo simulator for track structures and radiation chemistry at subcellular scale

Purpose

Track structure simulation codes can accurately reproduce the stochastic nature of particle–matter interactions in order to evaluate quantitatively radiation damage in biological cells such as DNA strand breaks and base damage. Such simulations handle large numbers of secondary charged particles and molecular species created in the irradiated medium. Every particle and molecular species are tracked step‐by‐step using a Monte Carlo method to calculate energy loss patterns and spatial distributions of molecular species inside a cell nucleus with high spatial accuracy. The Geant4‐DNA extension of the Geant4 general‐purpose Monte Carlo simulation toolkit allows for such track structure simulations and can be run on CPUs. However, long execution times have been observed for the simulation of DNA damage in cells. We present in this work an improvement of the computing performance of such simulations using ultraparallel processing on a graphical processing unit (GPU).

Methods

A new Monte Carlo simulator named MPEXS‐DNA, allowing high computing performance by using a GPU, has been developed for track structure and radiolysis simulations at the subcellular scale. MPEXS‐DNA physics and chemical processes are based on Geant4‐DNA processes available in Geant4 version 10.02 p03. We have reimplemented the Geant4‐DNA process codes of the physics stage (electromagnetic processes of charged particles) and the chemical stage (diffusion and chemical reactions for molecular species) for microdosimetry simulation by using the CUDA language. MPEXS‐DNA can calculate a distribution of energy loss in the irradiated medium caused by charged particles and also simulate production, diffusion, and chemical interactions of molecular species from water radiolysis to quantitatively assess initial damage to DNA. The validation of MPEXS‐DNA physics and chemical simulations was performed by comparing various types of distributions, namely the radial dose distributions for the physics stage, and the G‐value profiles for each chemical product and their linear energy transfer dependency for the chemical stage, to existing experimental data and simulation results obtained by other simulation codes, including PARTRAC.

Results

For physics validation, radial dose distributions calculated by MPEXS‐DNA are consistent with experimental data and numerical simulations. For chemistry validation, MPEXS‐DNA can also reproduce G‐value profiles for each molecular species with the same tendency as existing experimental data. MPEXS‐DNA also agrees with simulations by PARTRAC reasonably well. However, we have confirmed that there are slight discrepancies in G‐value profiles calculated by MPEXS‐DNA for molecular species such as H₂ and H₂O₂ when compared to experimental data and PARTRAC simulations. The differences in G‐value profiles between MPEXS‐DNA and PARTRAC are caused by the different chemical reactions considered. MPEXS‐DNA can drastically boost the computing performance of track structure and radiolysis simulations. By using NVIDIA's GPU devices adopting the Volta architecture, MPEXS‐DNA has achieved speedup factors up to 2900 against Geant4‐DNA simulations with a single CPU core.

Conclusion

The MPEXS‐DNA Monte Carlo simulation achieves similar accuracy to Monte Carlo simulations performed using other codes such as Geant4‐DNA and PARTRAC, and its predictions are consistent with experimental data. Notably, MPEXS‐DNA allows calculations that are, at maximum, 2900 times faster than conventional simulations using a CPU.

Investigating energy deposition in glandular tissues for mammography using multiscale Monte Carlo simulations

PURPOSE

To investigate energy deposition in glandular tissues of the breast on macro- and microscopic length scales in the context of mammography.

METHODS

Multiscale mammography models of breasts are developed, which include segmented, voxelized macroscopic tissue structure as well as nine regions of interest (ROIs) embedded throughout the breast tissue containing explicitly-modelled cells. Using a 30 kVp Mo/Mo spectrum, Monte Carlo (MC) techniques are used to calculate dose to ∼mm voxels containing glandular and/or adipose tissues, as well as energy deposition on cellular length scales. ROIs consist of at least 1000 mammary epithelial cells and ∼200 adipocytes; specific energy (energy imparted per unit mass; stochastic analogue of the absorbed dose) is calculated within mammary epithelial cell nuclei.

RESULTS

Macroscopic dose distributions within segmented breast tissue demonstrate considerable variation in energy deposition depending on depth and tissue structure. Doses to voxels containing glandular tissue vary between ∼0.1 and ∼4 times the mean glandular dose (MGD, averaged over the entire breast). Considering microscopic length scales, mean specific energies for mammary epithelial cell nuclei are ∼30% higher than the corresponding glandular voxel dose. Additionally, due to the stochastic nature of radiation, there is considerable variation in energy deposition throughout a cell population within a ROI: for a typical glandular voxel dose of 4 mGy, the standard deviation of the specific energy for mammary epithelial cell nuclei is 85% relative to the mean. Thus, for a glandular voxel dose of 4 mGy at the centre of the breast, corresponding mammary epithelial cell nuclei will receive specific energies up to ∼9 mGy (considering the upper end of the 1σ standard deviation of the specific energy), while a ROI located 2 cm closer to the radiation source will receive specific energies up to ∼40 mGy. Energy deposition within mammary epithelial cell nuclei is sensitive to cell model details including cellular elemental compositions and nucleus size, underlining the importance of realistic cellular models.

CONCLUSIONS

There is considerable variation in energy deposition on both macro- and microscopic length scales for mammography, with glandular voxel doses and corresponding cell nuclei specific energies many times higher than the MGD in parts of the breast. These results should be considered for radiation-induced cancer risk evaluation in mammography which has traditionally focused on a single metric such as the MGD.

Effects of mucus thickness and goblet cell hyperplasia on microdosimetric quantities characterizing the bronchial epithelium upon radon exposure

PURPOSE

The most exposed tissue upon radon exposure is the bronchial epithelium where goblet cells serve as responsive and adaptable front-line defenders. They can rapidly produce a vast amount of mucus, and can change in number, in response to airway insults. The objective of the present study is to quantify the effects of mucus discharge and goblet cell hyperplasia on the microscopic dose consequences of macroscopic radon exposures.

METHODS

For this purpose, computational models of the bronchial epithelium and alpha-particle transport have been prepared and applied to quantify the hits received and doses absorbed by cell nuclei in case of different mucus thicknesses and goblet cell number.

RESULTS AND CONCLUSIONS

Both mucus discharge and induction of goblet cell hyperplasia reduce radiation burden at the cellular level, and as such they both can be considered as radioadaptive responses to radon exposure. As compared to basal cell hyperplasia, goblet cell hyperplasia is more effective in reducing the microscopic dose consequences of a given macroscopic exposure. Such changes in exposure geometry highlight the need for improvements in the application of biokinetic and dosimetry models for incorporated radionuclides as well as the dose and dose rate effectiveness factor.