Monte Carlo Simulation of Base and Nucleotide Excision Repair of Clustered DNA Damage Sites. I. Model Properties and Predicted Trends

Modelling radiobiology

Radiotherapy has played an essential role in cancer treatment for over a century, and remains one of the best-studied methods of cancer treatment. Because of its close links with the physical sciences, it has been the subject of extensive quantitative mathematical modelling, but a complete understanding of the mechanisms of radiotherapy has remained elusive. In part this is because of the complexity and range of scales involved in radiotherapy-from physical radiation interactions occurring over nanometres to evolution of patient responses over months and years. This review presents the current status and ongoing research in modelling radiotherapy responses across these scales, including basic physical mechanisms of DNA damage, the immediate biological responses this triggers, and genetic- and patient-level determinants of response. Finally, some of the major challenges in this field and potential avenues for future improvements are also discussed.

A quantitative modelling approach for DNA repair on a population scale

The great advances of sequencing technologies allow the in vivo measurement of nuclear processes-such as DNA repair after UV exposure-over entire cell populations. However, data sets usually contain only a few samples over several hours, missing possibly important information in between time points. We developed a data-driven approach to analyse CPD repair kinetics over time in Saccharomyces cerevisiae. In contrast to other studies that consider sequencing signals as an average behaviour, we understand them as the superposition of signals from independent cells. By motivating repair as a stochastic process, we derive a minimal model for which the parameters can be conveniently estimated. We correlate repair parameters to a variety of genomic features that are assumed to influence repair, including transcription rate and nucleosome density. The clearest link was found for the transcription unit length, which has been unreported for budding yeast to our knowledge. The framework hence allows a comprehensive analysis of nuclear processes on a population scale.

Monte Carlo Simulation of Double-Strand Break Induction and Conversion after Ultrasoft X-rays Irradiation

This paper estimates the yields of DNA double-strand breaks (DSBs) induced by ultrasoft X-rays and uses the DSB yields and the repair outcomes to evaluate the relative biological effectiveness (RBE) of ultrasoft X-rays. We simulated the yields of DSB induction and predicted them in the presence and absence of oxygen, using a Monte Carlo damage simulation (MCDS) software, to calculate the RBE. Monte Carlo excision repair (MCER) simulations were also performed to calculate the repair outcomes (correct repairs, mutations, and DSB conversions). Compared to ⁶⁰Co γ-rays, the RBE values for ultrasoft X-rays (titanium K-shell, aluminum K-shell, copper L-shell, and carbon K-shell) for DSB induction were respectively 1.3, 1.9, 2.3, and 2.6 under aerobic conditions and 1.3, 2.1, 2.5, and 2.9 under a hypoxic condition (2% O₂). The RBE values for enzymatic DSBs were 1.6, 2.1, 2.3, and 2.4, respectively, indicating that the enzymatic DSB yields are comparable to the yields of DSB induction. The synergistic effects of DSB induction and enzymatic DSB formation further facilitate cell killing and the advantage in cancer treatment.

Impact of Hypoxia on Relative Biological Effectiveness and Oxygen Enhancement Ratio for a 62-MeV Therapeutic Proton Beam

This study uses the yields of double-strand breaks (DSBs) to determine the relative biological effectiveness (RBE) of proton beams, using cell survival as a biological endpoint. DSB induction is determined when cells locate at different depths (6 positions) along the track of 62 MeV proton beams. The DNA damage yields are estimated using Monte Carlo Damage Simulation (MCDS) software. The repair outcomes are estimated using Monte Carlo excision repair (MCER) simulations. The RBE for cell survival at different oxygen concentrations is calculated using the repair-misrepair-fixation (RMF) model. Using 60Co γ-rays (linear energy transfer (LET) = 2.4 keV/μm) as the reference radiation, the RBE for DSB induction and enzymatic DSB under aerobic condition (21% O2) are in the range 1.0-1.5 and 1.0-1.6 along the track depth, respectively. In accord with RBE obtained from experimental data, RMF model-derived RBE values for cell survival are in the range of 1.0-3.0. The oxygen enhancement ratio (OER) for cell survival (10%) decreases from 3.0 to 2.5 as LET increases from 1.1 to 22.6 keV/μm. The RBE values for severe hypoxia (0.1% O2) are in the range of 1.1-4.4 as LET increases, indicating greater contributions of direct effects for protons. Compared with photon therapy, the overall effect of 62 MeV proton beams results in greater cell death and is further intensified under hypoxic conditions.

The Effects of Dimethylsulfoxide and Oxygen on DNA Damage Induction and Repair Outcomes for Cells Irradiated by 62 MeV Proton and 3.31 MeV Helium Ions

Reactive oxygen species (ROS) play an essential role in radiation-induced indirect actions. In terms of DNA damage, double strand breaks (DSBs) have the greatest effects on the repair of DNA damage, cell survival and transformation. This study evaluated the biological effects of the presence of ROS and oxygen on DSB induction and mutation frequency. The relative biological effectiveness (RBE) and oxygen enhancement ratio (OER) of 62 MeV therapeutic proton beams and 3.31 MeV helium ions were calculated using Monte Carlo damage simulation (MCDS) software. Monte Carlo excision repair (MCER) simulations were used to calculate the repair outcomes (mutation frequency). The RBE values of proton beams decreased to 0.75 in the presence of 0.4 M dimethylsulfoxide (DMSO) and then increases to 0.9 in the presence of 2 M DMSO while the RBE values of 3.31 MeV helium ions increased from 2.9 to 5.7 (0–2 M). The mutation frequency of proton beams also decreased from 0.008–0.065 to 0.004–0.034 per cell per Gy by the addition of 2 M DMSO, indicating that ROS affects both DSB induction and repair outcomes. These results show that the combined use of DMSO in normal tissues and an increased dose in tumor regions increases treatment efficiency.

