Modelling the effect of incorporated halogenated pyrimidine on radiation-induced DNA strand breaks

Electronic properties of DNA-related molecules containing a bromine atom

PURPOSE

To clarify the radiosensitization mechanism masking the Auger effect of the cells possessing brominated DNA, the electronic properties of DNA-related molecules containing Br were investigated by X-ray spectroscopy and specific heat measurement.

MATERIALS AND METHODS

X-ray absorption near-edge structure (XANES) and X-ray photoemission spectroscopy (XPS) were used to measure the electronic properties of the nucleotides with and without Br. We determined the specific heat of 5-bromouracil crystals with thymine as a reference molecule at low temperatures of 3-48 K to calculate the microscopic state numbers.

RESULTS

Obtained XANES and XPS spectra indicated that both the lowest unoccupied molecular orbital (LUMO) and the core-levels were not affected by the Br incorporation. The state numbers of 5-bromouracil calculated from the specific heats obtained around 25 K was about 1.5 times larger than that for thymine below 20 K, although the numbers were almost the same below 5 K.

DISCUSSION

Our results suggest that the Br atom may not contribute substantially to the LUMO and core-level electronic states of the molecule, but rather to the microscopic states related to the excitation of lattice vibrations, which may be involved in valence electronic states.

Photoelectron-photoion(s) coincidence studies of molecules of biological interest

Photoelectron-photoion(s) coincidence, PEPICO, experiments with synchrotron radiation have become one of the most powerful tools to investigate dissociative photoionization thanks to their selectivity. In this paper their application to the study of molecular species of biological interest in the gas phase is reviewed. Some applications of PEPICO to the study of potential radiosensitizers, amino acids and small peptides and opportunities offered by the advent of novel methods for the production of beams of these molecules are discussed.

Electron-Induced Dissociation of the Potential Radiosensitizer 5-Selenocyanato-2'-deoxyuridine

5-Selenocyanato-2'-deoxyuridine (SeCNdU) is a recently proposed radiosensitizer based on 2'-deoxyuridine (dU) with the electron-affinic selenocyanato (-SeCN) side group attached at the C5 position of uracil. Since electron interaction processes may be an important source of DNA damage by ionizing radiation, we have studied low-energy dissociative electron attachment to SeCNdU in the gas phase. Negative ion formation has been obtained by means of mass spectrometry, where a rich fragmentation pattern is observed even at ∼0 eV. The reaction pathways exhibiting the highest ion yields are C₄N₂O₂H₂Se^•- and CN^-, both involving a cleavage of the Se-CN bond. The heaviest fragment anion observed is C₉N₂O₅H₁₀Se^•-, where besides the charged species, the hydrogen and cyano radicals are also formed. Further decomposition channels also yield the highly reactive hydroxyl radical, which possesses a high DNA damage potential. All observed channels have experimentally determined onsets at 0 eV, which are supported by calculations performed at the M06-2X/aug-cc-pVTZ level. The calculations comprise the thermochemical thresholds at standard and experimental (428.15 K, 3 × 10^-11 atm) conditions together with the adiabatic electron affinities. The present study shows that low-energy electrons very effectively decompose SeCNdU upon attachment of thermal electrons, producing a large variety of charged fragments and radicals.

Radiation damage to single stranded oligonucleotide trimers labelled with 5-iodopyrimidines

The radiolysis of deoxygenated aqueous solution containing trimeric oligonucleotides labelled with iodinated pyrimidines and Tris-HCl as the hydroxyl radical scavenger leads to electron attachment to the halogenated bases that mainly results in single strand breaks. The iodinated trimers are 2-fold more sensitive to solvated electrons than the brominated oligonucleotides, which is explained by the barrier-free dissociation of the iodinated base anions. The present study fills the literature gap concerning the chemistry triggered by ionizing radiation in the iodinated pyrimidines incorporated into DNA.

5-Selenocyanatouracil: A Potential Hypoxic Radiosensitizer. Electron Attachment Induced Formation of Selenium Centered Radical

