Radiobiological quantities in proton-therapy: Estimation and validation using Geant4-based Monte Carlo simulations

Monte Carlo simulations of cell survival in proton SOBP

Objective. The objective of this study is to develop a multi-scale modeling approach that accurately predicts radiation-induced DNA damage and survival fraction in specific cell lines.Approach. A Monte Carlo based simulation framework was employed to make the predictions. The FLUKA Monte Carlo code was utilized to estimate absorbed doses and fluence energy spectra, which were then used in the Monte Carlo Damage Simulation code to compute DNA damage yields in Chinese hamster V79 cell lines. The outputs were converted into cell survival fractions using a previously published theoretical model. To reduce the uncertainties of the predictions, new values for the parameters of the theoretical model were computed, expanding the database of experimental points considered in the previous estimation. Simulated results were validated against experimental data, confirming the applicability of the framework for proton beams up to 230 MeV. Additionally, the impact of secondary particles on cell survival was estimated.Main results. The simulated survival fraction versus depth in a glycerol phantom is reported for eighteen different configurations. Two proton spread out Bragg peaks at several doses were simulated and compared with experimental data. In all cases, the simulations follow the experimental trends, demonstrating the accuracy of the predictions up to 230 MeV.Significance. This study holds significant importance as it contributes to the advancement of models for predicting biological responses to radiation, ultimately contributing to more effective cancer treatment in proton therapy.

The Proton-Boron Reaction Increases the Radiobiological Effectiveness of Clinical Low- and High-Energy Proton Beams: Novel Experimental Evidence and Perspectives

Protontherapy is a rapidly expanding radiotherapy modality where accelerated proton beams are used to precisely deliver the dose to the tumor target but is generally considered ineffective against radioresistant tumors. Proton-Boron Capture Therapy (PBCT) is a novel approach aimed at enhancing proton biological effectiveness. PBCT exploits a nuclear fusion reaction between low-energy protons and ¹¹B atoms, i.e. p+¹¹B→ 3α (p-B), which is supposed to produce highly-DNA damaging α-particles exclusively across the tumor-conformed Spread-Out Bragg Peak (SOBP), without harming healthy tissues in the beam entrance channel. To confirm previous work on PBCT, here we report new in-vitro data obtained at the 62-MeV ocular melanoma-dedicated proton beamline of the INFN-Laboratori Nazionali del Sud (LNS), Catania, Italy. For the first time, we also tested PBCT at the 250-MeV proton beamline used for deep-seated cancers at the Centro Nazionale di Adroterapia Oncologica (CNAO), Pavia, Italy. We used Sodium Mercaptododecaborate (BSH) as ¹¹B carrier, DU145 prostate cancer cells to assess cell killing and non-cancer epithelial breast MCF-10A cells for quantifying chromosome aberrations (CAs) by FISH painting and DNA repair pathway protein expression by western blotting. Cells were exposed at various depths along the two clinical SOBPs. Compared to exposure in the absence of boron, proton irradiation in the presence of BSH significantly reduced DU145 clonogenic survival and increased both frequency and complexity of CAs in MCF-10A cells at the mid- and distal SOBP positions, but not at the beam entrance. BSH-mediated enhancement of DNA damage response was also found at mid-SOBP. These results corroborate PBCT as a strategy to render protontherapy amenable towards radiotherapy-resilient tumor. If coupled with emerging proton FLASH radiotherapy modalities, PBCT could thus widen the protontherapy therapeutic index.

Hypoxia Transcriptomic Modifications Induced by Proton Irradiation in U87 Glioblastoma Multiforme Cell Line

