Search

Computational models of tumor growth are valuable for simulating the dynamics of cancer progression and treatment responses. In particular, agent-based models (ABMs) tracking individual agents and their interactions are useful for their flexibility and ability to model complex behaviors. However, ABMs have often been confined to small domains or, when scaled up, have neglected crucial aspects like vasculature. Additionally, the integration into tumor ABMs of precise radiation dose calculations using gold-standard Monte Carlo (MC) methods, crucial in contemporary radiotherapy, has been lacking. Here, we introduce AMBER, an Agent-based fraMework for radioBiological Effects in Radiotherapy that computationally models tumor growth and radiation responses. AMBER is based on a voxelized geometry, enabling realistic simulations at relevant pre-clinical scales by tracking temporally discrete states stepwise. Its hybrid approach, combining traditional ABM techniques with continuous spatiotemporal fields of key microenvironmental factors such as oxygen and vascular endothelial growth factor, facilitates the generation of realistic tortuous vascular trees. Moreover, AMBER is integrated with TOPAS, an MC-based particle transport algorithm that simulates heterogeneous radiation doses. The impact of radiation on tumor dynamics considers the microenvironmental factors that alter radiosensitivity, such as oxygen availability, providing a full coupling between the biological and physical aspects. Our results show that simulations with AMBER yield accurate tumor evolution and radiation treatment outcomes, consistent with established volumetric growth laws and radiobiological understanding. Thus, AMBER emerges as a promising tool for replicating essential features of tumor growth and radiation response, offering a modular design for future expansions to incorporate specific biological traits., (© 2024. The Author(s), under exclusive licence to the Society for Mathematical Biology.)

Objective. To present and validate a method to simulate from first principles the effect of oxygen on radiation-induced double-strand breaks (DSBs) using the Monte Carlo Track-structure code TOPAS-nBio. Approach. Two chemical models based on the oxygen fixation hypothesis (OFH) were developed in TOPAS-nBio by considering an oxygen adduct state of DNA and creating a competition kinetic mechanism between oxygen and the radioprotective molecule WR-1065. We named these models 'simple' and 'detailed' due to the way they handle the hydrogen abstraction pathways. We used the simple model to obtain additional information for the •OH-DNA hydrogen abstraction pathway probability for the detailed model. These models were calibrated and compared with published experimental data of linear and supercoiling fractions obtained with R6K plasmids, suspended in dioxane as a hydroxyl scavenger, and irradiated with 137 Cs gamma-rays. The reaction rates for WR-1065 and O 2 with DNA were taken from experimental works. Single-Strand Breaks (SSBs) and DSBs as a function of the dose for a range of oxygen concentrations [O 2 ] (0.021%-21%) were obtained. Finally, the hypoxia reduction factor (HRF) was obtained from DSBs. Main Results. Validation results followed the trend of the experimental within 12% for the supercoiled and linear plasmid fractions for both models. The HRF agreed with measurements obtained with 137 Cs and 200-280 kVp x-ray within experimental uncertainties. However, the HRF at an oxygen concentration of 2.1% overestimated experimental results by a factor of 1.7 ± 0.1. Increasing the concentration of WR-1065 from 1 mM to 10-100 mM resulted in a HRF difference of 0.01, within the 8% statistical uncertainty between TOPAS-nBio and experimental data. This highlights the possibility of using these chemical models to recreate experimental HRF results. Significance. Results support the OFH as a leading cause of oxygen radio-sensitization effects given a competition between oxygen and chemical DNA repair molecules like WR-1065., (© 2024 Institute of Physics and Engineering in Medicine. All rights, including for text and data mining, AI training, and similar technologies, are reserved.)

Purpose: In radiotherapy (RT) for brain tumors, patient heterogeneity masks treatment effects, complicating the prediction and mitigation of radiation-induced brain necrosis. Therefore, understanding this heterogeneity is essential for improving outcome assessments and reducing toxicity., Experimental Design: We developed a clinically practical pipeline to clarify the relationship between dosimetric features and outcomes by identifying key variables. We processed data from a cohort of 130 patients treated with proton therapy for brain and head and neck tumors, utilizing an expert-augmented Bayesian network to understand variable interdependencies and assess structural dependencies. Critical evaluation involved a three-level grading system for each network connection and a Markov blanket analysis to identify variables directly impacting necrosis risk. Statistical assessments included log-likelihood ratio, integrated discrimination index, net reclassification index, and receiver operating characteristic (ROC)., Results: The analysis highlighted tumor location and proximity to critical structures such as white matter and ventricles as major determinants of necrosis risk. The majority of network connections were clinically supported, with quantitative measures confirming the significance of these variables in patient stratification (log-likelihood ratio = 12.17; P = 0.016; integrated discrimination index = 0.15; net reclassification index = 0.74). The ROC curve area was 0.66, emphasizing the discriminative value of nondosimetric variables., Conclusions: Key patient variables critical to understanding brain necrosis post-RT were identified, aiding the study of dosimetric impacts and providing treatment confounders and moderators. This pipeline aims to enhance outcome assessments by revealing at-risk patients, offering a versatile tool for broader applications in RT to improve treatment personalization in different disease sites., (©2024 American Association for Cancer Research.)

Radionuclides used for imaging and therapy can show high molecular specificity in the body with appropriate targeting ligands. We hypothesized that local energy delivered by molecularly targeted radionuclides could chemically activate prodrugs at disease sites while avoiding activation in off-target sites of toxicity. As proof-of-principle, we tested whether this strategy of " RA dionuclide i nduced D rug E ngagement for R elease" ( RAiDER ) could locally deliver combined radiation and chemotherapy to maximize tumor cytotoxicity while minimizing exposure to activated chemotherapy in off-target sites., Methods: We screened the ability of radionuclides to chemically activate a model radiation-activated prodrug consisting of the microtubule destabilizing monomethyl auristatin E caged by a radiation-responsive phenyl azide ("caged-MMAE") and interpreted experimental results using the radiobiology computational simulation suite TOPAS-nBio. RAiDER was evaluated in syngeneic mouse models of cancer using fibroblast activation protein inhibitor (FAPI) agents 99m Tc-FAPI-34 and 177 Lu-FAPI-04, the prostate-specific membrane antigen (PSMA) agent 177 Lu-PSMA-617, combined with caged-MMAE or caged-exatecan. Biodistribution in mice, combined with clinical dosimetry, estimated the relationship between radiopharmaceutical uptake in patients and anticipated concentrations of activated prodrug using RAiDER., Results: RAiDER efficiency varied by 250-fold across radionuclides ( 99m Tc> 177 Lu> 64 Cu> 68 Ga> 223 Ra> 18 F), yielding up to 1.22µM prodrug activation per Gy of exposure from 99m Tc. Computational simulations implicated low-energy electron-mediated free radical formation as driving prodrug activation. Clinically relevant radionuclide concentrations chemically activated caged-MMAE restored its ability to destabilize microtubules and increased its cytotoxicity by up to 600-fold compared to non-irradiated prodrug. Mice treated with 99m Tc-FAPI-34 and caged-MMAE accumulated up to 3000× greater concentrations of activated MMAE in tumors compared to other tissues. RAiDER with 99m Tc-FAPI-34 or 177 Lu-FAPI-04 delayed tumor growth, while monotherapies did not ( P <0.03). Clinically-guided dosimetry suggests sufficient radiation doses can be delivered to activate therapeutically meaningful levels of prodrug., Conclusion: This proof-of-concept study shows that RAiDER is compatible with multiple radionuclides commonly used in nuclear medicine and has the potential to improve the efficacy of radiopharmaceutical therapies to treat cancer safely. RAiDER thus shows promise as an effective strategy to treat disseminated malignancies and broadens the capability of radiopharmaceuticals to trigger diverse biological and therapeutic responses.

