Comprehensive track-structure based evaluation of DNA damage by light ions from radiotherapy-relevant energies down to stopping

W Friedland, E Schmitt, P Kundrát, M Dingfelder, G Baiocco, S Barbieri, A Ottolenghi, W Friedland, E Schmitt, P Kundrát, M Dingfelder, G Baiocco, S Barbieri, A Ottolenghi

Abstract

Track structures and resulting DNA damage in human cells have been simulated for hydrogen, helium, carbon, nitrogen, oxygen and neon ions with 0.25-256 MeV/u energy. The needed ion interaction cross sections have been scaled from those of hydrogen; Barkas scaling formula has been refined, extending its applicability down to about 10 keV/u, and validated against established stopping power data. Linear energy transfer (LET) has been scored from energy deposits in a cell nucleus; for very low-energy ions, it has been defined locally within thin slabs. The simulations show that protons and helium ions induce more DNA damage than heavier ions do at the same LET. With increasing LET, less DNA strand breaks are formed per unit dose, but due to their clustering the yields of double-strand breaks (DSB) increase, up to saturation around 300 keV/μm. Also individual DSB tend to cluster; DSB clusters peak around 500 keV/μm, while DSB multiplicities per cluster steadily increase with LET. Remarkably similar to patterns known from cell survival studies, LET-dependencies with pronounced maxima around 100-200 keV/μm occur on nanometre scale for sites that contain one or more DSB, and on micrometre scale for megabasepair-sized DNA fragments.

Conflict of interest statement

The authors declare no competing financial interests.

Figures

