PURPOSE: Simulation of indirect damage originating from the attack of free radical species produced by ionizing radiation on biological molecules based on the independent pair approximation is investigated in this work. In addition, a new approach, relying on the independent pair approximation that is at the origin of the independent reaction time (IRT) method, is proposed in the chemical stage of Geant4-DNA. METHODS: This new approach has been designed to respect the current Geant4-DNA chemistry framework while proposing a variant IRT method. Based on the synchronous algorithm, this implementation allows us to access the information concerning the position of radicals and may make it more convenient for biological damage simulations. Estimates of the evolution of free species as well as biological hits in a segment of DNA chromatin fiber in Geant4-DNA were compared for the dynamic time step approach of the step-by-step (SBS) method, currently used in Geant4-DNA, and this newly implemented IRT. RESULTS: Results show a gain in computation time of a factor of 30 for high LET particle tracks with a better than 10% agreement on the number of DNA hits between the value obtained with the IRT method as implemented in this work and the SBS method currently available in Geant4-DNA. CONCLUSION: Offering in Geant4-DNA more efficient methods for the chemical step based on the IRT method is a task in progress. For the calculation of biological damage, information on the position of chemical species is a crucial point. This can be achieved using the method presented in this paper.
- MeSH
- chromatin genetika MeSH
- DNA * genetika MeSH
- metoda Monte Carlo MeSH
- poškození DNA * MeSH
- reakční čas MeSH
- Publikační typ
- časopisecké články MeSH
PURPOSE: The complex relationship between linear energy transfer (LET) and cellular response to radiation is not yet fully elucidated. To better characterize DNA damage after irradiations with therapeutic protons, we monitored formation and disappearance of DNA double-strand breaks (DNA DSB) as a function of LET and time. Comparisons with conventional γ-rays and high LET carbon ions were also performed. MATERIALS AND METHODS: In the present work, we performed immunofluorescence-based assay to determine the amount of DNA DSB induced by different LET values along the 62 MeV therapeutic proton Spread out Bragg peak (SOBP) in three cancer cell lines, i.e. HTB140 melanoma, MCF-7 breast adenocarcinoma and HTB177 non-small lung cancer cells. Time dependence of foci formation was followed as well. To determine irradiation positions, corresponding to the desired LET values, numerical simulations were carried out using Geant4 toolkit. We compared γ-H2AX foci persistence after irradiations with protons to that of γ-rays and carbon ions. RESULTS: With the rise of LET values along the therapeutic proton SOBP, the increase of γ-H2AX foci number is detected in the three cell lines up to the distal end of the SOBP, while there is a decrease on its distal fall-off part. With the prolonged incubation time, the number of foci gradually drops tending to attain the residual level. For the maximum number of DNA DSB, irradiation with protons attain higher level than that of γ-rays. Carbon ions produce more DNA DSB than protons but not substantially. The number of residual foci produced by γ-rays is significantly lower than that of protons and particularly carbon ions. Carbon ions do not produce considerably higher number of foci than protons, as it could be expected due to their physical properties. CONCLUSIONS: In situ visualization of γ-H2AX foci reveal creation of more lesions in the three cell lines by clinically relevant proton SOBP than γ-rays. The lack of significant differences in the number of γ-H2AX foci between the proton and carbon ion-irradiated samples suggests an increased complexity of DNA lesions and slower repair kinetics after carbon ions compared to protons. For all three irradiation types, there is no major difference between the three cell lines shortly after irradiations, while later on, the formation of residual foci starts to express the inherent nature of tested cells, therefore increasing discrepancy between them.
- MeSH
- dvouřetězcové zlomy DNA účinky záření MeSH
- lidé MeSH
- lineární přenos energie * MeSH
- nádorové buněčné linie MeSH
- oprava DNA účinky záření MeSH
- protony * MeSH
- relativní biologická účinnost MeSH
- viabilita buněk účinky záření MeSH
- vztah dávky záření a odpovědi MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
PURPOSE: The simulation of individual particle tracks and the chemical stage following water radiolysis in biological tissue is an effective means of improving our knowledge of the physico-chemical contribution to the biological effect of ionizing radiation. However, the step-by-step simulation of the reaction kinetics of radiolytic species is the most time-consuming task in Monte Carlo track-structure simulations, with long simulation times that are an impediment to research. In this work, we present the implementation of the independent reaction times (IRT) method in Geant4-DNA Monte Carlo toolkit to improve the computational efficiency of calculating G-values, defined as the number of chemical species created or lost per 100 eV of deposited energy. METHODS: The computational efficiency of IRT, as implemented, is compared to that from available Geant4-DNA step-by-step simulations for electrons, protons and alpha particles covering a wide range of linear energy transfer (LET). The accuracy of both methods is verified using published measured data from fast electron irradiations for • OH and eaq- for time-dependent G-values. For IRT, simulations in the presence of scavengers irradiated by cobalt-60 γ-ray and 2 MeV protons are compared with measured data for different scavenging capacities. In addition, a qualitative assessment comparing measured LET-dependent G-values with Geant4-DNA calculations in pure liquid water is presented. RESULTS: The IRT improved the computational efficiency by three orders of magnitude relative to the step-by-step method while differences in G-values by 3.9% at 1 μs were found. At 7 ps, • OH and eaq- yields calculated with IRT differed from recent published measured data by 5% ± 4% and 2% ± 4%, respectively. At 1 μs, differences were 9% ± 5% and 6% ± 7% for • OH and eaq- , respectively. Uncertainties are one standard deviation. Finally, G-values at different scavenging capacities and LET-dependent G-values reproduced the behavior of measurements for all radiation qualities. CONCLUSION: The comprehensive validation of the Geant4-DNA capabilities to accurately simulate the chemistry following water radiolysis is an ongoing work. The implementation presented in this work is a necessary step to facilitate performing such a task.
PURPOSE: Geant4-DNA is used to calculate S-values for different subcellular distributions of low-energy electron sources in various cell geometries. METHOD: Calculations of cellular S-values for monoenergetic electron sources with energy from 1 to 100 keV and the Auger-electron emitting radionuclides Tc-99m, In-111, and I-125 have been made using the Geant4 Monte Carlo toolkit. The Geant4-DNA low-energy extension is employed for simulating collision-by-collision the complete slowing-down of electron tracks (down to 8 eV) in liquid water, used as a surrogate of human cells. The effect of cell geometry on S-values is examined by simulating electron tracks within different cell geometries, namely, a spherical, two ellipsoidal, and an irregular shape, all having equal cellular and nuclear volumes. Algorithms for randomly sampling the volume of the nucleus, cytoplasm, surface, and whole cell for each cell phantom are presented. RESULTS: Differences between Geant4-DNA and MIRD database up to 50% were found, although, for the present radionuclides, they mostly remain below 10%. For most source-target combinations the S-values for the spherical cell geometry were found to be within 20% of those for the ellipsoidal cell geometries, with a maximum deviation of 32%. Differences between the spherical and irregular geometries are generally larger reaching 100-300%. Most sensitive to the cell geometry is the absorbed dose to the nucleus when the source is localized on the cell surface. Interestingly, two published AAPM spectra for I-125 yield noticeable differences (up to 19%) in cellular S-values. CONCLUSION: Monte Carlo simulations of cellular S-values with Geant4-DNA reveal that, for the examined radionuclides, the widely used approximation of spherical cells is reasonably accurate (within 20-30%) even for ellipsoidal geometries. For irregular cell geometries the spherical approximation should be used with caution because, as in the present example, it may lead to erroneous results for the nuclear dose for the commonly encountered situation where the source is localized to the cell surface.
- MeSH
- absorpce radiace * MeSH
- biologické modely * MeSH
- dávka záření MeSH
- elektrony MeSH
- lidé MeSH
- metoda Monte Carlo MeSH
- počítačová simulace MeSH
- radiometrie metody MeSH
- statistické modely * MeSH
- velikost buňky * MeSH
- viabilita buněk účinky léků fyziologie MeSH
- vztah dávky záření a odpovědi MeSH
- Check Tag
- lidé MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
Clustered DNA damage induced by 10, 20 and 30 MeV protons in pBR322 plasmid DNA was investigated. Besides determination of strand breaks, additional lesions were detected using base excision repair enzymes. The plasmid was irradiated in dry form, where indirect radiation effects were almost fully suppressed, and in water solution containing only minimal residual radical scavenger. Simultaneous irradiation of the plasmid DNA in the dry form and in the solution demonstrated the contribution of the indirect effect as prevalent. The damage composition slightly differed when comparing the results for liquid and dry samples. The obtained data were also subjected to analysis concerning different methodological approaches, particularly the influence of irradiation geometry, models used for calculation of strand break yields and interpretation of the strand breaks detected with the enzymes. It was shown that these parameters strongly affect the results.
- MeSH
- biologické modely MeSH
- dvouřetězcové zlomy DNA MeSH
- elektroforéza v agarovém gelu MeSH
- enzymy opravy DNA metabolismus MeSH
- lineární přenos energie MeSH
- plazmidy metabolismus účinky záření MeSH
- poškození DNA * MeSH
- protony škodlivé účinky MeSH
- roztoky MeSH
- vztah dávky záření a odpovědi MeSH
- záření gama škodlivé účinky MeSH
- Publikační typ
- časopisecké články MeSH
- práce podpořená grantem MeSH
