The Monte Carlo track structure code PARTRAC has been further improved by implementing electron scattering cross-sections for liquid water and by explicitly modelling the interaction of water radicals with DNA. The model of the genome inside a human cell nucleus in its interphase is based on the atomic coordinates of the DNA double helix with an additional volume for the water shell. The DNA helix is wound around histone complexes, and these nucleosomes are folded into chromatin fibres and further to fibre loops, which are interconnected to build chromosomes with a territorial organisation. Simulations have been performed for the irradiation of human fibroblast cells with carbon K and aluminium K ultrasoft x-rays, 220 kVp x-rays and 60Co gamma-rays. The ratio single-strand breaks to double-strand breaks (ssb/dsb) for both types of ultrasoft x-rays is lower than for gamma-rays by a factor of 2. The contributions of direct and indirect effects to strand break induction are almost independent of photon energy. Strand break patterns from indirect effects reflect differences in the susceptibility of the DNA helix to OH* attack inside the chromatin fibre. Distributions of small DNA fragments (<3 kbp) are determined by the chromatin fibre structure irrespective of whether direct or indirect effects are causing the breaks. In the calculated fragment size distributions for larger DNA fragments (>30 kbp), a substantial deviation from random breakage is found only for carbon K irradiation, and is attributed to its inhomogeneous dose distribution inside the cell nucleus. For the other radiation qualities, the results for larger fragments can be approximated by random breakage distributions calculated for a yield of dsb which is about 10% lower than the average for the whole genome. The excess of DNA fragments detected experimentally in the 8-300 kbp region after x-ray irradiation is not seen in our simulation results.

The quality of DNA damage induced by protons and alpha-particles of various linear energy transfer (LET) was studied. The aim was to single out specific lesions in the DNA molecule that might lead to biological endpoints such as inactivation. A DNA model coupled with a track structure code (MOCA-15) were used to simulate the lesions induced on the two helixes. Four categories of DNA breaks were considered: single-strand breaks (ssb), blunt-ended double-strand breaks (dsb, with no or few overlapping bases), sticky-ended double-strand breaks (with cohesive free ends of many bases), and deletions (complex lesions which involve at least two dsb within a small number of base pairs). Calculations were carried out assuming various sets of parameters characterizing the production of these different DNA breaks. No large variations in the yields of ssb and blunt- or sticky-ended dsb were found in the LET range between 10 and 200 keV/mu m. On the other hand, the yield of deletions increases up to about 100 keV/mu m and seems to reach a plateau at higher LET values. In the LET interval from 30 to 60 keV/mu m, protons proved to be more efficient than alpha-particles in inducing deletions. The induction of these complex lesions is thus dependent not simply on LET but also on the characteristics of the track structure. Comparison with RBE values for cell killing shows that this special class of dsb might play an important role in radiation-induced cell inactivation.

The yield of DNA double-strand breaks (dsb) and DNA complex lesions induced by protons and alpha-particles of various energies was simulated using a Monte Carlo track structure code (MOCA15) and a simple model of the DNA molecule. DNA breaks of different complexity were analysed. The linear energy transfer (LET) and particle-type dependence of lesions of higher complexity seems to confirm the importance of clustered damage in DNA as a relevant step leading to biological endpoints such as cell inactivation. The detailed structure of proton and alpha-particle tracks was analysed to identify the main characteristics possibly responsible for such a dependence. The role of the primary ion and of its secondary electrons in inducing dsb and complex lesions is described, showing that the relative contribution of secondary electron tracks alone in inducing clustered lesions is almost negligible at high LET, but tends to dominate below = 10 keV/micron. This is consistent with the observed similar effectiveness of low-LET fast particle radiation and sparsely ionizing radiation such as x-rays. The dependence on LET and particle type is mainly due to energy deposition events of the primary ion together with short range electrons surrounding the ion track; the yield of complex lesions due to secondary electron tracks alone is substantially LET independent. The radial distributions of the energy contributing to the induction of complex lesions were analyzed and compared with the radial distributions of energy deposition of the full tracks. The results suggest that the stochastic behaviour (i.e. cluster properties) of the energy deposition pattern within a radius of a few nanometers around the ion track plays a relevant role in determining the biological radiation effectiveness.

A new physical module for the biophysical simulation code PARTRAC has recently been developed, based on newly derived electron inelastic-scattering cross-sections in liquid water. In the present work, two modules of PARTRAC describing the production, diffusion and interaction of chemical species were developed with the specific purpose of quantifying the role of the uncertainties in the parameters controlling the early stages of liquid water radiolysis. A set of values for such parameters was identified, and time-dependent yields and frequency distributions of chemical species produced by electrons of different energies were calculated. The calculated yields were in good agreement with available data and simulations, thus confirming the reliability of the code. As the primary-electron energy decreases down to 1 keV, the *OH decay kinetics were found to get faster, reflecting variations in the spatial distribution of the initial energy depositions. In agreement with analogous works, an opposite trend was found for energies of a few hundred eV, due to the very small number of species involved. The spreading effects shown at long times by *OH frequency distributions following 1 keV irradiation were found to be essentially due to stochastic aspects of the chemical stage, whereas for 1 MeV tracks the physical and pre-chemical stages also were found to play a significant role. Relevant differences in the calculated e(aq) -yields were found by coupling the physics of PARTRAC with descriptions of the pre-chemical and chemical stages adopted in different models. This indicates a strict interrelation of the various stages, and thus a strong dependence of the parameter values on the assumptions made for the preceding and subsequent stages of the process. Although equally acceptable results can be obtained starting from different assumptions, it is necessary to keep control of such uncertainties, since they can significantly influence the modeling of radical attack on DNA and, more generally, radiobiological damage estimation. This study confirms the need for new, independently derived data on specific steps of water radiolysis, to be included in comprehensive biophysical simulation codes.