About this Journal Submit a Manuscript Table of Contents
ISRN Biophysics
Volume 2012 (2012), Article ID 725350, 6 pages
http://dx.doi.org/10.5402/2012/725350
Research Article

Fragmentation of Plasmid DNA Produced by Gamma Radiation: A Theoretical Approach

1Physics Institute, Rua do Matão, Travessa R, No. 187, University of São Paulo, 05508-090 São Paulo, SP, Brazil
2FESP, São Paulo Engineering School, São Paulo, SP, Brazil
3Laboratory of Devices Research, University of Campinas, Campinas, SP, Brazil

Received 16 August 2011; Accepted 18 September 2011

Academic Editor: E. Dague

Copyright © 2012 R. A. S. Silva et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Breaks in DNA, resulting in fragmented parts, can be produced by ionizing radiation which, in turn, is the starting point in the search for novel physical aspects of DNA strands. Double-strand breaks in particular cause disruption of the DNA strand, splitting it into several fragments. In order to study effects produced by radiation in plasmid DNA, a new simple mechanical model for this molecule is proposed. In this model, a Morse-like potential and a high-LET component are used to describe the DNA-radiation interaction. Two power laws, used to fit results of the model, suggest that, firstly, distribution of fragment size is nonextensive and, secondly, that a transition phase is present in the DNA fragment distribution pattern.

1. Introduction

One of the most important molecules for living beings is DNA, and as a result it has been widely studied. All genetic information controlling cellular reproduction is in DNA, and this is the most sensitive molecule in cell nuclei [13]. The interaction of radiation with organic molecules and the consequences of this interaction for biological systems have been thoroughly investigated.

Ionizing radiations such as gammas, neutrons, or alpha particles can produce breaks in DNA resulting in fragmented parts that would be the starting point in the search for novel physical aspects of DNA strands. Among these, nonextensive properties (nonlinearities) at nanoscale as a likely possibility would be mentioned. In fact, some evidence in this direction could be inferred from [4].

It is commonly accepted that secondary electrons formed in the process of ionization are mostly responsible for DNA damage, either by directly breaking the DNA strands, or by reacting with water molecules, producing more secondary electrons and free radicals, which, in turn, could provoke additional damage to DNA. Many different kinds of lesions in DNA are caused by ionizing radiation, more especially base damage and strand breaks. The latter are particularly important because they are irreparable when produced in large and clustered quantities [5]. According to experimental evidence, initial-level breaks and residual/unrepaired-level breaks, can be correlated with the various biological end-points such as loss of reproductivity, cell death, capacity, mutations, or carcinogenesis [6]. Nowadays, DNA-radiation interaction is, however, not totally understood. A detailed study of damage produced by radiation is very important for maintenance of life. It is necessary to stress that DNA damage depends on the kind of radiation (alpha, beta, or gamma) and also that electrons are fundamental in this damage mechanism [7].

Theoretical studies in radiation biology with high-LET charged particles have been helpful in understanding some types of DNA damage as double-strand breaks [8]. High-LET radiation features a highly structured deposition of energy along the charged particle tracks [9]. Moreover, high-LET radiation is expected to induce more severe and complex lesions correlated along individual particle tracks because of energy deposition characteristics [911]. A high quantity of biological damage is induced at sites close to particle trajectory, and damage clustering is caused by this energy localization. Double-strand breaks in particular cause disruption of the DNA strand, separating it into several fragments. Size distribution pattern of these fragments could reveal important aspects of the breaking mechanisms as discussed in this paper.

Presently, DNA fragmentation can be easily quantified through atomic force microscopy (AFM). This technique has been largely used in biological science [12], combining high spatial resolution with the ability to image biological samples in their physiological environment. Radiobiological research uses this powerful AFM imaging capability [13]. In fact, the investigation of DSBs induced by X-rays, neutrons, electrons, and alpha particles in plasmid DNA has been performed with AFM [1416]. Results on damage induced in plasmid DNA by heavy ions have also been reported elsewhere [17].

