The sensitivity of radiobiological models in carbon ion radiotherapy (CIRT) and its consequences on the clinical treatment plan: Differences between LEM and MKM models.

PURPOSE: Carbon ion radiotherapy (CIRT) relies on relative biological effectiveness (RBE)-weighted dose calculations. Japanese clinics predominantly use the microdosimetric kinetic model (MKM), while European centers utilize the local effect model (LEM). Despite both models estimating RBE-distributions in tissue, their physical and mathematical assumptions differ, leading to significant disparities in RBE-weighted doses. Several European clinics adopted Japanese treatment schedules, necessitating adjustments in dose prescriptions and organ at risk (OAR) constraints. In the context of these two clinically used standards for RBE-weighted dose estimation, the objective of this study was to highlight specific scenarios for which the translations between models diverge, as shortcomings between them can influence clinical decisions. METHODS: Our aim was to discuss planning strategies minimizing those discrepancies, ultimately striving for more accurate and robust treatments. Evaluations were conducted in a virtual water phantom and patient CT-geometry, optimizing LEM RBE-weighted dose first and recomputing MKM thereafter. Dose-averaged linear energy transfer (LETd) distributions were also assessed. RESULTS: Results demonstrate how various parameters influence LEM/MKM translation. Similar LEM-dose distributions lead to markedly different MKM-dose distributions and variations in LETd. Generally, a homogeneous LEM RBE-weighted dose aligns with lower MKM values in most of the target volume. Nevertheless, paradoxical MKM hotspots may emerge (at the end of the range), potentially influencing clinical outcomes. Therefore, translation between models requires great caution. CONCLUSIONS: Understanding the relationship between these two clinical standards enables combining European and Japanese based experiences. The implementation of optimal planning strategies ensures the safety and acceptability of the clinical plan for both models and therefore enhances plan robustness from the RBE-weighted dose and LETd distribution point of view. This study emphasizes the importance of optimal planning strategies and the need for comprehensive CIRT plan quality assessment tools. In situations where simultaneous LEM and MKM computation capabilities are lacking, it can provide guidance in plan design, ultimately contributing to enhanced CIRT outcomes.

DNA-damage RBE assessment for combined boron and gadolinium neutron capture therapy.

PURPOSE: Neutron capture therapy (NCT) by 10B and 157Gd agents is a unique irradiation-based method which can be used to treat brain tumors. Current study aims to quantitatively evaluate the relative biological effectiveness (RBE) and dose distributions during the combined BNCT and GdNCT modalities through a hybrid Monte Carlo (MC) simulation approach. METHODS: Snyder head phantom as well as a cubic hypothetical tumor was at first modeled by Geant4 MC Code. Then, the energy spectra and dose distribution relevant to the released secondary particles during the combined Gd/BNCT were scored for different concentrations of 157Gd and 10B inside tumor volume. Finally, the scored energy spectra were imported to the MCDS code to estimate both RBESSB and RBEDSB values for different 157Gd concentrations. RESULTS: The results showed that combined Gd/BNCT increases the fluence-averaged RBESSB values by about 1.7 times when 157Gd concentration increments from 0 to 2000 µg/g for both considered cell oxygen levels (pO2 = 10% and 100%). Besides, a reduction of about 26% was found for fluence-averaged RBEDSB values with an increment of 157Gd concentration in tumor volume. CONCLUSION: From the results, it can be concluded that combined Gd/BNCT technique can improve tumor coverage with higher dose levels but in the expense of RBEDSB reduction which can affect the clinical efficacy of the NCT technique.

High-LET charged particles: radiobiology and application for new approaches in radiotherapy.

