A kinetic model of single-strand annealing for the repair of DNA double-strand breaks

A Boolean network model of the double-strand break repair pathway choice

Double strand break (DSB) repair is critical to maintaining the integrity of the genome. DSB repair deficiency underlies multiple pathologies, including cancer, chromosome instability syndromes, and, potentially, neurodevelopmental defects. DSB repair is mainly handled by two pathways: highly accurate homologous recombination (HR), which requires a sister chromatid for template-based repair, limited to S/G2 phases of the cell cycle, and canonical non-homologous end joining (c-NHEJ), available throughout the cell cycle in which minimum homology is sufficient for highly efficient yet error-prone repair. Some circumstances, such as cancer, require alternative highly mutagenic DSB repair pathways like microhomology-mediated end-joining (MMEJ) and single-strand annealing (SSA), which are triggered to attend to DNA damage. These non-canonical repair alternatives are emerging as prominent drivers of resistance in drug-based tumor therapies. Multiple DSB repair options require tight inter-pathway regulation to prevent unscheduled activities. In addition to this complexity, epigenetic modifications of the histones surrounding the DSB region are emerging as critical regulators of the DSB repair pathway choice. Modeling approaches to understanding DSBs repair pathway choice are advantageous to perform simulations and generate predictions on previously uncharacterized aspects of DSBs response. In this work, we present a Boolean network model of the DSB repair pathway choice that incorporates the knowledge, into a dynamic system, of the inter-pathways regulation involved in DSB repair, i.e., HR, c-NHEJ, SSA, and MMEJ. Our model recapitulates the well-characterized HR activity observed in wild-type cells in response to DSBs. It also recovers clinically relevant behaviors of BRCA1/FANCS mutants, and their corresponding drug resistance mechanisms ascribed to DNA repair gain-of-function pathogenic variants. Since epigenetic modifiers are dynamic and possible druggable targets, we incorporated them into our model to better characterize their involvement in DSB repair. Our model predicted that loss of the TIP60 complex and its corresponding histone acetylation activity leads to activation of SSA in response to DSBs. Our experimental validation showed that TIP60 effectively prevents activation of RAD52, a key SSA executor, and confirms the suitable use of Boolean network modeling for understanding DNA DSB repair.

Influence of Acute Exercise on DNA Repair and PARP Activity before and after Irradiation in Lymphocytes from Trained and Untrained Individuals

Several studies indicate that acute exercise induces DNA damage, whereas regular exercise increases DNA repair kinetics. Although the molecular mechanisms are not completely understood, the induction of endogenous reactive oxygen species (ROS) during acute exhaustive exercise due to metabolic processes might be responsible for the observed DNA damage, while an adaptive increase in antioxidant capacity due to regular physical activity seems to play an important protective role. However, the protective effect of physical activity on exogenously induced DNA damage in human immune cells has been poorly investigated. We asked the question whether individuals with a high aerobic capacity would have an enhanced response to radiation-induced DNA damage. Immune cells are highly sensitive to radiation and exercise affects lymphocyte dynamics and immune function. Therefore, we measured endogenous and radiation-induced DNA strand breaks and poly (ADP-ribose) polymerase-1 (PARP1) activity in peripheral blood mononuclear cells (PBMCs) from endurance-trained (maximum rate of oxygen consumption measured during incremental exercise V'O_2max > 55 mL/min/kg) and untrained (V'O_2max < 45 mL/min/kg) young healthy male volunteers before and after exhaustive exercise. Our results indicate that: (i) acute exercise induces DNA strand breaks in lymphocytes only in untrained individuals, (ii) following acute exercise, trained individuals repaired radiation-induced DNA strand breaks faster than untrained individuals, and (iii) trained subjects retained a higher level of radiation-induced PARP1 activity after acute exercise. The results of the present study indicate that increased aerobic fitness can protect immune cells against radiation-induced DNA strand breaks.

MODELLING DSB REPAIR KINETICS FOR DNA DAMAGE INDUCED BY PROTON AND CARBON IONS

Proton and carbon therapy are the main choices of particle therapy for cancer treatment. Particle dose distribution is superior to conventional photon therapy dose distribution due to Bragg peak. However, the basic biology of cellular damage and cell death is not well understood. The aim of this work is to present a mechanistic model of double strand break (DSB) repair that predicts the repair kinetics of damage induced by particles employed in cancer therapy. Monte Carlo Track Damage Simulation (MCDS) was employed to model DNA damage. The frequency of DSB and SSB was computed for proton and carbon ions. DSBs were subjected to repair model to calculate the repair kinetics. Two distinct DSB repair models dependent on the cell cycle were proposed. The DSB repair model contains non-homologous end joining (NHEJ), homologous recombination (HR) and back up non-homologous end joining (B-NHEJ) repair processes. The DSB complexity results in the switch in the repair pathway from NHEJ to a slower process that starts with DSB end resection. DSB end resection in early S and G1 phases of the cell cycle enhances the B-NHEJ repair pathway, while in late S and G2 phases of the cell cycle promotes HR repair pathway. The repair model was transformed to a set of nonlinear differential equations. The model calculates the overall repair kinetics and protein temporal repair activity at the site of damage. The damage and repair model provides a detailed mechanistic understanding of all processes that are involved in the damage induction and repair. The number of DSB and their complexity increase as the particle energy decreases due to the proximity of particle interactions in water. The repair kinetics show a biphasic behaviour that is due to the NHEJ fast repair of simple type DSB and HR slow repair of complex type DSB.