Effects of indirect actions and oxygen on relative biological effectiveness: estimate of DSB inductions and conversions induced by therapeutic proton beams

Purpose: This study evaluated the DNA double strand breaks (DSBs) induced by indirect actions and its misrepairs to estimate the relative biological effectiveness (RBE) of proton beams.Materials and methods: From experimental data, DSB induction was evaluated in cells irradiated by 62 MeV proton beams in the presence of dimethylsulphoxide (DMSO) and under hypoxic conditions. The DNA damage yields for calculating the RBE were estimated using Monte Carlo Damage Simulation (MCDS) software. The repair outcomes (correct repairs, mutations and DSB conversions) were estimated using Monte Carlo Excision Repair (MCER) simulations.Results: The values for RBE of 62 MeV protons (LET = 1.051 keV/μm) for DSB induction and enzymatic DSB under aerobic condition (21% O₂) was 1.02 and 0.94, respectively, as comparing to ⁶⁰Co γ-rays (LET = 2.4 keV/μm). DMSO mitigated the inference of indirect action and reduced DSB induction to a greater extent when damaged by protons rather than γ-rays, resulting in a decreased RBE of 0.86. DMSO also efficiently prevented enzymatic DSB yields triggered by proton irradiation and reduced the RBE to 0.83. However, hypoxia (2% O₂) produced a similar level of DSB induction with respect to the protons and γ-rays, with a comparable RBE of 1.02.Conclusions: The RBE values of proton beams estimated from DSB induction and enzymatic DSB decreased by 16% and 12%, respectively, in the presence of DMSO. Our findings indicate that the overall effects of DSB induction and enzymatic DSB could intensify the tumor killing, while alleviate normal tissue damage when indirect actions are effectively interrupted.

Proton transport modeling in a realistic biological environment by using TILDA-V

Whether it is in radiobiology to identify DNA lesions or in medicine to adapt the radiotherapeutic protocols, a detailed understanding of the radiation-induced interactions in living matter is required. Monte Carlo track-structure codes have been successfully developed to describe these interactions and predict the radiation-induced energy deposits at the nanoscale level in the medium of interest. In this work, the quantum-mechanically based Monte Carlo track-structure code TILDA-V has been used to compute the slowing-down of protons in water and DNA. Stopping power and range are then reported and compared with existing data. Then, a first application of TILDA-V to cellular irradiations is also reported in order to highlight the absolute necessity of taking into account a realistic description of the cellular environment in microdosimetry.

Petri net-based model of the human DNA base excision repair pathway

Cellular DNA is daily exposed to several damaging agents causing a plethora of DNA lesions. As a first aid to restore DNA integrity, several enzymes got specialized in damage recognition and lesion removal during the process called base excision repair (BER). A large number of DNA damage types and several different readers of nucleic acids lesions during BER pathway as well as two sub-pathways were considered in the definition of a model using the Petri net framework. The intuitive graphical representation in combination with precise mathematical analysis methods are the strong advantages of the Petri net-based representation of biological processes and make Petri nets a promising approach for modeling and analysis of human BER. The reported results provide new information that will aid efforts to characterize in silico knockouts as well as help to predict the sensitivity of the cell with inactivated repair proteins to different types of DNA damage. The results can also help in identifying the by-passing pathways that may lead to lack of pronounced phenotypes associated with mutations in some of the proteins. This knowledge is very useful when DNA damage-inducing drugs are introduced for cancer therapy, and lack of DNA repair is desirable for tumor cell death.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

A New Standard DNA Damage (SDD) Data Format

Our understanding of radiation-induced cellular damage has greatly improved over the past few decades. Despite this progress, there are still many obstacles to fully understand how radiation interacts with biologically relevant cellular components, such as DNA, to cause observable end points such as cell killing. Damage in DNA is identified as a major route of cell killing. One hurdle when modeling biological effects is the difficulty in directly comparing results generated by members of different research groups. Multiple Monte Carlo codes have been developed to simulate damage induction at the DNA scale, while at the same time various groups have developed models that describe DNA repair processes with varying levels of detail. These repair models are intrinsically linked to the damage model employed in their development, making it difficult to disentangle systematic effects in either part of the modeling chain. These modeling chains typically consist of track-structure Monte Carlo simulations of the physical interactions creating direct damages to DNA, followed by simulations of the production and initial reactions of chemical species causing so-called "indirect" damages. After the induction of DNA damage, DNA repair models combine the simulated damage patterns with biological models to determine the biological consequences of the damage. To date, the effect of the environment, such as molecular oxygen (normoxic vs. hypoxic), has been poorly considered. We propose a new standard DNA damage (SDD) data format to unify the interface between the simulation of damage induction in DNA and the biological modeling of DNA repair processes, and introduce the effect of the environment (molecular oxygen or other compounds) as a flexible parameter. Such a standard greatly facilitates inter-model comparisons, providing an ideal environment to tease out model assumptions and identify persistent, underlying mechanisms. Through inter-model comparisons, this unified standard has the potential to greatly advance our understanding of the underlying mechanisms of radiation-induced DNA damage and the resulting observable biological effects when radiation parameters and/or environmental conditions change.

Induction of DNA Damage by Light Ions Relative to 60Co γ-rays

