Microdosimetric Kinetic Model Research Articles

Biological dose optimization incorporating intra-tumoural cellular radiosensitivity heterogeneity in ion-beam therapy treatment planning

Objective. Treatment plans of ion-beam therapy have been made under an assumption that all cancer cells within a tumour equally respond to a given radiation dose. However, an intra-tumoural cellular radiosensitivity heterogeneity clearly exists, and it may lead to an overestimation of therapeutic effects of the radiation. The purpose of this study is to develop a biological model that can incorporate the radiosensitivity heterogeneity into biological optimization for ion-beam therapy treatment planning. Approach. The radiosensitivity heterogeneity was modeled as the variability of a cell-line specific parameter in the microdosimetric kinetic model following the gamma distribution. To validate the developed intra-tumoural-radiosensitivity-heterogeneity-incorporated microdosimetric kinetic (HMK) model, a treatment plan with H-ion beams was made for a chordoma case, assuming a radiosensitivity heterogeneous region within the tumour. To investigate the effects of the radiosensitivity heterogeneity on the biological effectiveness of H-, He-, C-, O-, and Ne-ion beams, the relative biological effectiveness (RBE)-weighted dose distributions were planned for a cuboid target with the stated ion beams without considering the heterogeneity. The planned dose distributions were then recalculated by taking the heterogeneity into account. Main results. The cell survival fraction and corresponding RBE-weighted dose were formulated based on the HMK model. The first derivative of the RBE-weighted dose distribution was also derived, which is needed for fast biological optimization. For the patient plan, the biological optimization increased the dose to the radiosensitivity heterogeneous region to compensate for the heterogeneity-induced reduction in biological effectiveness of the H-ion beams. The reduction in biological effectiveness due to the heterogeneity was pronounced for low linear energy transfer (LET) beams but moderate for high-LET beams. The RBE-weighted dose in the cuboid target decreased by 7.6% for the H-ion beam, while it decreased by just 1.4% for the Ne-ion beam. Significance. Optimal treatment plans that consider intra-tumoural cellular radiosensitivity heterogeneity can be devised using the HMK model.

Evaluation of Helium Ion Radiotherapy in Combination with Gemcitabine in Pancreatic Cancer In Vitro.

Pancreatic cancer is one of the most aggressive and lethal cancers. New treatment strategies are highly warranted. Particle radiotherapy could offer a way to overcome the radioresistant nature of pancreatic cancer because of its biological and physical characteristics. Within particles, helium ions represent an attractive therapy option to achieve the highest possible conformity while at the same time protecting the surrounding normal tissue. The aim of this study was to evaluate the cytotoxic efficacy of helium ion irradiation in pancreatic cancer in vitro. Human pancreatic cancer cell lines AsPC-1, BxPC-3 and Panc-1 were irradiated with photons and helium ions at various doses and treated with gemcitabine. Photon irradiation was performed with a biological cabin X-ray irradiator, and helium ion irradiation was performed with a spread-out Bragg peak using the raster scanning technique at the Heidelberg Ion Beam Therapy Center (HIT). The cytotoxic effect on pancreatic cancer cells was measured with clonogenic survival. The survival curves were compared to the predicted curves that were calculated via the modified microdosimetric kinetic model (mMKM). The experimental relative biological effectiveness (RBE) of helium ion irradiation ranged from 1.0 to 1.7. The predicted survival curves obtained via mMKM calculations matched the experimental survival curves. Mainly additive cytotoxic effects were observed for the cell lines AsPC-1, BxPC-3 and Panc-1. Our results demonstrate the cytotoxic efficacy of helium ion radiotherapy in pancreatic cancer in vitro as well as the capability of mMKM calculation and its value for biological plan optimization in helium ion therapy for pancreatic cancer. A combined treatment of helium irradiation and chemotherapy with gemcitabine leads to mainly additive cytotoxic effects in pancreatic cancer cell lines. The data generated in this study may serve as the radiobiological basis for future experimental and clinical works using helium ion radiotherapy in pancreatic cancer treatment.

Calculation of biological effectiveness of SOBP proton beams: a TOPAS Monte Carlo study