Modeling the interplay between DNA-PK, Artemis, and ATM in non-homologous end-joining repair in G1 phase of the cell cycle

Modeling a biological process equips us with more comprehensive insight into the process and a more advantageous experimental design. Non-homologous end joining (NHEJ) is a major double-strand break (DSB) repair pathway that occurs throughout the cell cycle. The objective of the current work is to model the fast and slow phases of NHEJ in G1 phase of the cell cycle following exposure to ionizing radiation (IR). The fast phase contains the major components of NHEJ; Ku70/80 complex, DNA-dependent protein kinase catalytic subunit (DNA-PKcs), and XLF/XRCC4/ligase IV complex (XXL). The slow phase in G1 phase of the cell cycle is associated with more complex lesions and involves ATM and Artemis proteins in addition to the major components. Parameters are mainly obtained from experimental data. The model is successful in predicting the kinetics of DSB foci in 13 normal, ATM-deficient, and Artemis-deficient mammalian fibroblast cell lines in G1 phase of the cell cycle after exposure to low doses of IR. The involvement of ATM provides the model with the potency to be connected to different signaling pathways. Ku70/80 concentration and DNA-binding rate as well as XXL concentration and enzymatic activity are introduced as the best targets for affecting NHEJ DSB repair process. On the basis of the current model, decreasing concentration and DNA binding rate of DNA-PKcs is more effective than inhibiting its activity towards the Artemis protein.

The effect of stochasticity on repair of DNA double strand breaks throughout non-homologous end joining pathway

DNA double strand breaks (DSBs) are the most lethal lesions of DNA induced by ionizing radiation, industrial chemicals and a wide variety of drugs used in chemotherapy. In the context of DNA damage response system modelling, uncertainty may arise in several ways such as number of induced DSBs, kinetic rates and measurement error in observable quantities. Therefore, using the stochastic approaches is imperative to gain further insight into the dynamic behaviour of DSBs repair process. In this article, a continuous-time Markov chain (CTMC) model of the non-homologous end joining (NHEJ) mechanism is formulated according to the DSB complexity. Additionally, a Metropolis Monte Carlo method is used to perform maximum likelihood estimation of the kinetic rate constants. Here, the effects of fluctuating kinetic rates and DSBs induction rate of the NHEJ mechanism are investigated. The stochastic realizations of the total yield of simple and complex DSBs ligation are simulated to compare their asymptotic dynamics. Furthermore, it has been proved that the total yield of DSBs has a normal distribution for sufficiently large number of DSBs. In order to estimate the expected duration of repairing DSBs, the probability distribution of DSBs lifetime is calculated based on the CTMC NHEJ model. Moreover, the variability of total yield of DSBs during constant low-dose radiation is evaluated in the presented model. The findings indicate that in stochastic NHEJ model, when there is no new DSBs induction through the repair process, all DSBs are eventually repaired. However, when DSBs are induced by constant low-dose radiation, a number of DSBs remains un-repaired.

Small de novo CNVs as biomarkers of parental exposure to low doses of ionizing radiation of caesium-137

The radiological accident in Goiania in 1987 caused a trail of human contamination, animal, plant and environmental by a radionuclide. Exposure to ionizing radiation results in different types of DNA lesions. The mutagenic effects of ionizing radiation on the germline are special concern because they can endures for several generations, leading to an increase in the rate of mutations in children of irradiated parents. Thus, to evaluate the biological mechanisms of ionizing radiation in somatic and germline cells, with consequent determination of the rate mutations, is extremely important for the estimation of genetic risks. Recently it was established that Chromosomal Microarray Analysis is an important tool for detecting wide spectra of gains or losses in the human genome. Here we present the results of the effect of accidental exposure to low doses of ionizing radiation on the formation of CNVs in the progeny of a human population accidentally exposed to Caesium-137 during the radiological accident in Goiânia, Brazil.

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.

Mechanistic Modelling and Bayesian Inference Elucidates the Variable Dynamics of Double-Strand Break Repair