Objective . To allow the estimation of secondary cancer risks from radiation therapy treatment plans in a comprehensive and user-friendly Monte Carlo (MC) framework. Method . Patient planning computed tomography scans were extended superior-inferior using the International Commission on Radiological Protection's Publication 145 computational mesh phantoms and skeletal matching. Dose distributions were calculated with the TOPAS MC system using novel mesh capabilities and the digital imaging and communications in medicine radiotherapy extension interface. Finally, in-field and out-of-field cancer risk was calculated using both sarcoma and carcinoma risk models with two alternative parameter sets. Result . The TOPAS MC framework was extended to facilitate epidemiological studies on radiation-induced cancer risk. The framework is efficient and allows automated analysis of large datasets. Out-of-field organ dose was small compared to in-field dose, but the risk estimates indicate a non-negligible contribution to the total radiation induced cancer risk. Significance . This work equips the TOPAS MC system with anatomical extension, mesh geometry, and cancer risk model capabilities that make state-of-the-art out-of-field dose calculation and risk estimation accessible to a large pool of users. Furthermore, these capabilities will facilitate further refinement of risk models and sensitivity analysis of patient specific treatment options., (© 2024 Institute of Physics and Engineering in Medicine.)

Objective. To investigate different dosimetric aspects of 90 Y-IsoPet™ intratumoral therapy in canine soft tissue sarcomas, model the spatial spread of the gel post-injection, evaluate absorbed dose to clinical target volumes, and assess dose distributions and treatment efficacy. Approach. Six canine cases treated with 90 Y-IsoPet™ for soft tissue sarcoma at the Veterinary Health Center, University of Missouri are analyzed in this retrospective study. The dogs received intratumoral IsoPet™ injections, following a grid pattern to achieve a near-uniform dose distribution in the clinical target volume. Two dosimetry methods were performed retrospectively using the Monte Carlo toolkit OpenTOPAS: imaging-based dosimetry obtained from post-injection PET/CT scans, and stylized phantom-based dosimetry modeled from the planned injection points to the gross tumor volume. For the latter, a Gaussian parameter with variable sigma was introduced to reflect the spatial spread of IsoPet™. The two methods were compared using dose-volume histograms (DVHs) and dose homogeneity, allowing an approximation of the closest sigma for the spatial spread of the gel post-injection. In addition, we compared Monte Carlo-based dosimetry with voxel S-value (VSV)-based dosimetry to investigate the dosimetric differences. Main results. Imaging-based dosimetry showed differences between Monte Carlo and VSV calculations in tumor high-density areas with higher self-absorption. Stylized phantom-based dosimetry indicated a more homogeneous target dose with increasing sigma. The sigma approximation of the 90 Y-IsoPet™ post-injection gel spread resulted in a median sigma of approximately 0.44 mm across all cases to reproduce the dose heterogeneity observed in Monte Carlo calculations. Significance. The results indicate that dose modeling based on planned injection points can serve as a first-order approximation for the delivered dose in 90 Y-IsoPet™ therapy for canine soft tissue sarcomas. The dosimetry evaluation highlights the non-uniformity of absorbed doses despite the gel spread, emphasizing the importance of considering tumor dose heterogeneity in treatment evaluation. Our findings suggest that using Monte Carlo for dose calculation seems more suitable for this type of tumor where high-density areas might play an important role in dosimetry., (© 2024 Institute of Physics and Engineering in Medicine.)

Numerous dose rate effects have been described over the past 6-7 decades in the radiation biology and radiation oncology literature depending on the dose rate range being discussed. This review focuses on the impact and understanding of altering dose rates in the context of radiation therapy, but does not discuss dose rate effects as relevant to radiation protection. The review starts with a short historic review of early studies on dose rate effects, considers mechanisms thought to underlie dose rate dependencies, then discusses some current issues in clinical findings with altered dose rates, the importance of dose rate in brachytherapy, and the current timely topic of the use of very high dose rates, so-called FLASH radiotherapy. The discussion includes dose rate effects in vitro in cultured cells, in in vivo experimental systems and in the clinic, including both tumors and normal tissues. Gaps in understanding dose rate effects are identified, as are opportunities for improving clinical use of dose rate modulation., (© 2024 by Radiation Research Society. All rights of reproduction in any form reserved.)

Objective. To present a new set of lithium-ion cross-sections for (i) ionization and excitation processes down to 700 eV, and (ii) charge-exchange processes down to 1 keV u -1 . To evaluate the impact of the use of these cross-sections on micro a nano dosimetric quantities in the context of boron neutron capture (BNC) applications/techniques. Approach. The Classical Trajectory Monte Carlo method was used to calculate Li ion charge-exchange cross sections in the energy range of 1 keV u -1 to 10 MeV u -1 . Partial Li ion charge states ionization and excitation cross-sections were calculated using a detailed charge screening factor. The cross-sections were implemented in Geant4-DNA v10.07 and simulations and verified using TOPAS-nBio by calculating stopping power and continuous slowing down approximation (CSDA) range against data from ICRU and SRIM. Further microdosimetric and nanodosimetric calculations were performed to quantify differences against other simulation approaches for low energy Li ions. These calculations were: lineal energy spectra ( yf ( y ) and yd ( y )), frequency mean lineal energyyF-, dose mean lineal energyyD-and ionization cluster size distribution analysis. Microdosimetric calculations were compared against a previous MC study that neglected charge-exchange and excitation processes. Nanodosimetric results were compared against pure ionization scaled cross-sections calculations. Main results. Calculated stopping power differences between ICRU and Geant4-DNA decreased from 33.78% to 6.9%. The CSDA range difference decreased from 621% to 34% when compared against SRIM calculations. Geant4-DNA/TOPAS calculated dose mean lineal energy differed by 128% from the previous Monte Carlo. Ionization cluster size frequency distributions for Li ions differed by 76%-344.11% for 21 keV and 2 MeV respectively. With a decrease in the N 1 within 9% at 10 keV and agreeing after the 100 keV. With the new set of cross-sections being able to better simulate low energy behaviors of Li ions. Significance. This work shows an increase in detail gained from the use of a more complete set of low energy cross-sections which include charge exchange processes. Significant differences to previous simulation results were found at the microdosimetric and nanodosimetric scales that suggest that Li ions cause less ionizations per path length traveled but with more energy deposits. Microdosimetry results suggest that the BNC's contribution to cellular death may be mainly due to alpha particle production when boron-based drugs are distributed in the cellular membrane and beyond and by Li when it is at the cell cytoplasm regions., (© 2024 Institute of Physics and Engineering in Medicine.)

