Comment on "Reproducibility study of Monte Carlo simulations for nanoparticle dose enhancement and biological modeling of cell survival curves" by Veltenet al[Biomed Phys Eng Express 2023;9:045004].

This comment highlights two methodological issues with the recent article by Velten et al [Biomed Phys Eng Express 2023;9:045004]: First, the approach taken in this work with a local effect model (LEM) in 2D leads to a significant overstimation of the number of radiation-induced lesions. This results in order of magnitude smaller predicted survival rates compared to the conventional LEM. Second, the dose without nanoparticles is used as the 'macroscopic dose' against which cell survival is plotted. However, for the considered gold concentrations, the average absorbed dose under secondary particle equilibrium is between 2 and 20 times higher with nanoparticles than without.

Correction: Belchior et al. Repair Kinetics of DSB-Foci Induced by Proton and α-Particle Microbeams of Different Energies. Life 2022, 12, 2040.

In the original publication [...].

Repair Kinetics of DSB-Foci Induced by Proton and α-Particle Microbeams of Different Energies.

In this work, the induction and repair of radiation-induced 53BP1 foci were studied in human umbilical vein endothelial cells irradiated at the PTB microbeam with protons and α-particles of different energies. The data were analyzed in terms of the mean number of 53BP1 foci induced by the different ion beams. The number of 53BP1 foci found at different times post-irradiation suggests that the disappearance of foci follows first order kinetics. The mean number of initially produced foci shows the expected increase with LET. The most interesting finding of this work is that the absolute number of persistent foci increases with LET but not their fraction. Furthermore, protons seem to produce more persistent foci as compared to α-particles of even higher LET. This may be seen as experimental evidence that protons may be more effective in producing severe DNA lesions, as was already shown in other work, and that LET may not be the best suited parameter to characterize radiation quality.

Investigation into the foundations of the track-event theory of cell survival and the radiation action model based on nanodosimetry.

This work aims at elaborating the basic assumptions behind the "track-event theory" (TET) and its derivate "radiation action model based on nanodosimetry" (RAMN) by clearly distinguishing between effects of tracks at the cellular level and the induction of lesions in subcellular targets. It is demonstrated that the model assumptions of Poisson distribution and statistical independence of the frequency of single and clustered DNA lesions are dispensable for multi-event distributions because they follow from the Poisson distribution of the number of tracks affecting the considered target volume. It is also shown that making these assumptions for the single-event distributions of the number of lethal and sublethal lesions within a cell would lead to an essentially exponential dose dependence of survival for practically relevant values of the absorbed dose. Furthermore, it is elucidated that the model equation used for consideration of repair within the TET is based on the assumption that DNA lesions induced by different tracks are repaired independently. Consequently, the model equation is presumably inconsistent with the model assumptions and requires an additional model parameter. Furthermore, the methodology for deriving model parameters from nanodosimetric properties of particle track structure is critically assessed. Based on data from proton track simulations it is shown that the assumption of statistically independent targets leads to the prediction of negligible frequency of clustered DNA damage. An approach is outlined how track structure could be considered in determining the model parameters, and the implications for TET and RAMN are discussed.

Water acting as a catalyst for electron-driven molecular break-up of tetrahydrofuran.

Low-energy electron-induced reactions in hydrated molecular complexes are important in various fields ranging from the Earth's environment to radiobiological processes including radiation therapy. Nevertheless, our understanding of the reaction mechanisms in particular in the condensed phase and the role of water in aqueous environments is incomplete. Here we use small hydrogen-bonded pure and mixed dimers of the heterocyclic molecule tetrahydrofuran (THF) and water as models for biochemically relevant systems. For electron-impact-induced ionization of these dimers, a molecular ring-break mechanism is observed, which is absent for the THF monomer. Employing coincident fragment ion mass and electron momentum spectroscopy, and theoretical calculations, we find that ionization of the outermost THF orbital initiates significant rearrangement of the dimer structure increasing the internal energy and leading to THF ring-break. These results demonstrate that the local environment in form of hydrogen-bonded molecules can considerably affect the stability of molecular covalent bonds.