The number of patients treated with charged-particle radiotherapy as well as the number of treatment centers is increasing worldwide, particularly regarding protons. However, high-linear energy transfer (LET) particles, mainly carbon ions, are of special interest for application in radiotherapy, as their special physical features result in high precision and hence lower toxicity, and at the same time in increased efficiency in cell inactivation in the target region, i.e., the tumor. The radiobiology of high-LET particles differs with respect to DNA damage repair, cytogenetic damage, and cell death type, and their increased LET can tackle cells' resistance to hypoxia. Recent developments and perspectives, e.g., the return of high-LET particle therapy to the US with a center planned at Mayo clinics, the application of carbon ion radiotherapy using cost-reducing cyclotrons and the application of helium is foreseen to increase the interest in this type of radiotherapy. However, further preclinical research is needed to better understand the differential radiobiological mechanisms as opposed to photon radiotherapy, which will help to guide future clinical studies for optimal exploitation of high-LET particle therapy, in particular related to new concepts and innovative approaches. Herein, we summarize the basics and recent progress in high-LET particle radiobiology with a focus on carbon ions and discuss the implications of current knowledge for charged-particle radiotherapy. We emphasize the potential of high-LET particles with respect to immunogenicity and especially their combination with immunotherapy.

Modeling of ionizing radiation-induced chromosome aberration and tumor prevalence based on two classes of DNA double-strand breaks clustering in chromatin domains.

There has been some controversy over the use of radiobiological models when modeling the dose-response curves of ionizing radiation (IR)-induced chromosome aberration and tumor prevalence, as those curves usually show obvious non-targeted effects (NTEs) at low doses of high linear energy transfer (LET) radiation. The lack of understanding the contribution of NTEs to IR-induced carcinogenesis can lead to distinct deviations of relative biological effectiveness (RBE) estimations of carcinogenic potential, which are widely used in radiation risk assessment and radiation protection. In this work, based on the initial pattern of two classes of IR-induced DNA double-strand breaks (DSBs) clustering in chromatin domains and the subsequent incorrect repair processes, we proposed a novel radiobiological model to describe the dose-response curves of two carcinogenic-related endpoints within the same theoretical framework. The representative experimental data was used to verify the consistency and validity of the present model. The fitting results indicated that, compared with targeted effect (TE) and NTE models, the current model has better fitting ability when dealing with the experimental data of chromosome aberration and tumor prevalence induced by multiple types of IR with different LETs. Notably, the present model without introducing an NTE term was adequate to describe the dose-response curves of IR-induced chromosome aberration and tumor prevalence with NTEs in low-dose regions. Based on the fitting parameters, the LET-dependent RBE values were calculated for three given low doses. Our results showed that the RBE values predicted by the current model gradually decrease with the increase of doses for the endpoints of chromosome aberration and tumor prevalence. In addition, the calculated RBE was also compared with those evaluated from other models. These analyses show that the proposed model can be used as an alternative tool to well describe dose-response curves of multiple carcinogenic-related endpoints and effectively estimate RBE in low-dose regions.

Calculating dose-averaged linear energy transfer in an analytical treatment planning system for carbon-ion radiotherapy.

BACKGROUND: Compelling evidence shows the association between the relative biological effectiveness (RBE) of carbon-ion radiotherapy (CIRT) and the dose averaged linear energy transfer (LETd). However, the ability to calculate the LETd in commercially available treatment planning systems (TPS) is lacking. PURPOSE: This study aims to develop a method of calculating the LETd of CIRT plans that could be robustly carried out in RayStation (V10B, Raysearch, Sweden). METHODS: The calculation used the fragment spectra in RayStation for the CIRT treatment planning. The dose-weighted averaging procedure was supported by the microdosimetric kinetic model (MKM). The MKM-based pencil beam dose engine (PBA, v4.2) for calculating RBE-weighted doses was reformulated to become a LET-weighted calculating engine. A separate module was then configured to inversely calculate the LETd from the absorbed dose of a plan and the associated fragment spectra. In this study, the ion and energy-specific LET table in the LETd module was further matched with the values decoded from the baseline data of the Syngo TPS (V13C, Siemens, Germany). The LETd distributions of several monoenergetic and modulated beams were calculated and validated against the values derived from the Syngo TPS and the published data. RESULTS: The differences in LETds of the monoenergetic beams between the new method and the traditional method were within 3% in the entrance and Bragg-peak regions. However, a larger difference was observed in the distal region. The results of the modulated beams were in good agreement with the works from the published literature. CONCLUSIONS: The method presented herein reformulates the MKM dose engine in the RayStation TPS to inversely calculate LETds. The robustness and accuracy were demonstrated.