Radiotherapy is commonly used to treat cancer, and localized energy deposited by radiotherapy has the potential to chemically uncage prodrugs; however, it has been challenging to demonstrate prodrug activation that is both sustained in vivo and truly localized to tumors without affecting off-target tissues. To address this, we developed a series of novel phenyl-azide-caged, radiation-activated chemotherapy drug-conjugates alongside a computational framework for understanding corresponding pharmacokinetic and pharmacodynamic (PK/PD) behaviors. We especially focused on an albumin-bound prodrug of monomethyl auristatin E (MMAE) and found it blocked tumor growth in mice, delivered a 130-fold greater amount of activated drug to irradiated tumor versus unirradiated tissue, was 7.5-fold more efficient than a non albumin-bound prodrug, and showed no appreciable toxicity compared to free or cathepsin-activatable drugs. These data guided computational modeling of drug action, which indicated that extended pharmacokinetics can improve localized and cumulative drug activation, especially for payloads with low vascular permeability and diffusivity and particularly in patients receiving daily treatments of conventional radiotherapy for weeks. This work thus offers a quantitative PK/PD framework and proof-of-principle experimental demonstration of how extending prodrug circulation can improve its localized activity in vivo ., Competing Interests: The authors declare the following competing financial interest(s): M.A.M. has received unrelated support from Genentech/Roche and Pfizer and research funding from Ionis Pharmaceuticals. R.W. has consulted for Boston Scientific, ModeRNA, Earli, and Accure Health, none of whom contributed to or were involved in this research. Patents are pending and/or awarded with the authors and Massachusetts General Hospital., (© 2024 The Authors. Published by American Chemical Society.)

Clinical treatment of cancer commonly incorporates X-ray radiation therapy (XRT), and developing spatially precise radiation-activatable drug delivery strategies may improve XRT efficacy while limiting off-target toxicities associated with systemically administered drugs. Nevertheless, achieving this has been challenging thus far because strategies typically rely on radical species with short lifespans, and the inherent nature of hypoxic and acidic tumor microenvironments may encourage spatially heterogeneous effects. It is hypothesized that the challenge could be bypassed by using scintillating nanoparticles that emit light upon X-ray absorption, locally forming therapeutic drug depots in tumor tissues. Thus a nanoparticle platform (Scintillating nanoparticle Drug Depot; SciDD) that enables the local release of cytotoxic payloads only after activation by XRT is developed, thereby limiting off-target toxicity. As a proof-of-principle, SciDD is used to deliver a microtubule-destabilizing payload MMAE (monomethyl auristatin E). With as little as a 2 Gy local irradiation to tumors, MMAE payloads are released effectively to kill tumor cells. XRT-mediated drug release is demonstrated in multiple mouse cancer models and showed efficacy over XRT alone (p < 0.0001). This work shows that SciDD can act as a local drug depot with spatiotemporally controlled release of cancer therapeutics., (© 2024 Wiley‐VCH GmbH.)

Objective . We provide optimal particle split numbers for speeding up TOPAS Monte Carlo simulations of linear accelerator (linac) treatment heads while maintaining accuracy. In addition, we provide a new TOPAS physics module for simulating photoneutron production and transport. Approach. TOPAS simulation of a Siemens Oncor linac was used to determine the optimal number of splits for directional bremsstrahlung splitting as a function of the field size for 6 MV and 18 MV x-ray beams. The linac simulation was validated against published data of lateral dose profiles and percentage depth-dose curves (PDD) for the largest square field (40 cm side). In separate simulations, neutron particle split and the custom TOPAS physics module was used to generate and transport photoneutrons, called 'TsPhotoNeutron'. Verification of accuracy was performed by comparing simulations with published measurements of: (1) neutron yields as a function of beam energy for thick targets of Al, Cu, Ta, W, Pb and concrete; and (2) photoneutron energy spectrum at 40 cm laterally from the isocenter of the Oncor linac from an 18 MV beam with closed jaws and MLC. Main results. The optimal number of splits obtained for directional bremsstrahlung splitting enhanced the computational efficiency by two orders of magnitude. The efficiency decreased with increasing beam energy and field size. Calculated lateral profiles in the central region agreed within 1 mm/2% from measured data, PDD curves within 1 mm/1%. For the TOPAS physics module, at a split number of 146, the efficiency of computing photoneutron yields was enhanced by a factor of 27.6, whereas it improved the accuracy over existing Geant4 physics modules. Significance. This work provides simulation parameters and a new TOPAS physics module to improve the efficiency and accuracy of TOPAS simulations that involve photonuclear processes occurring in high- Z materials found in linac components, patient devices, and treatment rooms, as well as to explore new therapeutic modalities such as very-high energy electron therapy., (© 2024 Institute of Physics and Engineering in Medicine.)

Objective . Clinical outcomes after proton therapy have shown some variability that is not fully understood. Different approaches have been suggested to explain the biological outcome, but none has yet provided a comprehensive and satisfactory rationale for observed toxicities. The relatively recent transition from passive scattering (PS) to pencil beam scanning (PBS) treatments has significantly increased the voxel-wise dose rate in proton therapy. In addition, the dose rate distribution is no longer uniform along the cross section of the target but rather highly heterogeneous, following the spot placement. We suggest investigating dose rate as potential contributor to a more complex proton RBE model. Approach . Due to the time structure of the PBS beam delivery the instantaneous dose rate is highly variable voxel by voxel. Several possible parameters to represent voxel-wise dose rate for a given clinical PBS treatment plan are detailed. These quantities were implemented in the scripting environment of our treatment planning system, and computations experimentally verified. Sample applications to treated patient plans are shown. Main results . Computed dose rates we experimentally confirmed. Dose rate maps vary depending on which method is used to represent them. Mainly, the underlying time and dose intervals chosen determine the topography of the resultant distributions. The maximum dose rates experienced by any target voxel in a given PBS treatment plan in our system range from ∼100 to ∼450 Gy(RBE)/min, a factor of 10-100 increase compared to PS. These dose rate distributions are very heterogeneous, with distinct hot spots. Significance . Voxel-wise dose rates for current clinical PBS treatment plans vary greatly from clinically established practice with PS. The exploration of different dose rate measures to evaluate potential correlations with observed clinical outcomes is suggested, potentially adding a missing component in the understanding of proton relative biological effectiveness (RBE)., (© 2024 Institute of Physics and Engineering in Medicine.)