Figure 1. Setup of irradiation simulations.
Figure 1. Setup of irradiation simulations.
A circular ion source (red) is located tangentially to the cell nucleus in a surrounding water shell. 5 particles per run start with a given initial energy from random positions on the source area in normal direction to it; arrows illustrate carbon ion tracks of 0.25 MeV/u initial energy (range ~6 μm). Determination of LET (Fig. 2a) is based on energy deposits including secondary electrons inside nuclear volume (orange arrows); energy deposited outside the nucleus (green parts of the arrows) is not considered. Local analysis of DNA damage by slow ions (Fig. 6) is based on DNA amount and particle stopping power (Fig. 2b) in 200-nm-slabs (thin black lines). Primary and secondary particles leaving the surrounding shell are no longer scored. Coloured structures of the cell nucleus refer to chromatin of different chromosomes.
Figure 2. LET in dependence on the…
Figure 2. LET in dependence on the initial particle energy and penetrated depth.
Panel a: LET calculated from energy deposited within the cell nucleus (symbols connected by solid lines to guide the eyes) and comparison to stopping powers from ICRU (dotted lines). Panel b: LET determined from energy deposited in 200 nm slabs perpendicular to particle transport direction (histograms) for ions with 0.5 MeV/u (thick lines) or 0.25 MeV/u initial energy (thin lines). Stopping power data from ICRU shown by symbols (0.5 MeV/u: filled squares, 0.25 MeV/u: empty squares).
Figure 3. Ion irradiation-induced DNA damage in…
Figure 3. Ion irradiation-induced DNA damage in dependence on LET.
In panel a–e, simulation results are shown by symbols connected by lines to guide the eyes; error bars in panel a–d denote two standard errors of the mean (~95% confidence intervals) of the performed simulations. Panel a: SB yields from direct effects (dashed lines), indirect effects (dash-dotted lines), and both effects (symbols and solid lines). Panel b: DSB yields, in which isolated ones as well as individual DSB in DSB clusters are scored. Panel c: Yields of DSB sites, which comprise isolated DSB and DSB clusters; a DSB cluster is scored as one DSB site. Panel d: Yields of DSB clusters; a cluster includes two or more DSB within not more than 25 bp distance. Panel e: Average multiplicity of DSB in a cluster. Panel f: DSB induction by 60Co γ-rays, He, B, N and Ne ions measured by the pulsed-field gel electrophoresis technique from DNA fragments sized between 5 kbp and 6 Mbp (symbols) and corresponding simulation results considering detectable DNA fragments (solid lines).
Figure 4. Ion irradiation-induced DNA damage in…
Figure 4. Ion irradiation-induced DNA damage in dependence on (Zeff/β)2.
Simulation results are shown by symbols connected by lines to guide the eyes; error bars denote two standard errors of the mean (~95% confidence intervals) of the performed simulations. Panel a: SSB, i.e. breaks in DSB are not included (symbols and solid lines), with prediction of MCDS (dashed line). Panel b: DSB yields (symbols and solid lines) with prediction of MCDS (dashed line). Panel c: Yields of DSB sites, which comprise isolated DSB and DSB clusters; a DSB cluster is scored as one DSB site. Panel d: Yields of DSB clusters; a cluster includes two or more DSB within not more than 25 bp distance.
Figure 5. Simulation results on the yields…
Figure 5. Simulation results on the yields of DNA fragments by single ion tracks.
Simulation results are shown by symbols, connected by lines to guide the eyes. Error bars denote two standard errors of the mean of the performed simulations. The fragment size intervals (panel a: 60–100 bp; panel b: 0.3–3 kbp; panel c: 10–100 kbp; panel d: 0.3–3 Mbp) correspond to fragments formed when an ion track induces a pair of DSB sites within a nucleosome, chromatin fiber, a chromatin fiber loop, and a chromatin domain (or a giant loop), respectively.
Figure 6. Simulated DNA damage induced in…
Figure 6. Simulated DNA damage induced in slabs of 200 nm thickness by ions with 0.25 MeV/u initial energy in dependence on the LET within that slab (Fig. 2b).
Symbols connected by solid lines to guide the eye refer to subsequent slabs; dotted lines show the total damage yields in dependence on the nucleus-averaged LET (presented in Fig. 3b–d). Panel a: SB yields from direct and indirect effects. Panel b: DSB yields, in which isolated ones as well as individual DSB in DSB clusters are scored. Panel c: Yields of DSB sites, which comprise isolated DSB and DSB clusters; a DSB cluster is scored as one DSB site. Panel d: Yields of DSB clusters; a cluster includes two or more DSB within not more than 25 bp distance. The LET and the yields of DNA damage in the first slab are indicated by arrows; when the ion slows down, the damage yields proceed along the solid lines.