The propensity of 5-selenocyanatouracil (SeCNU) to decomposition induced by attachment of electron was scrutinized with the G3B3 composite quantum-chemical method and radiolytic studies. Favorable thermodynamic (Gibbs free reaction energy of -13.65 kcal/mol) and kinetic (Gibbs free activation energy of 1.22 kcal/mol) characteristics revealed by the G3B3 free energy profile suggest SeCNU to be sensitive to electron attachment. The title compound was synthesized in the reaction between uracil and selenocyanogen chloride in acetic acid. Then, an aqueous and deoxygenated solution of the HPLC purified compound containing tert-butanol as a hydroxyl radical scavenger was irradiated with X-rays. SeCNU radio-degradation results in two major products: the U-Se-Se-U dimer and the adduct of the ^●OtBu radical to the U-Se^● radical, U-Se-OtBu. The effects of radiolysis as well as the results of G3B3 calculations point to U-Se^● as the primary product of dissociative electron attachment to SeCNU. The MTT test shows that SeCNU is nontoxic in vitro in concentrations equal to or lower than 10^-6 M. Ionizing radiation will probably induce cytotoxic intra- and interstrand DNA cross-links as well as protein-DNA cross-links in the genomic DNA labeled with SeCNU.

Fragmentation of pure and hydrated clusters of 5Br-uracil by low energy carbon ions: observation of hydrated fragments

The fragmentation of the isolated 5-bromouracil (5BrU) molecule and pure and nano-hydrated 5BrU clusters induced by low energy ¹²C⁴⁺ ions has been studied.

Radiation track, DNA damage and response-a review

The purpose of this paper has been to review the current status and progress of the field of radiation biophysics, and draw attention to the fact that physics, in general, and radiation physics in particular, with the aid of mathematical modeling, can help elucidate biological mechanisms and cancer therapies. We hypothesize that concepts of condensed-matter physics along with the new genomic knowledge and technologies and mechanistic mathematical modeling in conjunction with advances in experimental DNA (Deoxyrinonucleic acid molecule) repair and cell signaling have now provided us with unprecedented opportunities in radiation biophysics to address problems in targeted cancer therapy, and genetic risk estimation in humans. Obviously, one is not dealing with 'low-hanging fruit', but it will be a major scientific achievement if it becomes possible to state, in another decade or so, that we can link mechanistically the stages between the initial radiation-induced DNA damage; in particular, at doses of radiation less than 2 Gy and with structural changes in genomic DNA as a precursor to cell inactivation and/or mutations leading to genetic diseases. The paper presents recent development in the physics of radiation track structure contained in the computer code system KURBUC, in particular for low-energy electrons in the condensed phase of water for which we provide a comprehensive discussion of the dielectric response function approach. The state-of-the-art in the simulation of proton and carbon ion tracks in the Bragg peak region is also presented. The paper presents a critical discussion of the models used for elastic scattering, and the validity of the trajectory approach in low-electron transport. Brief discussions of mechanistic and quantitative aspects of microdosimetry, DNA damage and DNA repair are also included as developed by the authors' work.

Site- and state-selected photofragmentation of 2Br-pyrimidine

The fragmentation of the 2Br-pyrimidine molecule following direct valence photoionization or inner shell excitation has been studied by electron-ion coincidence experiments. 2Br-pyrimidine has been chosen as a model for the class of pyrimidinic building blocks of three nucleic acids and several radiosensitizers. It is known that the site- and state-localization of energy deposition, typical of inner shell excitation, results in the enhancement of the total ion yield as well as in changes in the relative intensity of the different fragmentation channels. Here we address the question of the origin of this selective fragmentation by using electron-ion coincidence techniques. The results show that the fragmentation is strongly selective in the final singly charged ion state, independently of the process that leads to the population of that state, and the dominant fragmentation patterns correlate with the nearest appearance potential.

5-Thiocyanato-2'-deoxyuridine as a possible radiosensitizer: electron-induced formation of uracil-C5-thiyl radical and its dimerization

