Parisi, Alessio

Abstract Purpose Modeling relative biological effectiveness (RBE) is central to carbon ion radiotherapy treatment planning. The modified Microdosimetric Kinetic Model (mMKM) is a clinically established RBE framework that has guided treatment protocols at many existing carbon centers, while the Mayo Clinic Florida Microdosimetric Kinetic Model (MCF MKM) is a recently developed alternative. This work aims to implement the MCF MKM in the open‐source treatment planning system matRad and to quantitatively compare its RBE‐weighted dose predictions with those of the clinically established mMKM using identical physical dose distributions across multiple disease sites. These findings will help assess their dosimetric equivalence and inform protocol development for carbon ion radiotherapy at MCF. Methods Monte Carlo simulations of the MCF carbon beamline were performed to generate physical (IDD, LET, lateral spread) and biological base data for integration of each RBE model into matRad . Treatment plans were generated for six patients, each corresponding to a different disease site, using clinical beam configurations with carbon‐reference dose prescriptions, and plans were optimized using the MCF MKM. To isolate differences attributable solely to the RBE model, the resulting physical dose distributions were held fixed and RBE‐weighted doses were recalculated using the mMKM. Dose volume histogram (DVH) metrics and spatial dose‐difference maps were used to compare target coverage and organ‐at‐risk doses between the two models. Results Across patient cases, RBE‐weighted dose distributions from MCF MKM and mMKM showed strong agreement. Differences in target coverage were small, with CTV D95% differing by less than 1.6% across all disease sites and maximum target dose differences not exceeding 0.88%. Organ‐at‐risk dose deviations were limited, with differences of 3.0% or less across evaluated DVH metrics. Spatial dose‐difference maps showed that the largest discrepancies occurred in regions of steep dose gradients near target to organ‐at‐risk interfaces, while overall dose conformity and plan quality remained comparable between the two models. Conclusions: This study served as the first systematic model comparison of the MCF MKM and mMKM within a treatment planning environment. These findings suggest that the MCF MKM and mMKM produce dosimetrically consistent RBE‐weighted dose predictions under realistic planning conditions using carbon‐reference parameters. Accordingly, the fractionation schemes developed from years of clinical experience with mMKM implementations may serve as a practical foundation for protocol development at MCF.

Abstract Background Toxicities to organs at risk adjacent to the target volume represent a critical limitation to dose escalation in radiotherapy. While external beam therapy with accelerated ions allows for a reduced dose to surrounding healthy organs, its biological effectiveness can be markedly higher than that of conventional x‐rays and should be considered during treatment planning. Current biophysical models possess limitations in computing in vivo normal tissue complication probability (NTCP) curves for clinically relevant ion beams. Purpose This study introduces and validates a novel computational framework capable of predicting NTCP curves for serial organs following ion exposures. Methods The proposed framework integrates the linear‐quadratic model, the critical volume NTCP model, a novel two‐stage fitting process of fractionated photon NTCP data, histological data, and the Mayo Clinic Florida Microdosimetric Kinetic Model (MCF MKM) of relative biological effectiveness. The model was applied to predict the NTCP for the spinal cord myelopathy of rats exposed at different depths along the spread‐out Bragg peaks (SOBPs) of protons, helium, carbon, and oxygen ions. Radiation transport simulations were conducted to assess the microdosimetric pattern of energy deposition at the subcellular scale. Results The two‐stage photon fitting process yielded a α photons /β photons ratio of 2.3 ± 0.6 Gy and NTCP model parameters for calculations considering either endothelial or glial cells as the primary targets for the spinal cord myelopathy. Without using any ion data for calibration, the predicted NTCP showed an overall good agreement with the corresponding animal data. The model successfully reproduced the experimentally observed trends in the NTCP dependence with respect to fractionation, ion type, and exposure position along the SOBPs. Conclusions The model's ability to predict the NTCP curves for ion exposures, based solely on photon in vivo data, provides mechanistic insights and can guide preclinical investigations into the biological effects of different ions and fractionation schemes.

