Low-energy electrons inside active DNA models: A tool to elucidate the radiation action mechanisms

To postulate radiation action mechanisms and to test them by Monte Carlo simulation, a complex computer model was developed consisting of major components for the generation of a radiation spectrum, biomolecular structures, and electron track structures in liquid water. As the radiation source¹²⁵I is employed here; it is an excellent test radiation due to its exactly localized position in the DNA molecule and high biological toxicity as a consequence of the emission of short-ranging Auger electrons. A linear DNA plasmid model (Pomplun 1991) which can actively respond to radical attack (Terrissol and Pomplun 1994) has been modified into a nucleosome model representing the double-helix of DNA with 146 basepairs and more than 9000 atoms surrounding the histones. The introduction of this new target structure allows a more realistic simulation of cellular conditions. Using the model's decay accumulation aspect, the situation of many break and survival experiments can be approximated and the influence of several cellular parameters tested. As a first step, a correlation between the size of energy depositions and strand-break patterns was sought.     

Model for radial dependence of frequency distributions for energy imparted in nanometer volumes from HZE particles

This paper develops a deterministic model of frequency distributions for energy imparted (total energy deposition) in small volumes similar to DNA molecules from high-energy ions of interest for space radiation protection and cancer therapy. Frequency distributions for energy imparted are useful for considering radiation quality and for modeling biological damage produced by ionizing radiation. For high-energy ions, secondary electron (delta-ray) tracks originating from a primary ion track make dominant contributions to energy deposition events in small volumes. Our method uses the distribution of electrons produced about an ion's path and incorporates results from Monte Carlo simulation of electron tracks to predict frequency distributions for ions, including their dependence on radial distance. The contribution from primary ion events is treated using an impact parameter formalism of spatially restricted linear energy transfer (LET) and energy-transfer straggling. We validate our model by comparing it directly to results from Monte Carlo simulations for proton and alpha-particle tracks. We show for the first time frequency distributions of energy imparted in DNA structures by several high-energy ions such as cosmic-ray iron ions. Our comparison with results from Monte Carlo simulations at low energies indicates the accuracy of the method.     

Monte Carlo simulation of DNA strand-break induction in supercoiled plasmid pBR322 DNA from indirect effects

We present a new Monte Carlo simulation code system (DBREAK) of the detailed events that occur when ionizing radiation interacts with water and DNA molecules. The model treats the initial energy deposition by radiation, the formation of chemically active species, subsequent diffusion-controlled chemical reactions, and induction of DNA strand breaks. DBREAK assumes one-hit single-strand break (SSB) and two-hit double-strand break (DSB) mechanisms. A high-resolution model of plasmid DNA structure has been introduced. The calculated results are compared with the results of previously performed experiments of the same type. Under aerobic conditions, 89.4% of the DNA damage was attributed to OH-radical and 10.5% and 0.1% to e^– _aq and H, respectively. We also compared the differences between liquid-water track structure and gas-phase-water track structure. The calculated yield of SSBs by liquid-water track structure exceeded that of gas-phase-water track structure by a factor of 1.2. Received: 13 February 1997 / Accepted in revised form: 26 August 1997     

Chromatin organization contributes to non-randomly distributed double-strand breaks after exposure to high-LET radiation

The influence of higher-order chromatin structure on the non-random distribution of DNA double-strand breaks induced by high-LET radiation was investigated. Five different chromatin structures (intact cells, condensed and decondensed chromatin, nucleoids and naked genomic DNA) from GM5758 cells or K562 cells were irradiated with (137)Cs gamma-ray photons and 125 keV/microm nitrogen ions (16-25 MeV/nucleon). DNA was purified with a modified lysis procedure to avoid release of heat-labile sites, and fragment size distributions and double-strand break yields were analyzed by different pulsed-field gel electrophoresis protocols. Whereas double-strand breaks in photon-irradiated cells were randomly distributed, irradiation of intact K562 cells with high-LET nitrogen ions produced an excess of non-randomly distributed DNA fragments 10 kb-1 Mbp in size. Complete removal of proteins eliminated this non-random component. There was a gradual increase in the yield of double-strand breaks for each chromatin decondensation step, and compared to intact cells, the yields for naked DNA (in buffer without scavengers) increased 83 and 25 times after photon and nitrogen-ion irradiation, respectively. The corresponding relative biological effectiveness decreased from 1.6-1.8 for intact cells to 0.49 for the naked DNA. We conclude that the organization of DNA into chromatin fiber and higher-order structures is responsible for the majority of non-randomly distributed double-strand breaks induced by high-LET radiation. However, our data suggest a complex interaction between track structure and chromatin organization over several levels.     

Investigation into the foundations of the track-event theory of cell survival and the radiation action model based on nanodosimetry