A methodology to abridge microdosimetric distributions without a significant loss of the spectral information needed for the RBE computation in carbon ion therapy.

BACKGROUND: In order to compute the relative biological effectiveness (RBE) of ion radiation therapy with the Mayo Clinic Florida microdosimetric kinetic model (MCF MKM), it is necessary to process entire microdosimetric distributions. Therefore, a posteriori RBE recalculations (i.e., for a different cell line or another biological endpoint) would require whole spectral information. It is currently not practical to compute and store all this data for each clinical voxel. PURPOSE: To develop a methodology that allows to store a limited amount of physical information without losing accuracy in the RBE calculations nor the possibility of a posteriori RBE recalculations. METHODS: Computer simulations for four monoenergetic 12 C ion beams and a 12 C ion spread-out Bragg peak (SOBP) were performed to assess lineal energy distributions as a function of the depth within a water phantom. These distributions were used in combination with the MCF MKM to compute the in vitro clonogenic survival RBE for human salivary gland tumor cells (HSG cell line) and human skin fibroblasts (NB1RGB cell line). The RBE values were also calculated with a new abridged microdosimetric distribution methodology (AMDM) and compared with the reference RBE calculations using the entire distributions. RESULTS: The maximum relative deviation between the RBE values computed using the entire distributions and the AMDM was 0.61% (monoenergetic beams) and 0.49% (SOBP) for the HSG cell line, while 0.45% (monoenergetic beams) and 0.26% (SOBP) for the NB1RGB cell line. CONCLUSION: The excellent agreement between the RBE values computed using the entire lineal energy distributions and the AMDM represents a milestone for the clinical implementation of the MCF MKM.

Combined RBE and OER optimization in proton therapy with FLUKA based on EF5-PET.

INTRODUCTION: Tumor hypoxia is associated with poor treatment outcome. Hypoxic regions are more radioresistant than well-oxygenated regions, as quantified by the oxygen enhancement ratio (OER). In optimization of proton therapy, including OER in addition to the relative biological effectiveness (RBE) could therefore be used to adapt to patient-specific radioresistance governed by intrinsic radiosensitivity and hypoxia. METHODS: A combined RBE and OER weighted dose (ROWD) calculation method was implemented in a FLUKA Monte Carlo (MC) based treatment planning tool. The method is based on the linear quadratic model, with α and ß parameters as a function of the OER, and therefore a function of the linear energy transfer (LET) and partial oxygen pressure (pO2 ). Proton therapy plans for two head and neck cancer (HNC) patients were optimized with pO2 estimated from [18 F]-EF5 positron emission tomography (PET) images. For the ROWD calculations, an RBE of 1.1 (RBE1.1,OER ) and two variable RBE models, Rørvik (ROR) and McNamara (MCN), were used, alongside a reference plan without incorporation of OER (RBE1.1 ). RESULTS: For the HNC patients, treatment plans in line with the prescription dose and with acceptable target ROWD could be generated with the established tool. The physical dose was the main factor modulated in the ROWD. The impact of incorporating OER during optimization of HNC patients was demonstrated by the substantial difference found between ROWD and physical dose in the hypoxic tumor region. The largest physical dose differences between the ROWD optimized plans and the reference plan was 12.2 Gy. CONCLUSION: The FLUKA MC based tool was able to optimize proton treatment plans taking the tumor pO2 distribution from hypoxia PET images into account. Independent of RBE-model, both elevated LET and physical dose were found in the hypoxic regions, which shows the potential to increase the tumor control compared to a conventional optimization approach.

Treatment-planning approaches to intensity modulated proton therapy and the impact on dose-weighted linear energy transfer.