In this work, we have synthesized 5-thiocyanato-2'-deoxyuridine (SCNdU) along with the C6-deuterated nucleobase 5-thiocyanatouracil (6-D-SCNU) and studied their reactions with radiation-produced electrons. ESR spectra in γ-irradiated nitrogen-saturated frozen homogeneous solutions (7.5 M LiCl in H2O or D2O) of these compounds show that electron-induced S-CN bond cleavage occurs to form a thiyl radical (dU-5-S˙ or 6-D-U-5-S˙) and CN(-)via the initial π-anion radical (SCNdU˙(-)) intermediate in which the excess electron is on the uracil base. HPLC and LC-MS/MS studies of γ-irradiated N2-saturated aqueous solutions of SCNdU in the presence of sodium formate as a OH-radical scavenger at ambient temperature show the formation of the dU-5S-5S-dU dimer in preference to dU by about 10 to 1 ratio. This shows that both possible routes of electron-induced bond cleavage (dUC5-SCN and S-CN) in SCNdU˙(-) and dU-5-S˙ formation are preferred for the production of the σ-type uracilyl radical (dU˙) by 10 fold. DFT/M06-2x/6-31++G(d,p) calculations employing the polarizable continuum model (PCM) for aqueous solutions show that dU-5-S˙ and CN(-) formation was thermodynamically favored by over 15 kcal mol(-1) (ΔG) compared to dU˙ and SCN(-) production. The activation barriers for C5-S and S-CN bond cleavage in SCNdU˙(-) amount to 8.7 and 4.0 kcal mol(-1), respectively, favoring dU-5-S˙ and CN(-) formation. These results support the experimental observation of S-CN bond cleavage by electron addition to SCNdU that results in the formation of dU-5-S˙ and the subsequent dU-5S-5S-dU dimer. This establishes SCNdU as a potential radiosensitizer that could cause intra- and inter-strand crosslinking as well as DNA-protein crosslinking via S-S dimer formation.

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.

DSB repair model for mammalian cells in early S and G1 phases of the cell cycle: application to damage induced by ionizing radiation of different quality

The purpose of this work is to test the hypothesis that kinetics of double strand breaks (DSB) repair is governed by complexity of DSB. To test the hypothesis we used our recent published mechanistic mathematical model of DSB repair for DSB induced by selected protons, deuterons, and helium ions of different energies representing radiations of different qualities. In light of recent advances in experimental and computational techniques, the most appropriate method to study cellular responses in radiation therapy, and exposures to low doses of ionizing radiations is using mechanistic approaches. To this end, we proposed a 'bottom-up' approach to study cellular response that starts with the DNA damage. Monte Carlo track structure method was employed to simulate initial damage induced in the genomic DNA by direct and indirect effects. Among the different types of DNA damage, DSB are known to be induced in simple and complex forms. The DSB repair model in G1 and early S phases of the cell cycle was employed to calculate the repair kinetics. The model considers the repair of simple and complex DSB, and the DSB produced in the heterochromatin. The inverse sampling method was used to calculate the repair kinetics for each individual DSB. The overall repair kinetics for 500 DSB induced by single tracks of the radiation under test were compared with experimental results. The results show that the model is capable of predicting the repair kinetics for the DSB induced by radiations of different qualities within an accepted range of uncertainty.

Radiation induced base excision repair (BER): a mechanistic mathematical approach

This paper presents a mechanistic model of base excision repair (BER) pathway for the repair of single-stand breaks (SSBs) and oxidized base lesions produced by ionizing radiation (IR). The model is based on law of mass action kinetics to translate the biochemical processes involved, step-by-step, in the BER pathway to translate into mathematical equations. The BER is divided into two subpathways, short-patch repair (SPR) and long-patch repair (LPR). SPR involves in replacement of single nucleotide via Pol β and ligation of the ends via XRCC1 and Ligase III, while LPR involves in replacement of multiple nucleotides via PCNA, Pol δ/ɛ and FEN 1, and ligation via Ligase I. A hallmark of IR is the production of closely spaced lesions within a turn of DNA helix (named complex lesions), which have been attributed to a slower repair process. The model presented considers fast and slow component of BER kinetics by assigning SPR for simple lesions and LPR for complex lesions. In the absence of in vivo reaction rate constants for the BER proteins, we have deduced a set of rate constants based on different published experimental measurements including accumulation kinetics obtained from UVA irradiation, overall SSB repair kinetic experiments, and overall BER kinetics from live-cell imaging experiments. The model was further used to calculate the repair kinetics of complex base lesions via the LPR subpathway and compared to foci kinetic experiments for cells irradiated with γ rays, Si, and Fe ions. The model calculation show good agreement with experimental measurements for both overall repair and repair of complex lesions. Furthermore, using the model we explored different mechanisms responsible for inhibition of repair when higher LET and HZE particles are used and concluded that increasing the damage complexity can inhibit initiation of LPR after the AP site removal step in BER.

Comparison of the bromodeoxyuridine-mediated sensitization effects between low-LET and high-LET ionizing radiation on DNA double-strand breaks