Roadmap for metal nanoparticles in radiation therapy: current status, translational challenges, and future directions.

This roadmap outlines the potential roles of metallic nanoparticles (MNPs) in the field of radiation therapy. MNPs made up of a wide range of materials (from Titanium, Z = 22, to Bismuth, Z = 83) and a similarly wide spectrum of potential clinical applications, including diagnostic, therapeutic (radiation dose enhancers, hyperthermia inducers, drug delivery vehicles, vaccine adjuvants, photosensitizers, enhancers of immunotherapy) and theranostic (combining both diagnostic and therapeutic), are being fabricated and evaluated. This roadmap covers contributions from experts in these topics summarizing their view of the current status and challenges, as well as expected advancements in technology to address these challenges.

Determining dose enhancement factors of high-Z nanoparticles from simulations where lateral secondary particle disequilibrium exists.

Nanoparticles (NPs) containing high atomic number (high-Z) materials have been shown to enhance the radiobiological effectiveness of ionizing radiation. This effect is often attributed to an enhancement of the absorbed dose in the vicinity of the NPs, based on Monte Carlo simulations that show a significant local enhancement of the energy deposition on the microscopic scale. The results of such simulations may be significantly biased and lead to a severe overestimation of the dose enhancement if the condition of secondary particle equilibrium is not met in the simulation setup. This current work shows an approach to estimate a 'realistic' dose enhancement from the results of such biased simulations which is based on published photon interaction data and provides a way for correcting biased results.

From Energy Deposition of Ionizing Radiation to Cell Damage Signaling: Benchmarking Simulations by Measured Yields of Initial DNA Damage after Ion Microbeam Irradiation.

Advances in accelerator technology, which have enabled conforming radiotherapy with charged hadronic species, have brought benefits as well as potential new risks to patients. To better understand the effects of ionizing radiation on tumor and surrounding tissue, it is important to investigate and quantify the relationship between energy deposition at the nanometric scale and the initial biological events. Monte Carlo track structure simulation codes provide a powerful tool for investigating this relationship; however, their success and reliability are dependent on their improvement and development accordingly to the dedicated biological data to which they are challenged. For this aim, a microbeam facility that allows for fluence control, down to one ion per cell nucleus, was used to evaluate relative frequencies of DNA damage after interaction between the incoming ion and DNA according to radiation quality. Primary human cells were exposed to alpha particles of three different energies with respective linear energy transfers (LETs) of approximately 36, 85 or 170 keV·µm-1 at the cells' center position, or to protons (19 keV·µm-1). Statistical evaluation of nuclear foci formation (53BP1/γ-H2AX), observed using immunofluorescence and related to a particle traversal, was undertaken in a large population of cell nuclei. The biological results were adjusted to consider the factors that drive the experimental uncertainties, then challenged with results using Geant4-DNA code modeling of the ionizing particle interactions on a virtual phantom of the cell nucleus with the same mean geometry and DNA density as the cells used in our experiments. Both results showed an increase of relative frequencies of foci (or simulated DNA damage) in cell nuclei as a function of increasing LET of the traversing particles, reaching a quasi-plateau when the LET exceeded 80-90 keV·µm-1. For the LET of an alpha particle ranging from 80-90 to 170 keV·µm-1, 10-30% of the particle hits did not lead to DNA damage inducing 53BP1 or γ-H2AX foci formation.

INVESTIGATION INTO THE PROBABILITY FOR MISCOUNTING IN FOCI-BASED ASSAYS.

When early radiation damage to biological systems is studied based on the formation of foci at the location of DNA double-strand breaks, the foci observed in irradiated cells either may be induced by ionizing radiation (IR) interactions or they may be due to other causes that lead to observation of foci also in unirradiated cells. Generally, to take account of the latter, additional samples are taken where the exposure to IR is skipped in the protocol. The data analysis relies on statistical independence of the frequency distributions of background and radiation-induced foci. In microscopy, however, the observed spatial patterns of foci are 2D projections of the spatial distributions of foci in the observed cell nuclei. This may lead to missing foci when scoring their number, particularly if projections of foci overlap or coincide. This paper investigates to what extent the statistical independence of the frequency distribution of the number of foci coming from IR interaction or other causes is compromised by foci overlapping.