PURPOSE: We quantified the effect of various forward-based treatment-planning strategies in proton therapy on dose-weighted linear energy transfer (LETd). By maintaining the dosimetric quality at a clinically acceptable level, we aimed to evaluate the differences in LETd among various treatment-planning approaches and their practicality in minimizing biologic uncertainties associated with LETd. METHOD: Eight treatment-planning strategies that are achievable in commercial treatment-planning systems were applied on a cylindrical water phantom and four pediatric brain tumor cases. Each planning strategy was compared to either an opposed lateral plan (phantom study) or original clinical plan (patient study). Deviations in mean and maximum LETd from clinically acceptable dose distributions were compared. RESULTS: In the phantom study, using a range shifter and altering the robust scenarios during optimization had the largest effect on the mean clinical target volume LETd, which was reduced from 4.5 to 3.9 keV/µm in both cases. Variations in the intersection angle between beams had the largest effect on LETd in a ring defined 3 to 5 mm outside the target. When beam intersection angles were reduced from opposed laterals (180°) to 120°, 90°, and 60°, corresponding maximum LETd increased from 7.9 to 8.9, 10.9, and 12.2 keV/µm, respectively. A clear trend in mean and maximum LETd variations in the clinical cases could not be established, though spatial distribution of LETd suggested a strong dependence on patient anatomy and treatment geometry. CONCLUSION: Changes in LETd from treatment-plan setup follow intuitive trends in a controlled phantom experiment. Anatomical and other patient-specific considerations, however, can preclude generalizable strategies in clinical cases. For pediatric cranial radiation therapy, we recommend using opposed lateral treatment fields to treat midline targets.

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/cm2 Al seems to be a better choice than 20 g/cm2 for astronaut protection against GCR.

In vitro dose effect relationships of actinium-225- and lutetium-177-labeled PSMA-I&T.

PURPOSE: Targeting the prostate-specific membrane antigen (PSMA) using lutetium-177-labeled PSMA-specific tracers has become a very promising novel therapy option for prostate cancer (PCa). The efficacy of this therapy might be further improved by replacing the ß-emitting lutetium-177 with the α-emitting actinium-225. Actinium-225 is thought to have a higher therapeutic efficacy due to the high linear energy transfer (LET) of the emitted α-particles, which can increase the amount and complexity of the therapy induced DNA double strand breaks (DSBs). Here we evaluated the relative biological effectiveness of [225Ac]Ac-PSMA-I&T and [177Lu]Lu-PSMA-I&T by assessing in vitro binding characteristics, dosimetry, and therapeutic efficacy. METHODS AND RESULTS: The PSMA-expressing PCa cell line PC3-PIP was used for all in vitro assays. First, binding and displacement assays were performed, which revealed similar binding characteristics between [225Ac]Ac-PSMA-I&T and [177Lu]Lu-PSMA-I&T. Next, the assessment of the number of 53BP1 foci, a marker for the number of DNA double strand breaks (DSBs), showed that cells treated with [225Ac]Ac-PSMA-I&T had slower DSB repair kinetics compared to cells treated with [177Lu]Lu-PSMA-I&T. Additionally, clonogenic survival assays showed that specific targeting with [225Ac]Ac-PSMA-I&T and [177Lu]Lu-PSMA-I&T caused a dose-dependent decrease in survival. Lastly, after dosimetric assessment, the relative biological effectiveness (RBE) of [225Ac]Ac-PSMA-I&T was found to be 4.2 times higher compared to [177Lu]Lu-PSMA-I&T. CONCLUSION: We found that labeling of PSMA-I&T with lutetium-177 or actinium-225 resulted in similar in vitro binding characteristics, indicating that the distinct biological effects observed in this study are not caused by a difference in uptake of the two tracers. The slower repair kinetics of [225Ac]Ac-PSMA-I&T compared to [177Lu]Lu-PSMA-I&T correlates to the assumption that irradiation with actinium-225 causes more complex, more difficult to repair DSBs compared to lutetium-177 irradiation. Furthermore, the higher RBE of [225Ac]Ac-PSMA-I&T compared to [177Lu]Lu-PSMA-I&T underlines the therapeutic potential for the treatment of PCa.