In Glioblastoma Multiforme (GBM), hypoxia is associated with radioresistance and poor prognosis. Since standard GBM treatments are not always effective, new strategies are needed to overcome resistance to therapeutic treatments, including radiotherapy (RT). Our study aims to shed light on the biomarker network involved in a hypoxic (0.2% oxygen) GBM cell line that is radioresistant after proton therapy (PT). For cultivating cells in acute hypoxia, GSI’s hypoxic chambers were used. Cells were irradiated in the middle of a spread-out Bragg peak with increasing PT doses to verify the greater radioresistance in hypoxic conditions. Whole-genome cDNA microarray gene expression analyses were performed for samples treated with 2 and 10 Gy to highlight biological processes activated in GBM following PT in the hypoxic condition. We describe cell survival response and significant deregulated pathways responsible for the cell death/survival balance and gene signatures linked to the PT/hypoxia configurations assayed. Highlighting the molecular pathways involved in GBM resistance following hypoxia and ionizing radiation (IR), this work could suggest new molecular targets, allowing the development of targeted drugs to be suggested in association with PT.

Report on G4-Med, a Geant4 benchmarking system for medical physics applications developed by the Geant4 Medical Simulation Benchmarking Group

BACKGROUND

Geant4 is a Monte Carlo code extensively used in medical physics for a wide range of applications, such as dosimetry, micro- and nanodosimetry, imaging, radiation protection, and nuclear medicine. Geant4 is continuously evolving, so it is crucial to have a system that benchmarks this Monte Carlo code for medical physics against reference data and to perform regression testing.

AIMS

To respond to these needs, we developed G4-Med, a benchmarking and regression testing system of Geant4 for medical physics.

MATERIALS AND METHODS

G4-Med currently includes 18 tests. They range from the benchmarking of fundamental physics quantities to the testing of Monte Carlo simulation setups typical of medical physics applications. Both electromagnetic and hadronic physics processes and models within the prebuilt Geant4 physics lists are tested. The tests included in G4-Med are executed on the CERN computing infrastructure via the use of the geant-val web application, developed at CERN for Geant4 testing. The physical observables can be compared to reference data for benchmarking and to results of previous Geant4 versions for regression testing purposes.

RESULTS

This paper describes the tests included in G4-Med and shows the results derived from the benchmarking of Geant4 10.5 against reference data.

DISCUSSION

Our results indicate that the Geant4 electromagnetic physics constructor G4EmStandardPhysics_option4 gives a good agreement with the reference data for all the tests. The QGSP_BIC_HP physics list provided an overall adequate description of the physics involved in hadron therapy, including proton and carbon ion therapy. New tests should be included in the next stage of the project to extend the benchmarking to other physical quantities and application scenarios of interest for medical physics.

CONCLUSION

The results presented and discussed in this paper will aid users in tailoring physics lists to their particular application.

RBE variation in prostate carcinoma cells in active scanning proton beams: In-vitro measurements in comparison with phenomenological models

PURPOSE

In-vitro radiobiological studies are essential for modelling the relative biological effectiveness (RBE) in proton therapy. The purpose of this study was to experimentally determine the RBE values in proton beams along the beam path for human prostate carcinoma cells (Du-145). RBE-dose and RBE-LET_d (dose-averaged linear energy transfer) dependencies were investigated and three phenomenological RBE models, i.e. McNamara, Rørvik and Wilkens were benchmarked for this cell line.

METHODS

Cells were placed at multiple positions along the beam path, employing an in-house developed solid phantom. The experimental setup reflected the clinical prostate treatment scenario in terms of field size, depth, and required proton energies (127.2-180.1 MeV) and the physical doses from 0.5 to 6 Gy were delivered. The reference irradiation was performed with 200 kV X-ray beams. Respective (α/β) values were determined using the linear quadratic model and LET_d was derived from the treatment planning system at the exact location of cells.

RESULTS AND CONCLUSION

Independent of the cell survival level, all experimental RBE values were consistently higher in the target than the generic clinical RBE value of 1.1; with the lowest RBE value of 1.28 obtained at the beginning of the SOBP. A systematic RBE decrease with increasing dose was observed for the investigated dose range. The RBE values from all three applied models were considerably smaller than the experimental values. A clear increase of experimental RBE values with LET_d parameter suggests that proton LET must be taken into consideration for this low (α/β) tissue.