Mathematical modelling of the automated FADU assay for the quantification of DNA strand breaks and their repair in human peripheral mononuclear blood cells
BMC Biophysics volume 7, Article number: 9 (2014)
Cells continuously undergo DNA damage from exogenous agents like irradiation or genotoxic chemicals or from endogenous radicals produced by normal cellular metabolic activities. DNA strand breaks are one of the most common genotoxic lesions and they can also arise as intermediates of DNA repair activity. Unrepaired DNA damage can lead to genomic instability, which can massively compromise the health status of organisms. Therefore it is important to measure and quantify DNA damage and its repair.
We have previously published an automated method for measuring DNA strand breaks based on fluorimetric detection of alkaline DNA unwinding , and here we present a mathematical model of the FADU assay, which enables to an analytic expression for the relation between measured fluorescence and the number of strand breaks.
Assessment of the formation and also the repair of DNA strand breaks is a crucial functional parameter to investigate genotoxicity in living cells. A reliable and convenient method to quantify DNA strand breakage is therefore of significant importance for a wide variety of scientific fields, e.g. toxicology, pharmacology, epidemiology and medical sciences.
The fact that the DNA denatures (unwinds) under alkaline conditions has been used to establish FADU assay for measuring DNA strand breaks . This method is based on the limited denaturation of the DNA under precisely controlled conditions of pH, temperature and time. Under these conditions only chromosome ends plus open sites in the DNA serve as starting points for the unwinding, which proceeds bidirectionally from internal damage sites. After stopping the unwinding process at a specified time, a dye is added that interacts with the remaining double stranded DNA and during contact with DNA emits a fluorescence signal whose intensity is inversely related to the number of DNA strand breaks originally present: the less fluorescence signal, the more DNA strand breaks.
When X-radiation is applied to living cells, the number of strand breaks is known to increase linearly with the dose applied ,. However, the fluorescence signal intensity depends on the dose in a nonlinear fashion and displays saturation. Previously, the percentage of double stranded DNA, remaining after the alkaline unwinding, relative to undenatured control or more precisely the decadic logarithm of the intensity ratio with and without unwinding has been used for that purpose, but this does not fully result in a linear function with radiation dose. In this article, we present a mathematical model that describes the effect of DNA unwinding on the resulting strength of the fluorescence signal and thereby captures the non-linear relationship more precisely.
The main assumptions of the model are large cell numbers, homogeneously distributed damage sites, and small variation of cell parameters like initial damage, susceptibility to applied damage, and potential unwinding distance from the site of damage. The resulting formula relates the relative fluorescence intensity to the applied X-ray dose and fits well with experiments. To obtain this result, it suffices to use elementary arguments from probability theory. In particular, a detailed biochemical model of the unwinding process as in  is not required.
We consider N cells of identical type, which are subjected to a controlled DNA damage by a specified dose of X-radiation. After the treatment, the DNA in cell i exhibits a certain number of (single or double) strand breaks. We denote this damage D i and split it into a pre-existing, initial contribution and an induced damage which is proportional to the applied dose d . In mathematical terms, the FADU assay assigns to each D i a corresponding fluorescence signal F i of the cell and the total fluorescence of all cells is measured. Ideally, the fluorescence F i is related, after completion of the FADU assay, to the length of single and double stranded DNA in cell i according to
HereL i is the relative length of double stranded DNA in the cell so that 1-L i is the relative length of DNA, which has been unwound. In particular, c is the fluorescence signal that is obtained in the case of no unwinding L i = 1, and cB the one after complete unwinding L i = 0 (we call B the relative background fluorescence). Normalizing the total fluorescence F with the one in the case of no unwinding F 0 = Nc, we obtain
In other words, the relative fluorescence intensity I depends on the average relative length of double stranded DNA and the relative background fluorescence B under total unwinding. Later, the fact that N is large allows us to apply the mathematical law of large numbers to safely replace the average by an expected value.
It remains to model the relative length L i of double stranded DNA after the unwinding process and its relation to the damage D i . Since we assume that the process happens identically and independently in each cell, we consider a generic situation and drop the cell index i in the following discussion.
Our basic assumption is that the break points x k are uniformly distributed along the DNA strands and that the DNA unwinds a distance δ to the left and to the right of x k unless there is interference with the unwinding process going on at a neighboring break.
Since unwinding also starts at the ends of the chromosomes, it is useful to unify the description by hypothetically stringing the various chromosomes together. As far as unwinding is concerned, this does not change the result as long as the points which mark the relative stringing positions are added to the list of break points. This trick leads to an elegant description of the unwinding process: We start from m internal strand breaking points 0 < z1 <...... < z m < 1 and associate to each z k the two neighboring points that mark the potential unwinding interval which is also the actual unwinding if there is no interference with neighboring intervals. Obviously, such an interference appears between z k and z k+1 if and only if . Altogether, the contribution to the relative length L of double stranded DNA between z k and z k+1 is where u + = max(u, 0) denotes the positive part of a number u. Introducing the auxiliary points z 0 = 0, z m + 1 = 1 and summing up all the contributions, we obtain
where Δ k = z k + 1 -z k is the distance between consecutive break points. For m we use the splitting m + 1 = D 0 + D x into the number of breaks D x = βd due to X-ray damage and the number D 0 accounting for breaks at zero damage, resulting from normal cell metabolism. Note that in our model, D 0 is at least the number of chromosomes because we have added their terminal points as artificial breaks. Due to the stochastic nature of the damaging process, the cell metabolism and the unwinding, the number of strand breaks per dose β, the zero dose value D 0 and the unwinding distance δ should be considered variables whose value may vary randomly from cell to cell.
To estimate the resulting distribution of the relative fluorescence signal ƒ with some manageable expression it seems worthwhile to simplify the model further, even if it requires additional assumptions. For example, the situation simplifies a lot if we assume that all z 1, , z m are independently and uniformly distributed in the interval [0,1], which means that the particular role of the original terminal chromosome points is dropped. Since in a normal diploid human cell, the number of chromosomes is only 46, this may be acceptable as long as m is comparably large.
According to , the uniform distribution of the points z k implies that the distances Δ k = z k+1 -z k , k = 0,…, m between the points are independent and identically distributed like where U is uniformly distributed on [0,1]. In other words, Δ (m) has a probability density m(1-x)m -1 for x ∈ [0, 1]. This result allows us to compute the conditional expectation
For a given value δ we have in the case 2δ < 1
while E(Δ (m)- 2δ)+ = 0 when 2δ≥1. Combining this result with the previous formula, we find the conditional expectation
In order to compute the unconditional expectation E(L), reasonable assumptions on the probability distributions of D 0, β and δ are required. Finally, E(L) depends in a complicated way on the dose d, the mean values , the corresponding variances and maybe other parameters of the distributions.
However, since little information on the distributions of D 0, β and δ is available, we continue in a different way. We assume that the unconditional expectation E(L) can be approximated with an expression of the form
which is a reasonable assumption when D 0, β and δ have negligible variance among all cells in the assay. Then, the effective values and δ eff are close to the mean values .
Replacing the average in the formula for relative intensity I by the approximate expectation, we obtain
More compactly, this relationship can be written in the form
where the relation between the fit parameters P 0,λ and the model parameters is
Once B, P 0,λ are determined, the equivalent dose value d corresponding to the relative intensity P can easily be determined by inversion
Unfortunately, an estimate of which is essentially the total number of strand breaks is not available if only B, P 0,λ are known. However, a related dose-equivalent value
is accessible. In fact, applying the logarithm to (P 0 - B)/(1-B) and using the relation for λ, we find
Together with the expression for the dose value d, we obtain a formula that relates P to the value D tot .
In order to test our model, we use experimental results from the experiment specified in Material and Methods. Here, three replicates of the relative intensities at seven dose values are available for eleven donors. We denote them P jkl where the index j = 0,..,6 labels the X-ray dose, k = 1,..., 11 represents the donor and l = 1,2,3 indicates the independent repetitions of the experiment.
While each donor may exhibit different parameters P 0k , λ k , the background intensity B can be assumed identical in all experiments. To determine the unknown coefficients, we want to minimize the objective function
which leads to the same result as the minimization of
based on the averages
To quantify the quality of the fit, we check how much of the total sum of squares
is explained by the fitting which leads to an R 2-value
where is the residual in case of the optimal parameters. The optimization is carried out with the Levenberg-Marquardt algorithm of the Matlab software package with a resulting value of R 2 = 0.99.
Confidence intervals for the fit parameters are computed with the bootstrap method  exploiting the information contained in the threefold repetition of each experiment. For each index combination (j, k, l), N = 10000 bootstrap samples are generated by drawing intensities randomly from the values P jk1, P jk2, P jk3 with replacement. Then the computation of the optimal parameters is repeated for each n = 1,..., N and the resulting distribution is used to construct 95%-confidence intervals.
The resulting value for B is
DNA strand breaks are one of the most common genotoxic lesions. One of the methods used for measuring DNA strand breaks and their repair is the FADU assay. Under the selected alkaline conditions a large number of base pairs will be unwound to the right and left starting from each DNA "open side". Due to this amplification it is possible to detect very low number of DNA strand breaks. The main advantage of the FADU assay is that it can be performed in a fully automated fashion. However there are few methods that can be performed in a semi-automated version such as the COMET and γH2AX assays. Even though these techniques provide useful tools for measuring DNA strand breaks in a high-throughput fashion, they have some disadvantages. A comprehensive comparison between these assays has been published before . Briefly, γH2AX assay is considered a very sensitive method that specifically detects double DNA strand breaks however in some cases the presence of γH2AX foci in the absence of DNA damage has been demonstrated. The comet assay protocol includes many different steps making difficult the development of a fully automated version. Furthermore these also affect both intra-assay variability and inter-assay reproducibility.
We have previously published a modified and automated version of the FADU assay ,,. Compared with the original FADU assay , the automated FADU assay shows comparable sensitivity, yet increased robustness and throughput, and decreased operator time. However a suitable mathematical model, which allows translating the fluorescence intensity in DNA damage increases accuracy. In order to test whether the model is capable of describing the experimentally observed dependence between applied damage and resulting fluorescence intensity, we have measured the DNA strand breaks in peripheral mononuclear blood cells of 11 individuals. It turns out that the model is able to fit the dose-effect relation very accurately.
As shown in Table1, the resulting model parameters λ and P 0, which describe the cell's susceptibility to DNA damage and the level of endogenous DNA strand breaks respectively, show a certain variation among the individuals which is stronger for the P 0 value. An individual endogenous level of DNA strand breaks is expected; in general the importance of endogenous DNA damage has been reported before . For example, it is known that DNA double-strand breaks (DSBs) accumulate in senescing human cell cultures and in ageing mice . This is in accordance with our recent data showing that senescing human T lymphocytes have acquired more DNA strand breaks in vivo and an impaired DNA repair capacity in vitro. Additionally, endogenous DNA strand breaks accumulation in human cells has been associated with nutrition , cancer , diabetes ,, Down's syndrome , rheumatoid arthritis  and psychological stress disorders , among others. In other words, differences in the level of endogenous DNA strand breaks might reflect aging processes or health status of individuals.
The susceptibility to DNA damage λ also shows variation between individuals. This finding is not surprising and has been described before. Odarigi and colleagues reported an inter-individual variance in chromosomal damage after x-irradiation . Since radiation induces free oxygen radical formation, an individual damage response to x-ray might result from the different capacity of cells to scavenge these free radicals.
The investigation of DNA strand breaks and their repair is of great interest to toxicology, pharmacology, epidemiology and medical sciences at large. In general impaired DNA repair can lead to embryonic lethality, impaired growth, accelerated aging, shortened life span, and increased incidence of a variety of diseases, including a pronounced manifestation of cancer -. Therefore robust methods for quantification of DNA damage and its repair are necessary.
Finally, in our model, induced DNA strand breaks are assumed to be uniformly distributed. However chromatin state could have an impact on the distribution of DNA strand breaks. Chromatin state differs depending on cell type and cell status (cell cycle phase and/or transcription) therefore it is important to take these aspects into consideration. A mathematical model of DNA unwinding considering the chromatin state would be an interesting task for future research, especially when comparing different cell types or cell status. A more refined model could be based on a probability distribution, which describes the susceptibility to breaking under X-ray exposure along the DNA strand.
We have applied the above mathematical model for the quantification of DNA strand breaks in irradiated human PBMCs. This model can be also applied to other cell types and other DNA damaging agents as long as the damaging agent leads to uniformly distributed DNA strand breaks. Furthermore measuring the residual DNA strand breaks in damaged cells after several periods of time allow us to calculate the rate of repaired DNA strand breaks. Last but not least the amount of endogenous DNA strand breaks provides important information about the individual "cell-damage-status" and fluorescence signal after induced damage serves as biomarker of cell sensitivity to DNA strand breaks inducers, therefore an accurate model for DNA breakage quantification could be used for diagnostic and therapy purposes towards a personalized medicine, e.g. in cancer therapy.
Selection of volunteers was carried out in accordance with the Declaration of Helsinki and ethical approval was obtained from the Ethics Committee of University of Constance. A signed Informed Consent was obtained from each subject. Blood from healthy volunteers aged 23-55 was collected in sodium citrat tubes (Sarstedt, Nümbrecht, Germany). The automated FADU assay has previously been described in detail . Briefly, peripheral blood mononuclear cells (PBMC) were isolated by Biocoll (Biochrom, Berlin, Germany) density gradient centrifugation. Cells were counted using a cell counting device (Casy® counter), pelleted (10 min, 200 g), and resuspended in isotonic buffer (0.25 M meso-inositol; 10 mM sodium phosphate, pH 7.4; 1 mM magnesium chloride) or RPMI-1640 medium (Invitrogen) containing 100 U/ml penicillin (Invitrogen), 100 Mg/ml streptomycin (Invitrogen) and 10% fetal calf serum (FCS) at 4 × 106 cells per ml. In order to induce DNA strand breaks several aliquots of 100 l cell suspension were irradiated on ice at a dose rate of 1.9 Gy/min using an X-ray generator (CHF Müller, Hamburg, Germany, 70 keV, 1 mM Al-filter). Qualitative behaviour has been modelled, but it still depends on few quantitative parameters, which where determined based on the results obtained from unwinding experiments after various (0.3, 1, 1.9, 3.8, 7.5 and 14.9 Gy) irradiation doses (Figure 1).
The detection of DNA strand breaks is based on progressive DNA unwinding (denaturation) under highly controlled conditions of alkaline pH, time and temperature. The starting points for the unwinding process are DNA "open sites" like replication forks or chromosome ends, but also DNA strand breaks. For monitoring DNA damage a commercially available fluorescent dye (SYBR Green®) is used as marker for double stranded DNA, and a decrease in the fluorescence intensity of SYBR Green® indicates an increase of DNA unwinding and consequently a greater number of DNA strand breaks.
T and P0 are controls to be run in parallel with the experimentally treated cells. In P0 samples alkaline unwinding is allowed and represents the DNA strand breaks under physiological conditions (i.e. without exogenous DNA damage). In T samples the neutralisation follows the lysis and therefore unwinding cannot take place; T-samples provide a measure of total DNA content and yield a fluorescence signal defined as 100%. Px samples (P1,P2,P3,P4....Px) are the different extent of damage to be measured.
The work presented here was carried out in collaboration between all authors. MMV defined the research theme and designed methods and experiments. JS carried out the laboratory experiments and analyzed the data. TS co-worked on data interpretation and discussed analyses. MJ developed the mathematical model and carried out the parameter fitting. AB co-designed experiments, discussed analyses, interpretation, and presentation. MMV and MJ wrote the paper. All authors have contributed to, seen and approved the manuscript.
Moreno-Villanueva M, Pfeiffer R, Sindlinger T, Leake A, Müller M, Kirkwood T, Bürkle A: A modified and automated version of the `Fluorimetric Detection of Alkaline DNA Unwinding' method to quantify formation and repair of DNA strand breaks. BMC Biotechnol. 2009, 9: 39-10.1186/1472-6750-9-39.
Birnboim HC, Jevcak JJ: Fluorometric method for rapid detection of DNA strand breaks in human white blood cells produced by low doses of radiation. Cancer Res. 1981, 41: 1889-1892.
Wood DK, Weingeist DM, Bhatia SN, Engelward BP: Single cell trapping and DNA damage analysis using microwell arrays. PNAS. 2010, 1;107 (22): 10008-10013. 10.1073/pnas.1004056107.
Olive PL, Banáth JP: The comet assay: a method to measure DNA damage in individual cells. Nat Protoc. 2006, 1 (1): 23-29. 10.1038/nprot.2006.5.
Gaver DP, Jacobs PA, Carpenter RL, Burkhart JG: A mathematical model for intracellular effects of toxins on DNA adduction and repair. Bull Math Biol. 1997, 59 (1): 89-106. 10.1007/BF02459472.
Sokhansanj BA, Rodrigue GR, Fitch JP, Wilson DM: A quantitative model of human DNA base excision repair. I. Mechanistic insights. Nucleic Acids Res. 2002, 30 (8): 1817-1825. 10.1093/nar/30.8.1817.
Kesseler KJ, Kaufmann WK, Reardon JT, Elston TC, Sancar A: A mathematical model for human nucleotide excision repair: damage recognition by random order assembly and kinetic proofreading. J Theor Biol. 2007, 249 (2): 361-375. 10.1016/j.jtbi.2007.07.025.
Crooke PS, Parl FF: A mathematical model for DNA damage and repair. J Nucleic Acids. 2010, 2010: 352603-10.4061/2010/352603.
Devroye L: Non-Uniform Random Variate Generation. 1987, Springer Verlag, New York
Efron B, Tibshirani RJ: An Introduction to the Bootstrap. Monographs on Statistics and Applied Probability 57. 1994, Chapman&Hall/CRC,
Moreno-Villanueva M, Bürkle A: High-throughput Assays to Quantify the Formation of DNA Strand Breaks. High Throughput Screening Methods in Toxicity Testing. Edited by: Steinberg P. 2013, John Wiley & Sons, Inc, Hoboken, NJ, USA, 285-294. 10.1002/9781118538203.ch14.
Moreno-Villanueva M, Eltze T, Dressler D, Bernhardt J, Hirsch C, Wick P, von Scheven G, Lex K, Bürkle A: The automated FADU-assay, a potential high-throughput in vitro method for early screening of DNA breakage. ALTEX. 2011, 28 (4): 295-303. 10.14573/altex.2011.4.295.
De Bont R, van Larebeke N: Endogenous DNA damage in humans: a review of quantitative data. Mutagenesis. 2004, 19 (3): 169-185. 10.1093/mutage/geh025.
Sedelnikova OA, Horikawa I, Zimonjic DB, Popescu NC, Bonner WM, Barrett JC: Senescing human cells and ageing mice accumulate DNA lesions with unrepairable double-strand breaks. Nat Cell Biol. 2004, 6 (2): 168-170. 10.1038/ncb1095.
Brunner S, Herndler-Brandstetter D, Arnold CR, Wiegers GJ, Villunger A, Hackl M, Grillari J, Moreno-Villanueva M, Bürkle A, Grubeck-Loebenstein B: Upregulation of miR-24 is associated with a decreased DNA damage response upon etoposide treatment in highly differentiated CD8(+) T cells sensitizing them to apoptotic cell death. Aging Cell. 2012, 11 (4): 579-587. 10.1111/j.1474-9726.2012.00819.x.
Krajcovicová-Kudlácková M, Dusinská M: Oxidative DNA damage in relation to nutrition. Neoplasma. 2004, 51 (1): 30-33.
Vilenchik MM, Knudson AG: Endogenous DNA double-strand breaks: production, fidelity of repair, and induction of cancer. Proc Natl Acad Sci. 2003, 100 (22): 12871-12876. 10.1073/pnas.2135498100.
Lorenzi M, Montisano DF, Toledo S, Wong HC: Increased single strand breaks in DNA of lymphocytes from diabetic subjects. J Clin Invest. 1987, 79 (2): 653-656. 10.1172/JCI112863.
Blasiak J, Arabski M, Krupa R, Wozniak K, Zadrozny M, Kasznicki J, Zurawska M, Drzewoski J: DNA damage and repair in type 2 diabetes mellitus. Mutat Res. 2004, 554 (1-2): 297-304. 10.1016/j.mrfmmm.2004.05.011.
Zana M, Szécsény A, Czibula A, Bjelik A, Juhász A, Rimanóczy A, Szabó K, Vetró A, Szucs P, Várkonyi A, Pákáski M, Boda K, Raskó I, Janka Z, Kálmán J: Age-dependent oxidative stress-induced DNA damage in down's lymphocytes. Biochem Biophys Res Commun. 2006, 345 (2): 726-733. 10.1016/j.bbrc.2006.04.167. Epub 2006 May 5
Bhusate LL, Herbert KE, Scott DL, Perrett D: Increased DNA strand breaks in mononuclear cells from patients with rheumatoid arthritis. Ann Rheum Dis. 1992, 51 (1): 8-12. 10.1136/ard.51.1.8.
Morath J, Moreno-Villanueva M, Gilava H, Kolassa S, Ruf-Leuschner M, Schauer M, Elbert T, Bürkle A, Kolassa I-T: Effects of psychotherapy on DNA strand break accumulation originating from traumatic stress. Psychother Psychosom. 2014, 83 (5): 289-297. 10.1159/000362739.
Odagiri Y, Uchida H, Shibazaki S: Interindividual variation in cytogenetic response to X-ray and colchicine measured with the cytokinesis-block micronucleus assay. Mutat Res. 1997, 381 (1): 1-13. 10.1016/S0027-5107(97)00138-3.
Hakem R: DNA-damage repair; the good, the bad, and the ugly. EMBO J. 2008, 27: 589-605. 10.1038/emboj.2008.15.
Diderich K, Alanazi M, Hoeijmakers JH: Premature aging and cancer in nucleotide excision repair-disorders. DNA Repair. 2011, 10 (7): 772-780. 10.1016/j.dnarep.2011.04.025.
Keimling M, Deniz M, Varga D, Stahl A, Schrezenmeier H, Kreienberg R, Hoffmann I, König J, WiesMüller L: The power of DNA double-strand break (DSB) repair testing to predict breast cancer susceptibility. FASEB J. 2012, 26 (5): 2094-2104. 10.1096/fj.11-200790.
Varga D, Deniz M, Schwentner L, WiesMüller L: Ovarian cancer: in search of better marker systems based on DNA repair defects. Int J Mol Sci. 2013, 14 (1): 640-673. 10.3390/ijms14010640.
Toulany M, Rodemann HP: Membrane receptor signaling and control of DNA repair after exposure to ionizing radiation. Nuklearmedizin. 2010, 2010 (49 Suppl 1): S26-S30.
This work was supported by the European Union FP7 Health research, grant number HEALTH-F4-2008-200880 [MARK-AGE]. The authors would like to acknowledge Ms. Gudrun von Scheven for organizing the blood collection.
The authors declare that they have no competing interests.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Junk, M., Salzwedel, J., Sindlinger, T. et al. Mathematical modelling of the automated FADU assay for the quantification of DNA strand breaks and their repair in human peripheral mononuclear blood cells. BMC Biophys 7, 9 (2014). https://doi.org/10.1186/s13628-014-0009-z
- DNA repair
- DNA strand breaks
- Mathematical model