Towards harmonizing clinical linear energy transfer (LET) reporting in proton radiotherapy: a European multi-centric study.

BACKGROUND: Clinical data suggest that the relative biological effectiveness (RBE) in proton therapy (PT) varies with linear energy transfer (LET). However, LET calculations are neither standardized nor available in clinical routine. Here, the status of LET calculations among European PT institutions and their comparability are assessed. MATERIALS AND METHODS: Eight European PT institutions used suitable treatment planning systems with their center-specific beam model to create treatment plans in a water phantom covering different field arrangements and fulfilling commonly agreed dose objectives. They employed their locally established LET simulation environments and procedures to determine the corresponding LET distributions. Dose distributions D1.1 and DRBE assuming constant and variable RBE, respectively, and LET were compared among the institutions. Inter-center variability was assessed based on dose- and LET-volume-histogram parameters. RESULTS: Treatment plans from six institutions fulfilled all clinical goals and were eligible for common analysis. D1.1 distributions in the target volume were comparable among PT institutions. However, corresponding LET values varied substantially between institutions for all field arrangements, primarily due to differences in LET averaging technique and considered secondary particle spectra. Consequently, DRBE using non-harmonized LET calculations increased inter-center dose variations substantially compared to D1.1 and significantly in mean dose to the target volume of perpendicular and opposing field arrangements (p < 0.05). Harmonizing LET reporting (dose-averaging, all protons, LET to water or to unit density tissue) reduced the inter-center variability in LET to the order of 10-15% within and outside the target volume for all beam arrangements. Consequentially, inter-institutional variability in DRBE decreased to that observed for D1.1. CONCLUSION: Harmonizing the reported LET among PT centers is feasible and allows for consistent multi-centric analysis and reporting of tumor control and toxicity in view of a variable RBE. It may serve as basis for harmonized variable RBE dose prescription in PT.

Impact of RBE variations on risk estimates of temporal lobe necrosis in patients treated with intensity-modulated proton therapy for head and neck cancer.

BACKGROUND: Temporal lobe necrosis (TLN) is a potential late effect after radiotherapy for skull base head and neck cancer (HNC). Several photon-derived dose constraints and normal tissue complication probability (NTCP) models have been proposed, however variation in relative biological effectiveness (RBE) may challenge the applicability of these dose constraints and models in proton therapy. The purpose of this study was therefore to investigate the influence of RBE variations on risk estimates of TLN after Intensity-Modulated Proton Therapy for HNC. MATERIAL AND METHODS: Seventy-five temporal lobes from 45 previously treated patients were included in the analysis. Sixteen temporal lobes had radiation associated Magnetic Resonance image changes (TLIC) suspected to be early signs of TLN. Fixed (RWDFix) and variable RBE-weighed doses (RWDVar) were calculated using RBE = 1.1 and two RBE models, respectively. RWDFix and RWDVar for temporal lobes were compared using Friedman's test. Based on RWDFix, six NTCP models were fitted and internally validated through bootstrapping. Estimated probabilities from RWDFix and RWDVar were compared using paired Wilcoxon test. Seven dose constraints were evaluated separately for RWDFix and RWDVar by calculating the observed proportion of TLIC in temporal lobes meeting the specific dose constraints. RESULTS: RWDVar were significantly higher than RWDFix (p < 0.01). NTCP model performance was good (AUC:0.79-0.84). The median difference in estimated probability between RWDFix and RWDVar ranged between 5.3% and 20.0% points (p < 0.01), with V60GyRBE and DMax at the smallest and largest differences, respectively. The proportion of TLIC was higher for RWDFix (4.0%-13.1%) versus RWDVar (1.3%-5.3%). For V65GyRBE ≤ 0.03 cc the proportion of TLIC was less than 5% for both RWDFix and RWDVar. CONCLUSION: NTCP estimates were significantly influenced by RBE variations. Dmax as model predictor resulted in the largest deviations in risk estimates between RWDFix and RWDVar. V65GyRBE ≤ 0.03 cc was the most consistent dose constraint for RWDFix and RWDVar.

