ISRN BiophysicsVolume 2012 (2012), Article ID 725350, 6 pagesdoi:10.5402/2012/725350

Research Article

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

R. A. S. Silva,¹ J. D. T. Arruda-Neto,1,2 and L. Nieto³

¹Physics Institute, Rua do Matão, Travessa R, No. 187, University of São Paulo, 05508-090 São Paulo, SP, Brazil²FESP, São Paulo Engineering School, São Paulo, SP, Brazil³Laboratory 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 [1–3]. 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 [9–11]. 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 [14–16]. 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 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 is introduced to consider interactions between charged particle approximation and DNA. In this case, , where is related to radiation dose, LET is the average linear energy transfer [24], and Å³ 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 where represents the kinetic energy of the system, and 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 These integrals for all momenta are merely Gaussian integrals, providing 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 where is given by

Integrals in can be calculated using eigenfunctions and eigenvalues of a transfer integral operator [26, 27]:

In the thermodynamic limit, , the partition function is reduced to where is the ground-state eigenvalue for a Schrödinger-like equation given by with

The eigenvalue and the normalized eigenfunction for the ground state are, respectively, with and the normalization constant. In the limit , the average stretching will be dependent only on the ground-state eigenfunction as

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  nm, the number of links involving the neighboring nucleotides is given by  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

3. Results and Discussion

The best values for the parameters obtained from experimental results for DNA presented in the literature as  Å and  eV Å⁻² [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 . 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.

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 ( Å), the values of tend to be larger than the long distances ( Å), 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.

Figure 2: Variation of the average stretching as function of radiation dose, for and  Å.

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.

(a)

(b)

(c)

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.

Table 1: Fit of data in Figure 3(a), using two power laws. Parameter values are, , , , , = 0.4 and = 0.1.

1000
0,17782

950
0,18279

900
0,18829

850
0,19442

800
0,20131

750
0,20913

700
0,21813

650
0,22865

600
0,2412

550
0,25654

500
0,27595

450

0,55703

400

0,64064

350

0,73663

300

0,84905

250

0,98427

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 and , 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

S. H. Strogatz, “Quoted by S. Nadis in “News feature: all together now”,” Nature, vol. 421, pp. 780–782, 2003.

C. Tsallis, “Possible generalization of Boltzmann-Gibbs statistics,” Journal of Statistical Physics, vol. 52, no. 1-2, pp. 479–487, 1988.

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.

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.

U. Amaldi and G. Kraft, “Radiotherapy with beams of carbon ions,” Reports on Progress in Physics, vol. 68, no. 8, pp. 1861–1882, 2005.

E. L. Alpen, Radiation Biophysics, Academic Press, San Diego, Calif, USA, 1998.

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.

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.

M. Krämer and G. Kraft, “Track structure and DNA damage,” Advances in Space Research, vol. 14, no. 10, pp. 151–159, 1994.

D. T. Goodhead, “The initial physical damage produced by ionizing radiations,” International Journal of Radiation Biology, vol. 56, no. 5, pp. 623–634, 1989.

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.

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.

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.

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.

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.

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.

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.

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.

T. C. Lim, “The relationship between Lennard-Jones (12-6) and Morse potencial function,” Journal of Mathematical Chemistry, vol. 33, p. 1, 2003.

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.

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.

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.

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.

E. L. Alpen, Radiation Biophysics, Prentice Hall, NJ, USA, 1990.

J. E. Brady and G. E. Humiston, Química Geral. 2, LTC, Rio de Janeiro, Brazil, 1986.

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.

M. Peyrard, “Nonlinear dynamics and statistical physics of DNA,” Nonlinearity, vol. 17, no. 2, pp. R1–R40, 2004.

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.