This work aims at elaborating the basic assumptions behind the “track-event theory” (TET) and its derivate “radiation action model based on nanodosimetry” (RAMN) by clearly distinguishing between effects of tracks at the cellular level and the induction of lesions in subcellular targets. It is demonstrated that the model assumptions of Poisson distribution and statistical independence of the frequency of single and clustered DNA lesions are dispensable for multi-event distributions because they follow from the Poisson distribution of the number of tracks affecting the considered target volume. It is also shown that making these assumptions for the single-event distributions of the number of lethal and sublethal lesions within a cell would lead to an essentially exponential dose dependence of survival for practically relevant values of the absorbed dose. Furthermore, it is elucidated that the model equation used for consideration of repair within the TET is based on the assumption that DNA lesions induced by different tracks are repaired independently. Consequently, the model equation is presumably inconsistent with the model assumptions and requires an additional model parameter. Furthermore, the methodology for deriving model parameters from nanodosimetric properties of particle track structure is critically assessed. Based on data from proton track simulations it is shown that the assumption of statistically independent targets leads to the prediction of negligible frequency of clustered DNA damage. An approach is outlined how track structure could be considered in determining the model parameters, and the implications for TET and RAMN are discussed.

     

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.     

High-LET irradiation of a DNA-binding protein: protein-protein and DNA-protein crosslinks

The chromosomal protein MC1 is a monomeric protein of 93 amino acids that is able to bind any DNA but has a slight preferential affinity for some sequences and structures, like cruciform and minicircles. The protein has been irradiated with 36Ar18+ ions of 95 MeV/nucleon. The LET of these particles in water is close to 270 keV/microm. We tested the activity of the protein by measuring its ability to form complexes with DNA. We tested the integrity of the protein by measuring the molecular weight of the species formed. Compared with gamma radiation, we observed for the same dose a less efficient inactivation of the protein, a greater protection of the protein by the bound DNA, a lower induction of chain breakage, and a greater production of protein-protein and DNA-protein crosslinks. The results are discussed in terms of the quantitative and the qualitative differences between the two types of radiation: The global radical yield is slightly higher with gamma rays, whereas the density of radicals produced along the particle track is considerably higher with argon ions.     

Trinucleotide repeat instability: a hairpin curve at the crossroads of replication, recombination, and repair

The trinucleotide repeats that expand to cause human disease form hairpin structures in vitro that are proposed to be the major source of their genetic instability in vivo. If a replication fork is a train speeding along a track of double-stranded DNA, the trinucleotide repeats are a hairpin curve in the track. Experiments have demonstrated that the train can become derailed at the hairpin curve, resulting in significant damage to the track. Repair of the track often results in contractions and expansions of track length. In this review we introduce the in vitro evidence for why CTG/CAG and CCG/CGG repeats are inherently unstable and discuss how experiments in model organisms have implicated the replication, recombination and repair machinery as contributors to trinucleotide repeat instability in vivo.     

Intrinsic radiation sensitivity: cellular signaling is the key

The concept that the balance between DNA damage and repair determines intrinsic radiation sensitivity has dominated radiobiology for several decades. There is undeniably a cause- effect relationship between radiation-induced molecular alterations in the genomic DNA and cellular consequences. In the last decade, however, it has become obvious that the chromatin context affects the fate of damaged DNA and that cellular signaling is an important factor in defining intrinsic radiation sensitivity. Damaged DNA is the site of signal generation; however, alternative signaling at the plasma membrane is triggered: Reactive oxygen species (ROS) inactivate phosphatases and consequently cause activation of kinases localized at the plasma membrane; this includes ligand-independent activation of receptor kinases. Cells with an apparently functional DNA repair system may show increased radiation sensitivity due to deficiencies in specific kinases essential for repair activation and checkpoint control. Other signals that determine intrinsic radiosensitivity may affect proneness to apoptosis, the balance between DNA damage fixation and repair, and the translocation of proteins participating in the response to ionizing radiation. Interplay between the various signals decides the extent to which the repair of radiation-inflicted damage is supported or limited; in some cell types, this includes DNA-damage-independent processes guided by plasma membrane-generated signaling. Cellular signaling in the context of specific subcellular structures is the key to understanding how the molecular effects of radiation are expressed as biological consequences in various cell types. A systems approach should bring us closer to this end.     

DNA complex lesions induced by protons and ⋅-particles: track structure characteristics determining linear energy transfer and particle type dependence

The yield of DNA double-strand breaks (dsb) and DNA complex lesions induced by protons and α-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 α-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/μm. 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. Received: 20 December 1996 / Accepted in revised form: 5 March 1997     

Differential oxidation of deoxyribose in DNA by gamma and alpha-particle radiation