Objective. This study aims to investigate the biological effectiveness of Spread-Out Bragg-Peak (SOBP) proton beams with initial kinetic energies 50–250 MeV at different depths in water using TOPAS Monte Carlo code. Approach. The study modelled SOBP proton beams using TOPAS time feature. Various LET-based models and Repair-Misrepair-Fixation model were employed to calculate Relative Biological Effectiveness (RBE) for V79 cell lines at different on-axis depths based on TOPAS. Microdosimetric Kinetic Model and biological weighting function-based models, which utilize microdosimetric distributions, were also used to estimate the RBE. A phase-space-based method was adopted for calculating microdosimetric distributions. Main results. The trend of variation of RBE with depth is similar in all the RBE models, but the absolute RBE values vary based on the calculation models. RBE sharply increases at the distal edge of SOBP proton beams. In the entrance region of all the proton beams, RBE values at 4 Gy i.e. RBE(4 Gy) resulting from different models are in the range of 1.04–1.07, comparable to clinically used generic RBE of 1.1. Moving from the proximal to distal end of the SOBP, RBE(4 Gy) is in the range of 1.15–1.33, 1.13–1.21, 1.11–1.17, 1.13–1.18 and 1.17–1.21, respectively for 50, 100, 150, 200 and 250 MeV SOBP beams, whereas at the distal dose fall-off region, these values are 1.68, 1.53, 1.44, 1.42 and 1.40, respectively. Significance. The study emphasises application of depth-, dose- and energy- dependent RBE values in clinical application of proton beams.

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

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. 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 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. 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.

Nanodosimetric quantity-weighted dose optimization for carbon-ion treatment planning.

Dose verification of treatment plans is an essential step in radiotherapy workflows. In this work, we propose a novel method of treatment planning based on nanodosimetric quantity-weighted dose (NQWD), which could realize biological representation using pure physical quantities for biological-oriented carbon ion-beam treatment plans and their direct verification. The relationship between nanodosimetric quantities and relative biological effectiveness (RBE) was studied with the linear least-squares method for carbon-ion radiation fields. Next, under the framework of the matRad treatment planning platform, NQWD was optimized using the existing RBE-weighted dose (RWD) optimization algorithm. The schemes of NQWD-based treatment planning were compared with the RWD treatment plans in term of the microdosimetric kinetic model (MKM). The results showed that the nanodosimetric quantity F3 - 10 had a good correlation with the radiobiological effect reflected by the relationship between RBE and F3 - 10. Moreover, the NQWD-based treatment plans reproduced the RWD plans generally. Therefore, F3 - 10 could be adopted as a radiation quality descriptor for carbon-ion treatment planning. The novel method proposed herein not only might be helpful for rapid physical verification of biological-oriented ion-beam treatment plans with the development of experimental nanodosimetry, but also makes the direct comparison of ion-beam treatment plans in different institutions possible. Thus, our proposed method might be potentially developed to be a new strategy for carbon-ion treatment planning and improve patient safety for carbon-ion radiotherapy.

Integrating microdosimetricin vitroRBE models for particle therapy into TOPAS MC using the MicrOdosimetry-based modeliNg for RBE ASsessment (MONAS) tool.

Objective.In this paper, we present MONAS (MicrOdosimetry-based modelliNg for relative biological effectiveness (RBE) ASsessment) toolkit. MONAS is a TOPAS Monte Carlo extension, that combines simulations of microdosimetric distributions with radiobiological microdosimetry-based models for predicting cell survival curves and dose-dependent RBE.Approach.MONAS expands TOPAS microdosimetric extension, by including novel specific energy scorers to calculate the single- and multi-event specific energy microdosimetric distributions at different micrometer scales. These spectra are used as physical input to three different formulations of themicrodosimetric kinetic model, and to thegeneralized stochastic microdosimetric model(GSM2), to predict dose-dependent cell survival fraction and RBE. MONAS predictions are then validated against experimental microdosimetric spectra andin vitrosurvival fraction data. To show the MONAS features, we present two different applications of the code: (i) the depth-RBE curve calculation from a passively scattered proton SOBP and monoenergetic12C-ion beam by using experimentally validated spectra as physical input, and (ii) the calculation of the 3D RBE distribution on a real head and neck patient geometry treated with protons.Main results.MONAS can estimate dose-dependent RBE and cell survival curves from experimentally validated microdosimetric spectra with four clinically relevant radiobiological models. From the radiobiological characterization of a proton SOBP and12C fields, we observe the well-known trend of increasing RBE values at the distal edge of the radiation field. The 3D RBE map calculated confirmed the trend observed in the analysis of the SOBP, with the highest RBE values found in the distal edge of the target.Significance.MONAS extension offers a comprehensive microdosimetry-based framework for assessing the biological effects of particle radiation in both research and clinical environments, pushing closer the experimental physics-based description to the biological damage assessment, contributing to bridging the gap between a microdosimetric description of the radiation field and its application in proton therapy treatment with variable RBE.

