A systematic analysis of the particle irradiation data ensemble in the key of the microdosimetric kinetic model: Should clonogenic data be used for clinical relative biological effectiveness?

PURPOSE: To perform a systematic analysis of the Particle Irradiation Data Ensemble (PIDE) database for clonogenic survival assays in the context of the Microdosimetric Kinetic Model (MKM). METHODS AND MATERIAL: Our study used data from the PIDE database containing data on various cell lines and radiation types. Two main parameters of the MKM were determined experiment-wise: the domain radius, which accounts for the increase of the linear parameter as a function of LET or lineal energy, and the nucleus radius, which accounts for the overkilling effect at LET high enough. We used experiments with LET less and more than 75 keV/µm to determine domain and nucleus radius, respectively. Experiments with cells in asynchronous phase of the cell cycle and monoenergetic beams were considered, and data from 294 out of 461 available experiments with protons, alpha, and carbon beams were used. RESULTS: Domain and nucleus radii were determined for 32 cell lines as the median among cell-specific experiments after filtering experiments using protons, α-particles, and carbon ions, including 28 human cells and 12 rodent cells. The median values found for domain radii were 380 nm for normal human cells, 390 nm for tumor human cells, 295 nm for normal rodent cells, and 525 nm for tumor rodent cells (only one experiment with rodent tumor cells) with large variability across cell lines and across experiments on each cell line. CONCLUSIONS: Large inter-experiment variabilities were found for the same cell lines, based on enormous experimental uncertainties and different experimental conditions. Our analysis raises questions about how convenient is to use clonogenic data to feed RBE models to be utilized in the clinical practice in particle therapy.

Preparation of a radiobiology beam line at the 18 MeV proton cyclotron facility at CNA.

Proton therapy has gained interest in recent years due to its excellent clinical outcomes. However, the lack of accurate biological data, especially in the Bragg peak region of clinical beams, makes it difficult to implement biophysically optimized treatment plans in clinical practice. In this context, low energy proton accelerator facilities provide the perfect environment to collect good radiobiological data, as they can produce high LET beams with narrow energy distributions. This study presents the radiobiology beam line that has been designed at the 18 MeV proton cyclotron facility at the National Centre of Accelerators (CNA, Seville, Spain), to perform irradiations of mono-layer cell cultures. To ensure that all the cells receive the same dose with a suitable dose rate, low beam intensities and broad and homogeneous beam profiles are necessary. To do so, at the CNA an unfocused beam has been used, broadened with a 500 µm thick aluminium scattering foil. Homogeneous dose profiles, with deviations lower than 10% have been obtained over a circular surface of 35 mm diameter for an incident average energy of 12.8 MeV. Further, a Monte Carlo simulation of the beam line has been developed with Geant4, and benchmarked towards experimental measurements, with differences generally below 1%. Once validated, the code has been used, together with an ionization chamber, for dosimetry studies, to characterize the beam and monitor the dose. Finally, cultures of Human Bone Osteosarcoma cells (U2OS) have been successfully irradiated at the radiobiology beam line, investigating the effects of radiation in terms of DNA damage induction.

Calculation of clinical dose distributions in proton therapy from microdosimetry.

PURPOSE: To introduce a new algorithm-MicroCalc-for dose calculation by modeling microdosimetric energy depositions and the spectral fluence at each point of a particle beam. Proton beams are considered as a particular case of the general methodology. By comparing the results obtained against Monte Carlo computations, we aim to validate the microdosimetric formalism presented here and in previous works. MATERIALS AND METHODS: In previous works, we developed a function on the energy for the average energy imparted to a microdosimetric site per event and a model to compute the energetic spectrum at each point of the patient image. The number of events in a voxel is estimated assuming a model in which the voxel is completely filled by microdosimetric sites. Then, dose at every voxel is computed by integrating the average energy imparted per event multiplied by the number of events per energy beam of the spectral distribution within the voxel. Our method is compared with the proton convolution superposition (PCS) algorithm implemented in Eclipse™ and the fast Monte Carlo code MCsquare, which is here considered the benchmark, for in-water calculations, using in both cases clinically validated beam data. Two clinical cases are considered: a brain and a prostate case. RESULTS: For a SOBP beam in water, the mean difference at the central axis found for MicroCalc is of 0.86% against 1.03% for PCS. Three-dimensional gamma analyses in the PTVs compared with MCsquare for criterion (3%, 3 mm) provide gamma index of 95.07% with MicroCalc vs 94.50% with PCS for the brain case and 99.90% vs 100.00%, respectively, for the prostate case. For selected organs at risk in each case (brainstem and rectum), mean and maximum difference with respect to MCsquare dose are analyzed. In the brainstem, mean differences are 0.25 Gy (MicroCalc) vs 0.56 Gy (PCS), whereas for the rectum, these values are 0.05 Gy (MicroCalc) vs 0.07 Gy (PCS). CONCLUSIONS: The accuracy of MicroCalc seems to be, at least, not inferior to PCS, showing similar or better agreement with MCsquare in the considered cases. Additionally, the algorithm enables simultaneous computation of other quantities of interest. These results seem to validate the microdosimetric methodology in which the algorithm is based on.

An implementation to read and write IAEA phase-space files in GEANT4-based simulations.

PURPOSE: To develop a stand-alone code to make any application coded with the GEANT4 (GEometry ANd Tracking, version 4) toolkit capable of reading and writing phase-space (phsp) files in the format created by the IAEA (International Atomic Energy Agency), so that the exchange of phsp files between other validated Monte Carlo (MC) codes and GEANT4 is possible. METHODS: We present a stand-alone code, written in C++ object-oriented language, developed in a way that ensures the compatibility with future versions of the IAEA phsp format. The aim of the reader part is to get the information from a given IAEA phsp file and create the primary particles in a GEANT4 user application. On the other hand, the writer part of the code is the responsible for writing the IAEA phsp files during a run of the GEANT4 application. RESULTS: A testing simulation was written with GEANT4 to verify the performance of this code, with satisfactory results. An example of use in a GEANT4 application which simulates the treatment head of a radiotherapy linear electron accelerator (linac) is also shown, comparing dose calculations with experimental data. CONCLUSIONS: This stand-alone package, which can be used in any GEANT4 application, allows the exchange of validated phsp files between different MC codes and the use of phsp data from many different accelerators and fields in dosimetry studies. Furthermore, it also offers additional utilities of interest in medical applications.