References

    1. Hall E. J & Giaccia A. J. Radiobiology for the Radiologist. (Lippincott Williams & Wilkins, Philadelphia, 2006).
    1. Haberer T. et al.. The Heidelberg Ion Therapy Center. Radiother Oncol 73 Suppl 2, S186−190 (2004).
    1. Brahme A. Recent advances in light ion radiation therapy. Int J Radiat Oncol Biol Phys 58, 603–616 (2004).
    1. El Naqa I., Pater P. & Seuntjens J. Monte Carlo role in radiobiological modelling of radiotherapy outcomes. Phys Med Biol 57, R75–97 (2012).
    1. Wang C. K. The progress of radiobiological models in modern radiotherapy with emphasis on the uncertainty issue. Mutat Res 704, 175–181 (2010).
    1. Dale R. G., Jones B. & Carabe-Fernandez A. Why more needs to be known about RBE effects in modern radiotherapy. Appl Radiat Isot 67, 387–392 (2009).
    1. Paganetti H. et al.. Relative biological effectiveness (RBE) values for proton beam therapy. Int J Radiat Oncol Biol Phys 53, 407–421 (2002).
    1. Friedrich T., Scholz U., Elsasser T., Durante M. & Scholz M. Calculation of the biological effects of ion beams based on the microscopic spatial damage distribution pattern. Int J Radiat Biol 88, 103–107 (2012).
    1. Nikjoo H., Uehara S., Emfietzoglou D. & Cucinotta F. A. Track-structure codes in radiation research. Radiat Meas 41, 1052–1074 (2006).
    1. Friedland W., Dingfelder M., Kundrat P. & Jacob P. Track structures, DNA targets and radiation effects in the biophysical Monte Carlo simulation code PARTRAC. Mutat Res 711, 28–40 (2011).
    1. Alloni D., Campa A., Friedland W., Mariotti L. & Ottolenghi A. Track structure, radiation quality and initial radiobiological events: considerations based on the PARTRAC code experience. Int J Radiat Biol 88, 77–86 (2012).
    1. Friedland W. & Kundrát P. Modeling of radiation effects in cells and tissues. In Comprehensive Biomedical Physics, Vol. 9 (ed. Brahme A.) 105–142 (Elsevier, Amsterdam, 2014).
    1. Nikjoo H. & Liamsuwan T. Biophysical Basis of Ionizing Comprehensive Biomedical Physics, Vol. 9 (ed. Brahme A.) 65–104 (Elsevier, Amsterdam, 2014).
    1. Huang Y. W. et al.. Monte Carlo simulations of the relative biological effectiveness for DNA double strand breaks from 300 MeV u(−1) carbon-ion beams. Phys Med Biol 60, 5995–6012 (2015).
    1. Semenenko V. A. & Stewart R. D. Fast Monte Carlo simulation of DNA damage formed by electrons and light ions. Phys Med Biol 51, 1693–1706 (2006).
    1. Stewart R. D. et al.. Effects of radiation quality and oxygen on clustered DNA lesions and cell death. Radiat Res 176, 587–602 (2011).
    1. Stewart R. D. et al.. Rapid MCNP simulation of DNA double strand break (DSB) relative biological effectiveness (RBE) for photons, neutrons, and light ions. Phys Med Biol 60, 8249–8274 (2015).
    1. Dingfelder M., Inokuti M. & Paretzke H. G. Inelastic-collision cross sections of liquid water for interactions of energetic protons. Radiat Phys Chem 59, 255–275 (2000).
    1. Garcia-Molina R. et al.. Calculated depth-dose distributions for H+ and He+ beams in liquid water. Nucl Instrum Meth B 267, 2647–2652 (2009).
    1. Abril I., Garcia-Molina R., Denton C. D., Kyriakou I. & Emfietzoglou D. Energy loss of hydrogen- and helium-ion beams in DNA: calculations based on a realistic energy-loss function of the target. Radiat Res 175, 247–255 (2011).
    1. Nikjoo H., Uehara S. & Emfietzoglou D. Interaction of radiation with matter(CRC Press, Boca Raton, 2012).
    1. Liamsuwan T. & Nikjoo H. Cross sections for bare and dressed carbon ions in water and neon. Phys Med Biol 58, 641–672 (2013).
    1. ICRU. Stopping of ions heavier than helium. ICRU Report 73 (Oxford University Press, Oxford, 2005).
    1. Siegmund P., Schinner A. & Paul H. Errata and Addenda for ICRU Report 73, Stopping of ions heavier than helium. J ICRU 5, 1 (2009).
    1. Ziegler J. F., Ziegler M. D. & Biersack J. P. SRIM - The stopping and range of ions in matter (2010). Nucl Instrum Meth B 268, 1818–1823 (2010).
    1. Plante I. & Cucinotta F. A. Monte-Carlo Simulation of Ionizing Radiation Tracks. In Applications of Monte Carlo methods in biology, medicine and other fields of science(ed. Mode C. J.) (InTech, Rijeka, 2011).
    1. Barkas H. W. Nuclear Research Emulsions(Academic Press, New York, 1963).
    1. Schmitt E., Friedland W., Kundrat P., Dingfelder M & Ottolenghi A. Cross-section scaling for track structure simulations of low-energy ions in liquid water. Radiat Prot Dosim 166, 15–18 (2015).
    1. Baiocco G., Alloni D., Babini G., Mariotti L. & Ottolenghi A. Reaction mechanism interplay in determining the biological effectiveness of neutrons as a function of energy. Radiat Prot Dosim 166, 316–319 (2015).
    1. Baiocco G. et al.. The origin of neutron biological effectiveness as a function of energy. Sci Rep 6, 34033 (2016).
    1. Friedland W., Dingfelder M., Jacob P. & Paretzke H. G. Calculated DNA double-strand break and fragmentation yields after irradiation with He ions. Radiat Phys Chem 72, 279–286 (2005).
    1. Bernhardt P., Friedland W., Jacob P. & Paretzke H. G. Modeling of ultrasoft X-ray induced DNA damage using structured higher order DNA targets. Int J Mass Spectrom 223–224, 579–597 (2003).
    1. Paul H. Stopping Power of Matter for Ions. (2016).
    1. Dingfelder M., Hantke D., Inokuti M. & Paretzke H. G. Electron inelastic-scattering cross sections in liquid water. Radiat Phys Chem 53, 1–18 (1998).
    1. Boudaiffa B., Cloutier P., Hunting D., Huels M. A. & Sanche L. Resonant formation of DNA strand breaks by low-energy (3 to 20 eV) electrons. Science 287, 1658–1660 (2000).
    1. Alizadeh E., Orlando T. M. & Sanche L. Biomolecular Damage Induced by Ionizing Radiation: The Direct and Indirect Effects of Low-Energy Electrons on DNA. Annu Rev Phys Chem 66, 379–398 (2015).
    1. Alizadeh E. & Sanche L. Low-energy-electron interactions with DNA: approaching cellular conditions with atmospheric experiments. Eur Phys J D 68, 97–109 (2014).
    1. Kreipl M. S., Friedland W. & Paretzke H. G. Time- and space-resolved Monte Carlo study of water radiolysis for photon, electron and ion irradiation. Radiat Environ Biophys 48, 11–20 (2009).
    1. Kreipl M. S., Friedland W. & Paretzke H. G. Interaction of ion tracks in spatial and temporal proximity. Radiat Environ Biophys 48, 349–359 (2009).
    1. Friedland W. et al.. First steps towards systems radiation biology studies concerned with DNA and chromosome structure within living cells. Radiat Environ Biophys 47, 49–61 (2008).
    1. Georgakilas A. G., O’Neill P. & Stewart R. D. Induction and repair of clustered DNA lesions: what do we know so far? Radiat Res 180, 100–109 (2013).
    1. Verkhovtsev A., Surdutovich E. & Solov’yov A. V. Multiscale approach predictions for biological outcomes in ion-beam cancer therapy. Sci Rep 6, 27654 (2016).
    1. Pater P., Seuntjens J., El Naqa I. & Bernal M. A. On the consistency of Monte Carlo track structure DNA damage simulations. Med Phys 41, 121708 (2014).
    1. Charlton D. E., Nikjoo H. & Humm J. L. Calculation of initial yields of single- and double-strand breaks in cell nuclei from electrons, protons and alpha particles. Int J Radiat Biol 56, 1–19 (1989).
    1. ICRU. Linear energy transfer. ICRU Report 16 (International Commission on Radiation Units and Measurements, Washington, D.C., 1970).
    1. ICRU. Fundamental quantities and units for ionizing radiation (Revised). ICRU Report 85. J ICRU 11, 1 (2011).
    1. ICRU. Microdosimetry. ICRU Report 36 (International Commission on Radiation Units and Measurements, Bethesda, MD, 1983).
    1. Emfietzoglou D., Karava K., Papamichael G. & Moscovitch M. Monte-Carlo calculations of radial dose and restricted-let for protons in water. Radiat Prot Dosim 110, 871–879 (2004).
    1. Curtis S. B. Introduction to track structure and z*2/β2. (2016).
    1. ICRU. Stopping powers and ranges for protons and alpha particles. ICRU Report 49 (International Commission on Radiation Units and Measurements, Bethesda, MD, 1993).
    1. Friedrich T., Durante M. & Scholz M. Particle species dependence of cell survival RBE: Evident and not negligible. Acta Oncol 52, 589–603 (2013).
    1. Höglund E. DNA Fragmentation in Cultured Cells Exposed to High Linear Energy Transfer Radiation. Acta Universitatis Upsaliensis, Comprehensive Summaries of Uppsala Dissertations from the Faculty of Medicine 960, Uppsala(2000).
    1. Höglund E., Blomquist E., Carlsson J. & Stenerlow B. DNA damage induced by radiation of different linear energy transfer: initial fragmentation. Int J Radiat Biol 76, 539–547 (2000).
    1. Friedrich T., Durante M. & Scholz M. Modeling cell survival after photon irradiation based on double-strand break clustering in megabase pair chromatin loops. Radiat Res 178, 385–394 (2012).
    1. Herr L., Friedrich T., Durante M. & Scholz M. A Comparison of Kinetic Photon Cell Survival Models. Radiat Res 184, 494–508 (2015).
    1. Katz R., Dunn D. E. & Sinclair G. L. Thindown in Radiobiology. Radiat Prot Dosim 13, 281–284 (1985).
    1. Kundrát P. A semi-analytical radiobiological model may assist treatment planning in light ion radiotherapy. Phys Med Biol 52, 6813–6830 (2007).
    1. Alloni D. et al.. Energy dependence of the complexity of DNA damage induced by carbon ions. Radiat Prot Dosim 166, 86–90 (2015).
    1. Krämer M. et al.. Treatment planning for heavy-ion radiotherapy: physical beam model and dose optimization. Phys Med Biol 45, 3299–3317 (2000).
    1. Davidkova M., Kundrat P., Stepan V., Palajova Z. & Judas L. Lethal events in V79 cells irradiated by low-energy protons and correlations with distribution patterns of energy deposition, radical concentration and DNA damage. Appl Radiat Isot 67, 454–459 (2009).
    1. Friedland W. et al.. Simulation of light ion induced DNA damage patterns. Radiat Prot Dosim 122, 116–120 (2006).
    1. Kundrát P. & Stewart R. D. On the biophysical interpretation of lethal DNA lesions induced by ionising radiation. Radiat Prot Dosim 122, 169–172 (2006).
    1. Schmid T. E. et al.. Sub-micrometer 20 MeV protons or 45 MeV lithium spot irradiation enhances yields of dicentric chromosomes due to clustering of DNA double-strand breaks. Mutat Res 793, 30–40 (2015).
    1. Schmid T. E. et al.. Low LET protons focused to submicrometer shows enhanced radiobiological effectiveness. Phys Med Biol 57, 5889–5907 (2012).
    1. Nikjoo H., O’Neill P., Wilson W. E. & Goodhead D. T. Computational approach for determining the spectrum of DNA damage induced by ionizing radiation. Radiat Res 156, 577–583 (2001).
    1. Liang Y., Yang G., Liu F. & Wang Y. Monte Carlo simulation of ionizing radiation induced DNA strand breaks utilizing coarse grained high-order chromatin structures. Phys Med Biol 61, 445–460 (2016).
    1. Pater P. et al.. Proton and light ion RBE for the induction of direct DNA double strand breaks. Med Phys 43, 2131–2140 (2016).
    1. Friedland W., Jacob P., Bernhardt P., Paretzke H. G. & Dingfelder M. Simulation of DNA damage after proton irradiation. Radiat Res 159, 401–410 (2003).
    1. Alloni D., Campa A., Friedland W., Mariotti L. & Ottolenghi A. Integration of Monte Carlo simulations with PFGE experimental data yields constant RBE of 2.3 for DNA double-strand break induction by nitrogen ions between 125 and 225 keV/μm LET. Radiat Res 179, 690–697 (2013).

Source: PubMed

3
Tilaa