The incorporation of halogenated pyrmidines such as bromo- and iodo-deoxyuridines (BrdU, IdU) into DNA as thymidine analogs enhances cellular radiosensitivity when high-linear energy transfer (LET) radiation is not used. Although it is known that high-LET ionizing radiation confers fewer biological effects resulting from halogenated pyrimidine incorporation, the exact mechanisms of reduced radiosensitivity with high-LET radiation are not clear. We investigated the radiosensitization effects of halogenated pyrimidines with high-LET radiation using accelerated carbon and iron ions. Cells synchronized into the G1 phase after unifilar (1 cell cycle) and bifilar (2 cell cycles) substitution with 10 µM BrdU were exposed to various degrees of LET with heavy ions and X-rays. We then carried out a colony formation assay to measure cell survival. The γ-H2AX focus formation assay provided a measure of DNA double-strand break (DSB) formation and repair kinetics. Chromosomal aberration formations for the first post-irradiation metaphase were also scored. For both low-LET X-rays and carbon ions (13 keV/µm), BrdU incorporation led to impaired DNA repair kinetics, a larger initial number of DNA DSBs more frequent chromosomal aberrations at the first post-irradiated metaphase, and increased radiosensitivity for cell lethality. The enhancement ratio was higher after bifilar substitution. In contrast, no such synergistic enhancements were observed after high-LET irradiation with carbon and iron ions (70 and 200 keV/µm, respectively), even after bifilar substitution. Our results suggest that BrdU substitution did not modify the number and quality of DNA DSBs produced by high-LET radiation. The incorporation of halogenated pyrimidines may produce more complex/clustered DNA damage along with radicals formed by low-LET ionizing radiation. In contrast, the severity of damage produced by high-LET radiation may undermine the effects of BrdU and account for the observed minimal radiosensitization effects.

Monte Carlo role in radiobiological modelling of radiotherapy outcomes

Radiobiological models are essential components of modern radiotherapy. They are increasingly applied to optimize and evaluate the quality of different treatment planning modalities. They are frequently used in designing new radiotherapy clinical trials by estimating the expected therapeutic ratio of new protocols. In radiobiology, the therapeutic ratio is estimated from the expected gain in tumour control probability (TCP) to the risk of normal tissue complication probability (NTCP). However, estimates of TCP/NTCP are currently based on the deterministic and simplistic linear-quadratic formalism with limited prediction power when applied prospectively. Given the complex and stochastic nature of the physical, chemical and biological interactions associated with spatial and temporal radiation induced effects in living tissues, it is conjectured that methods based on Monte Carlo (MC) analysis may provide better estimates of TCP/NTCP for radiotherapy treatment planning and trial design. Indeed, over the past few decades, methods based on MC have demonstrated superior performance for accurate simulation of radiation transport, tumour growth and particle track structures; however, successful application of modelling radiobiological response and outcomes in radiotherapy is still hampered with several challenges. In this review, we provide an overview of some of the main techniques used in radiobiological modelling for radiotherapy, with focus on the MC role as a promising computational vehicle. We highlight the current challenges, issues and future potentials of the MC approach towards a comprehensive systems-based framework in radiobiological modelling for radiotherapy.

Temporary anion states of pyrimidine and halopyrimidines

The empty-level electronic structures of pyrimidine and its 2-chloro, 2-bromo, and 5-bromo derivatives have been studied with electron transmission spectroscopy (ETS) and dissociative electron attachment spectroscopy (DEAS) in the 0-5 eV energy range. The spectral features were assigned to the corresponding anion states with the support of theoretical calculations at the ab initio and density functional theory levels. The empty orbital energies obtained by simple Koopmans' theorem calculations, scaled with empirical equations, quantitatively reproduced the energies of vertical electron attachment to π* and σ* empty orbitals measured in the ET spectra and predicted vertical electron affinities close to zero for the three halo derivatives. The total anion currents of the halo derivatives, measured at the walls of the collision chamber as a function of the impact electron energy, presented intense maxima below 0.5 eV. The mass-selected spectra showed that, in this energy, range the total anion current is essentially due to halide fragment anions. The DEA cross sections of the bromo derivatives were found to be about six times larger than that of the chloro derivative. The absolute cross sections at incident electron energies close to zero were evaluated to be 10(-16)-10(-15) cm(2).

Dosimetry on sub-cellular level for intracellular incorporated auger-electron-emitting radionuclides: a comparison of Monte Carlo simulations and analytic calculations