Development of modified microdosimetric kinetic model for relative biological effectiveness in proton therapy.

To predict the biological effects of ionising radiation, the quantity of biological dose is introduced instead of the physical absorbed dose. In proton therapy, a constant relative biological effectiveness (RBE) of 1.1 is usually applied clinically as recommended by the International Commission of Radiation Units and Measurements. This study presents a new model, based on the modified microdosimetric kinetic model (MMKM), for calculating variable RBE values based on experimental data on the induction of DNA double-strand breaks (DSBs) within cells. The MMKM was proposed based on experimental data for the yield of DSBs in mammalian cells, which allows modification of the yield of primary lesions in the MKM. In this approach, a unique function named f(LET), which describes the relation between RBE and linear energy transfer (LET), was considered for charged particles. In the presented model (DMMKM), the MMKM approach was developed further by considering different f(LET)s for different relevant ions involved in energy deposition events in proton therapy. Although experimental data represent the dependence of the yield of primary lesions on the ion species, the DSB yield (assumed as the main primary lesion) is assumed independent of the ion species in the MMKM. In the DMMKM, by considering the yield of primary lesions as a function of the ion species, the α and ß values are in better agreement with the experimental data as compared to those of the MKM and MMKM approaches. The biological dose in the DMMKM is predicted to be lower than that in the MMKM. Further, in the proposed model, the variation of the ß parameter is higher than the constant value assumed in the MKM, at the distal end of the spread-out Bragg peak (SOBP). Moreover, the level of cell death estimated by the MMKM at the SOBP region is higher than that obtained based on the DMMKM. It is concluded that considering modified f(LET)s in the model developed here is more consistent with experimental results than when MMKM and MKM approaches are considered. The DMMKM examines the biological effects with full detail and will, therefore, be effective in improving proton therapy.

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.

The Mayo Clinic Florida Microdosimetric Kinetic Model of Clonogenic Survival: Application to Various Repair-Competent Rodent and Human Cell Lines.

The relative biological effectiveness (RBE) calculations used during the planning of ion therapy treatments are generally based on the microdosimetric kinetic model (MKM) and the local effect model (LEM). The Mayo Clinic Florida MKM (MCF MKM) was recently developed to overcome the limitations of previous MKMs in reproducing the biological data and to eliminate the need for ion-exposed in vitro data as input for the model calculations. Since we are considering to implement the MCF MKM in clinic, this article presents (a) an extensive benchmark of the MCF MKM predictions against corresponding in vitro clonogenic survival data for 4 rodent and 10 cell lines exposed to ions from 1H to 238U, and (b) a systematic comparison with published results of the latest version of the LEM (LEM IV). Additionally, we introduce a novel approach to derive an approximate value of the MCF MKM model parameters by knowing only the animal species and the mean number of chromosomes. The overall good agreement between MCF MKM predictions and in vitro data suggests the MCF MKM can be reliably used for the RBE calculations. In most cases, a reasonable agreement was found between the MCF MKM and the LEM IV.

Assessing the DNA Damaging Effectiveness of Ionizing Radiation Using Plasmid DNA.

Plasmid DNA is useful for investigating the DNA damaging effects of ionizing radiation. In this study, we have explored the feasibility of plasmid DNA-based detectors to assess the DNA damaging effectiveness of two radiotherapy X-ray beam qualities after undergoing return shipment of ~8000 km between two institutions. The detectors consisted of 18 µL of pBR322 DNA enclosed with an aluminum seal in nine cylindrical cavities drilled into polycarbonate blocks. We shipped them to Toronto, Canada for irradiation with either 100 kVp or 6 MV X-ray beams to doses of 10, 20, and 30 Gy in triplicate before being shipped back to San Diego, USA. The Toronto return shipment also included non-irradiated controls and we kept a separate set of controls in San Diego. In San Diego, we quantified DNA single strand breaks (SSBs), double strand breaks (DSBs), and applied Nth and Fpg enzymes to quantify oxidized base damage. The rate of DSBs/Gy/plasmid was 2.8±0.7 greater for the 100 kVp than the 6 MV irradiation. The 100 kVp irradiation also resulted in 5±2 times more DSBs/SSB than the 6 MV beam, demonstrating that the detector is sensitive enough to quantify relative DNA damage effectiveness, even after shipment over thousands of kilometers.

