Artem L. Ponomarev

Image-Based Modeling Reveals Dynamic Redistribution of DNA Damage into Nuclear Sub-Domains

Several proteins involved in the response to DNA double strand breaks (DSB) form microscopically visible nuclear domains, or foci, after exposure to ionizing radiation. Radiation-induced foci (RIF) are believed to be located where DNA damage occurs. To test this assumption, we analyzed the spatial distribution of 53BP1, phosphorylated ATM, and γH2AX RIF in cells irradiated with high linear energy transfer (LET) radiation and low LET. Since energy is randomly deposited along high-LET particle paths, RIF along these paths should also be randomly distributed. The probability to induce DSB can be derived from DNA fragment data measured experimentally by pulsed-field gel electrophoresis. We used this probability in Monte Carlo simulations to predict DSB locations in synthetic nuclei geometrically described by a complete set of human chromosomes, taking into account microscope optics from real experiments. As expected, simulations produced DNA-weighted random (Poisson) distributions. In contrast, the distributions of RIF obtained as early as 5 min after exposure to high LET (1 GeV/amu Fe) were non-random. This deviation from the expected DNA-weighted random pattern can be further characterized by “relative DNA image measurements.” This novel imaging approach shows that RIF were located preferentially at the interface between high and low DNA density regions, and were more frequent than predicted in regions with lower DNA density. The same preferential nuclear location was also measured for RIF induced by 1 Gy of low-LET radiation. This deviation from random behavior was evident only 5 min after irradiation for phosphorylated ATM RIF, while γH2AX and 53BP1 RIF showed pronounced deviations up to 30 min after exposure. These data suggest that DNA damage–induced foci are restricted to certain regions of the nucleus of human epithelial cells. It is possible that DNA lesions are collected in these nuclear sub-domains for more efficient repair.

Locations of radiation-produced DNA double strand breaks along chromosomes: a stochastic cluster process formalism.

Ionizing radiation produces DNA double strand breaks (DSBs) in chromosomes. For densely ionizing radiation, the DSBs are not spaced randomly along a chromosome: recent data for size distributions of DNA fragments indicate break clustering on kbp-Mbp scales. Different DSB clusters on a chromosome are typically made by different, statistically independent, stochastically structured radiation tracks, and the average number of tracks involved can be small. We therefore model DSB positions along a chromosome as a stationary Poisson cluster process, i.e. a stochastic process consisting of secondary point processes whose locations are determined by a primary point process that is Poisson. Each secondary process represents a break cluster, typically consisting of 1-10 DSBs in a comparatively localized stochastic pattern determined by chromatin geometry and radiation track structure. Using this Poisson cluster process model, which we call the randomly located clusters (RLC) formalism, theorems are derived for how the DNA fragment-size distribution depends on radiation dose. The RLC dose-response relations become non-linear when the dose becomes so high that DSB clusters from different tracks overlap or adjoin closely. The RLC formalism generalizes previous models, fits current data adequately and facilitates mechanistically based extrapolations from high-dose experiments to the much lower doses of interest for most applications.

Chromatin loops are responsible for higher counts of small DNA fragments induced by high-LET radiation, while chromosomal domains do not affect the fragment sizes

Purpose: To apply a polymer model of DNA damage induced by high-LET (linear energy transfer) radiation and determine the influence of chromosomal domains and loops on fragment length distribution. Materials and methods: The yields of DSB (double-strand breaks) induced by high-LET radiation were calculated using a track structure model along with a polymer model of DNA packed in the cell nucleus. The cell nucleus was constructed to include the chromosomal domains and chromatin loops. The latter were generated by the random walk method. Results and conclusions: We present data for DSB yields per track per cell, DNA fragment sizes, the radial distribution of DSB with respect to the track center, and the distribution of 0, 1, 2, and more DSB from a single particle. Calculations were carried out for a range of particles including He (40 keV/μm), N (225 keV/μm), and Fe ions (150 keV/μm). Situations relevant to PFGE (pulsed-field gel electrophoresis) and microbeam experiments with direct irradiation of the cell nucleus were simulated to demonstrate the applicability of the model. Data show that chromosomal domains do not have a significant influence on fragment-size distribution, while the presence of DNA loops increases the frequencies of smaller fragments by nearly 30% for fragment sizes in the range from 2 kbp (bp = base pair) to 20 kbp.

NUCLEAR INTERACTIONS IN HEAVY ION TRANSPORT AND EVENT-BASED RISK MODELS

