ALBERT

Description: Purpose: To simulate radiation-induced chromosome aberrations in mammalian cells (e.g., rings, translocations, and dicentrics) and to calculate their frequency distributions following exposure to DNA double strand breaks (DSBs) produced by high-LET ions. Methods: The interphase genome was assumed to be comprised of a collection of 2 kbp rigid-block monomers following the random-walk geometry. Additional details for the modeling of chromosomal structure, such as chromosomal domains and chromosomal loops, were included. A radial energy profile for heavy ion tracks was used to simulate the high-LET pattern of induced DSBs. The induced DSB pattern depended on the ion charge and kinetic energy, but always corresponded to the DSB yield of 25 DSBs/cell/Gy. The sum of all energy contributions from Poisson-distributed particle tracks was taken to account for all possible one-track and multi-track effects. The relevant output of the model was DNA fragments produced by DSBs. The DSBs, or breakpoints, were defined by (x, y, z, l) positions, where x, y, z were the Euclidian coordinates of a DSB, and where l was the relative position along the genome. Results: The code was used to carry out Monte Carlo simulations for DSB rejoinings at low doses. The resulting fragments were analyzed to estimate the frequencies of specific types of chromosomal aberrations. Histograms for relative frequencies of chromosomal aberrations and P.D.F.s (probability density functions) of a given aberration type were produced. The relative frequency of dicentrics to rings was compared to empirical data to calibrate rejoining probabilities. Of particular interest was the predicted distribution of ring sizes, irrespective of their frequencies relative to other aberrations. Simulated ring sizes were . 4 kbp, which are far too small to be observed experimentally (i.e., by microscopy) but which, nevertheless, are conjectured to exist. Other aberrations, for example, inversions, translocations, as well as multi-centrics were also recorded. Conclusion: High-LET DNA damage affects the frequencies of chromosomal aberrations. The ratio of rings to dicentrics is correct for the genomic size cut-offs corresponding to available experimental data. The present work predicts a relative abundance of small rings following irradiation by heavy ions.

Description: A three-dimensional tissue model was incorporated into a new Monte Carlo algorithm that simulates passage of heavy ions in a tissue box . The tissue box was given as a realistic model of tissue based on confocal microscopy images. The action of heavy ions on the cellular matrix for 2- or 3-dimensional cases was simulated. Cells were modeled as a cell culture monolayer in one example, where the data were taken directly from microscopy (2-d cell matrix), and as a multi-layer obtained from confocal microscopy (3-d case). Image segmentation was used to identify cells with precise areas/volumes in an irradiated cell culture monolayer, and slices of tissue with many cell layers. The cells were then inserted into the model box of the simulated physical space pixel by pixel. In the case of modeled tissues (3-d), the tissue box had periodic boundary conditions imposed, which extrapolates the technique to macroscopic volumes of tissue. For the real tissue (3-d), specific spatial patterns for cell apoptosis and necrosis are expected. The cell patterns were modeled based on action cross sections for apoptosis and necrosis estimated from current experimental data. A spatial correlation function indicating a higher spatial concentration of damaged cells from heavy ions relative to the low-LET radiation cell damage pattern is presented. The spatial correlation effects among necrotic cells can help studying microlesions in organs, and probable effects of directionality of heavy ion radiation on epithelium and endothelium.

Description: We present new results obtained with our model based on the track structure and chromatin geometry that predicts the DSB spatial and genomic distributions in a cell nucleus with the full genome represented. The model generates stochastic patterns of DSBs in the physical space of the nucleus filled with the realistic configuration of human chromosomes. The model was re-used to find the distribution of DSBs in a physical volume corresponding to a visible phosphorylation focus believed to be associated with a DSB. The data shows whether there must more than one DSB per foci due to finite size of the visible focus, even if a single DSB is radiochemically responsible for the phosphorylation of DNA in its vicinity. The same model can predict patterns of closely located DSBs in a given gene, or in a DNA loop, one of the large-scale chromatin structures. We demonstrated for the example of the HPRT gene, how different sorts of radiation lead to proximity effect in DSB locations, which is important for modeling gene deletions. The spectrum of intron deletions and total gene deletions was simulated for the HPRT gene. The same proximity effect of DSBs in a loop can hinder DSB restitutions, as parts of the loop between DSBs is deleted with a higher likelihood. The distributions of DSBs and deletions of DNA in a loop are presented.