A quantitative dosimetric comparison was performed between Monte Carlo (MC) simulations and analytic calculations at the (sub) cellular level (V79 cells) for four nucleus-incorporated radiochemicals ((125)I/(123)I/(77)Br-UdR and A (125)IP) and two radiochemicals that localised mainly in the cytoplasm of cells ((125)I-dihydrorhodamine and Na(2)(51)CrO(4)). A microscopic investigation around the decay site of the three DNA-incorporated radionuclides ((125)I/(123)I/(77)Br-UdR) was also carried out. On the whole, deviations between MC and analytic calculations for the absorbed dose and dose rate to the cell nucleus were within ∼10%. The dose rate to the nucleus for the radiochemicals that mainly localised in the cytoplasm was greater than that for the nucleus-incorporated ones. Also evident was that the dose rate to the nucleus was approximately the same for the three DNA-incorporated radiochemicals. In contrast to the small differences found between MC and analytic calculations for the (average) absorbed dose to the nucleus, the dosimetric analysis at the microscopic level for the three DNA-incorporated radionuclides showed that the two computational approaches lead to a completely different energy deposition pattern around the decay site.

Molecular simulations of ultra-low-energy nitrogen ion bombardment of A-DNA in vacuum

For investigating mechanisms involved in low-energy ion beam induced mutation, besides experiments using low-energy and low-fluence ions to bombard naked DNA, molecular simulations were carried out as an effort towards the insight in molecular interactions between ions and DNA. In the current study, Monte Carlo (MC) and molecular dynamics (MD) simulations were applied. The results of MC simulations provide some clues about the interaction energies and sites of preference of N-ion bombardment on an A-DNA short duplex strand. MD simulations of a single N-ion moving towards the same DNA strand with different linear velocities corresponding to bombardment energies of 0.1, 1, 10 and 100 eV revealed information about changes in bond lengths and visibly distorted structures of bombarded nucleotides. The simulations demonstrated that ion-bombardment-induced DNA change in structure was not a random but preferential effect.

The Auger spectroscopy of pyrimidine and halogen-substituted pyrimidines

The C 1s and N 1s Auger spectra of pyrimidine, 2-chloropyrimidine, and 5-bromopyrimidine have been measured in an electron impact experiment at 1000 eV. In the case of the halogen-substituted pyrimidines, also the Cl 2p and Br 3d Auger spectra have been recorded. We have thoroughly analyzed and interpreted all the Auger spectra recorded here with the aid of accurate Green's function calculations with a large basis set. The spectra are extremely complex with thousands of states contributing and almost no single-state feature even near the double ionization threshold. Besides reproducing and explaining with great detail nearly all the main spectral features observed, the calculations have successfully unraveled the interplay among the different C 1s core hole chemical shifts in each molecule and how this affects some fingerprinting details in the composite C 1s Auger spectra.

Methoxyamine potentiates iododeoxyuridine-induced radiosensitization by altering cell cycle kinetics and enhancing senescence

We previously reported that methoxyamine (an inhibitor of base excision repair) potentiates iododeoxyuridine (IUdR)-induced radiosensitization in human tumor cells. In this study, we investigated the potential mechanisms of this enhanced cell death. Human colorectal carcinoma RKO cells were exposed to IUdR (3 micromol/L) and/or methoxyamine (3 mmol/L) for 48 hours before ionizing radiation (5 Gy). We found that IUdR/methoxyamine altered cell cycle kinetics and led to an increased G1 population but a decreased S population before ionizing radiation. Immediately following ionizing radiation (up to 6 hours), IUdR/methoxyamine-pretreated cells showed a stringent G1-S checkpoint but an insufficient G2-M checkpoint, whereas a prolonged G1 arrest, containing 2CG1 and 4CG1 cells, was found at later times up to 72 hours. Levels of cell cycle-specific markers [p21, p27, cyclin A, cyclin B1, and pcdc2(Y15)] and DNA damage signaling proteins [gammaH2AX, pChk1(S317), and pChk2(T68)] supported these altered cell cycle kinetics. Interestingly, we found that IUdR/methoxyamine pretreatment reduced ionizing radiation-induced apoptosis. Additionally, the extent of cell death through necrosis or autophagy seemed similar in all (IUdR +/- methoxyamine + ionizing radiation) treatment groups. However, a larger population of senescence-activated beta-galactosidase-positive cells was seen in IUdR/methoxyamine/ionizing radiation-treated cells, which was correlated with the increased activation of the senescence factors p53 and pRb. These data indicate that IUdR/methoxyamine pretreatment enhanced the effects of ionizing radiation by causing a prolonged G1 cell cycle arrest and by promoting stress-induced premature senescence. Thus, senescence, a novel ionizing radiation-induced tumor suppression pathway, may be effectively targeted by IUdR/methoxyamine pretreatment, resulting in an improved therapeutic gain for ionizing radiation.