CONSIDERATIONS ON THE TARGET SIZE IN A WALL-LESS NANODOSIMETER.

In nanodosimetry, the ionization component of charged particle track structure is characterized by measuring the frequency distribution of ionizations in target volumes that simulate nanometric sites in liquid water. For the Ion Counter nanodosimeter at PTB, the sensitive volume is defined by the electrical field and the extraction aperture. In this paper, a procedure is presented to define a cylindrical effective measurement target based on the second moments of the detection efficiency map. An analytical model of the efficiency map is developed to investigate the dependence of the simulated site size on the nanodosimeter's operating parameters. Within the limits of the simplifying assumptions, the model gives a reasonable approximation of the efficiency map.

PROSPECTS FOR METROLOGY RELATED TO BIOLOGICAL RADIATION EFFECTS OF ION BEAMS.

In recent years, several approaches have been proposed to provide an understanding of the enhanced relative biological effectiveness of ion beams based on multi-scale models of their radiation effects. Among these, the BioQuaRT project was the only one which focused on developing metrology for a multi-scale characterization of particle track structure. The progress made within the BioQuaRT project has motivated the formation of a department 'Radiation Effects' at PTB dedicated to metrological research on ionizing radiation effects. This paper gives an overview of the department's present research directions and shortly discusses ideas for the future development of metrology related to biological effects of ion beams that are based on a stakeholder consultation.

Analysis of Radiation-Induced Chromosomal Aberrations on a Cell-by-Cell Basis after Alpha-Particle Microbeam Irradiation: Experimental Data and Simulations.

There is a continued need for further clarification of various aspects of radiation-induced chromosomal aberration, including its correlation with radiation track structure. As part of the EMRP joint research project, Biologically Weighted Quantities in Radiotherapy (BioQuaRT), we performed experimental and theoretical analyses on chromosomal aberrations in Chinese hamster ovary cells (CHO-K1) exposed to α particles with final energies of 5.5 and 17.8 MeV (absorbed doses: ∼2.3 Gy and ∼1.9 Gy, respectively), which were generated by the microbeam at the Physikalisch-Technische Bundesanstalt (PTB) in Braunschweig, Germany. In line with the differences in linear energy transfer (approximately 85 keV/µm for 5.5 MeV and 36 keV/µm for 17.8 MeV α particles), the 5.5 MeV α particles were more effective than the 17.8 MeV α particles, both in terms of the percentage of aberrant cells (57% vs. 33%) and aberration frequency. The yield of total aberrations increased by a factor of ∼2, although the increase in dicentrics plus centric rings was less pronounced than in acentric fragments. The experimental data were compared with Monte Carlo simulations based on the BIophysical ANalysis of Cell death and chromosomal Aberrations model (BIANCA). This comparison allowed interpretation of the results in terms of critical DNA damage [cluster lesions (CLs)]. More specifically, the higher aberration yields observed for the 5.5 MeV α particles were explained by taking into account that, although the nucleus was traversed by fewer particles (nominally, 11 vs. 25), each particle was much more effective (by a factor of ∼3) at inducing CLs. This led to an increased yield of CLs per cell (by a factor of ∼1.4), consistent with the increased yield of total aberrations observed in the experiments.

Double differential cross sections for proton induced electron emission from molecular analogues of DNA constituents for energies in the Bragg peak region.

For track structure simulations in the Bragg peak region, measured electron emission cross sections of DNA constituents are required as input for developing parameterized model functions representing the scattering probabilities. In the present work, double differential cross sections were measured for the electron emission from vapor-phase pyrimidine, tetrahydrofuran, and trimethyl phosphate that are structural analogues to the base, the sugar, and the phosphate residue of the DNA, respectively. The range of proton energies was from 75 keV to 135 keV, the angles ranged from 15° to 135°, and the electron energies were measured from 10 eV to 200 eV. Single differential and total electron emission cross sections are derived by integration over angle and electron energy and compared to the semi-empirical Hansen-Kocbach-Stolterfoht (HKS) model and a quantum mechanical calculation employing the first Born approximation with corrected boundary conditions (CB1). The CB1 provides the best prediction of double and single differential cross section, while total cross sections can be fitted with semi-empirical models. The cross sections of the three samples are proportional to their total number of valence electrons.

