Monte Carlo simulation of DNA damage by low LET radiation using inhomogeneous higher order DNA targets

Rapid MCNP simulation of DNA double strand break (DSB) relative biological effectiveness (RBE) for photons, neutrons, and light ions

To account for particle interactions in the extracellular (physical) environment, information from the cell-level Monte Carlo damage simulation (MCDS) for DNA double strand break (DSB) induction has been integrated into the general purpose Monte Carlo N-particle (MCNP) radiation transport code system. The effort to integrate these models is motivated by the need for a computationally efficient model to accurately predict particle relative biological effectiveness (RBE) in cell cultures and in vivo. To illustrate the approach and highlight the impact of the larger scale physical environment (e.g. establishing charged particle equilibrium), we examined the RBE for DSB induction (RBEDSB) of x-rays, (137)Cs γ-rays, neutrons and light ions relative to γ-rays from (60)Co in monolayer cell cultures at various depths in water. Under normoxic conditions, we found that (137)Cs γ-rays are about 1.7% more effective at creating DSB than γ-rays from (60)Co (RBEDSB  =  1.017) whereas 60-250 kV x-rays are 1.1 to 1.25 times more efficient at creating DSB than (60)Co. Under anoxic conditions, kV x-rays may have an RBEDSB up to 1.51 times as large as (60)Co γ-rays. Fission neutrons passing through monolayer cell cultures have an RBEDSB that ranges from 2.6 to 3.0 in normoxic cells, but may be as large as 9.93 for anoxic cells. For proton pencil beams, Monte Carlo simulations suggest an RBEDSB of about 1.2 at the tip of the Bragg peak and up to 1.6 a few mm beyond the Bragg peak. Bragg peak RBEDSB increases with decreasing oxygen concentration, which may create opportunities to apply proton dose painting to help address tumor hypoxia. Modeling of the particle RBE for DSB induction across multiple physical and biological scales has the potential to aid in the interpretation of laboratory experiments and provide useful information to advance the safety and effectiveness of hadron therapy in the treatment of cancer.

Spectrum of Radiation-Induced Clustered Non-DSB Damage - A Monte Carlo Track Structure Modeling and Calculations

The aim of this report is to present the spectrum of initial radiation-induced cellular DNA damage [with particular focus on non-double-strand break (DSB) damage] generated by computer simulations. The radiation types modeled in this study were monoenergetic electrons (100 eV-1.5 keV), ultrasoft X-ray photons Ck, AlK and TiK, as well as some selected ions including 3.2 MeV/u proton; 0.74 and 2.4 MeV/u helium ions; 29 MeV/u nitrogen ions and 950 MeV/u iron ions. Monte Carlo track structure methods were used to simulate damage induction by these radiation types in a cell-mimetic condition from a single-track action. The simulations took into account the action of direct energy deposition events and the reaction of hydroxyl radicals on atomistic linear B-DNA segments of a few helical turns including the water of hydration. Our results permitted the following conclusions: a. The absolute levels of different types of damage [base damage, simple and complex single-strand breaks (SSBs) and DSBs] vary depending on the radiation type; b. Within each damage class, the relative proportions of simple and complex damage vary with radiation type, the latter being higher with high-LET radiations; c. Overall, for both low- and high-LET radiations, the ratios of the yields of base damage to SSBs are similar, being about 3.0 ± 0.2; d. Base damage contributes more to the complexity of both SSBs and DSBs, than additional SSB damage and this is true for both low- and high-LET radiations; and e. The average SSB/DSB ratio for low-LET radiations is about 18, which is about 5 times higher than that for high-LET radiations. The hypothesis that clustered DNA damage is more difficult for cells to repair has gained currency among radiobiologists. However, as yet, there is no direct in vivo experimental method to validate the dependence of kinetics of DNA repair on DNA damage complexity (both DSB and non-DSB types). The data on the detailed spectrum of DNA damage presented here, in particular the non-DSB type, provide a good basis for testing mechanistic models of DNA repair kinetics such as base excision repair.

Track structures, DNA targets and radiation effects in the biophysical Monte Carlo simulation code PARTRAC

This review describes the PARTRAC suite of comprehensive Monte Carlo simulation tools for calculations of track structures of a variety of ionizing radiation qualities and their biological effects. A multi-scale target model characterizes essential structures of the whole genomic DNA within human fibroblasts and lymphocytes in atomic resolution. Calculation methods and essential results are recapitulated regarding the physical, physico-chemical and chemical stage of track structure development of radiation damage induction. Recent model extension towards DNA repair processes extends the time dimension by about 12 orders of magnitude and paves the way for superior predictions of radiation risks.

Induction of strand breaks by low-energy electrons (8-68 eV) in a self-assembled monolayer of oligonucleotides: effective cross sections and attenuation lengths

