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

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.

Microdosimetric assessment about proton spread-out Bragg peak at different depths based on the normal human mesh-type cell population model

Objective.In clinical proton therapy, the spread-out Bragg peak (SOBP) is commonly used to fit the target shape. Dose depositions at microscopic sites vary, even with a consistent absorbed dose (D) in SOBP. In the present study, monolayer mesh-type cell population models were developed for microdosimetric assessment at different SOBP depths.Approach.Normal human bronchial epithelial (BEAS-2B) and hepatocytes (L-O2) mesh-type cell models were constructed based on fluorescence tomography images of normal human cells. Particle transport simulation in cell populations was performed coupled with Monte Carlo software PHITS. The relationship between microdosimetry and macrodosimetry of SOBP at different depths was described by analyzing the microdosimetric indicators such as specific energyz,specific energy distributionfz,D,and relative standard deviationσz/z¯within cells. Additionally, the microdosimetric distributions characteristics and their contributing factors were also discussed.Main results.The microscopic dose distribution is strongly influenced by cellular size, shape, and material. The mean specific energyz¯of nucleus and cytoplasm in the cell population is greater than the overall absorbed dose of the cell population model (Dp), with a maximumz¯/Dpof 1.1. The cellular dose distribution is different between the BEAS-2B mesh-type model and its concentric ellipsoid geometry-type model, which difference inz¯is about 10.3% for the nucleus and about 7.5% for the cytoplasm with the SOBP depth of 15 cm. WhenD= 2 Gy, the maximumzof L-O2 nucleus reaches 2.8 Gy andσz/z¯is 5.1% at the mid-depth SOBP (16-18 cm); while the maximumzof the BEAS-2B nucleus reaches 2.2 Gy with only 2.7% ofσz/z¯.Significance.The significant variation of microdosimetric distributions of SOBP different depths indicates the necessity to use mesh-type cell population models, which have the potential to be compared with biological results and build the bio-physical model.

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.

State-of-the-art and potential of experimental microdosimetry in ion-beam therapy

In radiotherapy, radiation-quality should be an expression of the biological and physical characteristics of ionizing radiation such as spatial distribution of ionization or energy deposition. Linear energy transfer (LET) and lineal energy (y) are two descriptors used to quantify the radiation quality. These two quantities are connected and exhibit similar features. In ion-beam therapy (IBT), lineal energy can be measured with microdosimeters, which are specifically designed to cope with the high fluence of particles in clinical beams, while the quantification of LET is generally based on calculations. In pre-clinical studies, microdosimetric spectra are used for the indirect determination of relative biological effectiveness (RBE), e.g., using the microdosimetric kinetic model (MKM) or biophysical response functions. In this context it is important to consider saturation effects, which occur when the highest values of y become less biologically relevant compared to the relative contribution they make to the physical dose. Recent clinical data suggests that local tumor control and normal tissue effects can be linked to macroscopic and microscopic dosimetry parameters. In particular, positive clinical outcomes have been correlated to the highest LET values in the density distribution, and there is no evident link to the saturation discussed above. A systematic collection of microdosimetric information in combination with clinical data in retrospective studies may clarify the role of radiation quality at the highest LET. In the clinical setting, microdosimetry is not widely used yet, despite its potential to be linked with LET by experimentally-determined y values. Through this connection, both play an important role in complex therapy techniques such as intensity modulated particle therapy (IMPT), LET-painting and multi-ion optimization. This review summarizes the current state of microdosimetry for IBT and its potential, as well as research and development needed to make experimental microdosimetry a mature procedure in a clinical context.

Biophysical characterization of collimated and uncollimated fields in pencil beam scanning proton therapy

Objective. The lateral dose fall-off in proton pencil beam scanning (PBS) technique remains the preferred choice for sparing adjacent organs at risk as opposed to the distal edge due to the proton range uncertainties and potentially high relative biological effectiveness. However, because of the substantial spot size along with the scattering in the air and in the patient, the lateral penumbra in PBS can be degraded. Combining PBS with an aperture can result in a sharper dose fall-off, particularly for shallow targets.Approach. The aim of this work was to characterize the radiation fields produced by collimated and uncollimated 100 and 140 MeV proton beams, using Monte Carlo simulations and measurements with a MiniPIX-Timepix detector. The dose and the linear energy transfer (LET) were then coupled with publishedin silicobiophysical models to elucidate the potential biological effects of collimated and uncollimated fields.Main results. Combining an aperture with PBS reduced the absorbed dose in the lateral fall-off and out-of-field by 60%. However, the results also showed that the absolute frequency-averaged LET (LET_F) values increased by a maximum of 3.5 keVμm^-1in collimated relative to uncollimated fields, while the dose-averaged LET (LET_D) increased by a maximum of 7 keVμm^-1. Despite the higher LET values produced by collimated fields, the predicted DNA damage yields remained lower, owing to the large dose reduction.Significance. This work demonstrated the dosimetric advantages of combining an aperture with PBS coupled with lower DNA damage induction. A methodology for calculating dose in water derived from measurements with a silicon-based detector was also presented. This work is the first to demonstrate experimentally the increase in LET caused by combining PBS with aperture, and to assess the potential DNA damage which is the initial step in the cascade of events leading to the majority of radiation-induced biological effects.