Description: A new Monte Carlo algorithm was developed that can model passage of heavy ions in a tissue, and their action on the cellular matrix for 2- or 3-dimensional cases. The build-up of secondaries such as projectile fragments, target fragments, other light fragments, and delta-rays was simulated. Cells were modeled as a cell culture monolayer in one example, where the data were taken directly from microscopy (2-d cell matrix). A simple model of tissue was given as abstract spheres with close approximation to real cell geometries (3-d cell matrix), as well as a realistic model of tissue was proposed based on microscopy images. Image segmentation was used to identify cells in an irradiated cell culture monolayer, or slices of tissue. The cells were then inserted into the model box pixel by pixel. In the case of cell monolayers (2-d), the image size may exceed the modeled box size. Such image was is moved with respect to the box in order to sample as many cells as possible. In the case of the simple tissue (3-d), the tissue box is modeled with periodic boundary conditions, which extrapolate the technique to macroscopic volumes of tissue. For real tissue, specific spatial patterns for cell apoptosis and necrosis are expected. The cell patterns were modeled based on action cross sections for apoptosis and necrosis estimated based on BNL data, and other experimental data.

Description: The code RITRACKS (Relativistic Ion Tracks) was developed to simulate detailed stochastic radiation track structures of ions of different types and energies. Many new capabilities were added to the code during the recent years. Several options were added to specify the times at which the tracks appear in the irradiated volume, allowing the simulation of dose-rate effects. The code has been used to simulate energy deposition in several targets: spherical, ellipsoidal and cylindrical. More recently, density changes as well as a spherical shell were implemented for spherical targets, in order to simulate energy deposition in walled tissue equivalent proportional counters. RITRACKS is used as a part of the new program BDSTracks (Biological Damage by Stochastic Tracks) to simulate several types of chromosome aberrations in various irradiation conditions. The simulation of damage to various DNA structures (linear and chromatin fiber) by direct and indirect effects has been improved and is ongoing. Many improvements were also made to the graphic user interface (GUI), including the addition of several labels allowing changes of units. A new GUI has been added to display the electron ejection vectors. The parallel calculation capabilities, notably the pre- and post-simulation processing on Windows and Linux machines have been reviewed to make them more portable between different systems. The calculation part is currently maintained in an Atlassian Stash repository for code tracking and possibly future collaboration.

Description: Exploration missions to Mars and other destinations raise many questions about the health of astronauts. The continuous exposure of astronauts to galactic cosmic rays is one of the main concerns for long-term missions. Cosmic ionizing radiations are composed of different ions of various charges and energies notably, highly charged energy (HZE) particles. The HZE particles have been shown to be more carcinogenic than low-LET radiation, suggesting the severity of chromosomal aberrations induced by HZE particles is one possible explanation. However, most mathematical models predicting cell death and mutation frequency are based on directly fitting various HZE dose response and are in essence empirical approaches. In this work, we assume a simple biological mechanism to model DNA repair and use it to simultaneously explain the low- and high-LET response using the exact same fitting parameters. Our work shows that the geometrical position of DNA repair along tracks of heavy ions are sufficient to explain why high-LET particles can induce more death and mutations. Our model is based on assuming DNA double strand breaks (DSBs) are repaired within repair domain, and that any DSBs located within the same repair domain cluster into one repair unit, facilitating chromosomal rearrangements and increasing the probability of cell death. We introduced this model in 2014 using simplified microdosimetry profiles to predict cell death. In this work, we collaborated with NASA Johnson Space Center to generate more accurate microdosimetry profiles derived by Monte Carlo techniques, taking into account track structure of HZE particles and simulating DSBs in realistic cell geometry. We simulated 224 data points (D, A, Z, E) with the BDSTRACKS model, leading to a large coverage of LET from ~10 to 2,400 keV/m. This model was used to generate theoretical RBE for various particles and energies for both cell death and mutation frequencies. The RBE LET dependence is in agreement with experimental data known in human and murine cells. It suggests that cell shape and its orientation with respect to the HZE particle beam can modify the biological response to radiation. Such discovery will be tested experimentally and, if proven accurate, will be another strong supporting evidence for DNA repair domains and their critical role in interpreting cosmic radiation sensitivity.