Effects of cellular radioresponse on therapeutic helium-, carbon-, oxygen-, and neon-ion beams: a simulation study

Objective. Helium, oxygen, and neon ions in addition to carbon ions will be used for hypofractionated multi-ion therapy to maximize the therapeutic effectiveness of charged-particle therapy. To use new ions in cancer treatments based on the dose-fractionation protocols established in carbon-ion therapy, this study examined the cell-line-specific radioresponse to therapeutic helium-, oxygen-, and neon-ion beams within wide dose ranges. Approach. Response of cells to ions was described by the stochastic microdosimetric kinetic model. First, simulations were made for the irradiation of one-field spread-out Bragg peak beams in water with helium, carbon, oxygen, and neon ions to achieve uniform survival fractions at 37%, 10%, and 1% for human salivary gland tumor (HSG) cells, the reference cell line for the Japanese relative biological effectiveness weighted dose system, within the target region defined at depths from 90 to 150 mm. The HSG cells were then replaced by other cell lines with different radioresponses to evaluate differences in the biological dose distributions of each ion beam with respect to those of carbon-ion beams. Main results. For oxygen- and neon-ion beams, the biological dose distributions within the target region were almost equivalent to those of carbon-ion beams, differing by less than 5% in most cases. In contrast, for helium-ion beams, the biological dose distributions within the target region were largely different from those of carbon-ion beams, more than 10% in several cases. Significance. From the standpoint of tumor control evaluated by the clonogenic cell survival, this study suggests that the dose-fractionation protocols established in carbon-ion therapy could be reasonably applied to oxygen- and neon-ion beams while some modifications in dose prescription would be needed when the protocols are applied to helium-ion beams. This study bridges the gap between carbon-ion therapy and hypofractionated multi-ion therapy.

Understanding Relative Biological Effectiveness and Clinical Outcome of Prostate Cancer Therapy Using Particle Irradiation: Analysis of Tumor Control Probability With the Modified Microdosimetric Kinetic Model

Purpose: Recent experimental studies and clinical trial results might indicate that - at least for some indications - continued use of the mechanistic model for relative biological effectiveness (RBE) applied at carbon ion therapy facilities in Europe for several decades (LEM-I) may be unwarranted. We present a novel clinical framework for prostate cancer treatment planning and tumor control probability (TCP) prediction based on the modified microdosimetric kinetic model (mMKM) for particle therapy. Methods: Treatment plans of 91 prostate tumor patients (proton: 46, carbon ions: 45) applying 66GyRBE [RBE=1.1 for protons and LEM-I, (α/β)x=2.0Gy, for carbon ions] in 20 fractions were recalculated using mMKM ((α/β)x=3.1Gy). Based solely on the response data of photon-irradiated patient groups stratified according to risk and usage of ADT, we derived parameters for an mMKM-based Poisson-TCP model. Subsequently, new carbon and helium ion plans, adhering to prescribed biological dose criteria, were generated. These were systematically compared to the clinical experience of Japanese centers employing an analogous fractionation scheme and existing proton plans. Results: mMKM predictions suggested significant biological dose deviation between proton and carbon ion arms. Patients irradiated with protons received 3.25±0.08GyRBEmMKM/Fx, while patients treated with carbon ions received 2.51±0.05GyRBEmMKM/Fx. TCP predictions were 86±3% for protons and 52±4% for carbon ions, matching the clinical outcome of 85% and 50%. Newly optimized carbon ion plans, guided by the mMKM/TCP model, effectively replicated clinical data from Japanese centers. Using mMKM, helium ions exhibited similar target coverage as proton and carbon ions and an improved rectum and bladder sparing compared to proton. Conclusions: Our mMKM-based model for prostate cancer treatment planning and TCP prediction was validated against clinical data for proton and carbon ion therapy and its application was extended to helium ion therapy. Based on the data presented in this work, mMKM seems to be a good candidate for clinical biological calculations in carbon ion therapy for prostate cancer.