DNA double-strand breaks are lesions that form during metabolism, DNA replication and exposure to mutagens. When a double-strand break occurs one of a number of repair mechanisms is recruited, all of which have differing propensities for mutational events. Despite DNA repair being of crucial importance, the relative contribution of these mechanisms and their regulatory interactions remain to be fully elucidated. Understanding these mutational processes will have a profound impact on our knowledge of genomic instability, with implications across health, disease and evolution. Here we present a new method to model the combined activation of non-homologous end joining, single strand annealing and alternative end joining, following exposure to ionising radiation. We use Bayesian statistics to integrate eight biological data sets of double-strand break repair curves under varying genetic knockouts and confirm that our model is predictive by re-simulating and comparing to additional data. Analysis of the model suggests that there are at least three disjoint modes of repair, which we assign as fast, slow and intermediate. Our results show that when multiple data sets are combined, the rate for intermediate repair is variable amongst genetic knockouts. Further analysis suggests that the ratio between slow and intermediate repair depends on the presence or absence of DNA-PKcs and Ku70, which implies that non-homologous end joining and alternative end joining are not independent. Finally, we consider the proportion of double-strand breaks within each mechanism as a time series and predict activity as a function of repair rate. We outline how our insights can be directly tested using imaging and sequencing techniques and conclude that there is evidence of variable dynamics in alternative repair pathways. Our approach is an important step towards providing a unifying theoretical framework for the dynamics of DNA repair processes.

Genome-based, mechanism-driven computational modeling of risks of ionizing radiation: The next frontier in genetic risk estimation?

Research activity in the field of estimation of genetic risks of ionizing radiation to human populations started in the late 1940s and now appears to be passing through a plateau phase. This paper provides a background to the concepts, findings and methods of risk estimation that guided the field through the period of its growth to the beginning of the 21st century. It draws attention to several key facts: (a) thus far, genetic risk estimates have been made indirectly using mutation data collected in mouse radiation studies; (b) important uncertainties and unsolved problems remain, one notable example being that we still do not know the sensitivity of human female germ cells to radiation-induced mutations; and (c) the concept that dominated the field thus far, namely, that radiation exposures to germ cells can result in single gene diseases in the descendants of those exposed has been replaced by the concept that radiation exposure can cause DNA deletions, often involving more than one gene. Genetic risk estimation now encompasses work devoted to studies on DNA deletions induced in human germ cells, their expected frequencies, and phenotypes and associated clinical consequences in the progeny. We argue that the time is ripe to embark on a human genome-based, mechanism-driven, computational modeling of genetic risks of ionizing radiation, and we present a provisional framework for catalyzing research in the field in the 21st century.

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.

Ionizing radiation and genetic risks. XVII. Formation mechanisms underlying naturally occurring DNA deletions in the human genome and their potential relevance for bridging the gap between induced DNA double-strand breaks and deletions in irradiated germ cells

While much is known about radiation-induced DNA double-strand breaks (DSBs) and their repair, the question of how deletions of different sizes arise as a result of the processing of DSBs by the cell's repair systems has not been fully answered. In order to bridge this gap between DSBs and deletions, we critically reviewed published data on mechanisms pertaining to: (a) repair of DNA DSBs (from basic studies in this area); (b) formation of naturally occurring structural variation (SV) - especially of deletions - in the human genome (from genomic studies) and (c) radiation-induced mutations and structural chromosomal aberrations in mammalian somatic cells (from radiation mutagenesis and radiation cytogenetic studies). The specific aim was to assess the relative importance of the postulated mechanisms in generating deletions in the human genome and examine whether empirical data on radiation-induced deletions in mouse germ cells are consistent with predictions of these mechanisms. The mechanisms include (a) NHEJ, a DSB repair process that does not require any homology and which functions in all stages of the cell cycle (and is of particular relevance in G0/G1); (b) MMEJ, also a DSB repair process but which requires microhomology and which presumably functions in all cell cycle stages; (c) NAHR, a recombination-based DSB repair mechanism which operates in prophase I of meiosis in germ cells; (d) MMBIR, a microhomology-mediated, replication-based mechanism which operates in the S phase of the cell cycle, and (e) strand slippage during replication (involved in the origin of small insertions and deletions (INDELs). Our analysis permits the inference that, between them, these five mechanisms can explain nearly all naturally occurring deletions of different sizes identified in the human genome, NAHR and MMBIR being potentially more versatile in this regard. With respect to radiation-induced deletions, the basic studies suggest that those arising as a result of the operation of NHEJ/MMEJ processes, as currently formulated, are expected to be relatively small. However, data on induced mutations in mouse spermatogonial stem cells (irradiation in G0/G1 phase of the cell cycle and DSB repair presumed to be via NHEJ predominantly) show that most are associated with deletions of different sizes, some in the megabase range. There is thus a 'discrepancy' between what the basic studies suggest and the empirical observations in mutagenesis studies. This discrepancy, however, is only an apparent but not a real one. It can be resolved by considering the issue of deletions in the broader context of and in conjunction with the organization of chromatin in chromosomes and nuclear architecture, the conceptual framework for which already exists in studies carried out during the past fifteen years or so. In this paper, we specifically hypothesize that repair of DSBs induced in chromatin loops may offer a basis to explain the induction of deletions of different sizes and suggest an approach to test the hypothesis. We emphasize that the bridging of the gap between induced DSB and resulting deletions of different sizes is critical for current efforts in computational modeling of genetic risks.

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.

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.