Description: The formation of DNA double-strand breaks (DSBs) and chromosome aberrations is an important consequence of ionizing radiation. To simulate DNA double-strand breaks and the formation of chromosome aberrations, we have recently merged the codes RITRACKS (Relativistic Ion Tracks) and NASARTI (NASA Radiation Track Image). The program RITRACKS is a stochastic code developed to simulate detailed event-by-event radiation track structure: [1] This code is used to calculate the dose in voxels of 20 nm, in a volume containing simulated chromosomes, [2] The number of tracks in the volume is calculated for each simulation by sampling a Poisson distribution, with the distribution parameter obtained from the irradiation dose, ion type and energy. The program NASARTI generates the chromosomes present in a cell nucleus by random walks of 20 nm, corresponding to the size of the dose voxels, [3] The generated chromosomes are located within domains which may intertwine, and [4] Each segment of the random walks corresponds to approx. 2,000 DNA base pairs. NASARTI uses pre-calculated dose at each voxel to calculate the probability of DNA damage at each random walk segment. Using the location of double-strand breaks, possible rejoining between damaged segments is evaluated. This yields various types of chromosomes aberrations, including deletions, inversions, exchanges, etc. By performing the calculations using various types of radiations, it will be possible to obtain relative biological effectiveness (RBE) values for several types of chromosome aberrations.

Description: We have developed a model that can simulate different types of radiation induced chromosomal aberrations (CA's) and can provide predictions on the frequency and size of chromosomes with terminal deletions. Chromosomes with terminal deletions lack telomeres and this can elicit sister chromatid unions and the prolonged breakage/fusion/bridge (B/F/B) cycles that have been observed in mammalian tumors. The loss of a single telomere has been shown to cause extensive genomic instability through the B/F/B cycle process. Our model uses a stochastic process of DNA broken end joining, in which a realistic spectrum of CA's is created from improperly joined DNA free ends formed by DNA double strand breaks (DSBs). The distribution of the DNA free ends is given by a mechanistic model that takes into account the chromatin structure and track structure for high-LET radiation. The model allows for DSB clustering from high-LET radiation and simulates the formation of CA's in stages that correspond to the actual time after radiation exposure. The time scale for CA formation is derived from experimental data on DSB repair kinetics. At any given time a nucleus may have intact chromosomes, CA's, and/or unrepaired fragments, some of which are defined as terminal deletions, if they are capped by one telomere. The model produces a spectrum of terminal deletions with their corresponding probabilities and size distributions for different heavy ions exposures for the first division after exposure. This data provides valuable information because there is limited experimental data available in the literature on the on the actual size of terminal deletions. We compare our model output to the available experimental data and make a reasonable extrapolation on the number of chromosomes lacking telomeres in human lymphocytes exposed to heavy ions. This model generates data which may lead to predictions on the rate of genomic instability in cells after exposure to high charge and energy nuclei affecting astronauts during space missions.

Description: New experimental data show how chromosomal aberrations for low- and high-LET radiation are dependent on DSB repair deficiencies in wild-type, AT and NBS cells. We simulated the development of chromosomal aberrations in these cells lines in a stochastic track-structure-dependent model, in which different cells have different kinetics of DSB repair. We updated a previously formulated model of chromosomal aberrations, which was based on a stochastic Monte Carlo approach, to consider the time-dependence of DSB rejoining. The previous version of the model had an assumption that all DSBs would rejoin, and therefore we called it a time-independent model. The chromosomal-aberrations model takes into account the DNA and track structure for low- and high-LET radiations, and provides an explanation and prediction of the statistics of rare and more complex aberrations. We compared the program-simulated kinetics of DSB rejoining to the experimentally-derived bimodal exponential curves of the DSB kinetics. We scored the formation of translocations, dicentrics, acentric and centric rings, deletions, and inversions. The fraction of DSBs participating in aberrations was studied in relation to the rejoining time. Comparisons of simulated dose dependence for simple aberrations to the experimental dose-dependence for HF19, AT and NBS cells will be made.