The specific types and numbers of clusters of DNA lesions, including both DNA double-strand breaks (DSBs) and non-DSB clusters, are widely considered 1 of the most important initiating events underlying the relative biological effectiveness (RBE) of the light ions of interest in the treatment of cancer related to megavoltage x-rays and 60Co γ-rays. This review summarizes the categorization of DNA damage, reviews the underlying mechanisms of action by ionizing radiation, and quantifies the general trends in DSB and non-DSB cluster formation by light ions under normoxic and anoxic conditions, as predicted by Monte Carlo simulations that reflect the accumulated evidence from decades of research on radiation damage to DNA. The significance of the absolute and relative numbers of clusters and the local complexity of DSB and non-DSB clusters are discussed in relation to the formation of chromosome aberrations and the loss of cell reproductive capacity. Clinical implications of the dependence of DSB induction on ionization density is reviewed with an eye towards increasing the therapeutic ratio of proton and carbon ion therapy through the explicit optimization of RBE-weighted dose.

Geant4-DNA example applications for track structure simulations in liquid water: A report from the Geant4-DNA Project

This Special Report presents a description of Geant4-DNA user applications dedicated to the simulation of track structures (TS) in liquid water and associated physical quantities (e.g., range, stopping power, mean free path…). These example applications are included in the Geant4 Monte Carlo toolkit and are available in open access. Each application is described and comparisons to recent international recommendations are shown (e.g., ICRU, MIRD), when available. The influence of physics models available in Geant4-DNA for the simulation of electron interactions in liquid water is discussed. Thanks to these applications, the authors show that the most recent sets of physics models available in Geant4-DNA (the so-called "option4" and "option 6" sets) enable more accurate simulation of stopping powers, dose point kernels, and W-values in liquid water, than the default set of models ("option 2") initially provided in Geant4-DNA. They also serve as reference applications for Geant4-DNA users interested in TS simulations.

Computational modeling in glioblastoma: from the prediction of blood-brain barrier permeability to the simulation of tumor behavior

The integrated in silico-in vitro-in vivo approaches have fostered the development of new treatment strategies for glioblastoma patients and improved diagnosis, establishing the bridge between biochemical research and clinical practice. These approaches have provided new insights on the identification of bioactive compounds and on the complex mechanisms underlying the interactions among glioblastoma cells, and the tumor microenvironment. This review focuses on the key advances pertaining to computational modeling in glioblastoma, including predictive data on drug permeability across the blood-brain barrier, tumor growth and treatment responses. Structure- and ligand-based methods have been widely adopted, enabling the study of dynamic and evolutionary aspects of glioblastoma. Their potential applications as predictive tools and the advantages over other well-known methodologies are outlined. Challenges regarding in silico approaches for predicting tumor properties are also discussed.

Computational modeling of radiobiological effects in bone metastases for different radionuclides

PURPOSE

Computational simulation is a simple and practical way to study and to compare a variety of radioisotopes for different medical applications, including the palliative treatment of bone metastases. This study aimed to evaluate and compare cellular effects modelled for different radioisotopes currently in use or under research for treatment of bone metastases using computational methods.

METHODS

Computational models were used to estimate the radiation-induced cellular effects (Virtual Cell Radiobiology algorithm) post-irradiation with selected particles emitted by Strontium-89 (⁸⁹Sr), Samarium-153 (¹⁵³Sm), Lutetium-177 (¹⁷⁷Lu), and Radium-223 (²²³Ra).

RESULTS

Cellular kinetics post-irradiation using ⁸⁹Sr β^- particles, ¹⁵³Sm β^-particles, ¹⁷⁷Lu β^-particles and ²²³Ra α particles showed that the cell response was dose- and radionuclide-dependent. ¹⁷⁷Lu beta minus particles and, in particular, ²²³Ra alpha particles, yielded the lowest survival fraction of all investigated particles.

CONCLUSIONS

²²³Ra alpha particles induced the highest cell death of all investigated particles on metastatic prostate cells in comparison to irradiation with β^-radionuclides, two of the most frequently used radionuclides in the palliative treatment of bone metastases in clinical routine practice. Moreover, the data obtained suggest that the used computational methods might provide some perception about cellular effects following irradiation with different radionuclides.

Analysis of the track- and dose-averaged LET and LET spectra in proton therapy using the geant4 Monte Carlo code

PURPOSE

The motivation of this study was to find and eliminate the cause of errors in dose-averaged linear energy transfer (LET) calculations from therapeutic protons in small targets, such as biological cell layers, calculated using the geant 4 Monte Carlo code. Furthermore, the purpose was also to provide a recommendation to select an appropriate LET quantity from geant 4 simulations to correlate with biological effectiveness of therapeutic protons.

METHODS

The authors developed a particle tracking step based strategy to calculate the average LET quantities (track-averaged LET, LETt and dose-averaged LET, LETd) using geant 4 for different tracking step size limits. A step size limit refers to the maximally allowable tracking step length. The authors investigated how the tracking step size limit influenced the calculated LETt and LETd of protons with six different step limits ranging from 1 to 500 μm in a water phantom irradiated by a 79.7-MeV clinical proton beam. In addition, the authors analyzed the detailed stochastic energy deposition information including fluence spectra and dose spectra of the energy-deposition-per-step of protons. As a reference, the authors also calculated the averaged LET and analyzed the LET spectra combining the Monte Carlo method and the deterministic method. Relative biological effectiveness (RBE) calculations were performed to illustrate the impact of different LET calculation methods on the RBE-weighted dose.

RESULTS

Simulation results showed that the step limit effect was small for LETt but significant for LETd. This resulted from differences in the energy-deposition-per-step between the fluence spectra and dose spectra at different depths in the phantom. Using the Monte Carlo particle tracking method in geant 4 can result in incorrect LETd calculation results in the dose plateau region for small step limits. The erroneous LETd results can be attributed to the algorithm to determine fluctuations in energy deposition along the tracking step in geant 4. The incorrect LETd values lead to substantial differences in the calculated RBE.

CONCLUSIONS