Background: Ultra-high dose rate (FLASH) radiation has been reported to efficiently suppress tumor growth while sparing normal tissue; however, the mechanism of the differential tissue sparing effect is still not known. Oxygen has long been known to profoundly impact radiobiological responses, and radiolytic oxygen depletion has been considered to be a possible cause or contributor to the FLASH phenomenon., Purpose: This work investigates the impact of tissue pO 2 profiles, oxygen depletion per unit dose (g), and the oxygen concentration yielding half-maximum radiosensitization (the average of its maximum value and one) (k) in tumor and normal tissue., Methods: We developed a model that considers the dependent relationship between oxygen depletion and change of radiosensitivity by FLASH irradiation. The model assumed that FLASH irradiation depletes intracellular oxygen more rapidly than it diffuses into the cell from the extracellular environment. Cell survival was calculated based on the linear quadratic-linear model and the radiosensitivity related parameters were adjusted in 1 Gy increments of the administered dose. The model reproduced published experimental data that were obtained with different cell lines and oxygen concentrations, and was used to analyze the impact of parameter uncertainties on the radiobiological responses. This study expands the oxygen depletion analysis of FLASH to normal human tissue and tumor based on clinically determined aggregate and individual patient pO 2 profiles., Results: The results show that the pO 2 profile is the most essential factor that affects biological response and analyses based on the median pO 2 rather than the full pO 2 profile can be unreliable and misleading. Additionally, the presence of a small fraction of cells on the threshold of radiobiologic hypoxia substantially alters biological response due to FLASH oxygen depletion. We found that an increment in the k value is generally more protective of tumor than normal tissue due to a higher frequency of lower pO 2 values in tumors. Variation in the g value affects the dose at which oxygen depletion impacts response, but does not alter the dose-dependent response trends, if the g value is identical in both tumor and normal tissue., Conclusions: The therapeutic efficacy of FLASH oxygen depletion is likely patient and tissue-dependent. For breast cancer, FLASH is beneficial in a minority of cases; however, in a subset of well oxygenated tumors, a therapeutic gain may be realized due to induced normal tissue hypoxia., (© 2023 American Association of Physicists in Medicine.)

The induction and repair of DNA double-strand breaks (DSBs) are critical factors in the treatment of cancer by radiotherapy. To investigate the relationship between incident radiation and cell death through DSB induction many in silico models have been developed. These models produce and use custom formats of data, specific to the investigative aims of the researchers, and often focus on particular pairings of damage and repair models. In this work we use a standard format for reporting DNA damage to evaluate combinations of different, independently developed, models. We demonstrate the capacity of such inter-comparison to determine the sensitivity of models to both known and implicit assumptions. Specifically, we report on the impact of differences in assumptions regarding patterns of DNA damage induction on predicted initial DSB yield, and the subsequent effects this has on derived DNA repair models. The observed differences highlight the importance of considering initial DNA damage on the scale of nanometres rather than micrometres. We show that the differences in DNA damage models result in subsequent repair models assuming significantly different rates of random DSB end diffusion to compensate. This in turn leads to disagreement on the mechanisms responsible for different biological endpoints, particularly when different damage and repair models are combined, demonstrating the importance of inter-model comparisons to explore underlying model assumptions., (©2023 by Radiation Research Society. All rights of reproduction in any form reserved.)

Ovarian cancer (OC) is the most lethal gynecologic malignancy (5-y overall survival rate, 46%). OC is generally detected when it has already spread to the peritoneal cavity (peritoneal carcinomatosis). This study investigated whether gadolinium-based nanoparticles (Gd-NPs) increase the efficacy of targeted radionuclide therapy using [ 177 Lu]Lu-DOTA-trastuzumab (an antibody against human epidermal growth factor receptor 2). Gd-NPs have radiosensitizing effects in conventional external-beam radiotherapy and have been tested in clinical phase II trials. Methods: First, the optimal activity of [ 177 Lu]Lu-DOTA-trastuzumab (10, 5, or 2.5 MBq) combined or not with 10 mg of Gd-NPs (single injection) was investigated in athymic mice bearing intraperitoneal OC cell (human epidermal growth factor receptor 2-positive) tumor xenografts. Next, the therapeutic efficacy and toxicity of 5 MBq of [ 177 Lu]Lu-DOTA-trastuzumab with Gd-NPs (3 administration regimens) were evaluated. NaCl, trastuzumab plus Gd-NPs, and [ 177 Lu]Lu-DOTA-trastuzumab alone were used as controls. Biodistribution and dosimetry were determined, and Monte Carlo simulation of energy deposits was performed. Lastly, Gd-NPs' subcellular localization and uptake, and the cytotoxic effects of the combination, were investigated in 3 cancer cell lines to obtain insights into the involved mechanisms. Results: The optimal [ 177 Lu]Lu-DOTA-trastuzumab activity when combined with Gd-NPs was 5 MBq. Moreover, compared with [ 177 Lu]Lu-DOTA-trastuzumab alone, the strongest therapeutic efficacy (tumor mass reduction) was obtained with 2 injections of 5 mg of Gd-NPs/d (separated by 6 h) at 24 and 72 h after injection of 5 MBq of [ 177 Lu]Lu-DOTA-trastuzumab. In vitro experiments showed that Gd-NPs colocalized with lysosomes and that their radiosensitizing effect was mediated by oxidative stress and inhibited by deferiprone, an iron chelator. Exposure of Gd-NPs to 177 Lu increased the Auger electron yield but not the absorbed dose. Conclusion: Targeted radionuclide therapy can be combined with Gd-NPs to increase the therapeutic effect and reduce the injected activities. As Gd-NPs are already used in the clinic, this combination could be a new therapeutic approach for patients with ovarian peritoneal carcinomatosis., (© 2023 by the Society of Nuclear Medicine and Molecular Imaging.)

Background and Objective: MHV370, a dual antagonist of human Toll-like receptors (TLR) 7 and 8, suppresses cytokines and interferon-stimulated genes in vitro and in vivo, and  has demonstrated efficacy in murine models of lupus. This first-in-human study aimed to evaluate the safety, tolerability, pharmacokinetics and pharmacodynamics of single and multiple doses of MHV370 in healthy adults, as well as the effects of food consumption on a single dose of MHV370., Methods: This was a phase 1, randomised, placebo-controlled study conducted in three parts. In part A, participants received (3:1) a single ascending dose (SAD) of 1, 3, 10, 20, 40, 80, 160, 320, 640 and 1000 mg MHV370 or placebo. In part B, participants received (3:1) multiple ascending doses (MAD) of 25, 50, 100, 200 and 400 mg MHV370 twice daily (b.i.d) or placebo for 14 days. In part C, participants received an open-label single dose of 200 mg MHV370 under fasted or fed conditions. Safety, pharmacokinetic and pharmacodynamic parameters were evaluated., Results: MHV370 was well tolerated, and no safety signal was observed in the study. No dose-limiting adverse events occurred across the dose range evaluated. Plasma concentrations of MHV370 increased with dose (mean [SD] maximum plasma concentrations ranged from 0.97 [0.48] to 1670 [861.0] ng/mL for SAD of 3-1000 mg, 29.5 [7.98] to 759 [325.0] ng/mL for MAD of 25-400 mg b.i.d. on day 1). The intake of food did not have a relevant impact on the pharmacokinetics of MHV370. Pharmacodynamic data indicated time- and dose-dependent inhibition of TLR7-mediated CD69 expression on B cells (100% inhibition at 24 h post-dose starting from SAD 160 mg and MAD 50 mg b.i.d.) and TLR8-mediated TNF release after ex vivo stimulation (>90% inhibition at 24 h post-dose starting from SAD 320 mg and MAD 100 mg b.i.d.)., Conclusion: The safety, pharmacokinetic and pharmacodynamic data support the further development of MHV370 in systemic autoimmune diseases driven by the overactivation of TLR7 and TLR8., (© 2023. The Author(s).)