The physical description of the passage of heavy ions in tissue and shielding materials is of interest in radiobiology, cancer therapy and space exploration, including a human mission to Mars. Galactic cosmic rays (GCRs) consist of a large number of ion types and energies. Energy loss processes occur continuously along the path of heavy ions and are well described by the linear energy transfer (LET), straggling and multiple scattering algorithms. Nuclear interactions lead to much larger energy deposition than atomic-molecular collisions and alter the composition of heavy ion beams while producing secondary nuclei often in high multiplicity events. The major nuclear interaction processes of importance for describing heavy ion beams was reviewed, including nuclear fragmentation, elastic scattering and knockout-cascade processes. The quantum multiple scattering fragmentation model is shown to be in excellent agreement with available experimental data for nuclear fragmentation cross sections and is studied for application to thick target experiments. A new computer model, which was developed for the description of biophysical events from heavy ion beams at the NASA Space Radiation Laboratory (NSRL), called the GCR Event Risk-Based Model (GERMcode) is described.

3D visualisation of the stochastic patterns of the radial dose in nano-volumes by a Monte Carlo simulation of HZE ion track structure

The description of energy deposition by high charge and energy (HZE) nuclei is of importance for space radiation risk assessment and due to their use in hadrontherapy. Such ions deposit a large fraction of their energy within the so-called core of the track and a smaller proportion in the penumbra (or track periphery). We study the stochastic patterns of the radial dependence of energy deposition using Monte Carlo track structure codes RITRACKS and RETRACKS, that were used to simulate HZE tracks and calculate energy deposition in voxels of 40 nm. The simulation of a (56)Fe(26+) ion of 1 GeV u(-1) revealed zones of high-energy deposition which maybe found as far as a few millimetres away from the track core in some simulations. The calculation also showed that ∼43 % of the energy was deposited in the penumbra. These 3D stochastic simulations combined with a visualisation interface are a powerful tool for biophysicists which may be used to study radiation-induced biological effects such as double strand breaks and oxidative damage and the subsequent cellular and tissue damage processing and signalling.

Stochastic properties of radiation-induced DSB: DSB distributions in large scale chromatin loops, the HPRT gene and within the visible volumes of DNA repair foci

Purpose: We computed probabilities to have multiple double-strand breaks (DSB), which are produced in DNA on a regional scale, and not in close vicinity, in volumes matching the size of DNA damage foci, of a large chromatin loop, and in the physical volume of DNA containing the HPRT (human hypoxanthine phosphoribosyltransferase) locus. Materials and methods: The model is based on a Monte Carlo description of DSB formation by heavy ions in the spatial context of the entire human genome contained within the cell nucleus, as well as at the gene sequence level. Results: We showed that a finite physical volume corresponding to a visible DNA repair focus, believed to be associated with one DSB, can contain multiple DSB due to heavy ion track structure and the DNA supercoiled topography. A corrective distribution was introduced, which was a conditional probability to have excess DSB in a focus volume, given that there was already one present. The corrective distribution was calculated for 19.5 MeV/amu N ions, 3.77 MeV/amu α-particles, 1000 MeV/amu Fe ions, and X-rays. The corrected initial DSB yield from the experimental data on DNA repair foci was calculated. The DSB yield based on the corrective function converts the focus yield into the DSB yield, which is comparable with the DSB yield based on the earlier PFGE experiments. The distribution of DSB within the physical limits of the HPRT gene was analyzed by a similar method as well. Conclusion: This corrective procedure shows the applicability of the model and empowers the researcher with a tool to better analyze focus statistics. The model enables researchers to analyze the DSB yield based on focus statistics in real experimental situations that lack one-to-one focus-to-DSB correspondance.

Polymer chromosome models and Monte Carlo simulations of radiation breaking DNA.

MOTIVATION Chromatin breakage by ionizing radiation is relevant to studies of carcinogenesis, tumor radiotherapy, biodosimetry and molecular biology. This article focuses on computer analysis of chromosome irradiation in mammlian cells. METHODS Polymer physics and Monte Carlo numerical methods are used to develop a coarse-grained computational approach. Chromatin is modeled as a random walk on a cubic lattice, and the radiation tracks hitting the chromatin are modeled as straight lines hitting lattice sites. Each track can make a cluster of DSBs on a chromosome. RESULTS The results obtained replace conjectured DNA fragment-size distribution functions in the recently developed RLC formalism by more mechanistically motivated distributions. The discrete lattice algorithm reproduces features of current radiation experiments relevant to chromatin on large scales. It approximates the continuous formalism and experimental data with adequate precision. It was also found that assuming either fixed chromatin with correlations among different clusters of DSBs or moving chromatin with no such correlations gives virtually identical numerical predictions.