When the geant 4 particle tracking method is used to calculate the average LET values within targets with a small step limit, such as smaller than 500 μm, the authors recommend the use of LETt in the dose plateau region and LETd around the Bragg peak. For a large step limit, i.e., 500 μm, LETd is recommended along the whole Bragg curve. The transition point depends on beam parameters and can be found by determining the location where the gradient of the ratio of LETd and LETt becomes positive.

Positive Cofactor 4 (PC4) is critical for DNA repair pathway re-routing in DT40 cells

PC4 is an abundant single-strand DNA binding protein that has been implicated in transcription and DNA repair. Here, we show that PC4 is involved in the cellular DNA damage response. To elucidate the role, we used the DT40 chicken B cell model, which produces clustered DNA lesions at Ig loci via the action of activation-induced deaminase. Our results help resolve key aspects of immunoglobulin diversification and suggest an essential role of PC4 in repair pathway choice. We show that PC4 ablation in gene conversion (GC)-active cells significantly disrupts GC but has little to no effect on targeted homologous recombination. In agreement, the global double-strand break repair response, as measured by γH2AX foci analysis, is unperturbed 16 hours post irradiation. In cells with the pseudo-genes removed (GC inactive), PC4 ablation reduced the overall mutation rate while simultaneously increasing the transversion mutation ratio. By tagging the N-terminus of PC4, gene conversion and somatic hypermutation are all but abolished even when native non-tagged PC4 is present, indicating a dominant negative effect. Our data point to a very early and deterministic role for PC4 in DNA repair pathway re-routing.

Rapid MCNP simulation of DNA double strand break (DSB) relative biological effectiveness (RBE) for photons, neutrons, and light ions

To account for particle interactions in the extracellular (physical) environment, information from the cell-level Monte Carlo damage simulation (MCDS) for DNA double strand break (DSB) induction has been integrated into the general purpose Monte Carlo N-particle (MCNP) radiation transport code system. The effort to integrate these models is motivated by the need for a computationally efficient model to accurately predict particle relative biological effectiveness (RBE) in cell cultures and in vivo. To illustrate the approach and highlight the impact of the larger scale physical environment (e.g. establishing charged particle equilibrium), we examined the RBE for DSB induction (RBEDSB) of x-rays, (137)Cs γ-rays, neutrons and light ions relative to γ-rays from (60)Co in monolayer cell cultures at various depths in water. Under normoxic conditions, we found that (137)Cs γ-rays are about 1.7% more effective at creating DSB than γ-rays from (60)Co (RBEDSB  =  1.017) whereas 60-250 kV x-rays are 1.1 to 1.25 times more efficient at creating DSB than (60)Co. Under anoxic conditions, kV x-rays may have an RBEDSB up to 1.51 times as large as (60)Co γ-rays. Fission neutrons passing through monolayer cell cultures have an RBEDSB that ranges from 2.6 to 3.0 in normoxic cells, but may be as large as 9.93 for anoxic cells. For proton pencil beams, Monte Carlo simulations suggest an RBEDSB of about 1.2 at the tip of the Bragg peak and up to 1.6 a few mm beyond the Bragg peak. Bragg peak RBEDSB increases with decreasing oxygen concentration, which may create opportunities to apply proton dose painting to help address tumor hypoxia. Modeling of the particle RBE for DSB induction across multiple physical and biological scales has the potential to aid in the interpretation of laboratory experiments and provide useful information to advance the safety and effectiveness of hadron therapy in the treatment of cancer.

Comparative analysis of 11 different radioisotopes for palliative treatment of bone metastases by computational methods

PURPOSE

Throughout the years, the palliative treatment of bone metastases using bone seeking radiotracers has been part of the therapeutic resources used in oncology, but the choice of which bone seeking agent to use is not consensual across sites and limited data are available comparing the characteristics of each radioisotope. Computational simulation is a simple and practical method to study and to compare a variety of radioisotopes for different medical applications, including the palliative treatment of bone metastases. This study aims to evaluate and compare 11 different radioisotopes currently in use or under research for the palliative treatment of bone metastases using computational methods.

METHODS

Computational models were used to estimate the percentage of deoxyribonucleic acid (DNA) damage (fast Monte Carlo damage algorithm), the probability of correct DNA repair (Monte Carlo excision repair algorithm), and the radiation-induced cellular effects (virtual cell radiobiology algorithm) post-irradiation with selected particles emitted by phosphorus-32 ((32)P), strontium-89 ((89)Sr), yttrium-90 ((90)Y ), tin-117 ((117m)Sn), samarium-153 ((153)Sm), holmium-166 ((166)Ho), thulium-170 ((170)Tm), lutetium-177 ((177)Lu), rhenium-186 ((186)Re), rhenium-188 ((188)Re), and radium-223 ((223)Ra).

RESULTS

(223)Ra alpha particles, (177)Lu beta minus particles, and (170)Tm beta minus particles induced the highest cell death of all investigated particles and radioisotopes. The cell survival fraction measured post-irradiation with beta minus particles emitted by (89)Sr and (153)Sm, two of the most frequently used radionuclides in the palliative treatment of bone metastases in clinical routine practice, was higher than (177)Lu beta minus particles and (223)Ra alpha particles.

CONCLUSIONS

(223)Ra and (177)Lu hold the highest potential for palliative treatment of bone metastases of all radioisotopes compared in this study. Data reported here may prompt future in vitro and in vivo experiments comparing different radionuclides for palliative treatment of bone metastases, raise the need for the careful rethinking of the current widespread clinical use of (89)Sr and (153)Sm, and perhaps strengthen the use of (223)Ra and (177)Lu in the palliative treatment of bone metastases.

Effects of indirect actions and oxygen on relative biological effectiveness: estimate of DSB induction and conversion induced by gamma rays and helium ions