FLASH radiation therapy (FLASH-RT), delivered with ultrahigh dose rate (UHDR), may allow patients to be treated with less normal tissue toxicity for a given tumor dose compared with currently used conventional dose rate. Clinical trials are being carried out and are needed to test whether this improved therapeutic ratio can be achieved clinically. During the clinical trials, quality assurance and credentialing of equipment and participating sites, particularly pertaining to UHDR-specific aspects, will be crucial for the validity of the outcomes of such trials. This report represents an initial framework proposed by the NRG Oncology Center for Innovation in Radiation Oncology FLASH working group on quality assurance of potential UHDR clinical trials and reviews current technology gaps to overcome. An important but separate consideration is the appropriate design of trials to most effectively answer clinical and scientific questions about FLASH. This paper begins with an overview of UHDR RT delivery methods. UHDR beam delivery parameters are then covered, with a focus on electron and proton modalities. The definition and control of safe UHDR beam delivery and current and needed dosimetry technologies are reviewed and discussed. System and site credentialing for large, multi-institution trials are reviewed. Quality assurance is then discussed, and new requirements are presented for treatment system standard analysis, patient positioning, and treatment planning. The tables and figures in this paper are meant to serve as reference points as we move toward FLASH-RT clinical trial performance. Some major questions regarding FLASH-RT are discussed, and next steps in this field are proposed. FLASH-RT has potential but is associated with significant risks and complexities. We need to redefine optimization to focus not only on the dose but also on the dose rate in a manner that is robust and understandable and that can be prescribed, validated, and confirmed in real time. Robust patient safety systems and access to treatment data will be critical as FLASH-RT moves into the clinical trials., (Copyright © 2023 The Authors. Published by Elsevier Inc. All rights reserved.)

Purpose: Radiation-induced lymphopenia has gained attention recently as the result of its correlation with survival in a range of indications, particularly when combining radiation therapy (RT) with immunotherapy. The purpose of this study is to use a dynamic blood circulation model combined with observed lymphocyte depletion in patients to derive the in vivo radiosensitivity of circulating lymphocytes and study the effect of RT delivery parameters., Methods and Materials: We assembled a cohort of 17 patients with hepatocellular carcinoma treated with proton RT alone in 15 fractions (fx) using conventional dose rates (beam-on time [BOT], 120 seconds) for whom weekly absolute lymphocyte counts (ALCs) during RT and follow-up were available. We used HEDOS, a time-dependent, whole-body, blood flow computational framework, in combination with explicit liver blood flow modeling, to calculate the dose volume histograms for circulating lymphocytes for changing BOTs (1 second-300 seconds) and fractionations (5 fx, 15 fx). From this, we used the linear cell survival model and an exponential model to determine patient-specific lymphocyte radiation sensitivity, α, and recovery, σ, respectively., Results: The in vivo-derived patient-specific α had a median of 0.65 Gy - 1 (range, 0.30-1.38). Decreasing BOT to 1 second led to an increased average end-of-treatment ALC of 27.5%, increasing to 60.3% when combined with the 5-fx regimen. Decreasing to 5 fx at the conventional dose rate led to an increase of 17.0% on average. The benefit of both increasing dose rate and reducing the number of fractions was patient specificࣧpatients with highly sensitive lymphocytes benefited most from decreasing BOT, whereas patients with slow lymphocyte recovery benefited most from the shorter fractionation regimen., Conclusions: We observed that increasing dose rate at the same fractionation reduced ALC depletion more significantly than reducing the number of fractions. High-dose-rates led to an increased sparing of lymphocytes when shortening the fractionation regimen, particularly for patients with radiosensitive lymphocytes at elevated risk., (Copyright © 2023 Elsevier Inc. All rights reserved.)

Objective . To commission a proton, double-scattering FLASH beamline by maximizing efficiency and field size, enabling higher-linear energy transfer FLASH radiotherapy to cells and small animals using a spread-out Bragg peak (SOBP) treatment configuration. We further aim to provide a configuration guide for the design of future FLASH proton double-scattering (DS) beamlines. Approach . Beam spot size and spread were measured with film and implemented into TOol for PArticle Simulation (TOPAS). Monte Carlo simulations were optimized to verify the ideal positioning, dimensions, and material of scattering foils, secondary scatterers, ridge filters, range compensators, and apertures. A ridge filter with three discrete heights was used to create a spread-out Bragg peak (SOBP) and was experimentally verified using our in-house experimental FLASH beamline. The increase in dose rate was compared to nominal shoot-through techniques. Results . The configuration and scatterer distance producing the largest field size of acceptable flatness, without drastically compromising dose rate was determined to be an elliptical field of 2 cm × 1.5 cm (25% larger than a previous configuration). SOBP testing yielded three distinct but connected spikes in dose with flatness under 5%. Reducing the thickness of the (first) scattering foil by a factor of two was found to increase efficiency by 50%. The new settings increased the field size, provided a Bragg peak treatment option, and increased the maximum available dose rate by 85%, as compared to the previous, shoot through method. Significance . Beam line updates established FLASH dose rates of over 135 Gy s -1 (potentially higher) at our double-scattering beamline, increased the efficiency and field size, and enabled SOBP treatments by incorporating an optimized ridge filter. Based on our simulations we provide parametric suggestions when commissioning a new proton DS beamline. This enhanced FLASH beamline for SOBP irradiations with higher dose rates and larger field sizes will enable a wider variety of experimentation in future studies., (© 2023 Institute of Physics and Engineering in Medicine.)

