1
|
Mietelska M, Pietrzak M, Bancer A, Ruciński A, Szefliński Z, Brzozowska B. Ionization Detail Parameters for DNA Damage Evaluation in Charged Particle Radiotherapy: Simulation Study Based on Cell Survival Database. Int J Mol Sci 2024; 25:5094. [PMID: 38791135 PMCID: PMC11121214 DOI: 10.3390/ijms25105094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2024] [Revised: 04/28/2024] [Accepted: 05/03/2024] [Indexed: 05/26/2024] Open
Abstract
Details of excitation and ionization acts hide a description of the biological effects of charged particle traversal through living tissue. Nanodosimetry enables the introduction of novel quantities that characterize and quantify the particle track structure while also serving as a foundation for assessing biological effects based on this quantification. This presents an opportunity to enhance the planning of charged particle radiotherapy by taking into account the ionization detail. This work uses Monte Carlo simulations with Geant4-DNA code for a wide variety of charged particles and their radiation qualities to analyze the distribution of ionization cluster sizes within nanometer-scale volumes, similar to DNA diameter. By correlating these results with biological parameters extracted from the PIDE database for the V79 cell line, a novel parameter R2 based on ionization details is proposed for the evaluation of radiation quality in terms of biological consequences, i.e., radiobiological cross section for inactivation. By incorporating the probability p of sub-lethal damage caused by a single ionization, we address limitations associated with the usually proposed nanodosimetric parameter Fk for characterizing the biological effects of radiation. We show that the new parameter R2 correlates well with radiobiological data and can be used to predict biological outcomes.
Collapse
Affiliation(s)
- Monika Mietelska
- Biomedical Physics Division, Institute of Experimental Physics, Faculty of Physics, University of Warsaw, 02-093 Warsaw, Poland;
- Radiological Metrology and Biomedical Physics Division, Nuclear Facilities Operations Department, National Centre for Nuclear Research, 05-400 Świerk, Poland; (M.P.); (A.B.)
| | - Marcin Pietrzak
- Radiological Metrology and Biomedical Physics Division, Nuclear Facilities Operations Department, National Centre for Nuclear Research, 05-400 Świerk, Poland; (M.P.); (A.B.)
- Laboratory of Translational Imaging in Oncology, Inserm, Institut Curie, Université Paris Saclay, 91401 Orsay, France
| | - Aleksandr Bancer
- Radiological Metrology and Biomedical Physics Division, Nuclear Facilities Operations Department, National Centre for Nuclear Research, 05-400 Świerk, Poland; (M.P.); (A.B.)
| | | | | | - Beata Brzozowska
- Biomedical Physics Division, Institute of Experimental Physics, Faculty of Physics, University of Warsaw, 02-093 Warsaw, Poland;
| |
Collapse
|
2
|
Ballisat L, De Sio C, Beck L, Guatelli S, Sakata D, Shi Y, Duan J, Velthuis J, Rosenfeld A. Dose and DNA damage modelling of diffusing alpha-emitters radiation therapy using Geant4. Phys Med 2024; 121:103367. [PMID: 38701625 DOI: 10.1016/j.ejmp.2024.103367] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/08/2023] [Revised: 04/02/2024] [Accepted: 04/19/2024] [Indexed: 05/05/2024] Open
Abstract
PURPOSE Diffusing alpha-emitters radiation therapy (DaRT) is a brachytherapy technique using α-particles to treat solid tumours. The high linear energy transfer (LET) and short range of α-particles make them good candidates for the targeted treatment of cancer. Treatment planning of DaRT requires a good understanding of the dose from α-particles and the other particles released in the 224Ra decay chain. METHODS The Geant4 Monte Carlo toolkit has been used to simulate a DaRT seed to better understand the dose contribution from all particles and simulate the DNA damage due to this treatment. RESULTS Close to the seed α-particles deliver the majority of dose, however at radial distances greater than 4 mm, the contribution of β-particles is greater. The RBE has been estimated as a function of number of double strand breaks (DSBs) and complex DSBs. A maximum seed spacing of 5.5 mm and 6.5 mm was found to deliver at least 20 Gy RBE weighted dose between the seeds for RBEDSB and RBEcDSB respectively. CONCLUSIONS The DNA damage changes with radial distance from the seed and has been found to become less complex with distance, which is potentially easier for the cell to repair. Close to the seed α-particles contribute the majority of dose, however the contribution from other particles cannot be neglected and may influence the choice of seed spacing.
Collapse
Affiliation(s)
| | - Chiara De Sio
- School of Physics, University of Bristol, Bristol, UK
| | - Lana Beck
- School of Physics, University of Bristol, Bristol, UK
| | - Susanna Guatelli
- Centre for Medical Radiation Physics (CMRP), University of Wollongong, NSW, Australia
| | - Dousatsu Sakata
- Division of Health Sciences, Osaka University, Osaka 565-0871, Japan
| | - Yuyao Shi
- School of Physics, University of Bristol, Bristol, UK
| | - Jinyan Duan
- School of Physics, University of Bristol, Bristol, UK
| | - Jaap Velthuis
- School of Physics, University of Bristol, Bristol, UK
| | - Anatoly Rosenfeld
- Centre for Medical Radiation Physics (CMRP), University of Wollongong, NSW, Australia
| |
Collapse
|
3
|
Bae C, Hernández Millares R, Ryu S, Moon H, Kim D, Lee G, Jiang Z, Park MH, Kim KH, Koom WS, Ye SJ, Lee K. Synergistic Effect of Ferroptosis-Inducing Nanoparticles and X-Ray Irradiation Combination Therapy. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2024; 20:e2310873. [PMID: 38279618 DOI: 10.1002/smll.202310873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/24/2023] [Revised: 12/18/2023] [Indexed: 01/28/2024]
Abstract
Ferroptosis, characterized by the induction of cell death via lipid peroxidation, has been actively studied over the last few years and has shown the potential to improve the efficacy of cancer nanomedicine in an iron-dependent manner. Radiation therapy, a common treatment method, has limitations as a stand-alone treatment due to radiation resistance and safety as it affects even normal tissues. Although ferroptosis-inducing drugs help alleviate radiation resistance, there are no safe ferroptosis-inducing drugs that can be considered for clinical application and are still in the research stage. Here, the effectiveness of combined treatment with radiotherapy with Fe and hyaluronic acid-based nanoparticles (FHA-NPs) to directly induce ferroptosis, considering the clinical applications is reported. Through the induction of ferroptosis by FHA-NPs and apoptosis by X-ray irradiation, the therapeutic efficiency of cancer is greatly improved both in vitro and in vivo. In addition, Monte Carlo simulations are performed to assess the physical interactions of the X-rays with the iron-oxide nanoparticle. The study provides a deeper understanding of the synergistic effect of ferroptosis and X-ray irradiation combination therapy. Furthermore, the study can serve as a valuable reference for elucidating the role and mechanisms of ferroptosis in radiation therapy.
Collapse
Affiliation(s)
- Chaewon Bae
- Program in Nanoscience and Technology, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
| | - Rodrigo Hernández Millares
- Program in Biomedical Radiation Sciences, Department of Transdisciplinary Studies, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
| | - Suhyun Ryu
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
| | - Hyowon Moon
- Department of Radiation Oncology, Yonsei Cancer Center, Heavy Ion Therapy Research Institute, Yonsei University College of Medicine, Seoul, 03722, Republic of Korea
| | - Dongwoo Kim
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
| | - Gyubok Lee
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
| | - Zhuomin Jiang
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
| | - Min Hee Park
- THEDONEE, 1208, 156, Gwanggyo-ro, Yeongtong-gu, Suwon-si, Gyeonggi-do, 16506, Republic of Korea
| | - Kyung Hwan Kim
- Department of Radiation Oncology, Yonsei Cancer Center, Heavy Ion Therapy Research Institute, Yonsei University College of Medicine, Seoul, 03722, Republic of Korea
| | - Woong Sub Koom
- Department of Radiation Oncology, Yonsei Cancer Center, Heavy Ion Therapy Research Institute, Yonsei University College of Medicine, Seoul, 03722, Republic of Korea
| | - Sung-Joon Ye
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
- Biomedical Research Institute, Seoul National University Hospital, Seoul, 03080, South Korea
- Advanced Institute of Convergence Technology, Seoul National University, Suwon, 16229, South Korea
- Research Institute for Convergence Science, Seoul National University, Seoul, 08826, Republic of Korea
| | - Kangwon Lee
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
- Research Institute for Convergence Science, Seoul National University, Seoul, 08826, Republic of Korea
| |
Collapse
|
4
|
Rafiepour P, Sina S, Amoli ZA, Shekarforoush SS, Farajzadeh E, Mortazavi SMJ. A mechanistic simulation of induced DNA damage in a bacterial cell by X- and gamma rays: a parameter study. Phys Eng Sci Med 2024:10.1007/s13246-024-01424-x. [PMID: 38652348 DOI: 10.1007/s13246-024-01424-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2023] [Accepted: 04/07/2024] [Indexed: 04/25/2024]
Abstract
Mechanistic Monte Carlo simulations calculating DNA damage caused by ionizing radiation are highly dependent on the simulation parameters. In the present study, using the Geant4-DNA toolkit, the impact of different parameters on DNA damage induced in a bacterial cell by X- and gamma-ray irradiation was investigated. Three geometry configurations, including the simple (without DNA details), the random (a random multiplication of identical DNA segments), and the fractal (a regular replication of DNA segments using fractal Hilbert curves), were simulated. Also, three physics constructors implemented in Geant4-DNA, i.e., G4EmDNAPhysics_option2, G4EmDNAPhysics_option4, and G4EmDNAPhysics_option6, with two energy thresholds of 17.5 eV and 5-37.5 eV were compared for direct DNA damage calculations. Finally, a previously developed mathematical model of cell repair called MEDRAS (Mechanistic DNA Repair and Survival) was employed to compare the impact of physics constructors on the cell survival curve. The simple geometry leads to undesirable results compared to the random and fractal ones, highlighting the importance of simulating complex DNA structures in mechanistic simulation studies. Under the same conditions, the DNA damage calculated in the fractal geometry was more consistent with the experimental data. All physics constructors can be used alternatively with the fractal geometry, provided that an energy threshold of 17.5 eV is considered for recording direct DNA damage. All physics constructors represent a similar behavior in generating cell survival curves, although the slopes of the curves are different. Since the inverse of the slope of a bacterial cell survival curve (i.e., the D10-value) is highly sensitive to the simulation parameters, it is not logical to determine an optimal set of parameters for calculating the D10-value by Monte Carlo simulation.
Collapse
Affiliation(s)
- Payman Rafiepour
- Department of Nuclear Engineering, School of Mechanical Engineering, Shiraz University, Shiraz, Iran
| | - Sedigheh Sina
- Department of Nuclear Engineering, School of Mechanical Engineering, Shiraz University, Shiraz, Iran.
- Radiation research center, School of Mechanical Engineering, Shiraz University, Shiraz, Iran.
| | - Zahra Alizadeh Amoli
- Department of Food Hygiene and Public Health, School of Veterinary Medicine, Shiraz University, Shiraz, Iran
| | - Seyed Shahram Shekarforoush
- Department of Food Hygiene and Public Health, School of Veterinary Medicine, Shiraz University, Shiraz, Iran
| | - Ebrahim Farajzadeh
- Secondary Standard Dosimetry Laboratory (SSDL), Pars Isotope Co, Karaj, Iran
| | - Seyed Mohammad Javad Mortazavi
- Ionizing and Non-ionizing Radiation Protection Research Center (INIRPRC), Shiraz University of Medical Sciences, Shiraz, Iran
| |
Collapse
|
5
|
Hernández Millares R, Bae C, Kim SJ, Kim T, Park SY, Lee K, Ye SJ. Clonogenic assay and computational modeling using real cell images to study physical enhancement and cellular sensitization induced by metal nanoparticles under MV and kV X-ray irradiation. NANOSCALE 2024; 16:7110-7122. [PMID: 38501279 DOI: 10.1039/d3nr06257k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
This study was initiated due to the physically unexplainable tumor controls resulting from metal nanoparticle (MNP) experiments even under MV X-ray irradiation. A more accurate explanation of the mechanism of radiosensitization induced by MNP is warranted, considering both its physical dose enhancement and biological sensitization, as related research is lacking. Thus, we aimed to examine the intricate dynamics involved in MNP-induced radiosensitization. We conducted specifically designed clonogenic assays for the A549 lung cancer cell line with MNP irradiated by 6 MV and 300 kVp X-rays. Two types of MNP were employed: one based on iron oxide, promoting ferroptosis, and the other on gold nanoparticles known for inducing a significant dose enhancement, particularly at low-energy X-rays. We introduced the lethality enhancement factor (LEF) as the fraction in the cell killing attributed to biological sensitization. Subsequently, Monte Carlo simulations were conducted to evaluate the radial dose profiles for each MNP, corresponding to the physical enhancement. Finally, the local effect model was applied to the clonogenic assay results on real cell images. The LEF and the dose enhancement in the cytoplasm were incorporated to increase the accuracy in the average lethal events and, consequently, in the survival fraction. The results reveal an increased cell killing for both of the MNP under MV and kV X-ray irradiation. In both types of MNP, the LEF reveals a biological sensitization evident. The sensitizer enhancement ratio, derived from the calculations, exhibited only 3% and 1% relative differences compared to the conventional linear-quadratic model for gold and ferroptosis inducer nanoparticles, respectively. These findings indicate that MNPs sensitize cells via radiation through mechanisms akin to ferroptosis inducers, not exclusively relying on a physical dose enhancement. Their own contributions to survival fractions were successfully integrated into computational modeling.
Collapse
Affiliation(s)
- Rodrigo Hernández Millares
- Program in Biomedical Radiation Sciences, Department of Transdisciplinary Studies, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
| | - Chaewon Bae
- Program in Nanoscience and Technology, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea
| | - Seok-Jin Kim
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea.
| | - Taewan Kim
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea.
| | - So-Yeon Park
- Department of Radiation Oncology, Veterans Health Service Medical Center, Seoul, 05368, Republic of Korea
- Institute of Radiation Medicine, Seoul National University Medical Research Center, Seoul, 03080, Republic of Korea
| | - Kangwon Lee
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea.
- Research Institute for Convergence Science, Seoul National University, Seoul, 08826, Republic of Korea
| | - Sung-Joon Ye
- Department of Applied Bioengineering, Graduate School of Convergence Science and Technology, Seoul National University, Seoul, 08826, Republic of Korea.
- Research Institute for Convergence Science, Seoul National University, Seoul, 08826, Republic of Korea
- Advanced Institute of Convergence Technology, Seoul National University, Suwon 16229, Korea
| |
Collapse
|
6
|
Liu Y, Zhu K, Peng X, Luo S, Zhu J, Xiao W, He L, Wang X. Proton relative biological effectiveness for the induction of DNA double strand breaks based on Geant4. Biomed Phys Eng Express 2024; 10:035018. [PMID: 38181453 DOI: 10.1088/2057-1976/ad1bb9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2023] [Accepted: 01/05/2024] [Indexed: 01/07/2024]
Abstract
Uncertainties in the relative biological effectiveness (RBE) of proton remains a major barrier to the biological optimization of proton therapy. A large amount of experimental data suggest that proton RBE is variable. As an evolving Monte Carlo code toolkit, Geant4-DNA is able to simulate the initial DNA damage caused by particle beams through physical and chemical interactions at the nanometer scale over a short period of time. This contributes to evaluating the radiobiological effects induced by ionizing radiation. Based on the Geant4-DNA toolkit, this study constructed a DNA geometric model containing 6.32Gbp, simulated the relationship between radiochemical yields (G-values) and their corresponding chemical constructors, and calculated a detailed calculation of the sources of damage and the complexity of damage in DNA strand breaks. The damage model constructed in this study can simulate the relative biological effectiveness (RBE) in the proton Bragg peak region. The results indicate that: (1) When the electron energy is below 400 keV, the yield of OH·account for 18.1% to 25.3% of the total water radiolysis yields. (2) Under the influence of histone clearance function, the yield of indirect damage account for over 72.93% of the yield of DNA strand breaks (SBs). When linear energy transfer (LET) increased from 29.79 (keV/μm) to 64.29 (keV/μm), the yield of double strand breaks (DSB) increased from 17.27% to 32.65%. (3) By investigating the effect of proton Bragg peak depth on the yield of direct DSB (DSBdirect) and total DSB (DSBtotal), theRBEDSBtotandRBEDSBdirlevels of cells show that the RBE value of protons reaches 2.2 in the Bragg peak region.
Collapse
Affiliation(s)
- Yuchen Liu
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, People's Republic of China
| | - Kun Zhu
- State Key Laboratory of Radiation Medicine and Protection, School of Radiation Medicine and Protection, Soochow University, Suzhou 215123, People's Republic of China
| | - Xiaoyu Peng
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, People's Republic of China
| | - Siyuan Luo
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, People's Republic of China
| | - Jin Zhu
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, People's Republic of China
| | - Wancheng Xiao
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, People's Republic of China
| | - Lie He
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, People's Republic of China
| | - Xiaodong Wang
- School of Nuclear Science and Technology, University of South China, Hengyang 421001, People's Republic of China
| |
Collapse
|
7
|
Taheri A, Khandaker MU, Moradi F, Bradley DA. A simulation study on the radiosensitization properties of gold nanorods. Phys Med Biol 2024; 69:045029. [PMID: 38286017 DOI: 10.1088/1361-6560/ad2380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2023] [Accepted: 01/29/2024] [Indexed: 01/31/2024]
Abstract
Objective. Gold nanorods (GNRs) have emerged as versatile nanoparticles with unique properties, holding promise in various modalities of cancer treatment through drug delivery and photothermal therapy. In the rapidly evolving field of nanoparticle radiosensitization (NPRS) for cancer therapy, this study assessed the potential of gold nanorods as radiosensitizing agents by quantifying the key features of NPRS, such as secondary electron emission and dose enhancement, using Monte Carlo simulations.Approach. Employing the TOPAS track structure code, we conducted a comprehensive evaluation of the radiosensitization behavior of spherical gold nanoparticles and gold nanorods. We systematically explored the impact of nanorod geometry (in particular size and aspect ratio) and orientation on secondary electron emission and deposited energy ratio, providing validated results against previously published simulations.Main results. Our findings demonstrate that gold nanorods exhibit comparable secondary electron emission to their spherical counterparts. Notably, nanorods with smaller surface-area-to-volume ratios (SA:V) and alignment with the incident photon beam proved to be more efficient radiosensitizing agents, showing superiority in emitted electron fluence. However, in the microscale, the deposited energy ratio (DER) was not markedly influenced by the SA:V of the nanorod. Additionally, our findings revealed that the geometry of gold nanoparticles has a more significant impact on the emission of M-shell Auger electrons (with energies below 3.5 keV) than on higher-energy electrons.Significance. This research investigated the radiosensitization properties of gold nanorods, positioning them as promising alternatives to the more conventionally studied spherical gold nanoparticles in the context of cancer research. With increasing interest in multimodal cancer therapy, our findings have the potential to contribute valuable insights into the perspective of gold nanorods as effective multipurpose agents for synergistic photothermal therapy and radiotherapy. Future directions may involve exploring alternative metallic nanorods as well as further optimizing the geometry and coating materials, opening new possibilities for more effective cancer treatments.
Collapse
Affiliation(s)
- Ali Taheri
- Applied Physics and Radiation Technologies Group, CCDCU, School of Engineering and Technology, Sunway University, 47500 Bandar Sunway, Selangor, Malaysia
| | - Mayeen Uddin Khandaker
- Applied Physics and Radiation Technologies Group, CCDCU, School of Engineering and Technology, Sunway University, 47500 Bandar Sunway, Selangor, Malaysia
- Faculty of Graduate Studies, Daffodil International University, Daffodil Smart City, Birulia, Savar, Dhaka 1216, Bangladesh
| | - Farhad Moradi
- Fibre Optics Research Centre, Faculty of Engineering, Multimedia University, Jalan Multimedia 63100, Cyberjaya, Malaysia
| | - David Andrew Bradley
- Applied Physics and Radiation Technologies Group, CCDCU, School of Engineering and Technology, Sunway University, 47500 Bandar Sunway, Selangor, Malaysia
- School of Mathematics and Physics, University of Surrey, Guildford, GU2 7XH, United Kingdom
| |
Collapse
|
8
|
Feng H, Li W, Zhang Y, Chang C, Hua L, Feng Y, Lai Y, Geng L. Mechanistic modelling of relative biological effectiveness of carbon ion beams and comparison with experiments. Phys Med Biol 2024; 69:035020. [PMID: 38157549 DOI: 10.1088/1361-6560/ad1998] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2023] [Accepted: 12/29/2023] [Indexed: 01/03/2024]
Abstract
Objective.Relative biological effectiveness (RBE) plays a vital role in carbon ion radiotherapy, which is a promising treatment method for reducing toxic effects on normal tissues and improving treatment efficacy. It is important to have an effective and precise way of obtaining RBE values to support clinical decisions. A method of calculating RBE from a mechanistic perspective is reported.Approach.Ratio of dose to obtain the same number of double strand breaks (DSBs) between different radiation types was used to evaluate RBE. Package gMicroMC was used to simulate DSB yields. The DSB inductions were then analyzed to calculate RBE. The RBE values were compared with experimental results.Main results.Furusawa's experiment yielded RBE values of 1.27, 2.22, 3.00 and 3.37 for carbon ion beam with dose-averaged LET of 30.3 keVμm-1, 54.5 keVμm-1, 88 keVμm-1and 137 keVμm-1, respectively. RBE values computed from gMicroMC simulations were 1.75, 2.22, 2.87 and 2.97. When it came to a more sophisticated carbon ion beam with 6 cm spread-out Bragg peak, RBE values were 1.61, 1.63, 2.19 and 2.36 for proximal, middle, distal and distal end part, respectively. Values simulated by gMicroMC were 1.50, 1.87, 2.19 and 2.34. The simulated results were in reasonable agreement with the experimental data.Significance.As a mechanistic way for the evaluation of RBE for carbon ion radiotherapy by combining the macroscopic simulation of energy spectrum and microscopic simulation of DNA damages, this work provides a promising tool for RBE calculation supporting clinical applications such as treatment planning.
Collapse
Affiliation(s)
- Haonan Feng
- School of Physics, Beihang University, Beijing 102206, People's Republic of China
- Department of Medical Management, Chinese Academy of Science Heavy Ion Medicine (CASHIM) Co. Ltd, Beijing 100083, People's Republic of China
| | - Weiguang Li
- School of Physics, Beihang University, Beijing 102206, People's Republic of China
- Department of Medical Management, Chinese Academy of Science Heavy Ion Medicine (CASHIM) Co. Ltd, Beijing 100083, People's Republic of China
| | - Yibao Zhang
- Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Radiation Oncology, Peking University Cancer Hospital & Institute, Beijing 100142, People's Republic of China
| | - Cheng Chang
- Department of Medical Management, Chinese Academy of Science Heavy Ion Medicine (CASHIM) Co. Ltd, Beijing 100083, People's Republic of China
| | - Ling Hua
- Key Laboratory of Carcinogenesis and Translational Research (Ministry of Education/Beijing), Department of Radiation Oncology, Peking University Cancer Hospital & Institute, Beijing 100142, People's Republic of China
| | - Yiwen Feng
- Beijing University of Posts and Telecommunications, Beijing 100876, People's Republic of China
| | - Youfang Lai
- Department of Medical Management, Chinese Academy of Science Heavy Ion Medicine (CASHIM) Co. Ltd, Beijing 100083, People's Republic of China
| | - LiSheng Geng
- School of Physics, Beihang University, Beijing 102206, People's Republic of China
- Peng Huanwu Collaborative Center for Research and Education, Beihang University, Beijing 100191, People's Republic of China
- Beijing Key Laboratory of Advanced Nuclear Materials and Physics, Beihang University, Beijing 102206, People's Republic of China
- Southern Center for Nuclear-Science Theory (SCNT), Institute of Modern Physics, Chinese Academy of Sciences, Huizhou 516000, Guangdong Province, People's Republic of China
| |
Collapse
|
9
|
Friedrich T, Pfuhl T, Scholz M. Spectral composition of secondary electrons based on the Kiefer-Straaten ion track structure model. Phys Med Biol 2024; 69:025013. [PMID: 38118162 DOI: 10.1088/1361-6560/ad1765] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Accepted: 12/20/2023] [Indexed: 12/22/2023]
Abstract
The major part of energy deposition of ionizing radiation is caused by secondary electrons, independent of the primary radiation type. However, their spatial concentration and their spectral properties strongly depend on the primary radiation type and finally determine the pattern of molecular damage e.g. to biological targets as the DNA, and thus the final effect of the radiation exposure. To describe the physical and to predict the biological consequences of charged ion irradiation, amorphous track structure approaches have proven to be pragmatic and helpful. There, the local dose deposition in the ion track is equated by considering the emission and slowing down of the secondary electrons from the primary particle track. In the present work we exploit the model of Kiefer and Straaten and derive the spectral composition of secondary electrons as function of the distance to the track center. The spectral composition indicates differences to spectra of low linear energy transfer (LET) photon radiation, which we confirm by a comparison with Monte Carlo studies. We demonstrate that the amorphous track structure approach provides a simple tool for evaluating the spectral electron properties within the track structure. Predictions of the LET of electrons across the track structure as well as the electronic dose build-up effect are derived. Implications for biological effects and corresponding predicting models based on amorphous track structure are discussed.
Collapse
Affiliation(s)
- T Friedrich
- GSI Helmholtzzentrum für Schwerionenforschung, Planckstraße 1, D-64291 Darmstadt, Germany
| | - T Pfuhl
- GSI Helmholtzzentrum für Schwerionenforschung, Planckstraße 1, D-64291 Darmstadt, Germany
| | - M Scholz
- GSI Helmholtzzentrum für Schwerionenforschung, Planckstraße 1, D-64291 Darmstadt, Germany
| |
Collapse
|
10
|
Warmenhoven JW, Henthorn NT, McNamara AL, Ingram SP, Merchant MJ, Kirkby KJ, Schuemann J, Paganetti H, Prise KM, McMahon SJ. Effects of Differing Underlying Assumptions in In Silico Models on Predictions of DNA Damage and Repair. Radiat Res 2023; 200:509-522. [PMID: 38014593 DOI: 10.1667/rade-21-00147.1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2021] [Accepted: 10/05/2023] [Indexed: 11/29/2023]
Abstract
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.
Collapse
Affiliation(s)
- John W Warmenhoven
- Division of Cancer Sciences, University of Manchester, Manchester, United Kingdom
- The Christie NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester, United Kingdom
| | - Nicholas T Henthorn
- Division of Cancer Sciences, University of Manchester, Manchester, United Kingdom
- The Christie NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester, United Kingdom
| | - Aimee L McNamara
- Physics Division, Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Massachusetts
| | - Samuel P Ingram
- Division of Cancer Sciences, University of Manchester, Manchester, United Kingdom
- Christie Medical Physics and Engineering, The Christie NHS Foundation Trust, Manchester, United Kingdom
| | - Michael J Merchant
- Division of Cancer Sciences, University of Manchester, Manchester, United Kingdom
- The Christie NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester, United Kingdom
| | - Karen J Kirkby
- Division of Cancer Sciences, University of Manchester, Manchester, United Kingdom
- The Christie NHS Foundation Trust, Manchester Academic Health Science Centre, Manchester, United Kingdom
| | - Jan Schuemann
- Physics Division, Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Massachusetts
| | - Harald Paganetti
- Physics Division, Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Massachusetts
| | - Kevin M Prise
- Patrick G Johnston Centre for Cancer Research, Queen's University Belfast, Belfast, United Kingdom
| | - Stephen J McMahon
- Patrick G Johnston Centre for Cancer Research, Queen's University Belfast, Belfast, United Kingdom
| |
Collapse
|
11
|
Thibaut Y, Gonon G, Martinez JS, Petit M, Babut R, Vaurijoux A, Gruel G, Villagrasa C, Incerti S, Perrot Y. Experimental validation in a neutron exposure frame of the MINAS TIRITH for cell damage simulation. Phys Med Biol 2023; 68:225008. [PMID: 37848039 DOI: 10.1088/1361-6560/ad043d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Accepted: 10/17/2023] [Indexed: 10/19/2023]
Abstract
In the domains of medicine and space exploration, refining risk assessment models for protecting healthy tissue from ionizing radiation is crucial. Understanding radiation-induced effects requires biological experimentations at the cellular population level and the cellular scale modeling using Monte Carlo track structure codes. We present MINAS TIRITH, a tool using Geant4-DNA Monte Carlo-generated databases to study DNA damage distribution at the cell population scale. It introduces a DNA damage location module and proposes a method to convert double-strand breaks (DSB) into DNA Damage Response foci. We evaluate damage location precision and DSB-foci conversion parameters. MINAS TIRITH's accuracy is validated againstγ-H2AX foci distribution from cell population exposed to monoenergetic neutron beams (2.5 or 15.1 MeV) under different configurations, yielding mixed radiation fields. Strong agreement between simulation and experimental results was found demonstrating MINAS TIRITH's predictive precision in radiation-induced DNA damage topology. Additionally, modeling intercellular damage variability within a population subjected to a specific macroscopic dose identifies subpopulations, enhancing realistic fate models. This approach advances our understanding of radiation-induced effects on cellular systems for risk assessment improvement.
Collapse
Affiliation(s)
- Y Thibaut
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - G Gonon
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - J S Martinez
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - M Petit
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - R Babut
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - A Vaurijoux
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - G Gruel
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - C Villagrasa
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - S Incerti
- Université de Bordeaux, CNRS/IN2P3, LP2i, UMR 5797, F-33170 Gradignan, France
| | - Y Perrot
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| |
Collapse
|
12
|
Lamghari Y, Lu H, Bentourkia M. DNA damage by radiation as a function of electron energy and interaction at the atomic level with Monte Carlo simulation. Z Med Phys 2023; 33:489-498. [PMID: 35973908 PMCID: PMC10751702 DOI: 10.1016/j.zemedi.2022.07.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/19/2022] [Revised: 07/08/2022] [Accepted: 07/11/2022] [Indexed: 11/23/2022]
Abstract
In radiotherapy, X-ray or heavy ion beams target tumors to cause damage to their cell DNA. This damage is mainly induced by secondary low energy electrons. In this paper, we report the DNA molecular breaks at the atomic level as a function of electron energy and types of electron interactions using of Monte Carlo simulation. The number of DNA single and double strand breaks are compared to those from experimental results based on electron energies. In recent years, DNA atomistic models were introduced but still the simulations consider energy deposition in volumes of DNA or water equivalent material. We simulated a model of atomistic B-DNA in vacuum, forming 1122 base pairs of 30 nm in length. Each atom has been represented by a sphere whose radius equals the radius of van der Waals. We repeatedly simulated 10 million electrons for each energy from 4 eV to 500 eV and counted each interaction type with its position x, y, z in the volume of DNA. Based on the number and types of interactions at the atomic level, the number of DNA single and double strand breaks were calculated. We found that the dissociative electron attachment has the dominant effect on DNA strand breaks at energies below 10 eV compared to excitation and ionization. In addition, it is straightforward with our simulation to discriminate the strand and base breaks as a function of radiation interaction type and energy. In conclusion, the knowledge of DNA damage at the atomic level helps design direct internal therapeutic agents of cancer treatment.
Collapse
Affiliation(s)
- Youssef Lamghari
- Department of Nuclear Medicine and Radiobiology, University of Sherbrooke, Sherbrooke, QC, Canada
| | - Huizhong Lu
- Department of Nuclear Medicine and Radiobiology, University of Sherbrooke, Sherbrooke, QC, Canada
| | - M'hamed Bentourkia
- Department of Nuclear Medicine and Radiobiology, University of Sherbrooke, Sherbrooke, QC, Canada.
| |
Collapse
|
13
|
Kolovi S, Fois GR, Lanouar S, Chardon P, Miallier D, Baker LA, Bailly C, Beauger A, Biron DG, David K, Montavon G, Pilleyre T, Schoefs B, Breton V, Maigne L. Assessing radiation dosimetry for microorganisms in naturally radioactive mineral springs using GATE and Geant4-DNA Monte Carlo simulations. PLoS One 2023; 18:e0292608. [PMID: 37824461 PMCID: PMC10569590 DOI: 10.1371/journal.pone.0292608] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2022] [Accepted: 09/25/2023] [Indexed: 10/14/2023] Open
Abstract
Mineral springs in Massif Central, France can be characterized by higher levels of natural radioactivity in comparison to the background. The biota in these waters is constantly under radiation exposure mainly from the α-emitters of the natural decay chains, with 226Ra in sediments ranging from 21 Bq/g to 43 Bq/g and 222Rn activity concentrations in water up to 4600 Bq/L. This study couples for the first time micro- and nanodosimetric approaches to radioecology by combining GATE and Geant4-DNA to assess the dose rates and DNA damages to microorganisms living in these naturally radioactive ecosystems. It focuses on unicellular eukaryotic microalgae (diatoms) which display an exceptional abundance of teratological forms in the most radioactive mineral springs in Auvergne. Using spherical geometries for the microorganisms and based on γ-spectrometric analyses, we evaluate the impact of the external exposure to 1000 Bq/L 222Rn dissolved in the water and 30 Bq/g 226Ra in the sediments. Our results show that the external dose rates for diatoms are significant (9.7 μGy/h) and comparable to the threshold (10 μGy/h) for the protection of the ecosystems suggested by the literature. In a first attempt of simulating the radiation induced DNA damage on this species, the rate of DNA Double Strand Breaks per day is estimated to 1.11E-04. Our study confirms the significant mutational pressure from natural radioactivity to which microbial biodiversity has been exposed since Earth origin in hydrothermal springs.
Collapse
Affiliation(s)
- Sofia Kolovi
- Laboratoire de Physique de Clermont (LPC) - UMR6533, CNRS/IN2P3 Université Clermont Auvergne, Aubière, France
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
| | - Giovanna-Rosa Fois
- Laboratoire de Physique de Clermont (LPC) - UMR6533, CNRS/IN2P3 Université Clermont Auvergne, Aubière, France
| | - Sarra Lanouar
- Laboratoire de Physique de Clermont (LPC) - UMR6533, CNRS/IN2P3 Université Clermont Auvergne, Aubière, France
| | - Patrick Chardon
- Laboratoire de Physique de Clermont (LPC) - UMR6533, CNRS/IN2P3 Université Clermont Auvergne, Aubière, France
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
| | - Didier Miallier
- Laboratoire de Physique de Clermont (LPC) - UMR6533, CNRS/IN2P3 Université Clermont Auvergne, Aubière, France
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
| | - Lory-Anne Baker
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
- Laboratoire Microorganismes: Génome Environnement (LMGE) - UMR6023, CNRS, Université Clermont Auvergne, Clermont–Ferrand, France
- Laboratoire de Géographie Physique et Environnementale (GEOLAB) - UMR6042, CNRS, Université Clermont Auvergne, Clermont-Ferrand, France
| | - Céline Bailly
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
- Laboratoire de Physique Subatomique et des Technologies Associées (SUBATECH) - UMR6457, CNRS/IN2P3/IMT Atlantique/Université de Nantes, Nantes, France
| | - Aude Beauger
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
- Laboratoire Microorganismes: Génome Environnement (LMGE) - UMR6023, CNRS, Université Clermont Auvergne, Clermont–Ferrand, France
| | - David G. Biron
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
- Laboratoire Microorganismes: Génome Environnement (LMGE) - UMR6023, CNRS, Université Clermont Auvergne, Clermont–Ferrand, France
| | - Karine David
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
- Laboratoire de Physique Subatomique et des Technologies Associées (SUBATECH) - UMR6457, CNRS/IN2P3/IMT Atlantique/Université de Nantes, Nantes, France
| | - Gilles Montavon
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
- Laboratoire de Physique Subatomique et des Technologies Associées (SUBATECH) - UMR6457, CNRS/IN2P3/IMT Atlantique/Université de Nantes, Nantes, France
| | - Thierry Pilleyre
- Laboratoire de Physique de Clermont (LPC) - UMR6533, CNRS/IN2P3 Université Clermont Auvergne, Aubière, France
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
| | - Benoît Schoefs
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
- Metabolism, Molecular Engineering of Microalgae and Applications, Laboratoire de Biologie des Organismes, Stress, Santé Environnement, IUML FR3473, CNRS, Le Mans University, Le Mans, France
| | - Vincent Breton
- Laboratoire de Physique de Clermont (LPC) - UMR6533, CNRS/IN2P3 Université Clermont Auvergne, Aubière, France
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
| | - Lydia Maigne
- Laboratoire de Physique de Clermont (LPC) - UMR6533, CNRS/IN2P3 Université Clermont Auvergne, Aubière, France
- LTSER “Zone Atelier Territoires Uranifères”, Clermont-Ferrand, France
| | | |
Collapse
|
14
|
Nicolanti F, Caccia B, Cartoni A, Emfietzoglou D, Faccini R, Incerti S, Kyriakou I, Satta M, Tran HN, Mancini-Terracciano C. Calculation of electron interaction models in N 2 and O 2. Phys Med 2023; 114:102661. [PMID: 37703804 DOI: 10.1016/j.ejmp.2023.102661] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 06/19/2023] [Accepted: 08/05/2023] [Indexed: 09/15/2023] Open
Abstract
Cosmic rays have the potential to significantly affect the atmospheric composition by increasing the rate and changing the types of chemical reactions through ion production. The amount and states of ionization, and the spatial distribution of ions produced are still open questions for atmospheric models. To precisely estimate these quantities, it is necessary to simulate particle-molecule interactions, down to very low energies. Models enabling such simulations require interaction probabilities over a broad energy range and for all energetically allowed scattering processes. In this paper, we focus on electron interaction with the two most abundant molecules in the atmosphere, i.e., N2 and O2, as an initial step. A set of elastic and inelastic cross section models for electron transportation in oxygen and nitrogen molecules valid in the energy range 10 eV - 1 MeV, is presented. Comparison is made with available theoretical and experimental data and a reasonable good agreement is observed. Stopping power is calculated and compared with published data to assess the general consistency and reliability of our results. Good overall agreement is observed, with relative differences lower than 6% with the ESTAR database.
Collapse
Affiliation(s)
- F Nicolanti
- Physics Dep., Sapienza U. of Rome, p.le Aldo Moro, 5, 00185, Rome, Italy; INFN, Sec. of Rome, p.le Aldo Moro, 2, 00185, Rome, Italy.
| | - B Caccia
- ISS(Italian National Institute of Health), V. Regina Elena, 299, 00161, Rome, Italy
| | - A Cartoni
- Chemistry Dep., Sapienza U. of Rome, p.le Aldo Moro, 5, 00185, Rome, Italy
| | - D Emfietzoglou
- Med. Phys. Lab., Dept of Medicine, University of Ioannina, 45110, Ioannina, Greece
| | - R Faccini
- Physics Dep., Sapienza U. of Rome, p.le Aldo Moro, 5, 00185, Rome, Italy; INFN, Sec. of Rome, p.le Aldo Moro, 2, 00185, Rome, Italy
| | - S Incerti
- Université de Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, Chemin du Solarium, 19, 33170, Gradignan, France
| | - I Kyriakou
- Med. Phys. Lab., Dept of Medicine, University of Ioannina, 45110, Ioannina, Greece
| | - M Satta
- Chemistry Dep., Sapienza U. of Rome, p.le Aldo Moro, 5, 00185, Rome, Italy; ISMN-CNR, p. Aldo Moro, 7, 00185, Rome, Italy
| | - H N Tran
- Université de Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, Chemin du Solarium, 19, 33170, Gradignan, France
| | - C Mancini-Terracciano
- Physics Dep., Sapienza U. of Rome, p.le Aldo Moro, 5, 00185, Rome, Italy; INFN, Sec. of Rome, p.le Aldo Moro, 2, 00185, Rome, Italy
| |
Collapse
|
15
|
Carrasco-Hernandez J, Ramos-Méndez J, Padilla-Rodal E, Avila-Rodriguez MA. Cellular lethal damage of 64Cu incorporated in mammalian genome evaluated with Monte Carlo methods. Front Med (Lausanne) 2023; 10:1253746. [PMID: 37841004 PMCID: PMC10575761 DOI: 10.3389/fmed.2023.1253746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2023] [Accepted: 09/13/2023] [Indexed: 10/17/2023] Open
Abstract
Purpose Targeted Radionuclide Therapy (TRT) with Auger Emitters (AE) is a technique that allows targeting specific sites on tumor cells using radionuclides. The toxicity of AE is critically dependent on its proximity to the DNA. The aim of this study is to quantify the DNA damage and radiotherapeutic potential of the promising AE radionuclide copper-64 (64Cu) incorporated into the DNA of mammalian cells using Monte Carlo track-structure simulations. Methods A mammalian cell nucleus model with a diameter of 9.3 μm available in TOPAS-nBio was used. The cellular nucleus consisted of double-helix DNA geometrical model of 2.3 nm diameter surrounded by a hydration shell with a thickness of 0.16 nm, organized in 46 chromosomes giving a total of 6.08 giga base-pairs (DNA density of 14.4 Mbp/μm3). The cellular nucleus was irradiated with monoenergetic electrons and radiation emissions from several radionuclides including 111In, 125I, 123I, and 99mTc in addition to 64Cu. For monoenergetic electrons, isotropic point sources randomly distributed within the nucleus were modeled. The radionuclides were incorporated in randomly chosen DNA base pairs at two positions near to the central axis of the double-helix DNA model at (1) 0.25 nm off the central axis and (2) at the periphery of the DNA (1.15 nm off the central axis). For all the radionuclides except for 99mTc, the complete physical decay process was explicitly simulated. For 99mTc only total electron spectrum from published data was used. The DNA Double Strand Breaks (DSB) yield per decay from direct and indirect actions were quantified. Results obtained for monoenergetic electrons and radionuclides 111In, 125I, 123I, and 99mTc were compared with measured and calculated data from the literature for verification purposes. The DSB yields per decay incorporated in DNA for 64Cu are first reported in this work. The therapeutic effect of 64Cu (activity that led 37% cell survival after two cell divisions) was determined in terms of the number of atoms incorporated into the nucleus that would lead to the same DSBs that 100 decays of 125I. Simulations were run until a 2% statistical uncertainty (1 standard deviation) was achieved. Results The behavior of DSBs as a function of the energy for monoenergetic electrons was consistent with published data, the DSBs increased with the energy until it reached a maximum value near 500 eV followed by a continuous decrement. For 64Cu, when incorporated in the genome at evaluated positions (1) and (2), the DSB were 0.171 ± 0.003 and 0.190 ± 0.003 DSB/decay, respectively. The number of initial atoms incorporated into the genome (per cell) for 64Cu that would cause a therapeutic effect was estimated as 3,107 ± 28, that corresponds to an initial activity of 47.1 ± 0.4 × 10-3 Bq. Conclusion Our results showed that TRT with 64Cu has comparable therapeutic effects in cells as that of TRT with radionuclides currently used in clinical practice.
Collapse
Affiliation(s)
- Jhonatan Carrasco-Hernandez
- Departamento de Estructura de la Materia, Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - José Ramos-Méndez
- Department of Radiation Oncology, University of California, San Francisco, San Francisco, CA, United States
| | - Elizabeth Padilla-Rodal
- Departamento de Estructura de la Materia, Instituto de Ciencias Nucleares, Universidad Nacional Autónoma de México, Mexico City, Mexico
| | - Miguel A. Avila-Rodriguez
- Unidad Radiofarmacia-Ciclotrón, Facultad de Medicina, Universidad Nacional Autónoma de México, Mexico City, Mexico
| |
Collapse
|
16
|
Rumiantcev M, Li WB, Lindner S, Liubchenko G, Resch S, Bartenstein P, Ziegler SI, Böning G, Delker A. Estimation of relative biological effectiveness of 225Ac compared to 177Lu during [ 225Ac]Ac-PSMA and [ 177Lu]Lu-PSMA radiopharmaceutical therapy using TOPAS/TOPAS-nBio/MEDRAS. EJNMMI Phys 2023; 10:53. [PMID: 37695374 PMCID: PMC10495309 DOI: 10.1186/s40658-023-00567-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2023] [Accepted: 08/07/2023] [Indexed: 09/12/2023] Open
Abstract
AIM Over recent years, [225Ac]Ac-PSMA and [177Lu]Lu-PSMA radiopharmaceutical therapy have evolved as a promising treatment option for advanced prostate cancer. Especially for alpha particle emitter treatments, there is still a need for improving dosimetry, which requires accurate values of relative biological effectiveness (RBE). To achieve that, consideration of DNA damages in the cell nucleus and knowledge of the energy deposition in the location of the DNA at the nanometer scale are required. Monte Carlo particle track structure simulations provide access to interactions at this level. The aim of this study was to estimate the RBE of 225Ac compared to 177Lu. The initial damage distribution after radionuclide decay and the residual damage after DNA repair were considered. METHODS This study employed the TOol for PArtcile Simulation (TOPAS) based on the Geant4 simulation toolkit. Simulation of the nuclear DNA and damage scoring were performed using the TOPAS-nBio extension of TOPAS. DNA repair was modeled utilizing the Python-based program MEDRAS (Mechanistic DNA Repair and Survival). Five different cell geometries of equal volume and two radionuclide internalization assumptions as well as two cell arrangement scenarios were investigated. The radionuclide activity (number of source points) was adopted based on SPECT images of patients undergoing the above-mentioned therapies. RESULTS Based on the simulated dose-effect curves, the RBE of 225Ac compared to 177Lu was determined in a wide range of absorbed doses to the nucleus. In the case of spherical geometry, 3D cell arrangement and full radionuclide internalization, the RBE based on the initial damage had a constant value of approximately 2.14. Accounting for damage repair resulted in RBE values ranging between 9.38 and 1.46 for 225Ac absorbed doses to the nucleus between 0 and 50 Gy, respectively. CONCLUSION In this work, the consideration of DNA repair of the damage from [225Ac]Ac-PSMA and [177Lu]Lu-PSMA revealed a dose dependency of the RBE. Hence, this work suggested that DNA repair is an important aspect to understand response to different radiation qualities.
Collapse
Affiliation(s)
- Mikhail Rumiantcev
- Department of Nuclear Medicine, LMU University Hospital, LMU Munich, Munich, Germany.
| | - Wei Bo Li
- Federal Office for Radiation Protection, Medical and Occupational Radiation Protection, Oberschleißheim, Germany
| | - Simon Lindner
- Department of Nuclear Medicine, LMU University Hospital, LMU Munich, Munich, Germany
| | - Grigory Liubchenko
- Department of Nuclear Medicine, LMU University Hospital, LMU Munich, Munich, Germany
| | - Sandra Resch
- Department of Nuclear Medicine, LMU University Hospital, LMU Munich, Munich, Germany
| | - Peter Bartenstein
- Department of Nuclear Medicine, LMU University Hospital, LMU Munich, Munich, Germany
| | - Sibylle I Ziegler
- Department of Nuclear Medicine, LMU University Hospital, LMU Munich, Munich, Germany
| | - Guido Böning
- Department of Nuclear Medicine, LMU University Hospital, LMU Munich, Munich, Germany
| | - Astrid Delker
- Department of Nuclear Medicine, LMU University Hospital, LMU Munich, Munich, Germany
| |
Collapse
|
17
|
Han Y, Geng C, Liu Y, Wu R, Li M, Yu C, Altieri S, Tang X. Calculation of the DNA damage yield and relative biological effectiveness in boron neutron capture therapy via the Monte Carlo track structure simulation. Phys Med Biol 2023; 68:175028. [PMID: 37524085 DOI: 10.1088/1361-6560/acec2a] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2023] [Accepted: 07/31/2023] [Indexed: 08/02/2023]
Abstract
Objective.Boron neutron capture therapy (BNCT) is an advanced cellular-level hadron therapy that has exhibited remarkable therapeutic efficacy in the treatment of locally invasive malignancies. Despite its clinical success, the intricate nature of relative biological effectiveness (RBE) and mechanisms responsible for DNA damage remains elusive. This work aims to quantify the RBE of compound particles (i.e. alpha and lithium) in BNCT based on the calculation of DNA damage yields via the Monte Carlo track structure (MCTS) simulation.Approach. The TOPAS-nBio toolkit was employed to conduct MCTS simulations. The calculations encompassed four steps: determination of the angle and energy spectra on the nuclear membrane, quantification of the database containing DNA damage yields for ions with specific angle and energy, accumulation of the database and spectra to obtain the DNA damage yields of compound particles, and calculation of the RBE by comparison yields of double-strand break (DSB) with the reference gamma-ray. Furthermore, the impact of cell size and microscopic boron distribution was thoroughly discussed.Main results. The DSB yields induced by compound particles in three types of spherical cells (radius equal to 10, 8, and 6μm) were found to be 13.28, 17.34, 22.15 Gy Gbp-1for boronophenylalanine (BPA), and 1.07, 3.45, 8.32 Gy Gbp-1for sodium borocaptate (BSH). The corresponding DSB-based RBE values were determined to be 1.90, 2.48, 3.16 for BPA and 0.15, 0.49, 1.19 for BSH. The calculated DSB-based RBE showed agreement with experimentally values of compound biological effectiveness for melanoma and gliosarcoma. Besides, the DNA damage yield and DSB-based RBE value exhibited an increasing trend as the cell radius decreased. The impact of the boron concentration ratio on RBE diminished once the drug enrichment surpasses a certain threshold.Significance. This work is potential to provide valuable guidance for accurate biological-weighted dose evaluation in BNCT.
Collapse
Affiliation(s)
- Yang Han
- Department of Nuclear Science and Technology, Nanjing University of Aeronautics and Astronautics, Nanjing, People's Republic of China
- Department of Physics, University of Pavia, Pavia, Italy
| | - Changran Geng
- Department of Nuclear Science and Technology, Nanjing University of Aeronautics and Astronautics, Nanjing, People's Republic of China
| | - Yuanhao Liu
- Department of Nuclear Science and Technology, Nanjing University of Aeronautics and Astronautics, Nanjing, People's Republic of China
- Neuboron Medtech. Ltd, Nanjing, People's Republic of China
| | - Renyao Wu
- Department of Nuclear Science and Technology, Nanjing University of Aeronautics and Astronautics, Nanjing, People's Republic of China
| | - Mingzhu Li
- Department of Nuclear Science and Technology, Nanjing University of Aeronautics and Astronautics, Nanjing, People's Republic of China
| | - Chenxi Yu
- Department of Nuclear Science and Technology, Nanjing University of Aeronautics and Astronautics, Nanjing, People's Republic of China
| | - Saverio Altieri
- Department of Physics, University of Pavia, Pavia, Italy
- Istituto Nazionale di Fisica Nucleare (INFN), the section of Pavia, Pavia, Italy
| | - Xiaobin Tang
- Department of Nuclear Science and Technology, Nanjing University of Aeronautics and Astronautics, Nanjing, People's Republic of China
| |
Collapse
|
18
|
Suárez-García D, Cortés-Giraldo MA, Bertolet A. 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? Radiother Oncol 2023; 185:109730. [PMID: 37301260 PMCID: PMC10528084 DOI: 10.1016/j.radonc.2023.109730] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2022] [Revised: 05/20/2023] [Accepted: 05/31/2023] [Indexed: 06/12/2023]
Abstract
PURPOSE To 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 MATERIAL Our 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. RESULTS Domain 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. CONCLUSIONS Large 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.
Collapse
Affiliation(s)
- Daniel Suárez-García
- Departamento de Física Nuclear, Atómica y Molecular, Universidad de Sevilla, Sevilla, Spain
| | | | - Alejandro Bertolet
- Department of Radiation Oncology, Massachusetts General Hospital and Harvard Medical School, Boston, MA, USA.
| |
Collapse
|
19
|
Han Y, Geng C, D-Kondo JN, Li M, Ramos-Méndez J, Altieri S, Liu Y, Tang X. Microdosimetric Analysis for Boron Neutron Capture Therapy via Monte Carlo Track Structure Simulation with Modified Lithium Cross-sections. Radiat Phys Chem Oxf Engl 1993 2023; 209:110956. [PMID: 37206625 PMCID: PMC10191410 DOI: 10.1016/j.radphyschem.2023.110956] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Boron neutron capture therapy (BNCT) is a cellular-level hadron therapy achieving therapeutic effects via the synergistic action of multiple particles, including Lithium, alpha, proton, and photon. However, evaluating the relative biological effectiveness (RBE) in BNCT remains challenging. In this research, we performed a microdosimetric calculation for BNCT using the Monte Carlo track structure (MCTS) simulation toolkit, TOPAS-nBio. This paper reports the first attempt to derive the ionization cross-sections of low-energy (>0.025 MeV/u) Lithium for MCTS simulation based on the effective charge cross-section scalation method and phenomenological double-parameter modification. The fitting parameters λ1=1.101,λ2=3.486 were determined to reproduce the range and stopping power data from the ICRU report 73. Besides, the lineal energy spectra of charged particles in BNCT were calculated, and the influence of sensitive volume (SV) size was discussed. Condensed history simulation obtained similar results with MCTS when using Micron-SV while overestimating the lineal energy when using Nano-SV. Furthermore, we found that the microscopic boron distribution can significantly affect the lineal energy for Lithium, while the effect for alpha is minimal. Similar results to the published data by PHITS simulation were observed for the compound particles and monoenergetic protons when using micron-SV. Spectra with nano-SV reflected that the different track densities and absorbed doses in the nucleus together result in the dramatic difference in the macroscopic biological response of BPA and BSH. This work and the developed methodology could impact the research fields in BNCT where understanding radiation effects is crucial, such as the treatment planning system, source evaluation, and new boron drug development.
Collapse
Affiliation(s)
- Yang Han
- Nanjing University of Aeronautics and Astronautics, Department of Nuclear Science and Technology, Nanjing, 210016, China
- Joint International Research Laboratory on Advanced Particle Therapy, Nanjing 211100, China
- University of Pavia, Department of Physics, Pavia, 27100, Italy
| | - Changran Geng
- Nanjing University of Aeronautics and Astronautics, Department of Nuclear Science and Technology, Nanjing, 210016, China
- Joint International Research Laboratory on Advanced Particle Therapy, Nanjing 211100, China
- (+86 13601582233), (+86 15380796769)
| | - J. Naoki D-Kondo
- University of California San Francisco, Department of Radiation Oncology, San Francisco, CA 94115, USA
| | - Mingzhu Li
- Nanjing University of Aeronautics and Astronautics, Department of Nuclear Science and Technology, Nanjing, 210016, China
- Joint International Research Laboratory on Advanced Particle Therapy, Nanjing 211100, China
| | - José Ramos-Méndez
- University of California San Francisco, Department of Radiation Oncology, San Francisco, CA 94115, USA
| | - Saverio Altieri
- University of Pavia, Department of Physics, Pavia, 27100, Italy
- Istituto Nazionale di Fisica Nucleare (INFN), the section of Pavia, Pavia, 27100, Italy
| | - Yuanhao Liu
- Nanjing University of Aeronautics and Astronautics, Department of Nuclear Science and Technology, Nanjing, 210016, China
- Joint International Research Laboratory on Advanced Particle Therapy, Nanjing 211100, China
| | - Xiaobin Tang
- Nanjing University of Aeronautics and Astronautics, Department of Nuclear Science and Technology, Nanjing, 210016, China
- Joint International Research Laboratory on Advanced Particle Therapy, Nanjing 211100, China
- (+86 13601582233), (+86 15380796769)
| |
Collapse
|
20
|
Ballisat L, Beck L, De Sio C, Guatelli S, Sakata D, Incerti S, Tran HN, Duan J, Maclean K, Shi Y, Velthuis J, Rosenfeld A. In-silico calculations of DNA damage induced by α-particles in the 224Ra DaRT decay chain for a better understanding of the radiobiological effectiveness of this treatment. Phys Med 2023; 112:102626. [PMID: 37393861 DOI: 10.1016/j.ejmp.2023.102626] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 06/07/2023] [Accepted: 06/13/2023] [Indexed: 07/04/2023] Open
Abstract
Diffusing alpha-emitters radiation Therapy (DaRT) is an interstitial brachytherapy technique using 224Ra seeds. For accurate treatment planning a good understanding of the early DNA damage due to α-particles is required. Geant4-DNA was used to calculate the initial DNA damage and radiobiological effectiveness due to α-particles with linear energy transfer (LET) values in the range 57.5-225.9 keV/μm from the 224Ra decay chain. The impact of DNA base pair density on DNA damage has been modelled, as this parameter varies between human cell lines. Results show that the quantity and complexity of DNA damage changes with LET as expected. Indirect damage, due to water radical reactions with the DNA, decreases and becomes less significant at higher LET values as shown in previous studies. As expected, the yield of complex double strand breaks (DSBs), which are harder for a cell to repair, increases approximately linearly with LET. The level of complexity of DSBs and radiobiological effectiveness have been found to increase with LET as expected. The quantity of DNA damage has been shown to increase for increased DNA density in the expected base pair density range of human cells. The change in damage yield as a function of base pair density is largest for higher LET α-particles, an increase of over 50% for individual strand breaks between 62.7 and 127.4 keV/μm. This change in yield shows that the DNA base pair density is an important parameter for modelling DNA damage particularly at higher LET where the DNA damage is greatest and most complex.
Collapse
Affiliation(s)
| | - Lana Beck
- School of Physics, University of Bristol, Bristol, UK
| | - Chiara De Sio
- School of Physics, University of Bristol, Bristol, UK
| | - Susanna Guatelli
- Centre for Medical Radiation Physics (CMRP), University of Wollongong, NSW, Australia
| | - Dousatsu Sakata
- Division of Health Sciences, Osaka University, Osaka 565-0871, Japan
| | - Sébastien Incerti
- University of Bordeaux, CNRS, LP2I, UMR 5797, F-33170 Gradignan, France
| | - Hoang Ngoc Tran
- University of Bordeaux, CNRS, LP2I, UMR 5797, F-33170 Gradignan, France
| | - Jinyan Duan
- School of Physics, University of Bristol, Bristol, UK
| | - Katie Maclean
- School of Physics, University of Bristol, Bristol, UK
| | - Yuyao Shi
- School of Physics, University of Bristol, Bristol, UK
| | - Jaap Velthuis
- School of Physics, University of Bristol, Bristol, UK
| | - Anatoly Rosenfeld
- Centre for Medical Radiation Physics (CMRP), University of Wollongong, NSW, Australia
| |
Collapse
|
21
|
Muggiolu G, Torfeh E, Simon M, Devès G, Seznec H, Barberet P. Recruitment Kinetics of XRCC1 and RNF8 Following MeV Proton and α-Particle Micro-Irradiation. BIOLOGY 2023; 12:921. [PMID: 37508352 PMCID: PMC10376363 DOI: 10.3390/biology12070921] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/03/2023] [Revised: 06/19/2023] [Accepted: 06/20/2023] [Indexed: 07/30/2023]
Abstract
Time-lapse fluorescence imaging coupled to micro-irradiation devices provides information on the kinetics of DNA repair protein accumulation, from a few seconds to several minutes after irradiation. Charged-particle microbeams are valuable tools for such studies since they provide a way to selectively irradiate micrometric areas within a cell nucleus, control the dose and the micro-dosimetric quantities by means of advanced detection systems and Monte Carlo simulations and monitor the early cell response by means of beamline microscopy. We used the charged-particle microbeam installed at the AIFIRA facility to perform micro-irradiation experiments and measure the recruitment kinetics of two proteins involved in DNA signaling and repair pathways following exposure to protons and α-particles. We developed and validated image acquisition and processing methods to enable a systematic study of the recruitment kinetics of GFP-XRCC1 and GFP-RNF8. We show that XRCC1 is recruited to DNA damage sites a few seconds after irradiation as a function of the total deposited energy and quite independently of the particle LET. RNF8 is recruited to DNA damage sites a few minutes after irradiation and its recruitment kinetics depends on the particle LET.
Collapse
Affiliation(s)
| | - Eva Torfeh
- University Bordeaux, CNRS, LP2I, UMR 5797, 33170 Gradignan, France
| | - Marina Simon
- University Bordeaux, CNRS, LP2I, UMR 5797, 33170 Gradignan, France
| | - Guillaume Devès
- University Bordeaux, CNRS, LP2I, UMR 5797, 33170 Gradignan, France
| | - Hervé Seznec
- University Bordeaux, CNRS, LP2I, UMR 5797, 33170 Gradignan, France
| | | |
Collapse
|
22
|
Bian J, Duran J, Shin WG, Ramos-Méndez J, Sankey JC, Childress L, Seuntjens J, Enger SA. GEANT4-DNA simulation of temperature-dependent and pH-dependent yields of chemical radiolytic species. Phys Med Biol 2023; 68:124002. [PMID: 37230081 DOI: 10.1088/1361-6560/acd90d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2022] [Accepted: 05/25/2023] [Indexed: 05/27/2023]
Abstract
Objective.GEANT4-DNA can simulate radiation chemical yield (G-value) for radiolytic species such as the hydrated electron (eaq-) with the independent reaction times (IRT) method, however, only at room temperature and neutral pH. This work aims to modify the GEANT4-DNA source code to enable the calculation ofG-values for radiolytic species at different temperatures and pH values.Approach.In the GEANT4-DNA source code, values of chemical parameters such as reaction rate constant, diffusion coefficient, Onsager radius, and water density were replaced by corresponding temperature-dependent polynomials. The initial concentration of hydrogen ion (H+)/hydronium ion (H3O+) was scaled for a desired pH using the relationship pH = -log10[H+]. To validate our modifications, two sets of simulations were performed. (A) A water cube with 1.0 km sides and a pH of 7 was irradiated with an isotropic electron source of 1 MeV. The end time was 1μs. The temperatures varied from 25 °C to 150 °C. (B) The same setup as (A) was used, however, the temperature was set to 25 °C while the pH varied from 5 to 9. The results were compared with published experimental and simulated work.Main results.The IRT method in GEANT4-DNA was successfully modified to simulateG-values for radiolytic species at different temperatures and pH values. Our temperature-dependent results agreed with experimental data within 0.64%-9.79%, and with simulated data within 3.52%-12.47%. The pH-dependent results agreed well with experimental data within 0.52% to 3.19% except at a pH of 5 (15.99%) and with simulated data within 4.40%-5.53%. The uncertainties were below ±0.20%. Overall our results agreed better with experimental than simulation data.Significance.Modifications in the GEANT4-DNA code enabled the calculation ofG-values for radiolytic species at different temperatures and pH values.
Collapse
Affiliation(s)
- Jingyi Bian
- Medical Physics Unit, Department of Oncology, Faculty of Medicine, McGill University, Montreal, Quebec, Canada
| | - Juan Duran
- Medical Physics Unit, Department of Oncology, Faculty of Medicine, McGill University, Montreal, Quebec, Canada
| | - Wook-Geun Shin
- Physics Division, Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA-02114, United States of America
| | - Jose Ramos-Méndez
- Department of Radiation Oncology, University of California San Francisco, CA, United States of America
| | - Jack C Sankey
- Department of Physics, McGill University, Montreal, Quebec, Canada
| | - Lilian Childress
- Department of Physics, McGill University, Montreal, Quebec, Canada
| | - Jan Seuntjens
- Medical Physics Unit, Department of Oncology, Faculty of Medicine, McGill University, Montreal, Quebec, Canada
- Research Institute of the McGill University Health Centre, Montreal, Quebec, Canada
| | - Shirin A Enger
- Medical Physics Unit, Department of Oncology, Faculty of Medicine, McGill University, Montreal, Quebec, Canada
- Research Institute of the McGill University Health Centre, Montreal, Quebec, Canada
- Lady Davis Institute for Medical Research, Jewish General Hospital, Montreal, Quebec, Canada
| |
Collapse
|
23
|
Chappuis F, Tran HN, Zein SA, Bailat C, Incerti S, Bochud F, Desorgher L. The general-purpose Geant4 Monte Carlo toolkit and its Geant4-DNA extension to investigate mechanisms underlying the FLASH effect in radiotherapy: Current status and challenges. Phys Med 2023; 110:102601. [PMID: 37201453 DOI: 10.1016/j.ejmp.2023.102601] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/17/2023] [Revised: 04/06/2023] [Accepted: 05/01/2023] [Indexed: 05/20/2023] Open
Abstract
FLASH radiotherapy is a promising approach to cancer treatment that offers several advantages over conventional radiotherapy. With this novel technique, high doses of radiation are delivered in a short period of time, inducing the so-called FLASH effect - a phenomenon characterized by healthy tissue sparing without alteration of tumor control. The mechanisms behind the FLASH effect remain unknown. One way to approach this problem is to gain insight into the initial parameters that can distinguish FLASH from conventional irradiation by simulating particle transport in aqueous media using the general-purpose Geant4 Monte Carlo toolkit and its Geant4-DNA extension. This review article discusses the current status of Geant4 and Geant4-DNA simulations to investigate mechanisms underlying the FLASH effect, as well as the challenges faced in this research field. One of the primary challenges is to accurately simulate the experimental irradiation parameters. Another challenge is the temporal extension of the simulations. This review also focuses on two hypotheses to explain the FLASH effect - namely the oxygen depletion hypothesis and the inter-track interactions hypothesis - and discusses how the Geant4 toolkit can be used to investigate them. The aim of this review is to provide an overview of Geant4 and Geant4-DNA simulations for FLASH radiotherapy and to highlight the challenges that need to be overcome in order to better study the FLASH effect.
Collapse
Affiliation(s)
- Flore Chappuis
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland
| | - Hoang Ngoc Tran
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | - Sara A Zein
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | - Claude Bailat
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland
| | - Sébastien Incerti
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | - François Bochud
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland
| | - Laurent Desorgher
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland.
| |
Collapse
|
24
|
Wanstall HC, Henthorn NT, Jones J, Santina E, Chadwick AL, Angal-Kalinin D, Morris G, Warmenhoven JW, Smith R, Mathisen S, Merchant MJ, Jones RM. Quantification of damage to plasmid DNA from 35 MeV electrons, 228 MeV protons and 300 kVp X-rays in varying hydroxyl radical scavenging environments. JOURNAL OF RADIATION RESEARCH 2023:7153712. [PMID: 37154587 DOI: 10.1093/jrr/rrad032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/29/2022] [Revised: 03/22/2023] [Indexed: 05/10/2023]
Abstract
The pBR322 plasmid DNA was irradiated with 35 MeV electrons, 228 MeV protons and 300 kVp X-rays to quantify DNA damage and make comparisons of DNA damage between radiation modalities. Plasmid was irradiated in a medium containing hydroxyl radical scavengers in varying concentrations. This altered the amount of indirect hydroxyl-mediated DNA damage, to create an environment that is more closely associated with a biological cell. We show that increasing hydroxyl scavenger concentration significantly reduced post-irradiation DNA damage to pBR322 plasmid DNA consistently and equally with three radiation modalities. At low scavenging capacities, irradiation with both 35 MeV electrons and 228 MeV protons resulted in increased DNA damage per dose compared with 300 kVp X-rays. We quantify both single-strand break (SSB) and double-strand break (DSB) induction between the modalities as a ratio of yields relative to X-rays, referred to as relative biological effectiveness (RBE). RBESSB values of 1.16 ± 0.15 and 1.18 ± 0.08 were calculated for protons and electrons, respectively, in a low hydroxyl scavenging environment containing 1 mM Tris-HCl for SSB induction. In higher hydroxyl scavenging capacity environments (above 1.1 × 106 s-1), no significant differences in DNA damage induction were found between radiation modalities when using SSB induction as a measure of RBE. Considering DSB induction, significant differences were only found between X-rays and 35 MeV electrons, with an RBEDSB of 1.72 ± 0.91 for 35 MeV electrons, indicating that electrons result in significantly more SSBs and DSBs per unit of dose than 300 kVp X-rays.
Collapse
Affiliation(s)
- Hannah C Wanstall
- Department of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
- Manchester Academic Health Science Centre, The Christie NHS Foundation Trust, Wilmslow Road, Manchester M20 4BX, UK
- The Cockcroft Institute, Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
| | - Nicholas T Henthorn
- Manchester Academic Health Science Centre, The Christie NHS Foundation Trust, Wilmslow Road, Manchester M20 4BX, UK
- Division of Cancer Sciences, Faculty of Biology, Medicine and Health, School of Medical Sciences, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - James Jones
- The Cockcroft Institute, Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
- ASTeC, STFC Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
| | - Elham Santina
- Division of Cancer Sciences, Faculty of Biology, Medicine and Health, School of Medical Sciences, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Amy L Chadwick
- Division of Cancer Sciences, Faculty of Biology, Medicine and Health, School of Medical Sciences, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Deepa Angal-Kalinin
- The Cockcroft Institute, Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
- ASTeC, STFC Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
| | - Geoffrey Morris
- The Cockcroft Institute, Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
- ASTeC, STFC Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
| | - John-William Warmenhoven
- Manchester Academic Health Science Centre, The Christie NHS Foundation Trust, Wilmslow Road, Manchester M20 4BX, UK
- Division of Cancer Sciences, Faculty of Biology, Medicine and Health, School of Medical Sciences, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Rob Smith
- The Cockcroft Institute, Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
- ASTeC, STFC Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
| | - Storm Mathisen
- The Cockcroft Institute, Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
- ASTeC, STFC Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
| | - Michael J Merchant
- Manchester Academic Health Science Centre, The Christie NHS Foundation Trust, Wilmslow Road, Manchester M20 4BX, UK
- Division of Cancer Sciences, Faculty of Biology, Medicine and Health, School of Medical Sciences, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Roger M Jones
- Department of Physics and Astronomy, Faculty of Science and Engineering, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
- The Cockcroft Institute, Daresbury Laboratory, Daresbury, Warrington WA4 4AD, UK
| |
Collapse
|
25
|
Margis S, Kyriakou I, Incerti S, Bordage MC, Emfietzoglou D. Sub-keV corrections to binary encounter cross section models for electron ionization of liquid water with application to the Geant4-DNA Monte Carlo code. Appl Radiat Isot 2023; 194:110693. [PMID: 36731390 DOI: 10.1016/j.apradiso.2023.110693] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2022] [Revised: 12/13/2022] [Accepted: 01/23/2023] [Indexed: 01/26/2023]
Abstract
INTRODUCTION The electron ionization cross section of water is one of the most important input in Monte Carlo studies of cellular radiobiological effects. Analytical cross section models of the binary-encounter type have the potential of reducing simulation time and facilitate application to a variety of biological materials (other than water). The Binary-Encounter-Bethe (BEB) and Binary-Encounter-Dipole (BED) models of NIST are perhaps the most popular of such models giving reliable results for atoms and molecules in the gas-phase over a wide energy range. However, the use of such models to sub-keV electron energies in liquid water raises concerns due to the neglect of condensed phase effects that leads to a significant overestimation when compared to medium-specific dielectric models. PURPOSE To modify the BEB and BED models towards better agreement with the recommended low-energy dielectric model of Geant4-DNA (Option 4). To implement the new modifications to the existing BEB model of the Option 6 physics constructor of Geant4-DNA and re-evaluate fundamental transport quantities for sub-keV electrons. METHODS In analogy to a Yukawa potential a simple, yet physically-motivated, modification of the Burgess correction term is proposed to account for the reduction of the Coulomb interaction due to the polarizability of the target. The magnitude of the correction is guided by the dielectric-based ionization cross section implemented in Option 4. RESULTS Differential, total and stopping ionization cross sections for low-energy electrons in liquid water are presented. When combined with the Vriens correction (which is not included in Option 6), the proposed modification to the BEB and BED models brings the ionization and stopping cross sections in much better agreement against those used in the Option 4 dielectric model of Geant4-DNA, with up to 30% and 10% deviation, respectively. Implementation of the new correction to the Option 6 constructor of Geant4-DNA and re-evaluation of fundamental transport quantities, such as electron penetration ranges and dose-point-kernels, reduced the discrepancies from Option 4 at sub-keV energies from 20 to 100% (or more) to well below 10% in most cases. CONCLUSIONS A simple modification to the BEB and BED analytic models was found to improve their performance for sub-keV electrons in liquid water medium. Implementation of the new modification to the Option 6 constructor of Geant4-DNA significantly improved the agreement with the recommended low-energy Option 4 constructor for a variety of fundamental quantities related to electron transport.
Collapse
Affiliation(s)
- Stefanos Margis
- Medical Physics Laboratory, University of Ioannina Medical School, 45110, Ioannina, Greece
| | - Ioanna Kyriakou
- Medical Physics Laboratory, University of Ioannina Medical School, 45110, Ioannina, Greece
| | - Sebastien Incerti
- Bordeaux University, CNRS/IN2P3, CENBG, UMR 5797, F-33170, Gradignan, France
| | | | - Dimitris Emfietzoglou
- Medical Physics Laboratory, University of Ioannina Medical School, 45110, Ioannina, Greece.
| |
Collapse
|
26
|
Chappuis F, Grilj V, Tran HN, Zein SA, Bochud F, Bailat C, Incerti S, Desorgher L. Modeling of scavenging systems in water radiolysis with Geant4-DNA. Phys Med 2023; 108:102549. [PMID: 36921424 DOI: 10.1016/j.ejmp.2023.102549] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 10/05/2022] [Revised: 01/11/2023] [Accepted: 02/13/2023] [Indexed: 03/16/2023] Open
Abstract
PURPOSE This paper presents the capabilities of the Geant4-DNA Monte Carlo toolkit to simulate water radiolysis with scavengers using the step-by-step (SBS) or the independent reaction times (IRT) methods. It features two examples of application areas: (1) computing the escape yield of H2O2 following a 60Co γ-irradiation and (2) computing the oxygen depletion in water irradiated with 1 MeV electrons. METHODS To ease the implementation of the chemical stage in Geant4-DNA, we developed a user interface that helps define the chemical reactions and set the concentration of scavengers. The first application area example required two computational steps to perform water radiolysis using NO2- and NO3- as scavengers and a 60Co irradiation. The oxygen depletion computation technique for the second application area example consisted of simulating track segments of 1 MeV electrons and determining the radio-induced loss and gain of oxygen molecules. RESULTS The production of H2O2 under variable scavenging levels is consistent with the literature; the mean relative difference between the SBS and IRT methods is 7.2 % ± 0.5 %. For the oxygen depletion 1 µs post-irradiation, the mean relative difference between both methods is equal to 9.8 % ± 0.3 %. The results in the microsecond scale depend on the initial partial pressure of oxygen in water. In addition, the computed oxygen depletions agree well with the literature. CONCLUSIONS The Geant4-DNA toolkit makes it possible to simulate water radiolysis in the presence of scavengers. This feature offers perspectives in radiobiology, with the possibility of simulating cell-relevant scavenging mechanisms.
Collapse
Affiliation(s)
- Flore Chappuis
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland
| | - Veljko Grilj
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland
| | - Hoang Ngoc Tran
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | - Sara A Zein
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | - François Bochud
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland
| | - Claude Bailat
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland
| | - Sébastien Incerti
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | - Laurent Desorgher
- Institute of Radiation Physics (IRA), Lausanne University Hospital and University of Lausanne, CH-1007 Lausanne, Switzerland.
| |
Collapse
|
27
|
Modelling tissue specific RBE for different radiation qualities based on a multiscale characterization of energy deposition. Radiother Oncol 2023; 182:109539. [PMID: 36806602 DOI: 10.1016/j.radonc.2023.109539] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2022] [Revised: 02/03/2023] [Accepted: 02/13/2023] [Indexed: 02/18/2023]
Abstract
PURPOSE We present the nanoCluE model, which uses nano- and microdosimetric quantities to model RBE for protons and carbon ions. Under the hypothesis that nano- and microdosimetric quantities correlates with the generation of complex DNA double strand breakes, we wish to investigate whether an improved accuracy in predicting LQ parameters may be achieved, compared to some of the published RBE models. METHODS The model is based on experimental LQ data for protons and carbon ions. We generated a database of track structure data for a number of proton and carbon ion kinetic energies with the Geant4-DNA Monte Carlo code. These data were used to obtain both a nanodosimetric quantity and a set of microdosimetric quantities. The latter were tested with different parameterizations versus experimental LQ-data to select the variable and parametrization that yielded the best fit. RESULTS For protons, the nanoCluE model yielded, for the ratio of the linear LQ term versus the test data, a root mean square error (RMSE) of 1.57 compared to 1.31 and 1.30 for two earlier other published proton models. For carbon ions the RMSE was 2.26 compared to 3.24 and 5.24 for earlier published carbon ion models. CONCLUSION These results demonstrate the feasibility of the nanoCluE RBE model for carbon ions and protons. The increased accuracy for carbon ions as compared to two other considered models warrants further investigation.
Collapse
|
28
|
Hallier DC, Smales GJ, Seitz H, Hahn MB. Bio-SAXS of single-stranded DNA-binding proteins: radiation protection by the compatible solute ectoine. Phys Chem Chem Phys 2023; 25:5372-5382. [PMID: 36637121 DOI: 10.1039/d2cp05053f] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Small-angle X-ray scattering (SAXS) can be used for structural determination of biological macromolecules and polymers in their native states (e.g. liquid phase). This means that the structural changes of (bio-)polymers, such as proteins and DNA, can be monitored in situ to understand their sensitivity to changes in chemical environments. In an attempt to improve the reliability of such experiments, the reduction of radiation damage occurring from exposure to X-rays is required. One such method, is to use scavenger molecules to protect macromolecules against radicals produced during radiation exposure, such as reactive oxygen species (ROS). In this study we investigate the feasibility of applying the compatible solute, osmolyte and radiation protector Ectoine (THP(B)), as a scavenger molecule during SAXS measurements of the single-stranded DNA-binding protein Gene-V Protein (G5P/GVP). In this case, we monitor the radiation induced changes of G5P during bio-SAXS measurments and the resulting microscopic energy-damage relation was determined from microdosimetric calculations by Monte-Carlo based particle scattering simulations with TOPAS/Geant4 and a custom target-model. This resulted in a median-lethal energy deposit of pure G5P at 4 mg mL-1 of E1/2 = 7 ± 5 eV, whereas a threefold increase of energy-deposit was needed under the presence of Ectoine to reach the same level of damage. This indicates that Ectoine increases the possible exposure time before radiation-damage to G5P is observed. Furthermore, the dominant type of damage shifted from aggregation in pure solutions towards a fragmentation for solutions containing Ectoine as a cosolute. These results are interpreted in terms of indirect radiation damage by reactive secondary species, as well as post-irradiation effects, related to preferential-exclusion of the cosolute from the protein surface. Hence, Ectoine is shown to provide a non-disturbing way to improve structure-determination of proteins via bio-SAXS in future studies.
Collapse
Affiliation(s)
- Dorothea C Hallier
- Universität Potsdam, Institut für Biochemie und Biologie, 14476 Potsdam, Germany.,Fraunhofer Institute for Cell Therapy and Immunology, Branch Bioanalytics and Bioprocesses (IZI-BB), 14476 Potsdam, Germany.,Bundesanstalt für Materialforschung und -prüfung (BAM), 12205 Berlin, Germany.
| | - Glen J Smales
- Bundesanstalt für Materialforschung und -prüfung (BAM), 12205 Berlin, Germany.
| | - Harald Seitz
- Universität Potsdam, Institut für Biochemie und Biologie, 14476 Potsdam, Germany.,Fraunhofer Institute for Cell Therapy and Immunology, Branch Bioanalytics and Bioprocesses (IZI-BB), 14476 Potsdam, Germany
| | - Marc Benjamin Hahn
- Bundesanstalt für Materialforschung und -prüfung (BAM), 12205 Berlin, Germany.
| |
Collapse
|
29
|
Chatzipapas KP, Tran NH, Dordevic M, Zivkovic S, Zein S, Shin W, Sakata D, Lampe N, Brown JMC, Ristic‐Fira A, Petrovic I, Kyriakou I, Emfietzoglou D, Guatelli S, Incerti S. Simulation of DNA damage using Geant4‐DNA: an overview of the “molecularDNA” example application. PRECISION RADIATION ONCOLOGY 2023. [DOI: 10.1002/pro6.1186] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023] Open
Affiliation(s)
| | - Ngoc Hoang Tran
- University of Bordeaux, CNRS, LP2I Bordeaux, UMR 5797 Gradignan France
| | - Milos Dordevic
- Vinca Institute of Nuclear Sciences, National Institute of the Republic of Serbia University of Belgrade, Vinca Belgrade Serbia
| | - Sara Zivkovic
- Vinca Institute of Nuclear Sciences, National Institute of the Republic of Serbia University of Belgrade, Vinca Belgrade Serbia
| | - Sara Zein
- University of Bordeaux, CNRS, LP2I Bordeaux, UMR 5797 Gradignan France
| | - Wook‐Geun Shin
- Physics Division, Department of Radiation Oncology Massachusetts General Hospital & Harvard Medical School Boston Massachusetts USA
| | | | | | - Jeremy M. C. Brown
- Department of Physics and Astronomy Swinburne University of Technology Melbourne Australia
| | - Aleksandra Ristic‐Fira
- Vinca Institute of Nuclear Sciences, National Institute of the Republic of Serbia University of Belgrade, Vinca Belgrade Serbia
| | - Ivan Petrovic
- Vinca Institute of Nuclear Sciences, National Institute of the Republic of Serbia University of Belgrade, Vinca Belgrade Serbia
| | - Ioanna Kyriakou
- Medical Physics Laboratory Department of Medicine University of Ioannina Ioannina Greece
| | - Dimitris Emfietzoglou
- Medical Physics Laboratory Department of Medicine University of Ioannina Ioannina Greece
| | - Susanna Guatelli
- Centre for Medical Radiation Physics University of Wollongong Wollongong New South Wales Australia
| | - Sébastien Incerti
- University of Bordeaux, CNRS, LP2I Bordeaux, UMR 5797 Gradignan France
| |
Collapse
|
30
|
Rädler M, Buizza G, Kawula M, Palaniappan P, Gianoli C, Baroni G, Paganelli C, Parodi K, Riboldi M. Impact of secondary particles on the magnetic field generated by a proton pencil beam: a finite-element analysis based on Geant4-DNA simulations. Med Phys 2023; 50:1000-1018. [PMID: 36346042 DOI: 10.1002/mp.16062] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Revised: 09/22/2022] [Accepted: 09/27/2022] [Indexed: 11/11/2022] Open
Abstract
PURPOSE To investigate the static magnetic field generated by a proton pencil beam as a candidate for range verification by means of Monte Carlo simulations, thereby improving upon existing analytical calculations. We focus on the impact of statistical current fluctuations and secondary protons and electrons. METHODS We considered a pulsed beam (10 μ ${\umu}$ s pulse duration) during the duty cycle with a peak beam current of 0.2 μ $\umu$ A and an initial energy of 100 MeV. We ran Geant4-DNA Monte Carlo simulations of a proton pencil beam in water and extracted independent particle phase spaces. We calculated longitudinal and radial current density of protons and electrons, serving as an input for a magnetic field estimation based on a finite element analysis in a cylindrical geometry. We made sure to allow for non-solenoidal current densities as is the case of a stopping proton beam. RESULTS The rising proton charge density toward the range is not perturbed by energy straggling and only lowered through nuclear reactions by up to 15%, leading to an approximately constant longitudinal current. Their relative low density however (at most 0.37 protons/mm3 for the 0.2 μ ${\umu}$ A current and a beam cross-section of 2.5 mm), gives rise to considerable current density fluctuations. The radial proton current resulting from lateral scattering and being two orders of magnitude weaker than the longitudinal current is subject to even stronger fluctuations. Secondary electrons with energies above 10 eV, that far outnumber the primary protons, reduce the primary proton current by only 10% due to their largely isotropic flow. A small fraction of electrons (<1%), undergoing head-on collisions, constitutes the relevant electron current. In the far-field, both contributions to the magnetic field strength (longitudinal and lateral) are independent of the beam spot size. We also find that the nuclear reaction-related losses cause a shift of 1.3 mm to the magnetic field profile relative to the actual range, which is further enlarged to 2.4 mm by the electron current (at a distance of ρ = 50 $\rho =50$ mm away from the central beam axis). For ρ > 45 $\rho >45$ mm, the shift increases linearly. While the current density variations cause significant magnetic field uncertainty close to the central beam axis with a relative standard deviation (RSD) close to 100%, they average out at a distance of 10 cm, where the RSD of the total magnetic field drops below 2%. CONCLUSIONS With the small influence of the secondary electrons together with the low RSD, our analysis encourages an experimental detection of the magnetic field through sensitive instrumentation, such as optical magnetometry or SQUIDs.
Collapse
Affiliation(s)
- Martin Rädler
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Munich, Germany
| | - Giulia Buizza
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci, Milan, Italy
| | - Maria Kawula
- Department of Radiation Oncology, LMU Hospital, Munich, Germany
| | - Prasannakumar Palaniappan
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Munich, Germany
| | - Chiara Gianoli
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Munich, Germany
| | - Guido Baroni
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci, Milan, Italy
| | - Chiara Paganelli
- Department of Electronics, Information and Bioengineering, Politecnico di Milano, Piazza Leonardo da Vinci, Milan, Italy
| | - Katia Parodi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Munich, Germany
| | - Marco Riboldi
- Department of Medical Physics, Faculty of Physics, Ludwig-Maximilians-Universität München, Munich, Germany
| |
Collapse
|
31
|
Thibaut Y, Gonon G, Martinez JS, Petit M, Vaurijoux A, Gruel G, Villagrasa C, Incerti S, Perrot Y. MINAS TIRITH: a new tool for simulating radiation-induced DNA damage at the cell population level. Phys Med Biol 2023; 68. [PMID: 36623319 DOI: 10.1088/1361-6560/acb196] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Accepted: 01/09/2023] [Indexed: 01/11/2023]
Abstract
Objective. The mechanisms of radiation-induced DNA damage can be understood via the fundamental acquisition of knowledge through a combination of experiments and modeling. Currently, most biological experiments are performed by irradiating an entire cell population, whereas modeling of radiation-induced effects is usually performed via Monte Carlo simulations with track structure codes coupled to realistic DNA geometries of a single-cell nucleus. However, the difference in scale between the two methods hinders a direct comparison because the dose distribution in the cell population is not necessarily uniform owing to the stochastic nature of the energy deposition. Thus, this study proposed the MINAS TIRITH tool to model the distribution of radiation-induced DNA damage in a cell population.Approach. The proposed method is based on precomputed databases of microdosimetric parameters and DNA damage distributions generated using the Geant4-DNA Monte Carlo Toolkit. First, a specific energyzwas assigned to each cell of an irradiated population for a particular absorbed doseDabs,following microdosimetric formalism. Then, each cell was assigned a realistic number of DNA damage events according to the specific energyz,respecting the stochastic character of its occurrence.Main results. This study validated the MINAS TIRITH tool by comparing its results with those obtained using the Geant4-DNA track structure code and a Geant4-DNA based simulation chain for DNA damage calculation. The different elements of comparison indicated consistency between MINAS TIRITH and the Monte Carlo simulation in case of the dose distribution in the population and the calculation of the amount of DNA damage.Significance. MINAS TIRITH is a new approach for the calculation of radiation-induced DNA damage at the cell population level that facilitates reasonable simulation times compared to those obtained with track structure codes. Moreover, this tool enables a more direct comparison between modeling and biological experimentation.
Collapse
Affiliation(s)
- Y Thibaut
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - G Gonon
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - J S Martinez
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - M Petit
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - A Vaurijoux
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - G Gruel
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - C Villagrasa
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| | - S Incerti
- Université de Bordeaux, CNRS/IN2P3, LP2i, UMR 5797, F-33170 Gradignan, France
| | - Y Perrot
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), PSE-SANTE/SDOS/LDRI, PSE-SANTE/SERAMED/LRAcc, PSE-SANTE/SDOS/LMDN, BP 17, F-92262 Fontenay-aux-Roses, France
| |
Collapse
|
32
|
Sakata D, Hirayama R, Shin WG, Belli M, Tabocchini MA, Stewart RD, Belov O, Bernal MA, Bordage MC, Brown JMC, Dordevic M, Emfietzoglou D, Francis Z, Guatelli S, Inaniwa T, Ivanchenko V, Karamitros M, Kyriakou I, Lampe N, Li Z, Meylan S, Michelet C, Nieminen P, Perrot Y, Petrovic I, Ramos-Mendez J, Ristic-Fira A, Santin G, Schuemann J, Tran HN, Villagrasa C, Incerti S. Prediction of DNA rejoining kinetics and cell survival after proton irradiation for V79 cells using Geant4-DNA. Phys Med 2023; 105:102508. [PMID: 36549067 DOI: 10.1016/j.ejmp.2022.11.012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/02/2021] [Revised: 11/28/2022] [Accepted: 11/30/2022] [Indexed: 12/24/2022] Open
Abstract
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.
Collapse
Affiliation(s)
- Dousatsu Sakata
- Department of Accelerator and Medical Physics, Institute for Quantum Medical Science, National Institutes for Quantum Science and Technology (QST), Chiba 263-8555, Japan; Division of Health Sciences, Osaka University, Osaka 565-0871, Japan.
| | - Ryoichi Hirayama
- Department of Charged Particle Therapy Research, Institute for Quantum Medical Science, National Institutes for Quantum Science and Technology (QST), Chiba 263-8555, Japan
| | - Wook-Geun Shin
- Department of Radiation Oncology, Seoul National University Hospital, Seoul 03080, Republic of Korea
| | | | | | - Robert D Stewart
- Department of Radiation Oncology, University of Washington, WA 98195-6043, USA
| | - Oleg Belov
- Veksler and Baldin Laboratory of High Energy Physics, Joint Institute for Nuclear Research, 141980 Dubna, Russia; Institute of System Analysis and Management, Dubna State University, 141980 Dubna, Russia
| | - Mario A Bernal
- Instituto de Fisica Gleb Wataghin, Universidade Estadual de Campinas, Campinas, SP, Brazil
| | - Marie-Claude Bordage
- INSERM, Université Paul Sabatier, UMR 1037, CRCT, Toulouse, France; Université Toulouse III-Paul Sabatier, UMR 1037, CRCT, Toulouse, France
| | - Jeremy M C Brown
- Department of Physics and Astronomy, Swinburne University of Technology, Hawthorn, Australia; Centre For Medical Radiation Physics, University of Wollongong, Wollongong, Australia; Department of Radiation Science and Technology, Delft University of Technology, The Netherlands
| | - Milos Dordevic
- Vinca Institute of Nuclear Sciences - National Institute of the Republic of Serbia, University of Belgrade, Belgrade, Serbia
| | - Dimitris Emfietzoglou
- Medical Physics Laboratory, Department of Medicine, University of Ioannina, GR 45110, Ioannina, Greece
| | - Ziad Francis
- Saint Joseph University of Beirut, UR Mathématiques et Modélisation, Beirut, Lebanon
| | - Susanna Guatelli
- Centre For Medical Radiation Physics, University of Wollongong, Wollongong, Australia
| | - Taku Inaniwa
- Department of Accelerator and Medical Physics, Institute for Quantum Medical Science, National Institutes for Quantum Science and Technology (QST), Chiba 263-8555, Japan
| | - Vladimir Ivanchenko
- Geant4 Associates International Ltd, Hebden Bridge, UK; Tomsk State University, Tomsk, Russia
| | | | - Ioanna Kyriakou
- Medical Physics Laboratory, Department of Medicine, University of Ioannina, GR 45110, Ioannina, Greece
| | | | - Zhuxin Li
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | | | - Claire Michelet
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | | | - Yann Perrot
- IRSN, Institut de Radioprotection et de Surete Nucleaire, 92262 Fontenay-aux-Roses, France
| | - Ivan Petrovic
- Vinca Institute of Nuclear Sciences - National Institute of the Republic of Serbia, University of Belgrade, Belgrade, Serbia
| | - Jose Ramos-Mendez
- Department of Radiation Oncology, University of California San Francisco, San Francisco 94143, CA, USA
| | - Aleksandra Ristic-Fira
- Vinca Institute of Nuclear Sciences - National Institute of the Republic of Serbia, University of Belgrade, Belgrade, Serbia
| | | | - Jan Schuemann
- Physics Division, Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA, USA
| | - Hoang N Tran
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| | - Carmen Villagrasa
- IRSN, Institut de Radioprotection et de Surete Nucleaire, 92262 Fontenay-aux-Roses, France
| | - Sebastien Incerti
- Univ. Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, F-33170 Gradignan, France
| |
Collapse
|
33
|
Sitmukhambetov S, Dinh B, Lai Y, Banigan EJ, Pan Z, Jia X, Chi Y. Development and implementation of a metaphase DNA model for ionizing radiation induced DNA damage calculation. Phys Med Biol 2022; 68:10.1088/1361-6560/aca5ea. [PMID: 36533598 PMCID: PMC9969557 DOI: 10.1088/1361-6560/aca5ea] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2022] [Accepted: 11/24/2022] [Indexed: 11/25/2022]
Abstract
Objective. To develop a metaphase chromosome model representing the complete genome of a human lymphocyte cell to support microscopic Monte Carlo (MMC) simulation-based radiation-induced DNA damage studies.Approach. We first employed coarse-grained polymer physics simulation to obtain a rod-shaped chromatid segment of 730 nm in diameter and 460 nm in height to match Hi-C data. We then voxelized the segment with a voxel size of 11 nm per side and connected the chromatid with 30 types of pre-constructed nucleosomes and 6 types of linker DNAs in base pair (bp) resolutions. Afterward, we piled different numbers of voxelized chromatid segments to create 23 pairs of chromosomes of 1-5μm long. Finally, we arranged the chromosomes at the cell metaphase plate of 5.5μm in radius to create the complete set of metaphase chromosomes. We implemented the model in gMicroMC simulation by denoting the DNA structure in a four-level hierarchical tree: nucleotide pairs, nucleosomes and linker DNAs, chromatid segments, and chromosomes. We applied the model to compute DNA damage under different radiation conditions and compared the results to those obtained with G0/G1 model and experimental measurements. We also performed uncertainty analysis for relevant simulation parameters.Main results. The chromatid segment was successfully voxelized and connected in bps resolution, containing 26.8 mega bps (Mbps) of DNA. With 466 segments, we obtained the metaphase chromosome containing 12.5 Gbps of DNA. Applying it to compute the radiation-induced DNA damage, the obtained results were self-consistent and agreed with experimental measurements. Through the parameter uncertainty study, we found that the DNA damage ratio between metaphase and G0/G1 phase models was not sensitive to the chemical simulation time. The damage was also not sensitive to the specific parameter settings in the polymer physics simulation, as long as the produced metaphase model followed a similar contact map distribution.Significance. Experimental data reveal that ionizing radiation induced DNA damage is cell cycle dependent. Yet, DNA chromosome models, except for the G0/G1 phase, are not available in the state-of-the-art MMC simulation. For the first time, we successfully built a metaphase chromosome model and implemented it into MMC simulation for radiation-induced DNA damage computation.
Collapse
Affiliation(s)
| | - Bryan Dinh
- Department of Physics, the University of Texas at Arlington, Arlington, TX 76019, USA
| | - Youfang Lai
- Department of Physics, the University of Texas at Arlington, Arlington, TX 76019, USA
| | - Edward J. Banigan
- Institute for Medical Engineering & Science and Department of Physics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
| | - Zui Pan
- Graduate Nursing, the University of Texas at Arlington, Arlington, TX 76019, USA
| | - Xun Jia
- Department of Radiation Oncology and Molecular Radiation Sciences, Johns Hopkins University, MD 21231, USA
| | - Yujie Chi
- Department of Physics, the University of Texas at Arlington, Arlington, TX 76019, USA
| |
Collapse
|
34
|
Electron tracks simulation in water: Performance comparison between GPU CPU and the EUMED grid installation. Phys Med 2022; 104:56-66. [PMID: 36368091 DOI: 10.1016/j.ejmp.2022.10.020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/06/2021] [Revised: 09/27/2022] [Accepted: 10/23/2022] [Indexed: 11/11/2022] Open
Abstract
PURPOSE We explored different technologies to minimize simulation time of the Monte-Carlo method for track generation following the Geant4-DNA processes for electrons in water. METHODS A GPU software tool is developed for electron track simulations. A similar CPU version is also developed using the same collision models. CPU simulations were carried out on a single user desktop computer and on the computing grid France Grilles using 10 and 100 computing nodes. Computing time results for CPU, GPU, and grid simulations are compared with those using Geant4-DNA processes. RESULTS The CPU simulations better performs when the number of electrons is less than 104 with 100 eV initial energy, this number decreases as the energy increases. The GPU simulations gives better results when the number of electrons is more than 104 with initial energy of 100 eV, this number decreases to 103 for electrons with 10KeV and increases back with higher energy. The use of the grid introduces an additional queuing time which slows down the overall simulation performance. Thus, the Grid gives better performance when the number of electrons is over 105 with initial energy of 10KeV, and this number decreases as the energy increases. CONCLUSIONS The CPU is best suited for small numbers of primary incident electrons. The GPU is best suited when the number of primary incident particles occupies sufficient resources on GPU card in order to get an important computing power. The grid is best suited for simulations with high number of primary incident electrons with high initial energy.
Collapse
|
35
|
Peter JS, Schuemann J, Held KD, McNamara AL. Nano-scale simulation of neuronal damage by galactic cosmic rays. Phys Med Biol 2022; 67:10.1088/1361-6560/ac95f4. [PMID: 36172820 PMCID: PMC9951267 DOI: 10.1088/1361-6560/ac95f4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2022] [Accepted: 09/28/2022] [Indexed: 11/11/2022]
Abstract
The effects of realistic, deep space radiation environments on neuronal function remain largely unexplored.In silicomodeling 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 ofin silicorealizations 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.
Collapse
Affiliation(s)
- Jonah S Peter
- Biophysics Program, Harvard University, Boston, MA 02115, United States of America
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA 02114, United States of America
| | - Jan Schuemann
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA 02114, United States of America
| | - Kathryn D Held
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA 02114, United States of America
| | - Aimee L McNamara
- Department of Radiation Oncology, Massachusetts General Hospital & Harvard Medical School, Boston, MA 02114, United States of America
| |
Collapse
|
36
|
Hilgers G, Braunroth T, Rabus H. Correlated ionisations in two spatially separated nanometric volumes within the track structure of 241Am alpha particles: comparison with Monte Carlo simulations. Radiat Phys Chem Oxf Engl 1993 2022. [DOI: 10.1016/j.radphyschem.2022.110488] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022]
|
37
|
Baiocco G, Bartzsch S, Conte V, Friedrich T, Jakob B, Tartas A, Villagrasa C, Prise KM. A matter of space: how the spatial heterogeneity in energy deposition determines the biological outcome of radiation exposure. RADIATION AND ENVIRONMENTAL BIOPHYSICS 2022; 61:545-559. [PMID: 36220965 PMCID: PMC9630194 DOI: 10.1007/s00411-022-00989-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Accepted: 08/03/2022] [Indexed: 05/10/2023]
Abstract
The outcome of the exposure of living organisms to ionizing radiation is determined by the distribution of the associated energy deposition at different spatial scales. Radiation proceeds through ionizations and excitations of hit molecules with an ~ nm spacing. Approaches such as nanodosimetry/microdosimetry and Monte Carlo track-structure simulations have been successfully adopted to investigate radiation quality effects: they allow to explore correlations between the spatial clustering of such energy depositions at the scales of DNA or chromosome domains and their biological consequences at the cellular level. Physical features alone, however, are not enough to assess the entity and complexity of radiation-induced DNA damage: this latter is the result of an interplay between radiation track structure and the spatial architecture of chromatin, and further depends on the chromatin dynamic response, affecting the activation and efficiency of the repair machinery. The heterogeneity of radiation energy depositions at the single-cell level affects the trade-off between cell inactivation and induction of viable mutations and hence influences radiation-induced carcinogenesis. In radiation therapy, where the goal is cancer cell inactivation, the delivery of a homogenous dose to the tumour has been the traditional approach in clinical practice. However, evidence is accumulating that introducing heterogeneity with spatially fractionated beams (mini- and microbeam therapy) can lead to significant advantages, particularly in sparing normal tissues. Such findings cannot be explained in merely physical terms, and their interpretation requires considering the scales at play in the underlying biological mechanisms, suggesting a systemic response to radiation.
Collapse
Affiliation(s)
- Giorgio Baiocco
- Radiation Biophysics and Radiobiology Group, Physics Department, University of Pavia, Pavia, Italy.
| | - Stefan Bartzsch
- Institute for Radiation Medicine, Helmholtz Centre Munich, Munich, Germany
- Department of Radiation Oncology, Technical University of Munich, Munich, Germany
| | - Valeria Conte
- Istituto Nazionale Di Fisica Nucleare INFN, Laboratori Nazionali Di Legnaro, Legnaro, Italy
| | - Thomas Friedrich
- Department of Biophysics, GSI Helmholtz Centre for Heavy Ion Research, Darmstadt, Germany
| | - Burkhard Jakob
- Department of Biophysics, GSI Helmholtz Centre for Heavy Ion Research, Darmstadt, Germany
| | - Adrianna Tartas
- Biomedical Physics Division, Institute of Experimental Physics, University of Warsaw, Warsaw, Poland
| | - Carmen Villagrasa
- IRSN, Institut de Radioprotection et de Sûreté Nucléaire, Fontenay aux Roses, France
| | - Kevin M Prise
- Patrick G Johnston Centre for Cancer Research, Queen's University Belfast, Belfast, UK
| |
Collapse
|
38
|
Tabbakh F, Hosmane NS, Tajudin SM, Ghorashi AH, Morshedian N. Using 157Gd doped carbon and 157GdF4 nanoparticles in proton-targeted therapy for effectiveness enhancement and thermal neutron reduction: a simulation study. Sci Rep 2022; 12:17404. [PMID: 36258012 PMCID: PMC9579128 DOI: 10.1038/s41598-022-22429-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2022] [Accepted: 10/14/2022] [Indexed: 01/10/2023] Open
Abstract
There are two major problems in proton therapy. (1) In comparison with the gamma-ray therapy, proton therapy has only ~ 10% greater biological effectiveness, and (2) the risk of the secondary neutrons in proton therapy is another unsolved problem. In this report, the increase of biological effectiveness in proton therapy has been evaluated with better performance than 11B in the presence of two proposed nanomaterials of 157GdF4 and 157Gd doped carbon with the thermal neutron reduction due to the presence of 157Gd isotope. The present study is based on the microanalysis calculations using GEANT4 Monte Carlo tool and GEANT4-DNA package for the strand breaks measurement. It was found that the proposed method will increase the effectiveness corresponding to the alpha particles by more than 100% and also, potentially will decrease the thermal neutrons fluence, significantly. Also, in this work, a discussion is presented on a significant contribution of the secondary alpha particles in total effectiveness in proton therapy.
Collapse
Affiliation(s)
- Farshid Tabbakh
- grid.459846.20000 0004 0611 7306Plasma and Nuclear Fusion Research School, Nuclear Science and Technology Research Institute, Tehran, 14155-1339 Iran
| | - Narayan S. Hosmane
- grid.261128.e0000 0000 9003 8934Department of Chemistry and Biochemistry, Northern Illinois University, DeKalb, IL 60115-2862 USA
| | - Suffian M. Tajudin
- grid.449643.80000 0000 9358 3479Faculty of Health Sciences, Universiti Sultan Zainal Abidin, Kuala Terengganu, Terengganu Malaysia
| | - Amir-Hossein Ghorashi
- grid.459846.20000 0004 0611 7306Plasma and Nuclear Fusion Research School, Nuclear Science and Technology Research Institute, Tehran, 14155-1339 Iran
| | - Nader Morshedian
- grid.459846.20000 0004 0611 7306Plasma and Nuclear Fusion Research School, Nuclear Science and Technology Research Institute, Tehran, 14155-1339 Iran
| |
Collapse
|
39
|
Quantifying Radiosensitization of PSMA-Targeted Gold Nanoparticles on Prostate Cancer Cells at Megavoltage Radiation Energies by Monte Carlo Simulation and Local Effect Model. Pharmaceutics 2022; 14:pharmaceutics14102205. [PMID: 36297640 PMCID: PMC9611822 DOI: 10.3390/pharmaceutics14102205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 10/06/2022] [Accepted: 10/11/2022] [Indexed: 11/17/2022] Open
Abstract
Active targeting gold nanoparticles (AuNPs) are a very promising avenue for cancer treatment with many publications on AuNP mediated radiosensitization at kilovoltage (kV) photon energies. However, uncertainty on the effectiveness of AuNPs under clinically relevant megavoltage (MV) radiation energies hinders the clinical translation of AuNP-assisted radiation therapy (RT) paradigm. The aim of this study was to investigate radiosensitization mediated by PSMA-targeted AuNPs irradiated by a 6 MV radiation beam at different depths to explore feasibility of AuNP-assisted prostate cancer RT under clinically relevant conditions. PSMA-targeted AuNPs (PSMA-AuNPs) were synthesized by conjugating PSMA antibodies onto PEGylated AuNPs through EDC/NHS chemistry. Confocal fluorescence microscopy was used to verify the active targeting of the developed PSMA-AuNPs. Transmission electron microscopy (TEM) was used to demonstrate the intracellular biodistribution of PSMA-AuNPs. LNCaP prostate cancer cells treated with PSMA-AuNPs were irradiated on a Varian 6 MV LINAC under varying depths (2.5 cm, 10 cm, 20 cm, 30 cm) of solid water. Clonogenic assays were carried out to determine the in vitro cell survival fractions. A Monte Carlo (MC) model developed on TOPAS platform was then employed to determine the nano-scale radial dose distribution around AuNPs, which was subsequently used to predict the radiation dose response of LNCaP cells treated with AuNPs. Two different cell models, with AuNPs located within the whole cell or only in the cytoplasm, were used to assess how the intracellular PSMA-AuNP biodistribution impacts the prostate cancer radiosensitization. Then, MC-based microdosimetry was combined with the local effect model (LEM) to calculate cell survival fraction, which was benchmarked against the in vitro clonogenic assays at different depths. In vitro clonogenic assay of LNCaP cells demonstrated the depth dependence of AuNP radiosensitization under clinical megavoltage beams, with sensitization enhancement ratio (SER) of 1.14 ± 0.03 and 1.55 ± 0.05 at 2.5 cm depth and 30 cm depth, respectively. The MC microdosimetry model showed the elevated percent of low-energy photons in the MV beams at greater depth, consequently resulting in increased dose enhancement ratio (DER) of AuNPs with depth. The AuNP-induced DER reached ~5.7 and ~8.1 at depths of 2.5 cm and 30 cm, respectively. Microdosimetry based LEM accurately predicted the cell survival under 6 MV beams at different depths, for the cell model with AuNPs placed only in the cell cytoplasm. TEM results demonstrated the distribution of PSMA-AuNPs in the cytoplasm, confirming the accuracy of MC microdosimetry based LEM with modelled AuNPs distributed within the cytoplasm. We conclude that AuNP radiosensitization can be achieved under megavoltage clinical radiotherapy energies with a dependence on tumor depth. Furthermore, the combination of Monte Carlo microdosimetry and LEM will be a valuable tool to assist with developing AuNP-aided radiotherapy paradigm and drive clinical translation.
Collapse
|
40
|
Pietrzak M, Nettelbeck H, Perrot Y, Villagrasa C, Bancer A, Bug M, Incerti S. Intercomparison of nanodosimetric distributions in nitrogen simulated with Geant4 and PTra track structure codes. Phys Med 2022; 102:103-109. [PMID: 36162229 DOI: 10.1016/j.ejmp.2022.09.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/06/2022] [Revised: 07/25/2022] [Accepted: 09/12/2022] [Indexed: 11/29/2022] Open
Abstract
To facilitate the use of Geant4-DNA for radiation transport simulations in micro- and nanodosimeters, which are physically operated with tissue-equivalent gases such as nitrogen (and propane), this work aims to extend the cross section data available in Geant4-DNA to include those of nitrogen for electron energies ranging from 1 MeV down to the ionisation threshold. To achieve this, interaction cross section data for nitrogen that have been used with the in-house PTB PTra track structure code have been implemented in the current state-of-the-art Geant4-DNA simulation toolkit. An intercomparison has been performed between the two codes to validate this implementation. To quantify the agreement between the cross section models for nitrogen adopted in PTra and those implemented in Geant4-DNA, the simulation results of both codes were analysed using three physical parameters describing the ionisation cluster size distribution (ICSD): mean ionisation cluster size, variance of the cluster size and the probability to obtain a single ionisation within the target. Statistical analysis of the results indicates that the interaction cross section models for nitrogen used in PTra (elastic scattering, impact ionisations and electronic excitations) have been successfully implemented in Geant4-DNA. In addition, simulated ICSDs were compared to those measured with the Jet Counter nanodosimeter for energies between 100 and 2000 eV. For greater energies, the ICRP data for LET and particle range were used as a reference. The modified Geant4-DNA code and data successfully passed all these benchmarks fulfilling the requirement for their public release in the next version of the Geant4 toolkit.
Collapse
Affiliation(s)
- Marcin Pietrzak
- National Centre for Nuclear Research (NCBJ), Andrzeja Sołtana 7, 05400 Otwock, Poland; European Radiation Dosimetry Group e.V. (Eurados), Ingolstädter Landstrasse 1, Neuherberg, 85764, Germany.
| | - Heidi Nettelbeck
- Physikalisch-Technische Bundesanstalt (PTB), Bundesallee 100, 38116 Braunschweig, Germany; European Radiation Dosimetry Group e.V. (Eurados), Ingolstädter Landstrasse 1, Neuherberg, 85764, Germany
| | - Yann Perrot
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), 31 Avenue de la Division Leclerc, 92260 Fontenay-Aux-Roses, France; European Radiation Dosimetry Group e.V. (Eurados), Ingolstädter Landstrasse 1, Neuherberg, 85764, Germany; Geant4-DNA Collaboration
| | - Carmen Villagrasa
- Institut de Radioprotection et de Sûreté Nucléaire (IRSN), 31 Avenue de la Division Leclerc, 92260 Fontenay-Aux-Roses, France; European Radiation Dosimetry Group e.V. (Eurados), Ingolstädter Landstrasse 1, Neuherberg, 85764, Germany; Geant4-DNA Collaboration
| | - Aleksandr Bancer
- National Centre for Nuclear Research (NCBJ), Andrzeja Sołtana 7, 05400 Otwock, Poland; European Radiation Dosimetry Group e.V. (Eurados), Ingolstädter Landstrasse 1, Neuherberg, 85764, Germany
| | - Marion Bug
- Physikalisch-Technische Bundesanstalt (PTB), Bundesallee 100, 38116 Braunschweig, Germany; European Radiation Dosimetry Group e.V. (Eurados), Ingolstädter Landstrasse 1, Neuherberg, 85764, Germany
| | - Sebastien Incerti
- Université de Bordeaux, CNRS, LP2I Bordeaux, UMR 5797, 19 Chemin du Solarium, 33170 Gradignan, France; Geant4-DNA Collaboration
| |
Collapse
|
41
|
Sarrut D, Arbor N, Baudier T, Borys D, Etxebeste A, Fuchs H, Gajewski J, Grevillot L, Jan S, Kagadis GC, Kang HG, Kirov A, Kochebina O, Krzemien W, Lomax A, Papadimitroulas P, Pommranz C, Roncali E, Rucinski A, Winterhalter C, Maigne L. The OpenGATE ecosystem for Monte Carlo simulation in medical physics. Phys Med Biol 2022; 67:184001. [PMID: 36001985 PMCID: PMC11149651 DOI: 10.1088/1361-6560/ac8c83] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/20/2022] [Accepted: 08/24/2022] [Indexed: 11/12/2022]
Abstract
This paper reviews the ecosystem of GATE, an open-source Monte Carlo toolkit for medical physics. Based on the shoulders of Geant4, the principal modules (geometry, physics, scorers) are described with brief descriptions of some key concepts (Volume, Actors, Digitizer). The main source code repositories are detailed together with the automated compilation and tests processes (Continuous Integration). We then described how the OpenGATE collaboration managed the collaborative development of about one hundred developers during almost 20 years. The impact of GATE on medical physics and cancer research is then summarized, and examples of a few key applications are given. Finally, future development perspectives are indicated.
Collapse
Affiliation(s)
- David Sarrut
- Université de Lyon; CREATIS; CNRS UMR5220; Inserm U1294; INSA-Lyon; Université Lyon 1, Léon Bérard cancer center, Lyon, France
| | - Nicolas Arbor
- Université de Strasbourg, IPHC, CNRS, UMR7178, F-67037 Strasbourg, France
| | - Thomas Baudier
- Université de Lyon; CREATIS; CNRS UMR5220; Inserm U1294; INSA-Lyon; Université Lyon 1, Léon Bérard cancer center, Lyon, France
| | - Damian Borys
- Department of Systems Biology and Engineering, Silesian University of Technology, Gliwice, Poland
| | - Ane Etxebeste
- Université de Lyon; CREATIS; CNRS UMR5220; Inserm U1294; INSA-Lyon; Université Lyon 1, Léon Bérard cancer center, Lyon, France
| | - Hermann Fuchs
- MedAustron Ion Therapy Center, Wiener Neustadt, Austria
- Medical University of Vienna, Department of Radiation Oncology, Vienna, Vienna, Währinger Gürtel 18-20, A-1090 Wien, Austria
| | - Jan Gajewski
- Institute of Nuclear Physics Polish Academy of Sciences, Krakow, Poland
| | | | - Sébastien Jan
- Université Paris-Saclay, Inserm, CNRS, CEA, Laboratoire d'Imagerie Biomédicale Multimodale (BioMaps), F-91401 Orsay, France
| | - George C Kagadis
- 3DMI Research Group, Department of Medical Physics, School of Medicine, University of Patras, Patras, Greece
| | - Han Gyu Kang
- National Institutes for Quantum Science and Technology (QST), 4-9-1 Anagawa, Inage-ku, Chiba 263-8555, Japan
| | - Assen Kirov
- Memorial Sloan Kettering Cancer, New York, NY 10021, United States of America
| | - Olga Kochebina
- Université Paris-Saclay, Inserm, CNRS, CEA, Laboratoire d'Imagerie Biomédicale Multimodale (BioMaps), F-91401 Orsay, France
| | - Wojciech Krzemien
- High Energy Physics Division, National Centre for Nuclear Research, Otwock-Świerk, Poland
- Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, S. Lojasiewicza 11, 30-348 Krakow, Poland
- Centre for Theranostics, Jagiellonian University, Kopernika 40 St, 31 501 Krakow, Poland
| | - Antony Lomax
- Center for Proton Therapy, PSI, Switzerland
- Department of Physics, ETH Zurich, Switzerland
| | | | - Christian Pommranz
- Werner Siemens Imaging Center, Department of Preclinical Imaging and Radiopharmacy, Eberhard Karls University Tuebingen, Roentgenweg 13, D-72076 Tuebingen, Germany
- Institute for Astronomy and Astrophysics, Eberhard Karls University Tuebingen, Sand 1, D-72076 Tuebingen, Germany
| | - Emilie Roncali
- University of California Davis, Departments of Biomedical Engineering and Radiology, Davis, CA 95616, United States of America
| | - Antoni Rucinski
- Institute of Nuclear Physics Polish Academy of Sciences, Krakow, Poland
| | - Carla Winterhalter
- Center for Proton Therapy, PSI, Switzerland
- Department of Physics, ETH Zurich, Switzerland
| | - Lydia Maigne
- Université Clermont Auvergne, Laboratoire de Physique de Clermont, CNRS, UMR 6533, F-63178 Aubière, France
| |
Collapse
|
42
|
Du C, Wang Y, Xue H, Gao H, Liu K, Kong X, Zhang W, Yin Y, Qiu D, Wang Y, Sun L. Research on the proximity functions of microdosimetry of low energy electrons in liquid water based on different Monte Carlo codes. Phys Med 2022; 101:120-128. [PMID: 35988482 DOI: 10.1016/j.ejmp.2022.08.006] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 02/15/2022] [Revised: 07/25/2022] [Accepted: 08/03/2022] [Indexed: 11/17/2022] Open
Abstract
PURPOSE The proximity function is an important index in microdosimetry for describing the spatial distribution of energy, which is closely related to the biological effects of organs or tissues in the target area. In this work, the impact of parameters, such as physic models, cut-off energy, and initial energy, on the proximity function are quantitated and compared. METHODS According to the track structure (TS) and condensed history (CH) low-energy electromagnetic models, this paper chooses a variety of Monte Carlo (Monte Carlo, MC) codes (Geant4-DNA, PHITS, and Penelope) to simulate the track structure of low-energy electrons in liquid water and evaluates the influence of the electron initial energy, cut-off energy, energy spectrum, and physical model factors on the differential proximity function. RESULTS The results show that the initial energy of electrons in the low-energy part (especially less than 1 keV) has a greater impact on the differential proximity function, and the choice of cut-off energy has a greater impact on the differential proximity function corresponding to small radius sites (generally less than 10 nm). The difference in the electronic energy spectrum has little effect on the result, and the proximity functions of different physics models show better consistency under large radius sites. CONCLUSIONS This work comprehensively compares the differential proximity functions under different codes by setting a variety of simulation conditions and has basic guiding significance for helping users simulate and analyze the deposition characteristics of microscale electrons according to the selection of an appropriate methodology and cut-off energy.
Collapse
Affiliation(s)
- ChuanSheng Du
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China
| | - YiDi Wang
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China
| | - HuiYuan Xue
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China
| | - Han Gao
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China
| | - Kun Liu
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China
| | - XiangHui Kong
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China
| | - WenYue Zhang
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China
| | - YuChen Yin
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China
| | - Dong Qiu
- State Key Laboratory of Radiation Medicine and Protection, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China; School of Public Health, Medical College of Soochow University, China
| | - YouYou Wang
- The Second Affiliated Hospital of Soochow University, China
| | - Liang Sun
- State Key Laboratory of Radiation Medicine and Protection, China; School of Radiation Medicine and Protection, Soochow University, China; Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, Suzhou 215123, China.
| |
Collapse
|
43
|
Ramos-Méndez J, García-García O, Domínguez-Kondo J, LaVerne JA, Schuemann J, Moreno-Barbosa E, Faddegon B. TOPAS-nBio simulation of temperature-dependent indirect DNA strand break yields. Phys Med Biol 2022; 67. [DOI: 10.1088/1361-6560/ac79f9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Accepted: 06/17/2022] [Indexed: 11/12/2022]
Abstract
Abstract
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.
Collapse
|
44
|
Perspectives of Gas Phase Ion Chemistry: Spectroscopy and Modeling. CONDENSED MATTER 2022. [DOI: 10.3390/condmat7030046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
The study of ions in the gas phase has a long history and has involved both chemists and physicists. The interplay of their competences with the use of very sophisticated commercial and/or homemade instrumentations and theoretical models has improved the knowledge of thermodynamics and kinetics of many chemical reactions, even if still many stages of these processes need to be fully understood. The new technologies and the novel free-electron laser facilities based on plasma acceleration open new opportunities to investigate the chemical reactions in some unrevealed fundamental aspects. The synchrotron light source can be put beside the FELs, and by mass spectrometric techniques and spectroscopies coupled with versatile ion sources it is possible to really change the state of the art of the ion chemistry in different areas such as atmospheric and astro chemistry, plasma chemistry, biophysics, and interstellar medium (ISM). In this manuscript we review the works performed by a joint combination of the experimental studies of ion–molecule reactions with synchrotron radiation and theoretical models adapted and developed to the experimental evidence. The review concludes with the perspectives of ion–molecule reactions by using FEL instrumentations as well as pump probe measurements and the initial attempt in the development of more realistic theoretical models for the prospective improvement of our predictive capability.
Collapse
|
45
|
Salim R, Taherparvar P. Dosimetry assessment of theranostic Auger-emitting radionuclides in a micron-sized multicellular cluster model: A Monte Carlo study using Geant4-DNA simulations. Appl Radiat Isot 2022; 188:110380. [PMID: 35868198 DOI: 10.1016/j.apradiso.2022.110380] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2022] [Revised: 06/28/2022] [Accepted: 07/07/2022] [Indexed: 11/02/2022]
Abstract
The present work is aimed at improving the multicellular dosimetry of several Auger radionuclides of interest for targeted cancer therapy, including 99mTc, 111In, 123I, 125I, and 201Tl. For this purpose, using the Geant4-DNA Monte Carlo code, a cluster of 13 similar spherical cells with a hexagonal packed arrangement was modeled, and the mean absorbed doses per unit cumulated activity (S-values) were calculated by considering two target←source configurations, cell←cell and nucleus←nucleus. The obtained ratios of cross-dose to self-dose S-value in terms of the distance between the source and target regions were evaluated and also compared to those estimated by the Medical Internal Radiation Dose (MIRD) method. Besides, the contribution of the Coster-Kronig, Auger and internal conversion electrons to the S-values was provided for each radionuclide. According to the results, it can be concluded that in contrast to self-absorption, the cross-absorption due to the Auger-emitters has not a significant role in the total energy deposition within a cell in the cluster.
Collapse
Affiliation(s)
- R Salim
- Department of Physics, Faculty of Science, University of Guilan, Rasht, Iran
| | - P Taherparvar
- Department of Physics, Faculty of Science, University of Guilan, Rasht, Iran.
| |
Collapse
|
46
|
Lillhök J, Billnert-Maróti R, Anastasiadis A. MCNP 6.2 simulations of energy deposition in low-density volumes corresponding to unit-density volumes on the nanometre level. RADIAT MEAS 2022. [DOI: 10.1016/j.radmeas.2022.106831] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
47
|
A model for Geant4-DNA to simulate ionization and excitation of liquid water by protons travelling above 100 MeV. Radiat Phys Chem Oxf Engl 1993 2022. [DOI: 10.1016/j.radphyschem.2022.110363] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
48
|
Villate JMZ, Rojas JV, Hahn MB, Puerta JA. Synthesis of 198Au nanoparticles sub 10 nm due optimization on local dose by Monte Carlo simulations for cancer treatment. J Radioanal Nucl Chem 2022. [DOI: 10.1007/s10967-022-08355-5] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
49
|
Pfuhl T, Friedrich T, Scholz M. A double-strand-break model for the relative biological effectiveness of electrons based on ionization clustering. Med Phys 2022; 49:5562-5575. [PMID: 35686448 DOI: 10.1002/mp.15796] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2022] [Revised: 05/31/2022] [Accepted: 06/01/2022] [Indexed: 11/09/2022] Open
Abstract
BACKGROUND The effectiveness of ionizing radiation regarding DNA damage induction depends on its spatial energy deposition pattern. For electrons an increased effectiveness is observed at low kinetic energies due to the enhanced density of energy deposition events at electron track ends. PURPOSE A model is presented, which enables the calculation of the double-strand-break (DSB) yield and the relative biological effectiveness (RBE) for DSB induction of electrons. METHODS The model applies the mean free path between two ionizations and the assumption that two ionizations within a certain threshold distance are necessary to potentially lead to a DSB. Next to an expression for the electron RBE according to its common definition, a local RBE is determined, which describes the electrons' local effectiveness at a defined point on their track. RESULTS This local RBE allows a better understanding of microscopic processes resulting from radiation and can be used, for instance, to describe the mean effectiveness of the mixed electron radiation field as a function of the radial distance to the center of an ion track. CONCLUSIONS The presented model reflects the experimentally observed increased effectiveness of low-energetic electrons. It will be used in a future work to improve RBE predictions for ions performed with the local effect model.
Collapse
Affiliation(s)
- Tabea Pfuhl
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | - Thomas Friedrich
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| | - Michael Scholz
- GSI Helmholtzzentrum für Schwerionenforschung, Darmstadt, Germany
| |
Collapse
|
50
|
Gerken LRH, Gogos A, Starsich FHL, David H, Gerdes ME, Schiefer H, Psoroulas S, Meer D, Plasswilm L, Weber DC, Herrmann IK. Catalytic activity imperative for nanoparticle dose enhancement in photon and proton therapy. Nat Commun 2022; 13:3248. [PMID: 35668122 PMCID: PMC9170699 DOI: 10.1038/s41467-022-30982-5] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2022] [Accepted: 05/24/2022] [Indexed: 12/19/2022] Open
Abstract
Nanoparticle-based radioenhancement is a promising strategy for extending the therapeutic ratio of radiotherapy. While (pre)clinical results are encouraging, sound mechanistic understanding of nanoparticle radioenhancement, especially the effects of nanomaterial selection and irradiation conditions, has yet to be achieved. Here, we investigate the radioenhancement mechanisms of selected metal oxide nanomaterials (including SiO2, TiO2, WO3 and HfO2), TiN and Au nanoparticles for radiotherapy utilizing photons (150 kVp and 6 MV) and 100 MeV protons. While Au nanoparticles show outstanding radioenhancement properties in kV irradiation settings, where the photoelectric effect is dominant, these properties are attenuated to baseline levels for clinically more relevant irradiation with MV photons and protons. In contrast, HfO2 nanoparticles retain some of their radioenhancement properties in MV photon and proton therapies. Interestingly, TiO2 nanoparticles, which have a comparatively low effective atomic number, show significant radioenhancement efficacies in all three irradiation settings, which can be attributed to the strong radiocatalytic activity of TiO2, leading to the formation of hydroxyl radicals, and nuclear interactions with protons. Taken together, our data enable the extraction of general design criteria for nanoparticle radioenhancers for different treatment modalities, paving the way to performance-optimized nanotherapeutics for precision radiotherapy. Nanoparticles have recently received attention in radiation therapy since they can act as radioenhancers. In this article, the authors report on the dose enhancement capabilities of a series of nanoparticles based on their metal core composition and beam characteristics, obtaining designing criteria for their optimal performance in specific radiotreatments.
Collapse
Affiliation(s)
- Lukas R H Gerken
- Nanoparticle Systems Engineering Laboratory, Institute of Energy and Process Engineering (IEPE), Department of Mechanical and Process Engineering (D-MAVT), ETH Zurich, Sonneggstrasse 3, 8092, Zurich, Switzerland.,Particles Biology Interactions Laboratory, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Lerchenfeldstrasse 5, 9014, St. Gallen, Switzerland
| | - Alexander Gogos
- Nanoparticle Systems Engineering Laboratory, Institute of Energy and Process Engineering (IEPE), Department of Mechanical and Process Engineering (D-MAVT), ETH Zurich, Sonneggstrasse 3, 8092, Zurich, Switzerland.,Particles Biology Interactions Laboratory, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Lerchenfeldstrasse 5, 9014, St. Gallen, Switzerland
| | - Fabian H L Starsich
- Nanoparticle Systems Engineering Laboratory, Institute of Energy and Process Engineering (IEPE), Department of Mechanical and Process Engineering (D-MAVT), ETH Zurich, Sonneggstrasse 3, 8092, Zurich, Switzerland.,Particles Biology Interactions Laboratory, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Lerchenfeldstrasse 5, 9014, St. Gallen, Switzerland
| | - Helena David
- Nanoparticle Systems Engineering Laboratory, Institute of Energy and Process Engineering (IEPE), Department of Mechanical and Process Engineering (D-MAVT), ETH Zurich, Sonneggstrasse 3, 8092, Zurich, Switzerland
| | - Maren E Gerdes
- Nanoparticle Systems Engineering Laboratory, Institute of Energy and Process Engineering (IEPE), Department of Mechanical and Process Engineering (D-MAVT), ETH Zurich, Sonneggstrasse 3, 8092, Zurich, Switzerland
| | - Hans Schiefer
- Department of Radiation Oncology, Cantonal Hospital St. Gallen (KSSG), Rorschacherstrasse 95, CH-9007, St. Gallen, Switzerland
| | - Serena Psoroulas
- Center for Proton Therapy, Paul Scherrer Institute (PSI), Forschungsstrasse 111, 5232, Villigen PSI, Switzerland
| | - David Meer
- Center for Proton Therapy, Paul Scherrer Institute (PSI), Forschungsstrasse 111, 5232, Villigen PSI, Switzerland
| | - Ludwig Plasswilm
- Department of Radiation Oncology, Cantonal Hospital St. Gallen (KSSG), Rorschacherstrasse 95, CH-9007, St. Gallen, Switzerland.,Department of Radiation Oncology, University Hospital Bern (Inselspital), 3010, Bern, Switzerland
| | - Damien C Weber
- Center for Proton Therapy, Paul Scherrer Institute (PSI), Forschungsstrasse 111, 5232, Villigen PSI, Switzerland.,Department of Radiation Oncology, University Hospital Bern (Inselspital), 3010, Bern, Switzerland.,Department of Radiation Oncology, University Hospital Zürich, 8091, Zürich, Switzerland
| | - Inge K Herrmann
- Nanoparticle Systems Engineering Laboratory, Institute of Energy and Process Engineering (IEPE), Department of Mechanical and Process Engineering (D-MAVT), ETH Zurich, Sonneggstrasse 3, 8092, Zurich, Switzerland. .,Particles Biology Interactions Laboratory, Department of Materials Meet Life, Swiss Federal Laboratories for Materials Science and Technology (Empa), Lerchenfeldstrasse 5, 9014, St. Gallen, Switzerland.
| |
Collapse
|