Exploring Helium Ions' Potential for Post-Mastectomy Left-Sided Breast Cancer Radiotherapy.

Proton therapy presents a promising modality for treating left-sided breast cancer due to its unique dose distribution. Helium ions provide increased conformality thanks to a reduced lateral scattering. Consequently, the potential clinical benefit of both techniques was explored. An explorative treatment planning study involving ten patients, previously treated with VMAT (Volumetric Modulated Arc Therapy) for 50 Gy in 25 fractions for locally advanced, node-positive breast cancer, was carried out using proton pencil beam therapy with a fixed relative biological effectiveness (RBE) of 1.1 and helium therapy with a variable RBE described by the mMKM (modified microdosimetric kinetic model). Results indicated that target coverage was improved with particle therapy for both the clinical target volume and especially the internal mammary lymph nodes compared to VMAT. Median dose value analysis revealed that proton and helium plans provided lower dose on the left anterior descending artery (LAD), heart, lungs and right breast than VMAT. Notably, helium therapy exhibited improved ipsilateral lung sparing over protons. Employing NTCP models as available in the literature, helium therapy showed a lower probability of grade ≤ 2 radiation pneumonitis (22% for photons, 5% for protons and 2% for helium ions), while both proton and helium ions reduce the probability of major coronary events with respect to VMAT.

Generality Assessment of a Model Considering Heterogeneous Cancer Cells for Predicting Tumor Control Probability for Stereotactic Body Radiation Therapy Against Non-Small Cell Lung Cancer

The generality of a model for predicting tumor control probability from in vitro clonogenic survival considering of cancer stem-like cells, the so-called integrated microdosimetric-kinetic model, is presented by comparing the model to public data on stereotactic body radiation therapy for non-small cell lung cancer cells.

Modeling for predicting survival fraction of cells after ultra-high dose rate irradiation

Objective. FLASH radiotherapy (FLASH-RT) with ultra-high dose rate (UHDR) irradiation (i.e. > 40 Gy s−1) spares the function of normal tissues while preserving antitumor efficacy, known as the FLASH effect. The biological effects after conventional dose rate-radiotherapy (CONV-RT) with ≤0.1 Gy s−1 have been well modeled by considering microdosimetry and DNA repair processes, meanwhile modeling of radiosensitivities under UHDR irradiation is insufficient. Here, we developed an integrated microdosimetric-kinetic (IMK) model for UHDR-irradiation enabling the prediction of surviving fraction after UHDR irradiation. Approach. The IMK model for UHDR-irradiation considers the initial DNA damage yields by the modification of indirect effects under UHDR compared to CONV dose rate. The developed model is based on the linear-quadratic (LQ) nature with the dose and dose square coefficients, considering the reduction of DNA damage yields as a function of dose rate. Main results. The estimate by the developed model could successfully reproduce the in vitro experimental dose–response curve for various cell line types and dose rates. Significance. The developed model would be useful for predicting the biological effects under the UHDR irradiation.

Microdosimetry-based investigation of biological effectiveness of 252Cf brachytherapy source: TOPAS Monte Carlo study

Objective. To investigate biological effectiveness of 252Cf brachytherapy source using Monte Carlo-calculated microdosimetric distributions. Approach. 252Cf source capsule was placed at the center of the spherical water phantom and phase-space data were scored as a function of radial distance in water (R = 1–5 cm) using TOPAS Monte Carlo code. The phase-space data were used to calculate microdosimetric distributions at 1 μm site size. Using these distributions, Relative Biological Effectiveness (RBE), mean quality factor () and Oxygen Enhancement Ratio (OER) were calculated as a function of R. Main results. The overall shapes of the microdosimetric distributions are comparable at all the radial distances in water. However, slight variation in the bin-wise yield is observed with R. RBE, and OER are insensitive to R over the range 1–5 cm. Microdosimetric kinetic model based RBE values are about 2.3 and 2.8 for HSG tumour cells and V79 cells, respectively, whereas biological weighting function-based RBE is about 2.8. ICRP 60 and ICRU 40 recommendation-based values are about 14.5 and 16, respectively. Theory of dual radiation action based RBE is 11.4. The calculated value of OER is 1.6. Significance. This study demonstrates the relative insensitivity of RBE, and OER radially away from the 252Cf source along the distances of 1–5 cm in water.