Introduction: DNA damage is the main predictor of response to radiation therapy for cancer. Its Q8 quantification and characterization are paramount for treatment optimization, particularly in advanced modalities such as proton and alpha-targeted therapy., Methods: We present a novel approach called the Microdosimetric Gamma Model (MGM) to address this important issue. The MGM uses the theory of microdosimetry, specifically the mean energy imparted to small sites, as a predictor of DNA damage properties. MGM provides the number of DNA damage sites and their complexity, which were determined using Monte Carlo simulations with the TOPAS-nBio toolkit for monoenergetic protons and alpha particles. Complexity was used together with a illustrative and simplistic repair model to depict the differences between high and low LET radiations., Results: DNA damage complexity distributions were were found to follow a Gamma distribution for all monoenergetic particles studied. The MGM functions allowed to predict number of DNA damage sites and their complexity for particles not simulated with microdosimetric measurements (yF) in the range of those studied., Discussion: Compared to current methods, MGM allows for the characterization of DNA damage induced by beams composed of multi-energy components distributed over any time configuration and spatial distribution. The output can be plugged into ad hoc repair models that can predict cell killing, protein recruitment at repair sites, chromosome aberrations, and other biological effects, as opposed to current models solely focusing on cell survival. These features are particularly important in targeted alpha-therapy, for which biological effects remain largely uncertain. The MGM provides a flexible framework to study the energy, time, and spatial aspects of ionizing radiation and offers an excellent tool for studying and optimizing the biological effects of these radiotherapy modalities., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Bertolet, Chamseddine, Paganetti and Schuemann.)

Objective . The TOPAS-nBio Monte Carlo track structure simulation code, a wrapper of Geant4-DNA, was extended for its use in pulsed and longtime homogeneous chemistry simulations using the Gillespie algorithm. Approach . Three different tests were used to assess the reliability of the implementation and its ability to accurately reproduce published experimental results: (1) a simple model with a known analytical solution, (2) the temporal evolution of chemical yields during the homogeneous chemistry stage, and (3) radiolysis simulations conducted in pure water with dissolved oxygen at concentrations ranging from 10 μ M to 1 mM with [H 2 O 2 ] yields calculated for 100 MeV protons at conventional and FLASH dose rates of 0.286 Gy s -1 and 500 Gy s -1 , respectively. Simulated chemical yield results were compared closely with data calculated using the Kinetiscope software which also employs the Gillespie algorithm. Main results . Validation results in the third test agreed with experimental data of similar dose rates and oxygen concentrations within one standard deviation, with a maximum of 1% difference for both conventional and FLASH dose rates. In conclusion, the new implementation of TOPAS-nBio for the homogeneous long time chemistry simulation was capable of recreating the chemical evolution of the reactive intermediates that follow water radiolysis. Significance . Thus, TOPAS-nBio provides a reliable all-in-one chemistry simulation of the physical, physico-chemical, non-homogeneous, and homogeneous chemistry and could be of use for the study of FLASH dose rate effects on radiation chemistry., (© 2023 Institute of Physics and Engineering in Medicine.)

Ultra-high dose rate irradiation has been reported to protect normal tissues more than conventional dose rate irradiation. This tissue sparing has been termed the FLASH effect. We investigated the FLASH effect of proton irradiation on the intestine as well as the hypothesis that lymphocyte depletion is a cause of the FLASH effect. A 16 × 12 mm 2 elliptical field with a dose rate of ~120 Gy/s was provided by a 228 MeV proton pencil beam. Partial abdominal irradiation was delivered to C57BL/6j and immunodeficient Rag1 -/- /C57 mice. Proliferating crypt cells were counted at 2 days post exposure, and the thickness of the muscularis externa was measured at 280 days following irradiation. FLASH irradiation did not reduce the morbidity or mortality of conventional irradiation in either strain of mice; in fact, a tendency for worse survival in FLASH-irradiated mice was observed. There were no significant differences in lymphocyte numbers between FLASH and conventional-dose-rate mice. A similar number of proliferating crypt cells and a similar thickness of the muscularis externa following FLASH and conventional dose rate irradiation were observed. Partial abdominal FLASH proton irradiation at 120 Gy/s did not spare normal intestinal tissue, and no difference in lymphocyte depletion was observed. This study suggests that the effect of FLASH irradiation may depend on multiple factors, and in some cases dose rates of over 100 Gy/s do not induce a FLASH effect and can even result in worse outcomes.

Objective . The TOol for PArticle Simulation (TOPAS) is a Geant4-based Monte Carlo software application that has been used for both research and clinical studies in medical physics. So far, most users of TOPAS have focused on radiotherapy-related studies, such as modeling radiation therapy delivery systems or patient dose calculation. Here, we present the first set of TOPAS extensions to make it easier for TOPAS users to model medical imaging systems. Approach . We used the extension system of TOPAS to implement pre-built, user-configurable geometry components such as detectors (e.g. flat-panel and multi-planar detectors) for various imaging modalities and pre-built, user-configurable scorers for medical imaging systems (e.g. digitizer chain). Main results . We developed a flexible set of extensions that can be adapted to solve research questions for a variety of imaging modalities. We then utilized these extensions to model specific examples of cone-beam CT (CBCT), positron emission tomography (PET), and prompt gamma (PG) systems. The first of these new geometry components, the FlatImager, was used to model example CBCT and PG systems. Detected signals were accumulated in each detector pixel to obtain the intensity of x-rays penetrating objects or prompt gammas from proton-nuclear interaction. The second of these new geometry components, the RingImager, was used to model an example PET system. Positron-electron annihilation signals were recorded in crystals of the RingImager and coincidences were detected. The simulated data were processed using corresponding post-processing algorithms for each modality and obtained results in good agreement with the expected true signals or experimental measurement. Significance . The newly developed extension is a first step to making it easier for TOPAS users to build and simulate medical imaging systems. Together with existing TOPAS tools, this extension can help integrate medical imaging systems with radiotherapy simulations for image-guided radiotherapy., (© 2023 Institute of Physics and Engineering in Medicine.)

Objective . Irradiation at FLASH dose rates (>40 Gy s -1 ) has received great attention due to its reported normal tissue sparing effect. The FLASH effect was originally observed in electron irradiations but has since been shown to also occur with both photon and proton beams. Several mechanisms have been proposed to explain the tissue sparing at high dose rates, including effects involving oxygen, such as depletion of oxygen within the irradiated cells. In this study, we investigated the protective role of FLASH proton irradiation on the skin when varying the oxygen concentration. Approach . Our double scattering proton system provided a 1.2 × 1.6 cm 2 elliptical field at a dose rate of ∼130 Gy s -1 . The conventional dose rate was ∼0.4 Gy s -1 . The legs of the FVB/N mice were marked with two tattooed dots and fixed in a holder for exposure. To alter the skin oxygen concentration, the mice were breathing pure oxygen or had their legs tied to restrict blood flow. The distance between the two dots was measured to analyze skin contraction over time. Main results . FLASH irradiation mitigated skin contraction by 15% compared to conventional dose rate irradiation. The epidermis thickness and collagen deposition at 75 d following 25 to 30 Gy exposure suggested a long-term protective function in the skin from FLASH irradiation. Providing the mice with oxygen or reducing the skin oxygen concentration removed the dose-rate-dependent difference in response. Significanc e. FLASH proton irradiation decreased skin contraction, epidermis thickness and collagen deposition compared to standard dose rate irradiations. The observed oxygen-dependence of the FLASH effect is consistent with, but not conclusive of, fast oxygen depletion during the exposure., (© 2023 Institute of Physics and Engineering in Medicine.)

