Nano-microdosimetric investigation at the therapeutic proton irradiation line of CATANA

Microdosimetric Analysis for Boron Neutron Capture Therapy via Monte Carlo Track Structure Simulation with Modified Lithium Cross-sections

Boron neutron capture therapy (BNCT) is a cellular-level hadron therapy achieving therapeutic effects via the synergistic action of multiple particles, including Lithium, alpha, proton, and photon. However, evaluating the relative biological effectiveness (RBE) in BNCT remains challenging. In this research, we performed a microdosimetric calculation for BNCT using the Monte Carlo track structure (MCTS) simulation toolkit, TOPAS-nBio. This paper reports the first attempt to derive the ionization cross-sections of low-energy (>0.025 MeV/u) Lithium for MCTS simulation based on the effective charge cross-section scalation method and phenomenological double-parameter modification. The fitting parameters λ 1 = 1.101 , λ 2 = 3.486 were determined to reproduce the range and stopping power data from the ICRU report 73. Besides, the lineal energy spectra of charged particles in BNCT were calculated, and the influence of sensitive volume (SV) size was discussed. Condensed history simulation obtained similar results with MCTS when using Micron-SV while overestimating the lineal energy when using Nano-SV. Furthermore, we found that the microscopic boron distribution can significantly affect the lineal energy for Lithium, while the effect for alpha is minimal. Similar results to the published data by PHITS simulation were observed for the compound particles and monoenergetic protons when using micron-SV. Spectra with nano-SV reflected that the different track densities and absorbed doses in the nucleus together result in the dramatic difference in the macroscopic biological response of BPA and BSH. This work and the developed methodology could impact the research fields in BNCT where understanding radiation effects is crucial, such as the treatment planning system, source evaluation, and new boron drug development.

A matter of space: how the spatial heterogeneity in energy deposition determines the biological outcome of radiation exposure

The outcome of the exposure of living organisms to ionizing radiation is determined by the distribution of the associated energy deposition at different spatial scales. Radiation proceeds through ionizations and excitations of hit molecules with an ~ nm spacing. Approaches such as nanodosimetry/microdosimetry and Monte Carlo track-structure simulations have been successfully adopted to investigate radiation quality effects: they allow to explore correlations between the spatial clustering of such energy depositions at the scales of DNA or chromosome domains and their biological consequences at the cellular level. Physical features alone, however, are not enough to assess the entity and complexity of radiation-induced DNA damage: this latter is the result of an interplay between radiation track structure and the spatial architecture of chromatin, and further depends on the chromatin dynamic response, affecting the activation and efficiency of the repair machinery. The heterogeneity of radiation energy depositions at the single-cell level affects the trade-off between cell inactivation and induction of viable mutations and hence influences radiation-induced carcinogenesis. In radiation therapy, where the goal is cancer cell inactivation, the delivery of a homogenous dose to the tumour has been the traditional approach in clinical practice. However, evidence is accumulating that introducing heterogeneity with spatially fractionated beams (mini- and microbeam therapy) can lead to significant advantages, particularly in sparing normal tissues. Such findings cannot be explained in merely physical terms, and their interpretation requires considering the scales at play in the underlying biological mechanisms, suggesting a systemic response to radiation.

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).

Microdosimetry of a therapeutic proton beam with a mini-TEPC and a MicroPlus-Bridge detector for RBE assessment

Proton beams are widely used worldwide to treat localized tumours, the lower entrance dose and no exit dose, thus sparing surrounding normal tissues, being the main advantage of this treatment modality compared to conventional photon techniques. Clinical proton beam therapy treatment planning is based on the use of a general relative biological effectiveness (RBE) of 1.1 along the whole beam penetration depth, without taking into account the documented increase in RBE at the end of the depth dose profile, in the Bragg peak and beyond. However, an inaccurate estimation of the RBE can cause both underdose or overdose, in particular it can cause the unfavourable situation of underdosing the tumour and overdosing the normal tissue just beyond the tumour, which limits the treatment success and increases the risk of complications. In view of a more precise dose delivery that takes into account the variation of RBE, experimental microdosimetry offers valuable tools for the quality assurance of LET or RBE-based treatment planning systems. The purpose of this work is to compare the response of two different microdosimetry systems: the mini-TEPC and the MicroPlus-Bridge detector. Microdosimetric spectra were measured across the 62 MeV spread out Bragg peak of CATANA with the mini-TEPC and with the Bridge microdosimeter. The frequency and dose distributions of lineal energy were compared and the different contributions to the spectra were analysed, discussing the effects of different site sizes and chord length distributions. The shape of the lineal energy distributions measured with the two detectors are markedly different, due to the different water-equivalent sizes of the sensitive volumes: 0.85 μm for the TEPC and 17.3 μm for the silicon detector. When the Loncol's biological weighting function is applied to calculate the microdosimetric assessment of the RBE, both detectors lead to results that are consistent with biological survival data for glioma U87 cells. Both the mini-TEPC and the MicroPlus-Bridge detector can be used to assess the RBE variation of a 62 MeV modulated proton beam along its penetration depth. The microdosimetric assessment of the RBE based on the Loncol's weighting function is in good agreement with radiobiological results when the 10% biological uncertainty is taken into account.

Development of a new microdosimetric biological weighting function for the RBE10 assessment in case of the V79 cell line exposed to ions from 1H to 238U

An improved biological weighting function (IBWF) is proposed to phenomenologically relate microdosimetric lineal energy probability density distributions with the relative biological effectiveness (RBE) for the in vitro clonogenic cell survival (surviving fraction = 10%) of the most commonly used mammalian cell line, i.e. the Chinese hamster lung fibroblasts (V79). The IBWF, intended as a simple and robust tool for a fast RBE assessment to compare different exposure conditions in particle therapy beams, was determined through an iterative global-fitting process aimed to minimize the average relative deviation between RBE calculations and literature in vitro data in case of exposure to various types of ions from ¹H to ²³⁸U. By using a single particle- and energy- independent function, it was possible to establish an univocal correlation between lineal energy and clonogenic cell survival for particles spanning over an unrestricted linear energy transfer range of almost five orders of magnitude (0.2 keV µm^-1 to 15 000 keV µm^-1 in liquid water). The average deviation between IBWF-derived RBE values and the published in vitro data was ∼14%. The IBWF results were also compared with corresponding calculations (in vitro RBE₁₀ for the V79 cell line) performed using the modified microdosimetric kinetic model (modified MKM). Furthermore, RBE values computed with the reference biological weighting function (BWF) for the in vivo early intestine tolerance in mice were included for comparison and to further explore potential correlations between the BWF results and the in vitro RBE as reported in previous studies. The results suggest that the modified MKM possess limitations in reproducing the experimental in vitro RBE₁₀ for the V79 cell line in case of ions heavier than ²⁰Ne. Furthermore, due to the different modelled endpoint, marked deviations were found between the RBE values assessed using the reference BWF and the IBWF for ions heavier than ²H. Finally, the IBWF was unchangingly applied to calculate RBE values by processing lineal energy density distributions experimentally measured with eight different microdosimeters in 19 ¹H and ¹²C beams at ten different facilities (eight clinical and two research ones). Despite the differences between the detectors, irradiation facilities, beam profiles (pristine or spread out Bragg peak), maximum beam energy, beam delivery (passive or active scanning), energy degradation system (water, PMMA, polyamide or low-density polyethylene), the obtained IBWF-based RBE trends were found to be in good agreement with the corresponding ones in case of computer-simulated microdosimetric spectra (average relative deviation equal to 0.8% and 5.7% for ¹H and ¹²C ions respectively).

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.