In Vivo Validation of the BIANCA Biophysical Model: Benchmarking against Rat Spinal Cord RBE Data

A method to predict space radiation biological effectiveness for non-cancer effects following intense Solar Particle Events

In addition to the continuous exposure to cosmic rays, astronauts in space are occasionally exposed to Solar Particle Events (SPE), which involve less energetic particles but can deliver much higher doses. The latter can exceed several Gy in a few hours for the most intense SPEs, for which non-stochastic effects are thus a major concern. To identify adequate shielding conditions that would allow respecting the dose limits established by the various space agencies, the absorbed dose in the considered organ/tissue must be multiplied by the corresponding Relative Biological Effectiveness (RBE), which is a complex quantity depending on several factors including particle type and energy, considered biological effect, level of effect (and thus absorbed dose), etc. While in several studies only the particle-type dependence of RBE is taken into account, in this work we developed and applied a new approach where, thanks to an interface between the FLUKA Monte Carlo transport code and the BIANCA biophysical model, the RBE dependence on particle energy and absorbed dose was also considered. Furthermore, we included in the considered SPE spectra primary particles heavier than protons, which in many studies are neglected. This approach was then applied to the October 2003 SPE (the most intense SPE of solar cycle 23, also known as "Halloween event") and the January 2005 event, which was characterized by a lower fluence but a harder spectrum, i.e., with higher-energy particles. The calculation outcomes were then discussed and compared with the current dose limits established for skin and blood forming organs in case of 30-days missions. This work showed that the BIANCA model, if interfaced to a radiation transport code, can be used to calculate the RBE values associated to Solar Particle Events. More generally, this work emphasizes the importance of taking into account the RBE dependence on particle energy and dose when calculating equivalent doses.

Computational approaches in the estimation of radiobiological damage for human-malignant cells irradiated with clinical proton and carbon beams

PURPOSE

The use of Monte Carlo (MC) simulations capable of reproducing radiobiological effects of ionising radiation on human cell lines is of great importance, especially for cases involving protons and heavier ion beams. In the latter, huge uncertainties can arise mainly related to the effects of the secondary particles produced in the beam-tissue interaction. This paper reports on a detailed MC study performed using Geant4-based approach on three cancer cell lines, the HTB-177, CRL-5876 and MCF-7, that were previously irradiated with therapeutic proton and carbon ion beams.

METHODS

A Geant4-based approach used jointly with analytical calculations has been developed to provide a more realistic estimation of the radiobiological damage produced by proton and carbon beams in tissues, reproducing available data obtained from in vitro cell irradiations. The MC "Hadrontherapy" Geant4 application and the Local Effect Model: LEM I, LEM II and LEM III coupled with the different numerical approaches: RapidRusso (RR) and RapidScholz (RS) were used in the study.

RESULTS

Experimental survival curves are compared with those evaluated using the highlighted Geant4 MC-based approach via chi-square statistical analysis, for the combinations of radiobiological models and numerical approaches, as outlined above.

CONCLUSION

This study has presented a comparison of the survival data from MC simulations to experimental survival data for three cancer cell lines. An overall best level of agreement was obtained for the HTB-177 cells.

Radiation Damage in Biomolecules and Cells 2.0

It is well known that ionizing radiation, when it hits living cells, causes a plethora of different damage types at different levels [...].

A Mission to Mars: Prediction of GCR Doses and Comparison with Astronaut Dose Limits

Long-term human space missions such as a future journey to Mars could be characterized by several hazards, among which radiation is one the highest-priority problems for astronaut health. In this work, exploiting a pre-existing interface between the BIANCA biophysical model and the FLUKA Monte Carlo transport code, a study was performed to calculate astronaut absorbed doses and equivalent doses following GCR exposure under different shielding conditions. More specifically, the interface with BIANCA allowed us to calculate both the RBE for cell survival, which is related to non-cancer effects, and that for chromosome aberrations, related to the induction of stochastic effects, including cancer. The results were then compared with cancer and non-cancer astronaut dose limits. Concerning the stochastic effects, the equivalent doses calculated by multiplying the absorbed dose by the RBE for chromosome aberrations ("high-dose method") were similar to those calculated using the Q-values recommended by ICRP. For a 650-day mission at solar minimum (representative of a possible Mars mission scenario), the obtained values are always lower than the career limit recommended by ICRP (1 Sv), but higher than the limit of 600 mSv recently adopted by NASA. The comparison with the JAXA limits is more complex, since they are age and sex dependent. Concerning the deterministic limits, even for a 650-day mission at solar minimum, the values obtained by multiplying the absorbed dose by the RBE for cell survival are largely below the limits established by the various space agencies. Following this work, BIANCA, interfaced with an MC transport code such as FLUKA, can now predict RBE values for cell death and chromosome aberrations following GCR exposure. More generally, both at solar minimum and at solar maximum, shielding of 10 g/cm² Al seems to be a better choice than 20 g/cm² for astronaut protection against GCR.

First application of the BIANCA biophysical model to carbon-ion patient cases