Purpose: Track structure Monte Carlo (MC) codes have achieved successful outcomes in the quantitative investigation of radiation-induced initial DNA damage. The aim of the present study is to extend a Geant4-DNA radiobiological application by incorporating a feature allowing for the prediction of DNA rejoining kinetics and corresponding cell surviving fraction along time after irradiation, for a Chinese hamster V79 cell line, which is one of the most popular and widely investigated cell lines in radiobiology., Methods: We implemented the Two-Lesion Kinetics (TLK) model, originally proposed by Stewart, which allows for simulations to calculate residual DNA damage and surviving fraction along time via the number of initial DNA damage and its complexity as inputs., Results: By optimizing the model parameters of the TLK model in accordance to the experimental data on V79, we were able to predict both DNA rejoining kinetics at low linear energy transfers (LET) and cell surviving fraction., Conclusion: This is the first study to demonstrate the implementation of both the cell surviving fraction and the DNA rejoining kinetics with the estimated initial DNA damage, in a realistic cell geometrical model simulated by full track structure MC simulations at DNA level and for various LET. These simulation and model make the link between mechanistic physical/chemical damage processes and these two specific biological endpoints., (Copyright © 2022 Associazione Italiana di Fisica Medica e Sanitaria. Published by Elsevier Ltd. All rights reserved.)

The effects of realistic, deep space radiation environments on neuronal function remain largely unexplored. In silico modeling studies of radiation-induced neuronal damage provide important quantitative information about physico-chemical processes that are not directly accessible through radiobiological experiments. Here, we present the first nano-scale computational analysis of broad-spectrum galactic cosmic ray irradiation in a realistic neuron geometry. We constructed thousands of in silico realizations of a CA1 pyramidal neuron, each with over 3500 stochastically generated dendritic spines. We simulated the entire 33 ion-energy beam spectrum currently in use at the NASA Space Radiation Laboratory galactic cosmic ray simulator (GCRSim) using the TOol for PArticle Simulation (TOPAS) and TOPAS-nBio Monte Carlo-based track structure simulation toolkits. We then assessed the resulting nano-scale dosimetry, physics processes, and fluence patterns. Additional comparisons were made to a simplified 6 ion-energy spectrum (SimGCRSim) also used in NASA experiments. For a neuronal absorbed dose of 0.5 Gy GCRSim, we report an average of 250 ± 10 ionizations per micrometer of dendritic length, and an additional 50 ± 10, 7 ± 2, and 4 ± 2 ionizations per mushroom, thin, and stubby spine, respectively. We show that neuronal energy deposition by proton andα-particle tracks declines approximately hyperbolically with increasing primary particle energy at mission-relevant energies. We demonstrate an inverted exponential relationship between dendritic segment irradiation probability and neuronal absorbed dose for each ion-energy beam. We also find that there are no significant differences in the average physical responses between the GCRSim and SimGCRSim spectra. To our knowledge, this is the first nano-scale simulation study of a realistic neuron geometry using the GCRSim and SimGCRSim spectra. These results may be used as inputs to theoretical models, aid in the interpretation of experimental results, and help guide future study designs., (© 2022 Institute of Physics and Engineering in Medicine.)

With increasing focus on the translation of the observed FLASH effect into clinical practice, this paper presents treatment planning considerations for its development using proton therapy. Potential requirements to induce a FLASH effect are discussed along with the properties of existing proton therapy delivery systems and the changes in planning and delivery approaches required to satisfy these prerequisites. For the exploration of treatment planning approaches for FLASH, developments in treatment planning systems are needed. Flexibility in adapting to new information will be important in such an evolving area. Variations in definitions, threshold values and assumptions can make it difficult to compare different published studies and to interpret previous studies in the context of new information. Together with the fact that much is left to be understood about the underlying mechanism behind the FLASH effect, a systematic and comprehensive approach to information storage is encouraged. Collecting and retaining more detailed information on planned and realised dose delivery as well as reporting the assumptions made in planning studies creates the potential for research to be revisited and re-evaluated in the light of future improvements in understanding. Forward thinking at the time of study development can help facilitate retrospective analysis. This, we hope, will increase the available evidence and accelerate the translation of the FLASH effect into clinical benefit., (Copyright © 2022 The Author(s). Published by Elsevier B.V. All rights reserved.)

Track structure Monte Carlo simulations are a useful tool to investigate the damage induced to DNA by ionizing radiation. These simulations usually rely on simplified geometrical representations of the DNA subcomponents. DNA damage is determined by the physical and physicochemical processes occurring within these volumes. In particular, damage to the DNA backbone is generally assumed to result in strand breaks. DNA damage can be categorized as direct (ionization of an atom part of the DNA molecule) or indirect (damage from reactive chemical species following water radiolysis). We also consider quasi-direct effects, i.e., damage originated by charge transfers after ionization of the hydration shell surrounding the DNA. DNA geometries are needed to account for the damage induced by ionizing radiation, and different geometry models can be used for speed or accuracy reasons. In this work, we use the Monte Carlo track structure tool TOPAS-nBio, built on top of Geant4-DNA, for simulation at the nanometer scale to evaluate differences among three DNA geometrical models in an entire cell nucleus, including a sphere/spheroid model specifically designed for this work. In addition to strand breaks, we explicitly consider the direct, quasi-direct, and indirect damage induced to DNA base moieties. We use results from the literature to determine the best values for the relevant parameters. For example, the proportion of hydroxyl radical reactions between base moieties was 80%, and between backbone, moieties was 20%, the proportion of radical attacks leading to a strand break was 11%, and the expected ratio of base damages and strand breaks was 2.5-3. Our results show that failure to update parameters for new geometric models can lead to significant differences in predicted damage yields., (©2022 by Radiation Research Society. All rights of reproduction in any form reserved.)

Objective. Monte Carlo (MC) codes are increasingly used for accurate radiotherapy dose calculation. In proton therapy, the accuracy of the dose calculation algorithm is expected to have a more significant impact than in photon therapy due to the depth-dose characteristics of proton beams. However, MC simulations come at a considerable computational cost to achieve statistically sufficient accuracy. There have been efforts to improve computational efficiency while maintaining sufficient accuracy. Among those, parallelizing particle transportation using graphic processing units (GPU) achieved significant improvements. Contrary to the central processing unit, a GPU has limited memory capacity and is not expandable. It is therefore challenging to score quantities with large dimensions requiring extensive memory. The objective of this study is to develop an open-source GPU-based MC package capable of scoring those quantities. Approach. We employed a hash-table, one of the key-value pair data structures, to efficiently utilize the limited memory of the GPU and score the quantities requiring a large amount of memory. With the hash table, only voxels interacting with particles will occupy memory, and we can search the data efficiently to determine their address. The hash-table was integrated with a novel GPU-based MC code, moqui. Main results. The developed code was validated against an MC code widely used in proton therapy, TOPAS, with homogeneous and heterogeneous phantoms. We also compared the dose calculation results of clinical treatment plans. The developed code agreed with TOPAS within 2%, except for the fall-off and regions, and the gamma pass rates of the results were >99% for all cases with a 2 mm/2% criteria. Significance. We can score dose-influence matrix and dose-rate on a GPU for a 3-field H&N case with 10 GB of memory using moqui, which would require more than 100 GB of memory with the conventionally used array data structure., (© 2022 Institute of Physics and Engineering in Medicine.)