Variation in relative biological effectiveness for cognitive structures in proton therapy of pediatric brain tumors.

BACKGROUND: Clinically, a constant value of 1.1 is used for the relative biological effectiveness (RBE) of protons, whereas in vitro the RBE has been shown to vary depending on physical dose, tissue type, and linear energy transfer (LET). As the LET increases at the distal end of the proton beam, concerns exist for an elevated RBE in normal tissues. The aim of this study was therefore to investigate the heterogeneity of RBE to brain structures associated with cognition (BSCs) in pediatric suprasellar tumors. MATERIAL AND METHODS: Intensity-modulated proton therapy (IMPT) plans for 10 pediatric craniopharyngioma patients were re-calculated using 11 phenomenological and two plan-based variable RBE models. Based on LET, tissue dependence and number of data points used to fit the models, the three RBE models considered the most relevant for the studied endpoint were selected. Thirty BSCs were investigated in terms of RBE and dose/volume parameters. RESULTS: For a representative patient, the median (range) dose-weighted mean RBE (RBEd) across all BSCs from the plan-based models was among the lowest (1.09 (1.02-1.52) vs. the phenomenological models at 1.21 (0.78-2.24)). Omitting tissue dependency resulted in RBEd at 1.21 (1.04-2.24). Across all patients, the narrower RBE model selection gave median RBEd values from 1.22 to 1.30. CONCLUSION: For all BSCs, there was a systematic model-dependent variation in RBEd, mirroring the uncertainty in biological effects of protons. According to a refined selection of in vitro models, the RBE variation across BSCs was in effect underestimated when using a fixed RBE of 1.1.

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.

A Mathematical Radiobiological Model (MRM) to Predict Complex DNA Damage and Cell Survival for Ionizing Particle Radiations of Varying Quality.

Predicting radiobiological effects is important in different areas of basic or clinical applications using ionizing radiation (IR); for example, towards optimizing radiation protection or radiation therapy protocols. In this case, we utilized as a basis the 'MultiScale Approach (MSA)' model and developed an integrated mathematical radiobiological model (MRM) with several modifications and improvements. Based on this new adaptation of the MSA model, we have predicted cell-specific levels of initial complex DNA damage and cell survival for irradiation with 11Β, 12C, 14Ν, 16Ο, 20Νe, 40Αr, 28Si and 56Fe ions by using only three input parameters (particle's LET and two cell-specific parameters: the cross sectional area of each cell nucleus and its genome size). The model-predicted survival curves are in good agreement with the experimental ones. The particle Relative Biological Effectiveness (RBE) and Oxygen Enhancement Ratio (OER) are also calculated in a very satisfactory way. The proposed integrated MRM model (within current limitations) can be a useful tool for the assessment of radiation biological damage for ions used in hadron-beam radiation therapy or radiation protection purposes.

Assessment of the Relative Biological Efficiency of Pencil Beam Scanning of Protons in Mice in Vivo.

The effect of proton pencil beam scanning in the dose range of 4.5-15 Gy on the radiosensitivity of mice under irradiation in two regions of the Bragg curve was studied according to the criteria of 30-day survival, dynamics of death, and average lifespan of mice. The relative biological effectiveness (RBE) value of protons relative to X-ray radiation before and at the Bragg peak determined by the LD50/30 index was 0.86 and 0.94, respectively, and by the criterion of 30-day survival at a dose of 6.5 Gy it was 0.83 and 0.84, respectively. With similar RBE values for protons in different regions of the Bragg curve, significant differences in the dynamics of the course of radiation sickness were revealed, which indicates different damage to critical systems and organs of animals and the induction of compensatory mechanisms involved in the formation of stress responses at the organismal level.