Abstract Purpose This study investigates a novel optimization framework that integrates the newly developed MCF MKM relative biological effectiveness (RBE) model into intensity modulated particle therapy (IMPT) planning for carbon ion therapy. Unlike conventional RBE models such as MKM or LEM, the MCF MKM provides detailed microdosimetric information, which improves predictive accuracy of cell survival but elevates both memory and time requirement in the optimization process. Furthermore, to the implicit nature of the MCF MKM, the gradient information, which is essential in inverse planning, is hardly accessible. In response, a randomized gradient‐free optimization method (RGFM) was developed to mitigate the computational challenge and achieve improved RBE dose distributions. Materials and Methods An RBE dose optimization problem was formulated based on the discretized MCF MKM, known as the Abridged Microdosimetry Distribution Model (AMDM). The RGFM was derived to find feasible and high‐quality settings of monitor unit for IMPT plans without the need for explicit gradient calculations. To verify the clinical plausibility of the proposed method, it was tested on two clinical scenarios—prostate and brain cases in comparison with benchmark plans derived solely from physical dose optimization. Computational times and dose‐volume histograms (DVHs) for both physical and RBE doses were evaluated. Results The RGFM algorithm successfully optimized the MCF MKM‐based IMPT plans. Compared to the benchmark, the RBE dose distributions showed improved coverage of the target and reduced hot spots. Such improvements were observed to coincide with trade‐offs in physical dose distributions. The complexity of the MCF MKM extended the computational time, yet high‐quality solutions were attainable in a clinically reasonable timeframe. Multi‐field optimization (MFO) further enhanced target coverage and critical organ sparing over single‐field optimization (SFO). Conclusions This proof‐of‐principle study confirmed the feasibility of using a randomized gradient‐free method to optimize RBE dose based on the MCF MKM in IMPT planning for carbon ion therapy. The approach provides a pathway for clinical adoption of advanced RBE models, potentially leading to more biologically informed and effective carbon ion therapy treatments. Future work involves refining algorithms, incorporating parallel computing, and performing additional case studies to improve computational efficiency and generalizability.

Microdosimetry plays a critical role in particle therapy by quantifying energy deposition within microscopic domains to assess biological effects. This study evaluates the influence of different microdosimetric functions (MFs) and domain geometries (DGs) on relative biological effectiveness (RBE) predictions in carbon ion radiotherapy. Specifically, we compare the analytical microdosimetric function (AMF), calculated for spherical domains and implemented in PHITS, with the Kiefer–Chatterjee (KC) track structure model, which is conventionally applied to cylindrical geometries. To enable a direct comparison, we also introduce a novel implementation of the KC model for spherical domains. Using both models, specific energy distributions were calculated across a range of domain sizes and geometries. These distributions were input into the modified microdosimetric kinetic model (mMKM) to calculate RBE for the HSG cell line and compared against published in vitro data. The results show that both microdosimetric function and domain geometry significantly affect microdosimetric spectra and the resulting RBE, with deviations exceeding 10% when fixed mMKM parameters are used. Parameter optimization within the mMKM enables alignment across models. Our findings emphasize that microdosimetric function and domain geometry selection must be explicitly accounted for in microdosimetry-based RBE modeling, and that model parameters must be tuned accordingly to ensure consistent and biologically accurate predictions.