Immunofluorescent tagging of DNA double-strand break (DSB) markers, such as γ-H2AX and other DSB repair proteins, are powerful tools in understanding biological consequences following irradiation. However, whilst the technique is widespread, there are many uncertainties related to its ability to resolve and reliably deduce the number of foci when counting using microscopy. We present a new tool for simulating radiation-induced foci in order to evaluate microscope performance within in silico immunofluorescent images. Simulations of the DSB distributions were generated using Monte Carlo track-structure simulation. For each DSB distribution, a corresponding DNA repair process was modelled and the un-repaired DSBs were recorded at several time points. Corresponding microscopy images for both a DSB and (γ-H2AX) fluorescent marker were generated and compared for different microscopes, radiation types and doses. Statistically significant differences in miscounting were found across most of the tested scenarios. These inconsistencies were propagated through to repair kinetics where there was a perceived change between radiation-types. These changes did not reflect the underlying repair rate and were caused by inconsistencies in foci counting. We conclude that these underlying uncertainties must be considered when analysing images of DNA damage markers to ensure differences observed are real and are not caused by non-systematic miscounting., (© 2022. The Author(s).)

Current Monte Carlo simulations of DNA damage have been reported only at ambient temperature. The aim of this work is to use TOPAS-nBio to simulate the yields of DNA single-strand breaks (SSBs) and double-strand breaks (DSBs) produced in plasmids under low-LET irradiation incorporating the effect of the temperature changes in the environment. A new feature was implemented in TOPAS-nBio to incorporate reaction rates used in the simulation of the chemical stage of water radiolysis as a function of temperature. The implemented feature was verified by simulating temperature-dependent G -values of chemical species in liquid water from 20 °C to 90 °C. For radiobiology applications, temperature dependent SSB and DSB yields were calculated from 0 °C to 42 °C, the range of available published measured data. For that, supercoiled DNA plasmids dissolved in aerated solutions containing EDTA irradiated by Cobalt-60 gamma-rays were simulated. TOPAS-nBio well reproduced published temperature-dependent G -values in liquid water and the yields of SSB and DSB for the temperature range considered. For strand break simulations, the model shows that the yield of SSB and DSB increased linearly with the temperature at a rate of (2.94 ± 0.17) × 10 -10 Gy -1 Da -1 °C -1 ( R 2  = 0.99) and (0.13 ± 0.01) × 10 -10 Gy -1 Da -1 °C -1 ( R 2  = 0.99), respectively. The extended capability of TOPAS-nBio is a complementary tool to simulate realistic conditions for a large range of environmental temperatures, allowing refined investigations of the biological effects of radiation., (© 2022 Institute of Physics and Engineering in Medicine.)

Research efforts in FLASH radiotherapy have increased at an accelerated pace recently. FLASH radiotherapy involves ultra-high dose rates and has shown to reduce toxicity to normal tissue while maintaining tumor response in pre-clinical studies when compared to conventional dose rate radiotherapy. The goal of this review is to summarize the studies performed to-date with proton, electron, and heavy ion FLASH radiotherapy, with particular emphasis on the physical aspects of each study and the advantages and disadvantages of each modality. Beam delivery parameters, experimental set-up, and the dosimetry tools used are described for each FLASH modality. In addition, modeling efforts and treatment planning for FLASH radiotherapy is discussed along with potential drawbacks when translated into the clinical setting. The final section concludes with further questions that have yet to be answered before safe clinical implementation of FLASH radiotherapy.

Objective . To evaluate the pre-treatment and post-treatment imaging-based dosimetry of patients treated with 90Y-microspheres, including accurate estimations of dose to tumor, healthy liver and lung. To do so, the Monte Carlo (MC) TOPAS platform is in this work extended towards its utilization in radionuclide therapy. Approach . Five patients treated at the Massachusetts General Hospital were selected for this study. All patients had data for both pre-treatment SPECT-CT imaging using 99mTc-MAA as a surrogate of the 90Y-microspheres treatment and SPECT-CT imaging immediately after the 90Y activity administration. Pre- and post-treatment doses were computed with TOPAS using the SPECT images to localize the source positions and the CT images to account for tissue inhomoegeneities. We compared our results with analytical calculations following the voxel-based MIRD scheme. Main results . TOPAS results largely agreed with the MIRD-based calculations in soft tissue regions: the average difference in mean dose to the liver was 0.14 Gy GBq -1 (2.6%). However, dose distributions in the lung differed considerably: absolute differences in mean doses to the lung ranged from 1.2 to 6.3 Gy GBq -1 and relative differences from 153% to 231%. We also found large differences in the intra-hepatic dose distributions between pre- and post-treatment imaging, but only limited differences in the pulmonary dose. Significance . Doses to lung were found to be higher using TOPAS with respect to analytical calculations which may significantly underestimate dose to the lung, suggesting the use of MC methods for 90Y dosimetry. According to our results, pre-treatment imaging may still be representative of dose to lung in these treatments., (© 2021 Institute of Physics and Engineering in Medicine.)

The chemical stage of the Monte Carlo track-structure (MCTS) code Geant4-DNA was extended for its use in DNA strand break (SB) simulations and compared against published experimental data. Geant4-DNA simulations were performed using pUC19 plasmids (2686 base pairs) in a buffered solution of DMSO irradiated by 60 Co or 137 Cs γ -rays. A comprehensive evaluation of SSB yields was performed considering DMSO, DNA concentration, dose and plasmid supercoiling. The latter was measured using the super helix density value used in a Brownian dynamics plasmid generation algorithm. The Geant4-DNA implementation of the independent reaction times method (IRT), developed to simulate the reaction kinetics of radiochemical species, allowed to score the fraction of supercoiled, relaxed and linearized plasmid fractions as a function of the absorbed dose. The percentage of the number of SB after •OH + DNA and H• + DNA reactions, referred as SSB efficiency, obtained using MCTS were 13.77% and 0.74% respectively. This is in reasonable agreement with published values of 12% and 0.8%. The SSB yields as a function of DMSO concentration, DNA concentration and super helix density recreated the expected published experimental behaviors within 5%, one standard deviation. The dose response of SSB and DSB yields agreed with published measurements within 5%, one standard deviation. We demonstrated that the developed extension of IRT in Geant4-DNA, facilitated the reproduction of experimental conditions. Furthermore, its calculations were strongly in agreement with experimental data. These two facts will facilitate the use of this extension in future radiobiological applications, aiding the study of DNA damage mechanisms with a high level of detail., (© 2021 Institute of Physics and Engineering in Medicine.)