Results for fragment size distribution functions with electrons and neutrons were obtained by Pang and collaborators [16]. As discussed in this paper, their results present leads for the occurrence of nonextensivity and crossover at nanoscale. In this regard, it is noted that there are several studies in the literature showing nonlinear aspects of DNA molecule [18]. For instance, studies of long-range correlations in DNA indicate that the length (𝐿) of homogeneous sequences in noncoding regions obeys a power-law-type distribution (a fingerprint for nonextensiveness), while in widely encoded regions these distributions are exponential [4].

Thus, here, a simple mechanic model for DNA is proposed, aiming at the understanding of DNA physical characteristics when exposed to different radiation intensities (doses), particularly aspects associated with its fragmentation pattern, for example, the occurrence, or not, of non-extensivity at molecular level. The model simulates DNA strands by spring-mass lines. These masses, representing nucleotides, are connected by two components of elastic potential—one is the conventional linear harmonic oscillator and the other the nonlinear Morse potential. The former holds masses together, and the latter takes into account the interaction between masses and radiation particles. The reason for choosing a Morse potential is because of its ability to simulate stretching and compression of covalent bonds [19, 20].

The results of our model calculations are compared with those from experiments carried out with gammas in our laboratory, and with electrons and neutrons from the literature [21].

2. The Model

In the present model, each DNA base pair is represented by a mass 𝑚. These masses are joined together along a strand by spring potentials with coupling constants 𝑘. Furthermore, there is a nonlinear potential relating the distance between masses and charged particles. Base pair displacements are only longitudinal; that is, they would be along the helix axis. Each mass (𝑚𝑖) has a linear displacement, denoted by 𝑥, while 𝑟𝑖 is the approximate distance of the radiation particles to 𝑚𝑖.

We propose the DNA-radiation interaction leading to the breaking of covalent bonds be given by a Morse-like potential 𝑉(𝑥𝑖,𝑟𝑖). This potential is generally chosen to simulate covalent [19, 20] bonds and hydrogen bonds in small molecules [22, 23].