The effect of fitting the reference photon dose-response on the clonogenic survival predicted with the Mayo Clinic Florida microdosimetric kinetic model in case of accelerated ions.

The Mayo Clinic Florida microdosimetric kinetic model (MCF MKM) is a recently developed clonogenic survival model. Since the MCF MKM relies on novel strategies to a priori determine the cell-specific model parameters, the only experiment-specific input values are the α and ß terms of the linear-quadratic model (LQM) of clonogenic survival for the reference photon exposure. Because the two LQM terms are anti-correlated, the fitting process of the reference photon survival curve was found to significantly influence the MCF MKM calculations. This article reports this effect for two clinically relevant cell lines (human brain glioblastoma A-172, human healthy foreskin fibroblasts AG01522) and ions (1H and 12C ions).

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 with in vitro data 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 with in vitro data 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.

Carbon-ion radiotherapy (CIRT) as treatment of pancreatic cancer at HIT: initial radiation plan analysis of the prospective phase II PACK-study

PurposeTo analyze the dose objectives and constraints applied at the prospective phase II PACK-study at Heidelberg ion therapy center (HIT) for different radiobiological models. MethodsTreatment plans of 14 patients from the PACK-study were analyzed and recomputed in terms of physical, biological dose and dose-averaged linear energy transfer (LETd). Both LEM-I (local effect model 1) and the adapted NIRS-MKM (microdosimetric kinetic model), were used for relative biological effectiveness (RBE)-weighted dose calculations (DBio|HIT and DBio|NIRS). A new constraint to the gastrointestinal (GI) tract was derived from the National Institute of Radiological Science (NIRS) clinical experience and considered for plan reoptimization (DBio|NIRS-const_48Gy and DBio|NIRS-const_50.4Gy). The Lyman-Kutcher-Burman (LKB) model of Normal Tissue Complication Probability (NTCP) for GI toxicity endpoints was computed. Furthermore, the computed LETd distribution was evaluated and correlated with Local Control (LC). ResultsOnly two patients showed a LETd98% in the GTV greater than 44 keV/μm. A HIT-dose constraint to the GI of D2cm3=48.6GyRBEHIT was derived from the NIRS experience, in alternative to the standard at HIT Dmax = 45.6 GyRBEHIT. In comparison with the original DBio|HIT,DBio|NIRS-const_48GyandDBio|NIRS-const_50.4Gy resulted in an increase in the ITV’s D98% of 8.7% and 11.3%. The NTCP calculation resulted in a probability for gastrointestinal bleeding of 4.5%, 12.3% and 13.0%, for DBio|NIRS, DBio|NIRS-const_48Gy and DBio|NIRS-const_50.4Gy, respectively. ConclusionThe results indicate that the current standards applied at HIT for CIRT closely align with the Japanese experience. However, to enhance tumor coverage, a more relaxed constraint on the GI tract may be considered. As the PACK-trial progresses, further analyses of various clinical endpoints are anticipated.

Improvement of the hybrid approach between Monte Carlo simulation and analytical function for calculating microdosimetric probability densities in macroscopic matter

Objective. Estimation of the probability density of the microdosimetric quantities in macroscopic matter is indispensable for applying the concept of microdosimetry to medical physics and radiological protection. The Particle and Heavy Ion Transport code System (PHITS) enables estimating the microdosimetric probability densities due to its unique hybrid modality between the Monte Carlo and analytical approaches called the microdosimetric function. It can convert the deposition energies calculated by the macroscopic Monte Carlo radiation transport simulation to microdosimetric probability densities in water using an analytical function based on the track-structure simulations. Approach. In this study, we improved this function using the latest track-structure simulation codes implemented in PHITS. The improved function is capable of calculating the probability densities of not only the conventional microdosimetric quantities such as lineal energy but also the number of ionization events occurring in a target site, the so-called ionization cluster size distribution, for arbitrary site diameters from 3 nm to 1 μm. Main results. The accuracy of the improved function was well verified by comparing the microdosimetric probability densities measured by tissue-equivalent proportional counters with the corresponding data calculated in this study. Test calculations for clonogenic cell survival using the improved function coupled with the modified microdosimetric kinetic model suggested a slight increase of its relative biological effectiveness compared with our previous estimations. As a new application of the improved function, we calculated the relative biological effectiveness of the single-strand break and double-strand break yields for proton irradiations using the updated PHITS coupled with the simplified DNA damage estimation model, and confirmed its equivalence in accuracy and its superiority in computational time compared to our previously proposed method based on the track-structure simulation. Significance. From these features, we concluded that the improved function could expand the application fields of PHITS by bridging the gap between microdosimetry and macrodosimetry.