High-resolution (e, 2e + ion) study of electron-impact ionization and fragmentation of methane.

The ionization and fragmentation of methane induced by low-energy (E0 = 66 eV) electron-impact is investigated using a reaction microscope. The momentum vectors of all three charged final state particles, two outgoing electrons, and one fragment ion, are detected in coincidence. Compared to the earlier study [Xu et al., J. Chem. Phys. 138, 134307 (2013)], considerable improvements to the instrumental mass and energy resolutions have been achieved. The fragment products CH4 (+), CH3 (+), CH2 (+), CH(+), and C(+) are clearly resolved. The binding energy resolution of ΔE = 2.0 eV is a factor of three better than in the earlier measurements. The fragmentation channels are investigated by measuring the ion kinetic energy distributions and the binding energy spectra. While being mostly in consistence with existing photoionization studies the results show differences including missing fragmentation channels and previously unseen channels.

An (e, 2e + ion) study of low-energy electron-impact ionization and fragmentation of tetrahydrofuran with high mass and energy resolutions.

We study the low-energy (E0 = 26 eV) electron-impact induced ionization and fragmentation of tetrahydrofuran using a reaction microscope. All three final-state charged particles, i.e., two outgoing electrons and one fragment ion, are detected in triple coincidence such that the momentum vectors and, consequently, the kinetic energies for charged reaction products are determined. The ionic fragments are clearly identified in the experiment with a mass resolution of 1 amu. The fragmentation pathways of tetrahydrofuran are investigated by measuring the ion kinetic energy spectra and the binding energy spectra where an energy resolution of 1.5 eV has been achieved using the recently developed photoemission electron source. Here, we will discuss the fragmentation reactions for the cations C4H8O(+), C4H7O(+), C2H3O(+), C3H6(+), C3H5(+), C3H3(+), CH3O(+), CHO(+), and C2H3(+).

Ionization cross section data of nitrogen, methane, and propane for light ions and electrons and their suitability for use in track structure simulations.

Track structure Monte Carlo simulations are frequently applied in micro- and nanodosimetry to calculate the radiation transport in detail. The use of a well-validated set of cross section data in such simulation codes ensures accurate calculations of transport parameters, such as ionization yields. These cross section data are, however, scarce and often discrepant when measured by different groups. This work surveys literature data on ionization and charge-transfer cross sections of nitrogen, methane, and propane for electrons, protons, and helium particles, focusing on the energy range between 100 keV and 20 MeV. Based on the evaluated data, different models for the parametrization of the cross section data are implemented in the code ptra, developed for simulating proton and alpha particle transport in an ion-counting nanodosimeter. The suitability of the cross section data is investigated by comparing the calculated mean ionization cluster size and energy loss with experimental results in either nitrogen or propane. For protons, generally good agreement between measured and simulated data is found when the Rudd model is used in ptra. For alpha particles, however, a considerable influence of different parametrizations of cross sections for ionization and charge transfer is observed. The ptra code using the charge-transfer data is, nevertheless, successfully benchmarked by the experimental data for the calculation of nanodosimetric quantities, but remaining discrepancies still have to be further investigated (up to 13% lower energy loss and 19% lower mean ionization cluster size than in the experiment). A continuation of this work should investigate data for the energy loss per interaction as well as differential cross section data of nitrogen and propane. Interpolation models for ionization and charge-transfer data are proposed. The Barkas model, frequently used for a determination of the effective charge in the ionization cross section, significantly underestimates both the energy loss (by up to 19%) and the mean ionization cluster size (up to 65%) for alpha particles. It is, therefore, not recommended for particle-track simulations.