Objective. The main objective of this work consists of applying, for the first time, the BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) biophysical model to the RBE calculation for C-ion cancer patients, and comparing the outcomes with those obtained by the LEM I model, which is applied in clinics. Indeed, the continuous development of heavy-ion cancer therapy requires modelling of biological effects of ion beams on tumours and normal tissues. The relative biological effectiveness (RBE) of heavy ions is higher than that of protons, with a significant variation along the beam path. Therefore, it requires a precise modelling, especially for the pencil-beam scanning technique. Currently, two radiobiological models, LEM I and MKM, are in use for heavy ions in scanned pencil-beam facilities. Approach. Utilizing an interface with the FLUKA Particle Therapy Tool, BIANCA was applied to re-calculate the RBE-weighted dose distribution for carbon-ion treatment of three patients (chordoma, head-and-neck and prostate) previously irradiated at CNAO, where radiobiological optimization was based on LEM I. The predictions obtained by BIANCA were based either on chordoma cell survival (RBE surv ), or on dicentric aberrations in peripheral blood lymphocytes (RBE ab ), which are indicators of late normal tissue damage, including secondary tumours. The simulation outcomes were then compared with those provided by LEM I. Main results. While in the target and in the entrance channel BIANCA predictions were lower than those obtained by LEM I, the two models provided very similar results in the considered OAR. The observed differences between RBE surv and RBE ab (which were also dependent on fractional dose and LET) suggest that in normal tissues the information on cell survival should be integrated by information more closely related to the induction of late damage, such as chromosome aberrations. Significance. This work showed that BIANCA is suitable for treatment plan optimization in ion-beam therapy, especially considering that it can predict both cell survival and chromosome aberrations and has previously shown good agreement with carbon-ion experimental data.

Impact of DNA Repair Kinetics and Dose Rate on RBE Predictions in the UNIVERSE

Accurate knowledge of the relative biological effectiveness (RBE) and its dependencies is crucial to support modern ion beam therapy and its further development. However, the influence of different dose rates of the reference radiation and ion beam are rarely considered. The ion beam RBE-model within our "UNIfied and VERSatile bio response Engine" (UNIVERSE) is extended by including DNA damage repair kinetics to investigate the impact of dose-rate effects on the predicted RBE. It was found that dose-rate effects increase with dose and biological effects saturate at high dose-rates, which is consistent with data- and model-based studies in the literature. In a comparison with RBE measurements from a high dose in-vivo study, the predictions of the presented modification were found to be improved in comparison to the previous version of UNIVERSE and existing clinical approaches that disregard dose-rate effects. Consequently, DNA repair kinetics and the different dose rates applied by the reference and ion beams might need to be considered in biophysical models to accurately predict the RBE. Additionally, this study marks an important step in the further development of UNIVERSE, extending its capabilities in giving theoretical guidance to support progress in ion beam therapy.

Radiobiological damage by space radiation: extension of the BIANCA model to heavy ions up to iron, and pilot application to cosmic ray exposure

Space research seems to be object of a renewed interest, also considering that human missions to the Moon, and possibly Mars, are being planned. Among the risks affecting such missions, astronauts' exposure to space radiation is a major concern. In this work, the question of the evaluation of biological damage by Galactic Cosmic Rays (GCR) was addressed by a biophysical model called BIophysical ANalysis of Cell death and chromosome Aberrations (BIANCA), which simulates the induction of cell death and chromosome aberrations by different ions. While previously BIANCA has been validated for calculating cell death along hadrontherapy beams up to oxygen, herein the approach was extended up to Fe ions. Specifically, experimental survival curves available in literature for V79 cells irradiated by Si-, Ne-, Ar- and Fe-ions were reproduced, and a reference radiobiological database describing V79 cell survival as a function of ion type (1 ⩽Z⩽ 26), energy and dose was constructed. Analogous databases were generated for Chinese hamster ovary hamster cells and human skin fibroblasts, finding good agreement between simulations and data. Concerning chromosome aberrations, which are regarded as radiation risk biomarkers, dicentric data in human lymphocytes irradiated by heavy ions up to iron were reproduced, and a radiobiological database allowing calculation of lymphocyte dicentric yields as a function of dose, ion type (1 ⩽Z⩽ 26) and energy was constructed. Following interface between BIANCA and the FLUKA Monte Carlo transport code, a feasibility study was performed to calculate the relative biological effectiveness (RBE) of different GCR spectrum components, for both dicentrics and cell death. Fe-ions, although representing only 10% of the total absorbed dose, were found to be responsible for about 35%-40% of the RBE-weighted dose. Interestingly, the RBE for dicentrics was higher than that for cell survival. More generally, this work shows that BIANCA can calculate RBE values for cell death and lymphocyte dicentrics not only for ion therapy, but also for space radiation.

Extension of the BIANCA biophysical model up to Fe-ions and applications for space radiation research

BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) is a biophysical model, implemented as a Monte Carlo code, which simulates the induction of chromosome aberrations and cell death by different monochromatic ion beams (i.e., different ion types and energy values), as well as photons. In previous works, the model predictions for cell survival and lymphocyte dicentrics along therapeutic-like ion beams have been successfully benchmarked against experimental data. With the aim of evaluating the biological damage induced by Galactic Cosmic Rays (GCR), in this study BIANCA was extended up to Fe-ions. A radiobiological database describing human skin fibroblast cell survival and lymphocyte dicentrics as a function of ion type (1≤Z≤26) and energy, as well as dose, was constructed. Afterwards, interfacing BIANCA with the FLUKA Monte Carlo transport code, a feasibility study was performed to calculate the Relative Biological Effectiveness (RBE) of Galactic Cosmic Rays, both for dicentrics and for cell survival; the results were discussed with respect to available space radiation protection limits. Following this work, BIANCA can now provide RBE predictions of cell killing, which can be related to deterministic effects, and lymphocyte dicentrics, more related to stochastic effects, for space radiation exposure.

Healthy Tissue Damage Following Cancer Ion Therapy: A Radiobiological Database Predicting Lymphocyte Chromosome Aberrations Based on the BIANCA Biophysical Model

Chromosome aberrations are widely considered among the best biomarkers of radiation health risk due to their relationship with late cancer incidence. In particular, aberrations in peripheral blood lymphocytes (PBL) can be regarded as indicators of hematologic toxicity, which is a major limiting factor of radiotherapy total dose. In this framework, a radiobiological database describing the induction of PBL dicentrics as a function of ion type and energy was developed by means of the BIANCA (BIophysical ANalysis of Cell death and chromosome Aberrations) biophysical model, which has been previously applied to predict the effectiveness of therapeutic-like ion beams at killing tumour cells. This database was then read by the FLUKA Monte Carlo transport code, thus allowing us to calculate the Relative Biological Effectiveness (RBE) for dicentric induction along therapeutic C-ion beams. A comparison with previous results showed that, while in the higher-dose regions (e.g., the Spread-Out Bragg Peak, SOBP), the RBE for dicentrics was lower than that for cell survival. In the lower-dose regions (e.g., the fragmentation tail), the opposite trend was observed. This work suggests that, at least for some irradiation scenarios, calculating the biological effectiveness of a hadrontherapy beam solely based on the RBE for cell survival may lead to an underestimation of the risk of (late) damage to healthy tissues. More generally, following this work, BIANCA has gained the capability of providing RBE predictions not only for cell killing, but also for healthy tissue damage.

Biological effectiveness of He-3 and He-4 ion beams for cancer hadrontherapy: a study based on the BIANCA biophysical model

While cancer therapy with protons and C-ions is continuously spreading, in the near future patients will be also treated with He-ions which, in comparison to photons, combine the higher precision of protons with the higher relative biological effectiveness (RBE) of C-ions. Similarly to C-ions, also for He-ions the RBE variation along the beam must be known as precisely as possible, especially for active beam delivery systems. In this framework the BIANCA biophysical model, which has already been applied to calculate the RBE along proton and C-ion beams, was extended to⁴He-ions and, following interface with the FLUKA code, was benchmarked against cell survival data on CHO normal cells and Renca tumour cells irradiated at different positions along therapeutic-like⁴He-ion beams at the Heidelberg Ion-beam Therapy centre, where the first He-ion patient will be treated soon. Very good agreement between simulations and data was obtained, showing that BIANCA can now be used to predict RBE following irradiation with all ion types that are currently used, or will be used soon, for hadrontherapy. Thanks to the development of a reference simulation database describing V79 cell survival for ion and photon irradiation, these predictions can be cell-type specific because analogous databases can be produced, in principle, for any cell line. Furthermore, survival data on CHO cells irradiated by a He-3 beam were reproduced to compare the biophysical properties of He-4 and He-3 beams, which is currently an open question. This comparison showed that, at the same depth, He-4 beams tend to have a higher RBE with respect to He-3 beams, and that this difference is also modulated by the considered physical dose, as well as the cell radiosensitivity. However, at least for the considered cases, no significant difference was found for the ratio between the RBE-weighted dose in the SOBP and that in the entrance plateau.

The RBE in ion beam radiotherapy: In vivo studies and clinical application

Ion beams used for radiotherapy exhibit an increased relative biological effectiveness (RBE), which depends on several physical treatment parameters as well as on biological factors of the irradiated tissues. While the RBE is an experimentally well-defined quantity, translation to patients is complex and requires radiobiological studies, dedicated models to calculate the RBE in treatment planning as well as strategies for dose prescription. Preclinical in vivo studies and analysis of clinical outcome are important to validate and refine RBE-models. This review describes the concept of the experimental and clinical RBE and explains the fundamental dependencies of the RBE based on in vitro experiments. The available preclinical in vivo studies on normal tissue and tumor RBE for ions heavier than protons are reviewed in the context of the historical and present development of ion beam radiotherapy. In addition, the role of in vivo RBE-values in the development and benchmarking of RBE-models as well as the transition of these models to clinical application are described. Finally, limitations in the translation of experimental RBE-values into clinical application and the direction of future research are discussed.