Event-by-event approach to the oxygen-effect-incorporated stochastic microdosimetric kinetic model for hypofractionated multi-ion therapy.

An oxygen-effect-incorporated stochastic microdosimetric kinetic (OSMK) model was previously developed to estimate the survival fraction of cells exposed to charged-particle beams with wide dose and linear energy transfer (LET) ranges under various oxygen conditions. In the model, hypoxia-induced radioresistance was formulated based on the dose-averaged radiation quality. This approximation may cause inaccuracy in the estimation of the biological effectiveness of the radiation with wide variation in energy deposited to a sensitive volume per event, such as spread-out Bragg peak (SOBP) beams. The purpose of this study was to apply an alternative approach so as to consider the energy depositions on an event-by-event basis. The production probability of radiation-induced lesions per energy was formulated with oxygen partial pressure to account for the hypoxia-induced radioresistance. The reduction in the oxygen enhancement ratio for high-LET radiations was modeled by reducing the sensitive-volume size and increasing the saturation energy in microdosimetry. The modified OSMK model was tested against the reported survival data of three cell lines exposed to six species of ions with wide dose and LET ranges under aerobic and hypoxic conditions. The model reasonably reproduced the reported cell survival data. To evaluate the event-by-event approach, survival distributions of Chinese hamster ovary cells exposed to SOBP beams were estimated using the original and modified OSMK models. The differences in the estimated survival distributions between the models were marginal even under extreme hypoxia. The event-by-event approach improved the theoretical validity of the OSMK model. However, the original OSMK model can still provide an accurate estimation of the biological effectiveness of therapeutic radiations.

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?

PurposeTo 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 materialOur 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. ResultsDomain 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. ConclusionsLarge 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.

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

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. 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. 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. 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. 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.

Validation of the relative biological effectiveness of active-energy scanning carbon-ion radiotherapy on a commercial treatment planning system with a microdosimetic kinetic model

BackgroundThe study objective was to validate the relative biological effectiveness (RBE) calculated by the modified microdosimetric kinetic model in RayStation (Ray-MKM) for active-energy scanning carbon-ion radiotherapy.MethodsThe Ray-MKM was benchmarked using a spread-out Bragg-peak (SOBP) plan, which was suggested in literature from the National Institute of Radiobiological Science (NIRS) in Japan. The residual RBE differences from the MKM at NIRS (NIRS-MKM) were derived using several SOBP plans with different ranges, SOBP widths, and prescriptions. To investigate the origins of the differences, we compared the saturation-corrected dose-mean specific energy Z_{1D}^{*} of the aforementioned SOBPs. Furthermore, we converted the RBE-weighted doses with the Ray-MKM to those with local effect model I (LEM doses). The purpose was to investigate whether the Ray-MKM could reproduce the RBE-weighted conversion study.ResultsThe benchmark determined the value of the clinical dose scaling factor, F_{clin}, as 2.40. The target mean RBE deviations between the Ray-MKM and NIRS-MKM were median: 0.6 (minimum: 0.0 to maximum: 1.69) %. The Z_{1D}^{*} difference in-depth led to the RBE difference in-depth and was remarkable at the distal end. The converted LEM doses from the Ray-MKM doses were comparable (the deviation being − 1.8–0.7%) to existing literature.ConclusionThis study validated the Ray-MKM based on our active-energy scanning carbon-ion beam via phantom studies. The Ray-MKM could generate similar RBEs as the NIRS-MKM after benchmarking. Analysis based on Z_{1D}^{*} indicated that the different beam qualities and fragment spectra caused the RBE differences. Since the absolute dose differences at the distal end were small, we neglected them. Furthermore, each centre may determine its centre-specific F_{clin} based on this approach.