Emerging evidence points to the importance of deoxyribose oxidation in the toxicity of oxidative DNA damage, including the formation of protein-DNA crosslinks and base adducts. With the goal of understanding the differences in deoxyribose oxidation chemistry known to occur with different oxidants, we have compared the formation of one product of 3'-oxidation of deoxyribose in DNA, 3'-phosphoglycolaldehyde (PGA) residues, in isolated DNA and cells exposed to ionizing radiations. A recently developed gas chromatography/negative chemical ionization mass spectrometry method was used to quantify PGA residues in purified DNA and in human TK6 lymphoblastoid cells exposed to gamma radiation (60Co) and alpha particles (241Am). The level of PGA residues was then correlated with the total quantity of deoxyribose oxidation determined by plasmid topoisomer analysis. Alpha-particle irradiation (0-100 Gy) of purified DNA in 50 mM potassium phosphate (pH 7.4) produced a linear dose response of 0.13 PGA residues per 10(6) nucleotides per gray. When normalized to an estimate of the total number of deoxyribose oxidation events (2.0 per 10(6) nucleotides per gray), PGA formation occurred in 7% (+/-0.5) of deoxyribose oxidation events produced by alpha-particle radiation. In contrast, the efficiency of PGA formation in gamma-irradiated DNA was found to be 1% (+/-0.02), which indicates a shift in the chemistry of deoxyribose oxidation, possibly as a result of the different track structures of the two types of ionizing radiation. Studies with gamma radiation were extended to TK6 cells, in which it was observed that gamma radiation produced a linear dose response of 0.0019 PGA residues per 10(6) nucleotides per gray. This is consistent with an approximately 1000-fold quenching effect in cells, similar to the results of other published studies of oxidative DNA damage in vivo.     

Biophysical modeling of fragment length distributions of DNA plasmids after X and heavy-ion irradiation analyzed by atomic force microscopy

The investigation of fragment length distributions of plasmid DNA gives insight into the influence of localized energy distribution on the induction of DNA damage, particularly the clustering of double-strand breaks. We present an approach that determines the fragment length distributions of plasmid DNA after heavy-ion irradiation by using the Local Effect Model. We find a good agreement of our simulations with experimental fragment distributions derived from atomic force microscopy (AFM) studies by including experimental constraints typical for AFM. Our calculations reveal that by comparing the fragmentation in terms of fluence, we can uniquely distinguish the effect of different radiation qualities. For very high-LET irradiation using nickel or uranium ions, no difference between their fragment distributions can be expected for the same dose level. However, for carbon ions with an intermediate LET, the fragmentation pattern differs from the distribution for very high-LET particles. The results of the model calculations can be used to determine the optimal experimental parameters for a demonstration of the influence of track structure on primary radiation damage. Additionally, we compare the results of our model for two different plasmid geometries.     

Palindrome regeneration by template strand-switching mechanism at the origin of DNA replication of porcine circovirus via the rolling-circle melting-pot replication model

Palindromic sequences (inverted repeats) flanking the origin of DNA replication with the potential of forming single-stranded stem-loop cruciform structures have been reported to be essential for replication of the circular genomes of many prokaryotic and eukaryotic systems. In this study, mutant genomes of porcine circovirus with deletions in the origin-flanking palindrome and incapable of forming any cruciform structures invariably yielded progeny viruses containing longer and more stable palindromes. These results suggest that origin-flanking palindromes are essential for termination but not for initiation of DNA replication. Detection of template strand switching in the middle of an inverted repeat strand among the progeny viruses demonstrated that both the minus genome and a corresponding palindromic strand served as templates simultaneously during DNA biosynthesis and supports the recently proposed rolling-circle "melting-pot" replication model. The genome configuration presented by this model, a four-stranded tertiary structure, provides insights into the mechanisms of DNA replication, inverted repeat correction (or conversion), and illegitimate recombination of any circular DNA molecule with an origin-flanking palindrome.     

离子束生物效应及应用中值得研究的几个问题

重离子束生物效应及重离子束在生命科学中的应用研究,在国内外物理学与生命科学领域中得到了广泛的开展,但对出现的一些现象还没有深刻地揭露其本质,作出机理性解释。为了深入研究,本文提出一些值得研究的问题供参考,如:重离子径迹结构及能量沉积分布模型,DNA辐射敏感位点,质量沉积-分子改造,直接作用与间接作用,放射性核束的应用等。     

Induction and rejoining of DNA double-strand breaks in normal human skin fibroblasts after exposure to radiation of different linear energy transfer: possible roles of track structure and chromatin organization