Self-assembled monolayers of 5'-32P-labeled 3'-thiolated oligonucleotides chemisorbed on gold were bombarded by low-energy electrons (LEE) of 8-68 eV. Shorter 5'-32P-oligonucleotides produced by LEE-induced strand breaks were separated with denaturing polyacrylamide gel electrophoresis and quantified by phosphor imaging. The yields of short oligonucleotides (y) decrease exponentially with their length (n), following the equation y=ae-bn, where a and b are constants, which are related to the average effective cross section per nucleotide for DNA strand break (sigmaeff) and the attenuation length (AL=1b) of LEE, respectively. The AL decreases with LEE energies from 2.5+/-0.6 nm at 8 eV to 0.8+/-0.1 nm at 68 eV, whereas sigmaeff increases from (3+/-1)x10(-18) to (5.1+/-1.6)x10(-17) cm2 within the same energy range. The energy dependence of sigmaeff shows a resonance peak of (2.8+/-0.9)x10(-17) cm2 at 18 eV superimposed on a monotonically rising curve. Transient electron attachment to a sigma* anion state of the deoxyribose group, followed by dipolar dissociation into H- and the corresponding positive-ion radical, leading to C-O bond cleavage, is proposed to account for this maximum.

Clusters of DNA Double-Strand Breaks Induced by Different Doses of Nitrogen Ions for Various LETs: Experimental Measurements and Theoretical Analyses

The yields and clustering of DNA double-strand breaks (DSBs) were investigated in normal human skin fibroblasts exposed to gamma rays or to a wide range of doses of nitrogen ions with various linear energy transfers (LETs). Data obtained by pulsed-field gel electrophoresis on the dose and LET dependence of DNA fragmentation were analyzed with the randomly located clusters (RLC) formalism. The formalism considers stochastic clustering of DSBs along a chromosome due to chromatin structure, particle track structure, and multitrack action. The relative biological effectiveness (RBE) for the total DSB yield did not depend strongly on LET, but particles with higher LET produced higher fractions of small DNA fragments, corresponding in the formalism to an increase in the average number of DSBs per DSB cluster. The results are consistent with the idea that DSB clustering along chromosomes is what leads to large RBEs of high-LET radiations for major biological end points. At a given dose, large fragments are less affected by the variability in LET than small fragments, suggesting that the two free ends in large fragments are often produced by two different tracks. The formalism successfully described an extra increase in small DNA fragments as dose increases and a related decrease in large fragments, mainly due to interlacing of DSB clusters produced along a chromosome by different tracks, since interlacing cuts larger DNA fragments into smaller ones.

A model for the induction of DNA damages and their evolution into cell clonogenic inactivation

The dependence of the initial production of DNA damages on radiation quality was examined by using a proposed new model on the basis of target theory. For the estimation of DNA damage-production by different radiation qualities, five possible modes of radiation action, including both direct and indirect effects, were assumed inside a target the molecular structure of which was defined to consist of 10 base-pairs of DNA surrounded by water molecules. The induction of DNA damage was modeled on the basis of comparisons between the primary ionization mean free path and the distance between pairs of ionized atoms, such distance being characteristic on the mode of radiation action. The OH radicals per average energy to produce an ion pair on the nanosecond time scale was estimated and used for indirect action. Assuming a relation between estimated yields of DNA damages and experimental inactivation cross sections for AT-cells, the present model enabled the quantitative reproduction of experimental results for AT-cell killing under aerobic or hypoxic conditions. The results suggest a higher order organization of DNA in a way that there will be at least two types of water environment, one filling half the space surrounding DNA with a depth of 3.7-4.3 nm and the other filling all space with a depth 4.6-4.9 nm.

Evaluation of lesion clustering in irradiated plasmid DNA

PURPOSE

To measure the yield of DNA strand breaks and clustered lesions in plasmid DNA irradiated with protons, helium nuclei, and y-rays.

MATERIALS AND METHODS

Plasmid DNA was irradiated with 1.03, 19.3 and 249 MeV protons (linear energy transfer = 25.5, 2.7, and 0.39 keV microm(-1) respectively), 26 MeV helium nuclei (25.5 keV microm) and gamma-rays (137Cs or 60Co) in phosphate buffer containing 2 mM or 200 mM glycerol. Single-and double-strand breaks (SSB and DSB) were measured by gel electrophoresis, and clustered lesions containing base lesions were quantified by converting them into irreparable DSB in transformed bacteria.

RESULTS

For protons, SSB yield decreased with increasing LET (linear energy transfer). The yield of DSB and all clustered lesions seemed to reach a minimum around 3 keV microm(-1). There was a higher yield of SSB, DSB and total clustered lesions for protons compared to helium nuclei at 25.5 keV microm(-1). A difference in the yields between 137Cs and 60Co gamma-rays was also observed, especially for SSB.

CONCLUSION

In this work we have demonstrated the complex LET dependence of clustered-lesion yields, governed by interplay of the radical recombination and change in track structure. As expected, there was also a significant difference in clustered lesion yields between various radiation fields, having the same or similar LET values, but differing in nanometric track structure.