Extrapolation of the DNA Fragment-Size Distribution after High-Dose Irradiation to Predict Effects at Low Doses

Abstract Ponomarev, A. L., Cucinotta, F. A., Sachs, R. K., Brenner, D. J. and Peterson, L. E. Extrapolation of the DNA Fragment-Size Distribution after High-Dose Irradiation to Predict Effects at Low Doses. Radiat. Res. 156, 594–597 (2001). The patterns of DSBs induced in the genome are different for sparsely and densely ionizing radiations: In the former case, the patterns are well described by a random-breakage model; in the latter, a more sophisticated tool is needed. We used a Monte Carlo algorithm with a random-walk geometry of chromatin, and a track structure defined by the radial distribution of energy deposition from an incident ion, to fit the PFGE data for fragment-size distribution after high-dose irradiation. These fits determined the unknown parameters of the model, enabling the extrapolation of data for high-dose irradiation to the low doses that are relevant for NASA space radiation research. The randomly-located-clusters formalism was used to speed the simulations. It was shown that only one adjustable parameter, Q, the track efficiency parameter, was necessary to predict DNA fragment sizes for wide ranges of doses. This parameter was determined for a variety of radiations and LETs and was used to predict the DSB patterns at the HPRT locus of the human X chromosome after low-dose irradiation. It was found that high-LET radiation would be more likely than low-LET radiation to induce additional DSBs within the HPRT gene if this gene already contained one DSB.

Analysis of Flow Cytometry DNA Damage Response Protein Activation Kinetics after Exposure to X Rays and High-Energy Iron Nuclei

Abstract We developed a mathematical method to analyze flow cytometry data to describe the kinetics of γ-H2AX and pATF2 phosphorylation in normal human fibroblast cells after exposure to various qualities of low-dose radiation. Previously reported flow cytometry kinetics for these DSB repair phospho-proteins revealed that distributions of intensity were highly skewed, severely limiting the detection of differences in the very low-dose range. Distributional analysis revealed significant differences between control and low-dose samples when distributions were compared using the Kolmogorov-Smirnov test. Differences in radiation quality were found in the distribution shapes and when a nonlinear model was used to relate dose and time to the decay of the mean ratio of phospho-protein intensities of irradiated samples to controls. We analyzed cell cycle phase- and radiation quality-dependent characteristic repair times and residual phospho-protein levels with these methods. Characteristic repair times for γ-H2AX were higher after exposure to iron nuclei compared to X rays in G1 cells and in S/G2 cells. The RBE in G1 cells for iron nuclei relative to X rays for γ-H2AX was 2.1 ± 0.6 and 5.0 ± 3.5 at 2 and 24 h after irradiation, respectively. For pATF2, a saturation effect was observed with reduced expression at high doses, especially for iron nuclei, with much slower characteristic repair times (>7 h) compared to X rays. RBEs for pATF2 were 0.7 ± 0.1 and 1.7 ± 0.5 at 2 and 24 h, respectively. Significant differences in γ-H2AX and pATF2 levels when irradiated samples were compared to controls were noted even at the lowest dose analyzed (0.05 Gy). These results show that mathematical models can be applied to flow cytometry data to identify important and subtle differences after exposure to various qualities of low-dose radiation.

Computational Model of Chromosome Aberration Yield Induced by High- and Low-LET Radiation Exposures

We present a computational model for calculating the yield of radiation-induced chromosomal aberrations in human cells based on a stochastic Monte Carlo approach and calibrated using the relative frequencies and distributions of chromosomal aberrations reported in the literature. A previously developed DNA-fragmentation model for high- and low-LET radiation called the NASARadiationTrackImage model was enhanced to simulate a stochastic process of the formation of chromosomal aberrations from DNA fragments. The current version of the model gives predictions of the yields and sizes of translocations, dicentrics, rings, and more complex-type aberrations formed in the G0/G1 cell cycle phase during the first cell division after irradiation. As the model can predict smaller-sized deletions and rings (<3 Mbp) that are below the resolution limits of current cytogenetic analysis techniques, we present predictions of hypothesized small deletions that may be produced as a byproduct of properly repaired DNA double-strand breaks (DSB) by nonhomologous end-joining. Additionally, the model was used to scale chromosomal exchanges in two or three chromosomes that were obtained from whole-chromosome FISH painting analysis techniques to whole-genome equivalent values.