DNA double-strand breaks are nonrandomly induced by high-LET radiation. Differences in the induction and rejoining of DSBs after irradiation with ions having different LET were detected by fragment analysis. The data obtained indicate that the track structure of the traversing particle and its interaction with the different chromatin structures of the cellular DNA influence the yield as well as the distribution of the induced damage. The induction and rejoining of clustered DSBs induced by the same nitrogen ion fluence at LETs of 80-225 keV/microm were investigated by a detailed analysis of the DNA fragmentation patterns in normal human fibroblasts. The DSBs in the cells were allowed to rejoin during incubations for 0-20 h. Two separate pulsed-field gel electrophoresis protocols were used, optimized for separation of fragments in the size ranges 1-6 Mbp and 5 kbp-1.5 Mbp. A strong influence of LET on the level of DSB induction was evident. The DSB yield increased from 4.5 +/- 0.2 to 10.0 +/- 0.3 DSBs per particle traversal through the cell nucleus when LET increased from 80 to 225 keV/microm. Further, the size distribution of the DNA fragments showed a significant dependence on radiation quality, with an excess of fragments at 50-200 kbp and around 1 Mbp. Differences in repair kinetics were also evident, with slower rejoining for increasing LET, and the initial nonrandom fragment distributions were still present after 1 h of repair.     

DNA Topoisomerase Inhibitors: Trapping a DNA-Cleaving Machine in Motion

Type II topoisomerases regulate DNA topology by making a double-stranded break in one DNA duplex, transporting another DNA segment through this break and then resealing it. Bacterial type IIA topoisomerase inhibitors, such as fluoroquinolones and novel bacterial topoisomerase inhibitors, can trap DNA cleavage complexes with double- or single-stranded cleaved DNA. To study the mode of action of such compounds, 21 crystal structures of a “gyrase^CORE” fusion truncate of Staphyloccocus aureus DNA gyrase complexed with DNA and diverse inhibitors have been published, as well as 4 structures lacking inhibitors. These structures have the DNA in various cleavage states and appear to track trajectories along the catalytic paths of the DNA cleavage/religation steps. The various conformations sampled by these multiple “gyrase^CORE” structures show rigid body movements of the catalytic GyrA WHD and GyrB TOPRIM domains across the dimer interface. Conformational changes common to all compound-bound structures suggest common mechanisms for DNA cleavage-stabilizing compounds. The structures suggest that S. aureus gyrase uses a single moving-metal ion for cleavage and that the central four base pairs need to be stretched between the two catalytic sites, in order for a scissile phosphate to attract a metal ion to the A-site to catalyze cleavage, after which it is “stored” in another coordination configuration (B-site) in the vicinity. We present a simplified model for the catalytic cycle in which capture of the transported DNA segment causes conformational changes in the ATPase domain that push the DNA gate open, resulting in stretching and cleaving the gate-DNA in two steps.     

Nuclear dynamics of RAD52 group homologous recombination proteins in response to DNA damage

Recombination between homologous DNA molecules is essential for the proper maintenance and duplication of the genome, and for the repair of exogenously induced DNA damage such as double-strand breaks. Homologous recombination requires the RAD52 group proteins, including Rad51, Rad52 and Rad54. Upon treatment of mammalian cells with ionizing radiation, these proteins accumulate into foci at sites of DNA damage induction. We show that these foci are dynamic structures of which Rad51 is a stably associated core component, whereas Rad52 and Rad54 rapidly and reversibly interact with the structure. Furthermore, we show that the majority of the proteins are not part of the same multi-protein complex in the absence of DNA damage. Executing DNA transactions through dynamic multi-protein complexes, rather than stable holo-complexes, allows flexibility. In the case of DNA repair, for example, it will facilitate cross-talk between different DNA repair pathways and coupling to other DNA transactions, such as replication.     

Simulation of DNA damage after proton irradiation

The biophysical radiation track simulation model PARTRAC was improved by implementing new interaction cross sections for protons in water. Computer-simulated tracks of energy deposition events from protons and their secondary electrons were superimposed on a higher-order DNA target model describing the spatial coordinates of the whole genome inside a human cell. Induction of DNA double-strand breaks was simulated for proton irradiation with LET values between 1.6 and 70 keV/microm and various reference radiation qualities. The yield of DSBs after proton irradiation was found to rise continuously with increasing LET up to about 20 DSBs per Gbp and Gy, corresponding to an RBE up to 2.2. About half of this increase resulted from a higher yield of DSB clusters associated with small fragments below 10 kbp. Exclusion of experimentally unresolved multiple DSBs reduced the maximum DSB yield by 30% and shifted it to an LET of about 40 keV/microm. Simulated fragment size distributions deviated significantly from random breakage distributions over the whole size range after irradiation with protons with an LET above 10 keV/microm. Determination of DSB yields using equations derived for random breakage resulted in an underestimation by up to 20%. The inclusion of background fragments had only a minor influence on the distribution of the DNA fragments induced by radiation. Despite limited numerical agreement, the simulations reproduced the trends in proton-induced DNA DSBs and fragment induction found in recent experiments.