Clustered DNA damage other than double-strand breaks (DSBs) can be detrimental to cells and can lead to mutagenesis or cell death. In addition to DSBs induced by ionizing radiation, misrepair of non-DSB clustered damage contributes extra DSBs converted from DNA misrepair via pathways for base excision repair and nucleotide excision repair. This study aimed to quantify the relative biological effectiveness (RBE) when DSB induction and conversion from non-DSB clustered damage misrepair were used as biological endpoints. The results showed that both linear energy transfer (LET) and indirect action had a strong impact on the yields for DSB induction and conversion. RBE values for DSB induction and maximum DSB conversion of helium ions (LET = 120 keV/μm) to (60)Co gamma rays were 3.0 and 3.2, respectively. These RBE values increased to 5.8 and 5.6 in the absence of interference of indirect action initiated by addition of 2-M dimethylsulfoxide. DSB conversion was ∼1-4% of the total non-DSB damage due to gamma rays, which was lower than the 10% estimate by experimental measurement. Five to twenty percent of total non-DSB damage due to helium ions was converted into DSBs. Hence, it may be possible to increase the yields of DSBs in cancerous cells through DNA repair pathways, ultimately enhancing cell killing.

Effects of (5'S)-5',8-cyclo-2'-deoxyadenosine on the base excision repair of oxidatively generated clustered DNA damage. A biochemical and theoretical study

The presence of 5',8-cyclo-2'-deoxyadenosine (5'S)-cdA induces modifications in the geometry of the DNA duplex in the 5'-end direction of the strand and in the 3'-end direction of the complementary strand. As a consequence, the enzymes are probably not able to adjust their active sites in this rigid structure. Additionally, clustered DNA damage sites, a signature of ionising radiation, pose a severe challenge to a cell's repair machinery, particularly base excision repair (BER). To date, clusters containing a DNA base lesion, (5'S)-cdA, which is repaired by nucleotide excision repair, have not been explored. We have therefore investigated whether bistranded clusters containing (5'S)-cdA influence the repairability of an opposed AP site lesion, which is repaired by BER. Using synthetic oligonucleotides containing a bistranded cluster with (5'S)-cdA and an AP site at different interlesion separations, we have shown that in the presence of (5'S)-cdA on the 5'-end side, repair of the AP site by the BER machinery is retarded when the AP site is ≤8 bases from the (5'S)-cdA. However, if (5'S)-cdA is located on the 3'-end side with respect to the AP site, the effect on its repair is much weaker and totally disappears for distances ≥8 bases.

The relative contributions of DNA strand breaks, base damage and clustered lesions to the loss of DNA functionality induced by ionizing radiation

The majority of studies on lethal radiobiological damage have focused on double-strand breaks (DSBs), a type of clustered DNA damage and the evaluation of their toxicity, while other types of clustered DNA damage have received much less attention. The main purpose of this study is to evaluate the contribution of different lesions induced by ionizing radiation to the loss of plasmid DNA functionality. We employed a simple model system comprising E. coli transformed with an irradiated plasmid [pGEM-3Zf (-)] to determine the effect of DSBs and other lesions including base damage and clustered lesions on the functionality ("viability") of the plasmid. The yields of γ-radiation-induced single-strand breaks (SSBs) and DSBs were measured by gel electrophoresis. We found that the transformation efficiency decreases with radiation dose, but this decrease cannot be explained by the formation of DSBs. For example, at doses of 500 and 700 Gy, the relative transformation efficiency falls from 100% to 53% and 26%, respectively, while only 5.7% and 9.1% of the plasmids contain a DSB. In addition, it is also unlikely that randomly distributed base lesions could explain the loss of functionality of the plasmid, since cells can repair them efficiently. However, clustered lesions other than DSBs, which are difficult to repair and result in the loss of information on both DNA strands, have the potential to induce the loss of plasmid functionality. We therefore measured the yields of γ-radiation-induced base lesions and cluster damage, which are respectively converted into SSBs and DSBs by the base excision repair enzymes endonuclease III (Nth) and formamidopyrimidine-DNA glycosylase (Fpg). Our data demonstrate that the yield of cluster damage (i.e., lesions that yield DSBs following digestion) is 31 times higher than that of frank DSBs. This finding suggests that frank DSBs make a relatively minor contribution to the loss of DNA functionality induced by ionizing radiation, while other toxic lesions formed at a much higher frequencies than DSBs must be responsible for the loss of plasmid functionality. These lesions may be clustered lesions/locally multiply damaged sites (LMDS), including base damage, SSBs and/or intrastrand and interstrand crosslinks, leading to the loss of vital information in the DNA. Using a mathematical model, we estimate that at least three toxic lesions are required for the inactivation of plasmid functionality, in part because even these complex lesions can be repaired.

Developing an in silico model of the modulation of base excision repair using methoxyamine for more targeted cancer therapeutics

Base excision repair (BER) is a major DNA repair pathway involved in the processing of exogenous non-bulky base damages from certain classes of cancer chemotherapy drugs as well as ionising radiation (IR). Methoxyamine (MX) is a small molecule chemical inhibitor of BER that is shown to enhance chemotherapy and/or IR cytotoxicity in human cancers. In this study, the authors have analysed the inhibitory effect of MX on the BER pathway kinetics using a computational model of the repair pathway. The inhibitory effect of MX depends on the BER efficiency. The authors have generated variable efficiency groups using different sets of protein concentrations generated by Latin hypercube sampling, and they have clustered simulation results into high, medium and low efficiency repair groups. From analysis of the inhibitory effect of MX on each of the three groups, it is found that the inhibition is most effective for high efficiency BER, and least effective for low efficiency repair.

DNA and cellular effects of charged particles

Development of new radiotherapy strategies based on the use of hadrons, as well as reduction of uncertainties associated with radiation health risk during long-term space flights, requires increasing knowledge of the mechanisms underlying the biological effects of charged particles. It is well known that charged particles are more effective in damaging biological systems than photons. This capability has been related to the production of spatially correlated and/or clustered DNA damage, in particular two or more double-strand breaks (DSB) in close proximity or DSB associated with other lesions within a localized DNA region. These kinds of complex damages are rarely induced by photons. They are difficult to repair accurately and are therefore expected to produce severe consequences at the cellular level. This paper provides a review of radiation-induced cellular effects and will discuss the dependence of cell death and mutation induction on the linear energy transfer of various light and heavy ions. This paper will show the inadequacy of a single physical parameter for describing radiation quality, underlining the importance of the characteristics of the track structure at the submicrometer level to determine the biological effects. This paper will give a description of the physical properties of the track structure that can explain the differences in the spatial distributions of DNA damage, in particular DSB, induced by radiation of different qualities. In addition, this paper will show how a combined experimental and theoretical approach based on Monte Carlo simulations can be useful for providing information on the damage distribution at the nanoscale level. It will also emphasize the importance, especially for DNA damage evaluation at low doses, of the more recent functional approaches based on the use of fluorescent antibodies against proteins involved in the cellular processing of DNA damage. Advantages and limitations of the different experimental techniques will be discussed with particular emphasis on the still unsolved problem of the clustered DNA damage resolution. Development of biophysical models aimed to describe the kinetics of the DNA repair process is underway, and it is expected to support the experimental investigation of the mechanisms underlying the cellular radiation response.

Multiscale QM/MM molecular dynamics study on the first steps of guanine damage by free hydroxyl radicals in solution

Understanding the damage of DNA bases from hydrogen abstraction by free OH radicals is of particular importance to understanding the indirect effect of ionizing radiation. Previous studies address the problem with truncated DNA bases as ab initio quantum simulations required to study such electronic-spin-dependent processes are computationally expensive. Here, for the first time, we employ a multiscale and hybrid quantum mechanical-molecular mechanical simulation to study the interaction of OH radicals with a guanine-deoxyribose-phosphate DNA molecular unit in the presence of water, where all of the water molecules and the deoxyribose-phosphate fragment are treated with the simplistic classical molecular mechanical scheme. Our result illustrates that the presence of water strongly alters the hydrogen-abstraction reaction as the hydrogen bonding of OH radicals with water restricts the relative orientation of the OH radicals with respect to the DNA base (here, guanine). This results in an angular anisotropy in the chemical pathway and a lower efficiency in the hydrogen-abstraction mechanisms than previously anticipated for identical systems in vacuum. The method can easily be extended to single- and double-stranded DNA without any appreciable computational cost as these molecular units can be treated in the classical subsystem, as has been demonstrated here.

Computational modeling of cellular effects post-irradiation with low- and high-let particles and different absorbed doses

The use of computational methods to improve the understanding of biological responses to various types of radiation is an approach where multiple parameters can be modelled and a variety of data is generated. This study compares cellular effects modelled for low absorbed doses against high absorbed doses. The authors hypothesized that low and high absorbed doses would contribute to cell killing via different mechanisms, potentially impacting on targeted tumour radiotherapy outcomes. Cellular kinetics following irradiation with selective low- and high-linear energy transfer (LET) particles were investigated using the Virtual Cell (VC) radiobiology algorithm. Two different cell types were assessed using the VC radiobiology algorithm: human fibroblasts and human crypt cells. The results showed that at lower doses (0.01 to 0.2 Gy), all radiation sources used were equally able to induce cell death (p>0.05, ANOVA). On the other hand, at higher doses (1.0 to 8.0 Gy), the radiation response was LET and dose dependent (p<0.05, ANOVA). The data obtained suggests that the computational methods used might provide some insight into the cellular effects following irradiation. The results also suggest that it may be necessary to re-evaluate cellular radiation-induced effects, particularly at low doses that could affect therapeutic effectiveness.

Reactive molecular dynamics study on the first steps of DNA damage by free hydroxyl radicals

We employ a large scale molecular simulation based on bond-order ReaxFF to simulate the chemical reaction and study the damage to a large fragment of DNA molecule in the solution by ionizing radiation. We illustrate that the randomly distributed clusters of diatomic OH radicals that are primary products of megavoltage ionizing radiation in water-based systems are the main source of hydrogen abstraction as well as formation of carbonyl and hydroxyl groups in the sugar moiety that create holes in the sugar rings. These holes grow up slowly between DNA bases and DNA backbone, and the damage collectively propagates to a DNA single and double strand break.

RBE of low energy electrons and photons

Relative biological effectiveness (RBE) compares the severity of damage induced by a radiation under test at a dose D relative to the reference radiation D(x) for the same biological endpoint. RBE is an important parameter in estimation of risk from exposure to ionizing radiation (IR). The present work provides a review of the recently published data and the knowledge of the RBE of low energy electrons and photons. The review presents RBE values derived from experimental data and model calculations including cell inactivation, chromosome aberration, cell transformation, micronuclei formation and induction of double-strand breaks. Biophysical models, including physical features of radiation track, and microdosimetry parameters are presented, analysed and compared with experimental data. The biological effects of low energy electrons and photons are of particular interest in radiation biology as these are strongly absorbed in micrometer and sub-micrometer layers of tissue. RBE values not only depend on the electron and photon energies but also on the irradiation condition, cell type and experimental conditions.

A new calculation on spectrum of direct DNA damage induced by low-energy electrons

In this work, direct DNA damage induced by low-energy electrons (<5 keV) is simulated using Monte Carlo methods, and the resulting yield of various strand breaks and base damages in cellular environment is presented. The simulation is based on a new inelastic cross section for the production of electron track structure in liquid water, and on ionization cross sections of DNA bases to generate base radical. Especially, a systematic approach of simulating detailed base damage is suggested. This approach includes improvement of a volume model of DNA, generation of the DNA base sequence, conversion of ionization events in liquid water at hit site to the ionization interaction of electrons with DNA bases and development of an algorithm to convert a base radical to a damage. The results obtained in terms of strand breaks are compared with those of experiments and other theoretical calculations, and good agreement was obtained. The yield of detailed base damages and clustered DNA damages caused by the combination of various strand breaks and base damages is presented, and the corresponding distribution characteristics are analyzed. The influence of the relative content of base pairs A-T and G-C in a DNA segment on the yield of both strand breaks and base damages is also explored. The present work provides fundamental information on DNA damage and represents the first effort toward the goal of obtaining the spectrum of clustered DNA damage including detailed base damages, for the mechanistic interpretation and prediction of radiation effects.

A nanodosimetric model of radiation-induced clustered DNA damage yields

We present a nanodosimetric model for predicting the yield of double strand breaks (DSBs) and non-DSB clustered damages induced in irradiated DNA. The model uses experimental ionization cluster size distributions measured in a gas model by an ion counting nanodosimeter or, alternatively, distributions simulated by a Monte Carlo track structure code developed to simulate the nanodosimeter. The model is based on a straightforward combinatorial approach translating ionizations, as measured or simulated in a sensitive gas volume, to lesions in a DNA segment of one-two helical turns considered equivalent to the sensitive volume of the nanodosimeter. The two model parameters, corresponding to the probability that a single ion detected by the nanodosimeter corresponds to a single strand break or a single lesion (strand break or base damage) in the equivalent DNA segment, were tuned by fitting the model-predicted yields to previously measured double-strand break and double-strand lesion yields in plasmid DNA irradiated with protons and helium nuclei. Model predictions were also compared to both yield data simulated by the PARTRAC code for protons of a wide range of different energies and experimental DSB and non-DSB clustered DNA damage yield data from the literature. The applicability and limitations of this model in predicting the LET dependence of clustered DNA damage yields are discussed.

3-Methylindole is mutagenic and a possible pulmonary carcinogen

Previous work has shown that bioactivation of the cigarette smoke pneumotoxicant 3-methylindole (3MI) by pulmonary cytochrome P450 enzymes is directly associated with formation of DNA adducts. Here, we present evidence that normal human lung epithelial cells, exposed to low micromolar concentrations of 3MI, showed extensive DNA damage, as measured by the comet assay, with similar potency to the prototypical genotoxic agents, doxorubicin and irinotecan. The DNA damage caused by 3MI was predominantly caused by single-strand breaks. Furthermore, we show that this damage decreased with time, given a subtoxic concentration, with detectable DNA fragmentation peaking 4 h after exposure and diminishing to untreated levels within 24 h. Pretreatment with an inhibitor of poly(ADP-ribose) polymerase 1 (PARP1), NU1025, nearly doubled the DNA damage produced by 5 microM 3MI, implying that PARP1, which among other activities, functions to repair single-strand breaks in DNA, repaired at least some of the 3MI-induced DNA fragmentation. A key cellular response to DNA damage, phosphorylation, and nuclear localization of p53 was seen at subtoxic levels of 3MI exposure. 3MI was highly mutagenic, with essentially the same potency as the prototype carcinogen, benzo[a]pyrene, only when a lung-expressed CYP2F3 enzyme was used to dehydrogenate 3MI to its putative DNA-alkylating intermediate. Conversely, a rat liver S9 metabolic system did not bioactivate 3MI to its mutagenic intermediate(s). Concentrations higher than 25 microM caused apoptosis, which became extensive at 100 microM, similar to the response seen with 10 microM doxorubicin. Our findings indicate that there is a low concentration window in which 3MI can cause extensive DNA damage and mutation, without triggering apoptotic defenses, reinforcing the hypothesis that inhaled 3MI from cigarette smoke may be a potent lung-selective carcinogen.

Detection of complex DNA damage in gamma-irradiated acute lymphoblastic leukemia Pre-b NALM-6 cells

Bistranded complex DNA damage, i.e., double-strand breaks (DSBs) and non-DSB oxidative clustered DNA lesions, is hypothesized to challenge the repair mechanisms of the cell and consequently the genomic integrity. The oxidative clustered DNA lesions may be persistent and may accumulate in human cancer cells for long times after irradiation. To evaluate the detection and possible accumulation of oxidative clustered DNA lesions in leukemia cells exposed to doses equivalent to those used in radiotherapy, we measured the induction of DSBs and three different types of oxidative clustered DNA lesions in NALM-6 cells, a human acute lymphoblastic leukemia (ALL) pre-B cell line, after exposure to (137)Cs gamma rays. For the detection and measurement of DSBs and oxidative clustered DNA lesions, we used an adaptation of the neutral comet assay (single-cell gel electrophoresis) using E. coli repair enzymes (Endo IV, Fpg and Endo III) as enzymatic probes. We found a linear dose response for the induction of DSBs and oxidative clustered DNA lesions. Clustered DNA lesions were more prevalent than prompt DSBs. For each DSB induced by radiation, approximately 2.5 oxidative clustered DNA lesions were detected. To our knowledge, this is the first study to demonstrate the detection and linear induction of oxidative clustered DNA lesions with radiation dose in an ALL cell line. These results point to the biological significance of clustered DNA lesions.

BGRT: biologically guided radiation therapy-the future is fast approaching!

Rapid advances in functional and biological imaging, predictive assays, and our understanding of the molecular and cellular responses underpinning treatment outcomes herald the coming of the long-sought goal of implementing patient-specific biologically guided radiation therapy (BGRT) in the clinic. Biological imaging and predictive assays have the potential to provide patient-specific, three-dimensional information to characterize the radiation response characteristics of tumor and normal structures. Within the next decade, it will be possible to combine such information with advanced delivery technologies to design and deliver biologically conformed, individualized therapies in the clinic. The full implementation of BGRT in the clinic will require new technologies and additional research. However, even the partial implementation of BGRT treatment planning may have the potential to substantially impact clinical outcomes.

Induction and processing of oxidative clustered DNA lesions in 56Fe-ion-irradiated human monocytes

Space and cosmic radiation is characterized by energetic heavy ions of high linear energy transfer (LET). Although both low- and high-LET radiations can create oxidative clustered DNA lesions and double-strand breaks (DSBs), the local complexity of oxidative clustered DNA lesions tends to increase with increasing LET. We irradiated 28SC human monocytes with doses from 0-10 Gy of (56)Fe ions (1.046 GeV/ nucleon, LET = 148 keV/microm) and determined the induction and processing of prompt DSBs and oxidative clustered DNA lesions using pulsed-field gel electrophoresis (PFGE) and Number Average Length Analysis (NALA). The (56)Fe ions produced decreased yields of DSBs (10.9 DSB Gy(-1) Gbp(-1)) and clusters (1 DSB: approximately 0.8 Fpg clusters: approximately 0.7 Endo III clusters: approximately 0.5 Endo IV clusters) compared to previous results with (137)Cs gamma rays. The difference in the relative biological effectiveness (RBE) of the measured and predicted DSB yields may be due to the formation of spatially correlated DSBs (regionally multiply damaged sites) which result in small DNA fragments that are difficult to detect with the PFGE assay. The processing data suggest enhanced difficulty compared with gamma rays in the processing of DSBs but not clusters. At the same time, apoptosis is increased compared to that seen with gamma rays. The enhanced levels of apoptosis observed after exposure to (56)Fe ions may be due to the elimination of cells carrying high levels of persistent DNA clusters that are removed only by cell death and/or "splitting" during DNA replication.

Fast Monte Carlo simulation of DNA damage formed by electrons and light ions

The passage of ionizing radiation through living organisms initiates physical and chemical processes that create clusters of damaged nucleotides within one or two turns of the DNA. These clusters are widely considered an important initiating event for the induction of other biological endpoints, including cell killing and neoplastic transformation. Monte Carlo simulations of the DNA damage formation process are a useful adjunct to experiments because they provide additional information about the spatial configuration of damage within a cluster. In this paper, the fast Monte Carlo damage simulation (MCDS) algorithm is re-parameterized so that yields of double-strand breaks, single-strand breaks and sites of multiple base damage can be simulated for electrons, protons and alpha particles with kinetic energies on the order of GeV. The MCDS algorithm provides a useful, quasi-phenomenological scheme to interpolate damage yields from computationally expensive, but more detailed, track-structure simulations. The predicted characteristics of various classes of damage produced by electrons, protons and alpha particles, such as average number of lesions per DNA damage cluster and cluster length in base pairs, are presented. A study examining the effects on damage complexity of an extrinsic free radical scavenger, dimethyl sulfoxide, is also presented. The reported studies provide new information that will aid efforts to characterize the relative biological effectiveness of high-energy protons and other light ions, which are sometimes used in particle therapy for the treatment of cancer.

On the biophysical interpretation of lethal DNA lesions induced by ionising radiation

Although DNA damage is widely viewed as a critical target for the induction of cell killing by ionising radiation, the exact nature of DNA damage responsible for these effects is unknown. To address this issue, the probability of forming lethal damage by single proton tracks, derived from published survival data for Chinese hamster V79 cells irradiated by protons with energies from 0.57 to 5.01 MeV, has been compared to estimated yields of clustered DNA lesions and repair outcomes calculated with Monte Carlo models. The reported studies provide new information about the potential relationship between the induction and repair of clustered DNA damage and trends in the expected number of lethal events for protons with increasing linear energy transfer (LET). A good correlation was found between the number of lethal events in V79 cells and the induction of double-strand breaks (DSBs) consisting of three or more elementary DNA lesions. For the yields of other types of DNA damage, as well as point mutations formed through the misrepair of base damage and single-strand breaks, observed trends with increasing LET are not consistent with trends in the yields of lethal events. This observation suggests that the relative biological effectiveness (RBE) of protons of varying quality may be more closely related to the induction of complex DSBs rather than other forms of damage.

Monte Carlo Simulation of Base and Nucleotide Excision Repair of Clustered DNA Damage Sites. II. Comparisons of Model Predictions to Measured Data

Clustered damage sites other than double-strand breaks (DSBs) have the potential to contribute to deleterious effects of ionizing radiation, such as cell killing and mutagenesis. In the companion article (Semenenko et al., Radiat. Res. 164, 180-193, 2005), a general Monte Carlo framework to simulate key steps in the base and nucleotide excision repair of DNA damage other than DSBs is proposed. In this article, model predictions are compared to measured data for selected low-and high-LET radiations. The Monte Carlo model reproduces experimental observations for the formation of enzymatic DSBs in Escherichia coli and cells of two Chinese hamster cell lines (V79 and xrs5). Comparisons of model predictions with experimental values for low-LET radiation suggest that an inhibition of DNA backbone incision at the sites of base damage by opposing strand breaks is active over longer distances between the damaged base and the strand break in hamster cells (8 bp) compared to E. coli (3 bp). Model estimates for the induction of point mutations in the human hypoxanthine guanine phosphoribosyl transferase (HPRT) gene by ionizing radiation are of the same order of magnitude as the measured mutation frequencies. Trends in the mutation frequency for low- and high-LET radiation are predicted correctly by the model. The agreement between selected experimental data sets and simulation results provides some confidence in postulated mechanisms for excision repair of DNA damage other than DSBs and suggests that the proposed Monte Carlo scheme is useful for predicting repair outcomes.