Abstract Background Projectile and target fragmentation are nuclear phenomena that can influence the computation of the linear energy transfer (LET) and the relative biological effectiveness (RBE) in external radiotherapy with accelerated ions. Correlations between these two quantities are routinely established during radiobiological experiments to interpret the results and to develop and calibrate RBE models. Purpose This study systematically evaluates the impact of secondary fragments and irradiation scenarios on the dose‐averaged LET, the RBE, and their correlation in the case of exposures to clinically relevant ion beams. Methods 57 600 microdosimetric lineal energy spectra and corresponding LET distributions were simulated with the Monte Carlo code PHITS across different scenarios, including track segment calculations, pristine, and spread‐out Bragg peaks of 1 H, 4 He, 12 C, 16 O, and 20 Ne ions within water phantoms. The LET distributions were analyzed to calculate the dose‐average LET, both including or excluding the contribution of secondary ions of an element different from the primary beam. Similarly, the lineal energy distributions were processed in conjunction with the Mayo Clinic Florida microdosimetric kinetic model to compute the RBE for two theoretical cell lines ( α / β = 2 and 10 Gy in the case of 6 MV x‐rays). Correlations between the RBE and dose‐averaged LET were established by analyzing the simulation results within water phantoms and then compared to the corresponding trends from the track segment calculations. Results Excluding secondary fragments had a pronounced impact on the dose‐averaged LET and the RBE, particularly in the entrance regions of proton beams and close to the distal edge of heavier ions. The correlations between the RBE and the dose‐averaged LET were not universal, but highly dependent on the irradiation scenario. For proton beams only, the dose‐averaged LET of hydrogen ions served as a practical first‐order descriptor of the RBE. Track segment simulations, commonly used for calibrating and benchmarking RBE models, provided a reasonable approximation for low‐energy beams but failed to fully capture the complexity of polyenergetic radiation fields. Conclusions Secondary fragments can substantially affect the dose‐averaged LET and the RBE, even in proton beams. The dose‐averaged LET, including or not the contributions from secondary fragments, was generally unable to adequately capture RBE variations across different scenarios. A more comprehensive approach, considering microdosimetric distributions, is necessary to accurately describe RBE variations in ion therapy.

Ion radiotherapy requires accurate relative biological effectiveness (RBE) calculations to account for the markedly different biological effects of ions compared to photons. Microdosimetric RBE models rely on descriptions of the energy deposition at the microscopic scale, either through radial dose distributions (RDDs) or microdosimetric probability density distributions. While RDD approaches focus on the theoretical description of the energy deposition around the ion track, microdosimetric distributions offer the advantage of being experimentally measurable, which is crucial for quality assurance programs. As the results of microdosimetric RBE models depend on whether RDD or microdosimetric distributions are used, the model parameters are not interchangeable between these approaches. This study presents and validates a method to reproduce the published reference biological and clinical dose calculations at NIRS-QST for only carbon ion radiotherapy by using the modified microdosimetric kinetic model (MKM) alongside microdosimetric distributions instead of the reference RDD approach. To achieve this, Monte Carlo simulations were performed to estimate the variation of the radiation quality within and outside the field of pristine and spread-out Bragg peaks. By appropriately optimizing the modified MKM parameters for microdosimetric distributions assessed within water spheres, we successfully reproduced the results of calculations using the reference NIRS-QST RDD, generally within 2%.

Abstract The repair of DNA double-strand breaks is a crucial yet delicate process which is affected by a multitude of factors. In this study, our goal is to analyse the influence of the linear energy transfer (LET) on the DNA repair kinetics. By utilizing the database of repair of DNA and aggregating the results of 84 experiments, we conduct various model fits to evaluate and compare different hypothesis regarding the effect of LET on the rejoining of DNA ends. Despite the considerable research efforts dedicated to this topic over the past decades, our findings underscore the complexity of the relationship between LET and DNA repair kinetics. This study leverages big data analysis to capture overall trends that single experimental studies might miss, providing a valuable model for understanding how radiation quality impacts DNA damage and subsequent biological effects. Our results highlight the gaps in our current understanding, emphasizing the pressing need for further investigation into this phenomenon.