For our purposes, it is assumed that the potential depends on a combination of variables 𝑥𝑖 and 𝑟𝑖, that is, 𝑉(𝑥𝑖+𝑔(𝑟𝑖)) thus, the Morse-like potential is given by𝑉𝑥𝑖𝑟+g𝑖=𝐷𝑀𝑒𝑎(𝑥𝑖+𝑔(𝑟𝑖))12,(1) where 𝐷𝑀 is the dissociation energy of a base pair in a given strand and 𝑎 is a parameter with dimension of inverted length. 𝑔(𝑟), a function defined by 𝐾𝑔(𝑟)=𝑟,(2) is introduced to consider interactions between charged particle approximation and DNA. In this case, 𝐾=𝐷LET𝐶, where 𝐷 is related to radiation dose, LET is the average linear energy transfer [24], and 𝐶=(3.43/7.504) Å3 eV−1 (3.4 Å is the average distance between the adjacent bases in the same DNA strand, and 7.504 eV is the necessary energy to break a covalent bond [25].

The Hamiltonian for the system described above can be written as𝐻=𝐻𝑝+𝐻𝑣,(3) where𝐻𝑝=𝑖𝑝2𝑥𝑖2𝑚(4) represents the kinetic energy of the system, and𝐻𝑉=𝑖𝑘2𝑥𝑖𝑥𝑖12𝑥+𝑉𝑖,𝑟𝑖(5) is the potential energy. In (2) 𝑝𝑥 is the linear momentum. For the sake of simplicity, the system is considered homogenous, meaning that masses 𝑚 and strength constants 𝑘 are the same for all nucleotides. The ionizing radiations (alpha, beta, and gamma) can be defined by a LET component [24] in the Morse potential.

For a chain with 𝑁 base pairs, the classical partition function is given in terms of the Hamiltonian introduced in (1) and may be factored as𝑍=𝑁𝑛=1𝑑𝑥𝑛𝑑𝑝𝑛𝑒𝐻/𝑘𝐵𝑇=𝑍𝑥𝑍𝑝.(6) These integrals for all momenta are merely Gaussian integrals, providing𝑍𝑝=2𝜋𝑚𝑘𝐵𝑇𝑁/2,(7) where 𝑘𝐵 is the Boltzmann constant. The term in the 𝑥 variables, 𝑍𝑥, comes from a harmonic chain of oscillators with an on-site nonlinear potential for each oscillator. Therefore, partition function 𝑍𝑥 is written as𝑍𝑥=𝑁𝑛=1𝑑𝑥𝑛𝑒𝑓(𝑥𝑖,𝑥𝑖1,𝑟𝑖)/𝑘𝐵𝑇,(8) where 𝑓 is given by𝑓𝑥𝑖,𝑥𝑖1,𝑟𝑖=𝑘2𝑥𝑖𝑥𝑖12𝑥+𝑉𝑖,𝑟𝑖.(9)

Integrals in 𝑍𝑥 can be calculated using eigenfunctions and eigenvalues of a transfer integral operator [26, 27]:𝑑𝑥1𝑒𝑓(𝑥𝑖,𝑥𝑖1,𝑟𝑖)/𝑘𝐵𝑇𝜓𝑖𝑥𝑖1=𝑒𝜀𝑖/𝑘𝐵𝑇𝜓𝑖𝑥𝑖.(10)

In the thermodynamic limit, 𝑁, the partition function 𝑍𝑥 is reduced to𝑍𝑥=𝑒𝑁𝜀0/𝑘𝐵𝑇,(11) where 𝜀0 is the ground-state eigenvalue for a Schrödinger-like equation given by𝑘𝐵𝑇2𝑑2𝑘2𝑑𝑥2𝑖𝑥+𝑉𝑖𝑟,𝑔𝑖𝜓0𝑥𝑖=𝜀0𝜎0𝜓0𝑥𝑖(12) with𝜎0=12𝑘𝐵𝑘𝑇ln2𝜋𝑘𝐵𝑇.(13)

The eigenvalue and the normalized eigenfunction for the ground state are, respectively,𝜀0=𝜎0+𝑘𝐵𝑇𝑎𝐷𝑘2𝑘𝐵𝑇𝑎28𝑘,(14)𝜓0𝑥𝑖=𝐶𝑒[𝛿𝑒𝑎(𝑥+𝑔(𝑟))𝑎𝑥(𝛿0.5)](15) with 𝛿=(2𝑘𝐷)1/2/(𝑘𝐵𝑇𝑎) and 𝐶 the normalization constant. In the limit 𝑁, the average stretching 𝑥 will be dependent only on the ground-state eigenfunction as𝑥=𝜓02(𝑥)𝑥𝑑𝑥.(16)

It is suggested that the mean 𝑥 could be related to the number of fragments produced by radiation. In the absence of any external interference, the average distance between neighboring nucleotides in the same DNA strand is 3.4 Å, and for a linearized DNA with 𝐿=1000 nm, the number 𝑁 of links involving the neighboring nucleotides is given by 𝑁=1000 nm/(3.4 Å) = 2941. It is considered that the breaking probability of DNA is given by 𝑥/𝑟. Therefore, we can relate this rate to the value of 𝑁 and 𝐿 given above to obtain the distribution of the number of fragments and their size according to radiation dose. In this case, the number of fragments and length distribution are given, respectively, by𝑁=𝑁𝑥𝑟𝐿,𝐿=𝑥/𝑟.(17)

3. Results and Discussion

The best values for the parameters obtained from experimental results for DNA presented in the literature as 𝑎=2.81 Å and 𝑘=0.06 eV Å−2 [28]. Covalent bonds are broken in radiation-DNA interaction, and energy values for these breaks are already known. For example, the energy to break a double covalent bond between oxygen and carbon atoms is around 724 kJ/mol [brad](or 7504 eV/mol). Therefore, for a preliminary theoretical investigation, it is assumed that 𝐷𝑀=7504eV. The value of LET used in this work is 19.6 keV/μm for gamma radiation [24]. The qualitative behavior of the average stretching 𝑥 as function of the approximation distance 𝑟 is obtained from (16) and shown in Figure 1.

725350.fig.001
Figure 1: Average stretching between neighboring bases pairs as function of 𝑟 for three radiation dose: low (continuous curve), moderate (+) and high (square).

It is observed that (1) the low value for 𝐷 (continuous curve in Figure 1) corresponds to the lowest damage level of gamma radiation. The average values for 𝑥 (𝑥) are very low suggesting that only a few breaks were produced along the DNA chain; (2) the curve drawn with + symbols corresponds to a higher dose relative to the previous continuous curve, and a faster increase of 𝑥 is observed, indicating occurrence of a higher damage level; (3) in the third case (squared symbol curve), factor 𝐷 was considerably increased resulting in a steeper increase of 𝑥, which could be associated with longer exposure to radiation. Therefore, large values 𝑥 could indicate several breaks in the DNA sequence.

The role played by radiation dose in DNA strand breakage may also be seen in Figure 2 which shows the dependence of 𝑥 on radiation dose. In the case of high doses of radiation, with short interaction distance (𝑟=5 Å), the values of 𝑥 tend to be larger than the long distances (𝑟=10 Å), indicating a higher frequency of DNA strand breakage. Moreover, only interactions of short distance (𝑟 much less than 5 Å) produce high values of 𝑥 at low dose.

725350.fig.002
Figure 2: Variation of the average stretching as function of radiation dose, for 𝑟=5 and 𝑟=10 Å.

Figure 3 shows the distribution of DNA fragments as a function of radiation dose. The changes are clearly visible since the larger amount of smaller fragments is produced with longer exposure to radiation, and an increase in the number of small fragments with increasing radiation dose is observed. What is more, the results in Figure 3 show a behavior qualitatively like that of a system with damped oscillations, characteristic of a fully homogeneous system.

fig3
Figure 3: Number of fragments versus their length. Graphs (a), (b), and (c) refer to the low, medium, and high dose, respectively.

Very similar behavior observed for the results shown in Figure 4 was obtained by Pang and collaborators for distribution functions of fragment size using neutrons and electrons [21]. The curves shown were obtained by a simple fitting of power laws where it was found that the exponents are fractional. Two functions of power law type to describe the data were necessary, strongly suggesting the occurrence of a crossover in nanoscale. Two functions used to fit the data in Figure 3(a), and their values are shown in Table 1.

tab1
Table 1: Fit of data in Figure 3(a), using two power laws. Parameter values are, 𝑎=25, 𝑏=0.1, 𝑐=6.5, 𝑑=13, 𝛼1 = 0.4 and 𝛼2 = 0.1.
725350.fig.004
Figure 4: Size distribution of DNA fragments. Full circles refer to the histogram (a) of Figure 3; triangles and squares were calculated from the expressions 𝑝(𝐿)=𝑐+𝑑𝐿𝛼1 and 𝑝(𝐿)=(𝑎+𝑏𝐿)𝛼2, respectively. The lines (dotted and continuous) are to guide the eyes only.

4. Final Remarks and Conclusions

The development of biophysical models requires a simplified description of structural properties of the elements involved, as carried out in the model developed in this paper. With the help of statistical mechanics tools, the model used here to simulate the DNA-radiation interaction has yielded results qualitatively consistent with those obtained experimentally: the higher number of smaller fragments is produced with the longer exposure to radiation, that is, with higher radiation doses.

The DNA strand behaved as a set of coupled oscillators subject to a purely harmonic interaction between each oscillator and its first neighbors. The average stretch 𝑥 between the base pairs would be the elongation of oscillation and determines the greater or lesser propensity of the system to break.

The shape of the curve representing the number of fragment size as a function of radiation dose, as shown by the model developed here (Figure 3), is characteristic of a system with damping. Two power laws used to fit results of the model suggest a nonextensive distribution of fragment size, this being a precondition for fractalization. Furthermore, obtaining two power laws for the same dose strongly suggests a transition phase in the DNA fragment distribution pattern.

It is pertinent to stress that the model does not take into account the inhomogeneity of the DNA molecule, and this property can affect the process of fragmentation; that is, regions rich in AT (base adenine attached to base thymine) could break more easily than GC-rich regions (base guanine connected to base cytosine) since the former are linked by two hydrogen bonds while GC have three hydrogen bonds. These aspects are to be considered in a forthcoming study.

Acknowledgment

This work was supported by FAPESP and CNPq (Brazilian agencies).

References

  1. S. H. Strogatz, “Quoted by S. Nadis in “News feature: all together now”,” Nature, vol. 421, pp. 780–782, 2003.
  2. C. Tsallis, “Possible generalization of Boltzmann-Gibbs statistics,” Journal of Statistical Physics, vol. 52, no. 1-2, pp. 479–487, 1988. View at Publisher · View at Google Scholar · View at Scopus
  3. C. Tsallis, S. V. F. Levy, A. M. C. Souza, and R. Maynard, “Statistical-mechanical foundation of the ubiquity of lévy distributions in nature,” Physical Review Letters, vol. 75, no. 20, pp. 3589–3593, 1995. View at Publisher · View at Google Scholar · View at Scopus
  4. P. Bernaola-Galván, R. Román-Roldán, and J. L. Oliver, “Compositional segmentation and long-range fractal correlations in DNA sequences,” Physical Review E, vol. 53, no. 5, pp. 5181–5189, 1996. View at Scopus
  5. U. Amaldi and G. Kraft, “Radiotherapy with beams of carbon ions,” Reports on Progress in Physics, vol. 68, no. 8, pp. 1861–1882, 2005. View at Publisher · View at Google Scholar · View at Scopus
  6. E. L. Alpen, Radiation Biophysics, Academic Press, San Diego, Calif, USA, 1998.
  7. Z. G. Yu and X. Song, “Variable range hopping and electrical conductivity along the DNA double helix,” Physical Review Letters, vol. 86, no. 26 I, pp. 6018–6021, 2001. View at Publisher · View at Google Scholar · View at Scopus
  8. A. Chatterjee, “Radiobiological effects of high-let particles: DNA strand breaks,” Nuclear Instruments and Methods in Physics Research A, vol. 280, no. 2-3, pp. 439–448, 1989.
  9. M. Krämer and G. Kraft, “Track structure and DNA damage,” Advances in Space Research, vol. 14, no. 10, pp. 151–159, 1994. View at Scopus
  10. D. T. Goodhead, “The initial physical damage produced by ionizing radiations,” International Journal of Radiation Biology, vol. 56, no. 5, pp. 623–634, 1989. View at Scopus
  11. H. Nikjoo, P. O'Neill, M. Terrissol, and D. T. Goodhead, “Quantitative modelling of DNA damage using Monte Carlo track structure method,” Radiation and Environmental Biophysics, vol. 38, no. 1, pp. 31–38, 1999. View at Publisher · View at Google Scholar · View at Scopus
  12. H. G. Hansma and J. H. Hoh, “Biomolecular imaging with the atomic force microscope,” Annual Review of Biophysics and Biomolecular Structure, vol. 23, pp. 115–139, 1994.
  13. D. Pang, et al., “Atomic force microscopy imaging of DNA and DNA repair proteins: applications in radiobiological research,” Radiation Oncology Investigations, vol. 5, pp. 163–169, 1997.
  14. S. Boichot, M. Fromm, S. Cunniffe et al., “Investigation of radiation damage in DNA by using atomic force microscopy,” Radiation Protection Dosimetry, vol. 99, no. 1–4, pp. 143–145, 2002. View at Scopus
  15. D. Pang, B. L. Berman, S. Chasovskikh, J. E. Rodgers, and A. Dritschilo, “Investigation of neutron-induced damage in DNA by atomic force microscopy: experimental evidence of clustered DNA lesions,” Radiation Research, vol. 150, no. 6, pp. 612–618, 1998. View at Scopus
  16. D. Pang, J. E. Rodgers, B. L. Berman, S. Chasovskikh, and A. Dritschilo, “Spatial distribution of radiation-induced double-strand breaks in plasmid DNA as resolved by atomic force microscopy,” Radiation Research, vol. 164, no. 6, pp. 755–765, 2005. View at Publisher · View at Google Scholar · View at Scopus
  17. K. Psonka, S. Brons, M. Heiss, E. Gudowska-Nowak, and G. Taucher-Scholz, “Induction of DNA damage by heavy ions measured by atomic force microscopy,” Journal of Physics Condensed Matter, vol. 17, no. 18, pp. S1443–S1446, 2005. View at Publisher · View at Google Scholar · View at Scopus
  18. M. Gerard and A. R. Bishop, “Statistical mechanics of a nonlinear model for DNA denaturation,” Physical Review Letters, vol. 62, no. 23, pp. 2755–2758, 1989. View at Publisher · View at Google Scholar · View at Scopus
  19. T. C. Lim, “The relationship between Lennard-Jones (12-6) and Morse potencial function,” Journal of Mathematical Chemistry, vol. 33, p. 1, 2003.
  20. J. A. Greathouse, J. S. Durkin, J. P. Larentzos, and R. T. Cygan, “Implementation of a Morse potential to model hydroxyl behavior in phyllosilicates,” Journal of Chemical Physics, vol. 130, no. 13, Article ID 134713, 2009. View at Publisher · View at Google Scholar · View at Scopus
  21. D. Pang, et al., “Investigation of neutron-induced damage in DNA by atomic force microscopy: experimental evidence of clustered DNA lesions,” Journal of The Radiation Research Society, vol. 150, no. 6, pp. 612–618, 1998.
  22. Y. Gao, K. V. Devi-Prasad, and E. W. Prohofsky, “A self-consistent microscopic theory of hydrogen bond melting with application to poly (dG)·poly(dC),” The Journal of Chemical Physics, vol. 80, no. 12, pp. 6291–6298, 1983.
  23. S. Zdravkovic and M. V. Sataric, “Single-molecule unzippering experiments on DNA and Peyrard-Bishop-Dauxois model,” Physical Review E, vol. 73, Article ID 021905, pp. 1–11, 2006.
  24. E. L. Alpen, Radiation Biophysics, Prentice Hall, NJ, USA, 1990.
  25. J. E. Brady and G. E. Humiston, Química Geral. 2, LTC, Rio de Janeiro, Brazil, 1986.
  26. D. J. Scalapino, M. Sears, and R. A. Ferrell, “Statistical mechanics of one-dimensional ginzburg-landau fields,” Physical Review B, vol. 6, no. 9, pp. 3409–3416, 1972. View at Publisher · View at Google Scholar · View at Scopus
  27. M. Peyrard, “Nonlinear dynamics and statistical physics of DNA,” Nonlinearity, vol. 17, no. 2, pp. R1–R40, 2004. View at Scopus
  28. R. A. S. Silva, E. D. Filho, and J. R. Ruggiero, “A model coupling vibrational and rotational motion for the DNA molecule,” Journal of Biological Physics, vol. 34, no. 5, pp. 511–519, 2008. View at Publisher · View at Google Scholar