AbstractBackgroundAdaptive radiotherapy (ART) workflows have been increasingly adopted to achieve dose escalation and tissue sparing under shifting anatomic conditions, but the necessity of recontouring and the associated time burden hinders a real‐time or online ART workflow. In response to this challenge, approaches to auto‐segmentation involving deformable image registration, atlas‐based segmentation, and deep learning‐based segmentation (DLS) have been developed. Despite the particular promise shown by DLS methods, implementing these approaches in a clinical setting remains a challenge, namely due to the difficulty of curating a data set of sufficient size and quality so as to achieve generalizability in a trained model.PurposeTo address this challenge, we have developed an intentional deep overfit learning (IDOL) framework tailored to the auto‐segmentation task. However, certain limitations were identified, particularly the insufficiency of the personalized dataset to effectively overfit the model. In this study, we introduce a personalized hyperspace learning (PHL)‐IDOL segmentation framework capable of generating datasets that induce the model to overfit specific patient characteristics for medical image segmentation.MethodsThe PHL‐IDOL model is trained in two stages. In the first, a conventional, general model is trained with a diverse set of patient data (n = 100 patients) consisting of CT images and clinical contours. Following this, the general model is tuned with a data set consisting of two components: (a) selection of a subset of the patient data (m < n) using the similarity metrics (mean square error (MSE), peak signal‐to‐noise ratio (PSNR), structural similarity index (SSIM), and the universal quality image index (UQI) values); (b) adjust the CT and the clinical contours using a deformed vector generated from the reference patient and the selected patients using (a). After training, the general model, the continual model, the conventional IDOL model, and the proposed PHL‐IDOL model were evaluated using the volumetric dice similarity coefficient (VDSC) and the Hausdorff distance 95% (HD95%) computed for 18 structures in 20 test patients.ResultsImplementing the PHL‐IDOL framework resulted in improved segmentation performance for each patient. The Dice scores increased from 0.810.05 with the general model, 0.83 for the continual model, 0.83 for the conventional IDOL model to an average of 0.87 with the PHL‐IDOL model. Similarly, the Hausdorff distance decreased from 3.06 with the general model, 2.84 for the continual model, 2.79 for the conventional IDOL model and 2.36 for the PHL‐IDOL model. All the standard deviations were decreased by nearly half of the values comparing the general model and the PHL‐IDOL model.ConclusionThe PHL‐IDOL framework applied to the auto‐segmentation task achieves improved performance compared to the general DLS approach, demonstrating the promise of leveraging patient‐specific prior information in a task central to online ART workflows.

AbstractBackgroundPhenomenological relative biological effectiveness (RBE) models for proton therapy, based on the dose‐averaged linear energy transfer (LET), have been developed to address the apparent RBE increase towards the end of the proton range. The results of these phenomenological models substantially differ due to varying empirical assumptions and fitting functions. In contrast, more theory‐based approaches are used in carbon ion radiotherapy, such as the microdosimetric kinetic model (MKM). However, implementing microdosimetry‐based models in LET‐based proton therapy treatment planning systems poses challenges.PurposeThis work presents a LET‐based version of the MKM that is practical for clinical use in proton radiotherapy.MethodsAt first, we derived an approximation of the Mayo Clinic Florida (MCF) MKM for relatively‐sparsely ionizing radiation such as protons. The mathematical formalism of the proposed model is equivalent to the original MKM, but it maintains some key features of the MCF MKM, such as the determination of model parameters from measurable cell characteristics. Subsequently, we carried out Monte Carlo calculations with PHITS in different simulated scenarios to establish a heuristic correlation between microdosimetric quantities and the dose averaged LET of protons.ResultsA simple allometric function was found able to describe the relationship between the dose‐averaged LET of protons and the dose‐mean lineal energy, which includes the contributions of secondary particles. The LET‐based MKM was used to model the in vitro clonogenic survival RBE of five human and rodent cell lines (A549, AG01522, CHO, T98G, and U87) exposed to pristine and spread‐out Bragg peak (SOBP) proton beams. The results of the LET‐based MKM agree well with the biological data in a comparable or better way with respect to the other models included in the study. A sensitivity analysis on the model results was also performed.ConclusionsThe LET‐based MKM integrates the predictive theoretical framework of the MCF MKM with a straightforward mathematical description of the RBE based on the dose‐averaged LET, a physical quantity readily available in modern treatment planning systems for proton therapy.