Research Article Volume 7 Issue 3
1Institute of Cell Biophysics RAS, Russia
1Bashkir State University, Russia
Correspondence: Farit K Zakiryanov, Institute of Physics and Technology, Bashkir State University, Zaki Validi 32, Ufa, Russia, Tel 791-7406-5235
Received: May 04, 2018 | Published: June 6, 2018
Citation: Yakushevich LV, Balashova VN, Zakiryanov FK. Is it possible to govern the movement of the transcription bubbles of DNA by constant and periodic external fields?. Biom Biostat Int J. 2018;7(3):218–224. DOI: 10.15406/bbij.2018.07.00212
The influence of constant and periodic external fields on the living systems can lead to various changes in their functional properties. At the molecular level, these changes can be modeled as the changes in the dynamic properties of biomolecules. In this paper we investigate the effect of constant and periodic fields on the movement of transcription bubbles being small unwound regions (~10bp) in the DNA double helix, which result from the interaction of the RNA polymerase with the DNA promoter region at the initiation stage of the process of transcription. The bubbles are modeled mathematically by kinks being one–soliton solutions of the modified sine–Gordon equation. To obtain these solutions and calculate the dynamic characteristics of the bubbles including their velocity, coordinate, phase portrait, total energy and size, we apply the method of McLaughlin and Scott and quasi–homogeneous approximation. For definiteness all of the calculations were made for the pTTQ18 sequence. The results were used to analyze the bubbles dynamic behavior and to answer the question: is it possible to govern the movement of the transcription bubbles of DNA by constant and periodic external fields? Possible consequences on the gene expression and biological activity of cells are discussed.
Keywords: transcription bubble, DNA torque, periodic fields, sine–Gordon equation, method of McLaughlin and Scott, quasi–homogeneous approximation, nonlinear DNA dynamics
DNA, deoxyribo nucleic acid; AOT, angular optical trap; MTT, magnetic tweezers technique
The action of constant and periodic fields on the living systems is one of the hotly discussed themes of modern biophysics because of the growing amount and variety of electronic devices and their influence on the basic, fundamental life processes and, consequently, on human health. It is assumed, in particular, that external periodic fields with the frequency of terahertz diapason can cause the changes in the physicochemical properties of DNA, which in turn can lead to the changes in the transcription and thus to the changes in the gene expression and cell differentiation.1,2 However, until now the mechanisms of the action of these fields remain unclear. Special attention is also paid to the influence of constant fields because of significant progress in the experimental studies of the dynamics of single molecules of DNA and by the development of new techniques that allows to measure directly the DNA torque.3 The method of the angular optical trap (AOT)4,5 and the method of magnetic tweezers technique (MTT)6–8 are among them. However, these and other currently existing experimental methods do not allow measure the DNA torque directly in the biological processes in which the DNA molecule is involved.
In this paper we apply the methods of mathematical modeling to investigate the action of constant and periodic fields on the movement of transcription bubbles of DNA which are small unwound regions (~10bp) in the DNA double helix resulting from the interaction of the RNA polymerase with the DNA promoter region at the initiation stage of the process of transcription (Figure 1).9 Our aim is to clarify the mechanisms of the influence of the fields on the bubbles dynamics, to find the relationship between the DNA torque and the bubble velocity, to estimate the value of the DNA torque necessary for the moving of the bubbles with the velocity of the process of transcription, to find the relation between the dynamic behavior of the bubbles and their initial velocities, and to give an answer to the question is it possible to govern the movement of the transcription bubbles of DNA by constant and periodic external fields (Figure 1).
In physics, the DNA molecule is considered as a complex dynamic system consisting of a large number of coupled atoms and atomic groups which are arranged in a certain way in space. Such a system is not static, but movable. It has a large amount of internal motions caused by the effect of temperature, collisions with the molecules of the solution, and interactions with proteins. Among the variety of the internal motions of DNA there are the translational movements of transcription bubbles that can be considered as translational movements of quasi–particles in the potential field of the DNA.
Mathematically the transcription bubbles can be described as the soliton–like solutions of the nonlinear differential equations imitating the internal DNA mobility. In the case of homogeneous synthetic DNA having the sequence of identical bases, Englander and co–authors10 showed that the bubble movement can be modeled by the sine–Gordon equation with constant coefficients, this equation having the exact one–soliton solutions in the form of kinks. In the work of Englander and co–authors10 as well as in the works of other authors developing the Englander’s idea,11–20 just these solutions were used to simulate the DNA open states or bubbles. To take into account effects of dissipation and the action of some external field, they modified the sine–Gordon equation by adding the following two terms: one to model effect of dissipation and the other to model influence of an external field. With the help of the method of McLaughlin and Scott21 the kink–like solutions of the modified sine–Gordon equation were found in the two particular cases: when the external field is constant and when the external field is periodic.22,23
In the case of inhomogeneous DNA, the coefficients of the modified sine–Gordon equation are no longer constants, but depend on the sequence of bases (Table 1). However, in the quasi–homogeneous approximation, this equation can be reduced to the homogeneous one but with the changed coefficients the values of which depend on the concentrations of different types of bases. In this paper we apply the quasi–homogeneous approximation to solve the equation of that type and to calculate the main dynamic characteristics of the bubbles necessary to analyze the influence of constant and periodical fields on the bubble behavior.
For definiteness we take a sequence of small circular DNA molecule–plasmid pTTQ18 (Figure 2) which is widely used in genetic engineering to transfer genetic information and for genetic manipulations.24 The plasmid sequence contains four functionally important regions: promoter (54–333), terminator (611–8011), and two coding regions CDS–1 (1585–2544) and CDS–2 (2762–3622). The sequence length is 4563 bases. Of these, the sequence has 1105 adenines, 1090 thymines, 1193 guanines and 1175 cytosines. In the second section, we describe the model and methods used. In the further two sections the results on the velocity, coordinate, phase portrait, total energy and size of the bubbles moving under the action of the constant and periodic fields are presented. In the final section the obtained results are discussed and the main conclusions are made (Figure 2).
Let us begin with the discrete version of the modified sine–Gordon equation:25
Ind2φn(t)dt2+αndφn(t)dt−KRn(Rn+1φn+1(t)−2Rnφn(t)+Rn−1φn−1(t))+Vnsinφn(t)==M0+M1cos(Ωt). (1)
Here φn(t) is the angular displacement of the n–th base, In is the moment of inertia of the n–th base, K is the stiffness of the sugar–phosphate chain,ais the distance between adjacent base pairs, Vn is the factor characterizing the interaction between the complementary bases inside the n–th pair, n = 1, 2, ... N, N is the number of bases in the sequence, αn=R2nλ, λ is the dissipation factor, M0 is the constant torque, M1cos(Ωt)is a periodic field.
The values of the coefficients in the left side of the system of equations (1) are shown in Table 1. The values of the parameters of the external fields (M0,M1 and Ω ) are yet arbitrary.
Type of the |
In |
K |
Rn |
Vn |
a |
αn |
adenine |
7.61 |
6.75 |
5.8 |
2.09 |
3.4 |
4.25 |
thymine |
4.86 |
6.75 |
4.8 |
1.43 |
3.4 |
2.91 |
guanine |
8.22 |
6.75 |
5.7 |
3.12 |
3.4 |
4.10 |
cytosine |
4.11 |
6.75 |
4.7 |
2.12 |
3.4 |
2.79 |
Let us assume that the desired solutions of Eqs. (1) are sufficiently smooth functions. Then, we can apply the continuum approximation:
a→0, zn=na→z,
ϕn(t)≡ϕ(zn,t)→ϕ(z,t),
In≡I(zn)→I(z), Vn≡V(zn)→V(z),
Rn≡R(zn)→R(z), αn≡α(zn)→α(z).
As a result Eqs. (1) are transformed to:
I(z)d2φ(z,t)dt2+α(z)dφ(z,t)dt−KR(z)a2d2[R(z)φ(z,t)]dz2+V(z)sinφ(z,t)==M0+M1cos(Ωt). (2)
To simplify further calculations, we use the quasi–homogeneous approximation according to which we can replace the coefficients of the left–hand side of Eq. (2) by the averaged values:27
I(z)→ˉI=IANAN+ITNTN+IGNGN+ICNCN,
R(z)→ˉR=RANAN+RTNTN+RGNGN+RCNCN, (3)
V(z)→ˉV=VANAN+VTNTN+VGNGN+VCNCN,
α(z)→ˉα=αANAN+αTNTN+αGNGN+αCNCN,
where NAis the number of adenines, NT is the number of thymines, NGis the number of guanines, NC is the number of cytosines, and N=(NA+NT+NG+NC) is the total number of bases in the sequence.
After the averaging procedure, Eq. (2) takes the form similar to the sine–Gordon equation, but with the coefficients recalculated by formulas (3):
ˉI∂2φ(z,t)∂t2+ˉα∂φ(z,t)∂t−ˉK'a2∂2φ(z,t)∂z2+ˉVsinφ(z,t)=M0+M1cos(Ωt), (4)
where ˉK'=KˉR2. Numerical values of the coefficients that are averaged over the entire sequence of plasmid pTTQ18 are given in Table 2.
Type of the sequence |
ˉI |
ˉK' |
ˉV |
ˉα |
pTTQ18 |
6.21 |
1.88 |
2.21 |
3.51 |
Table 2 Coefficients of Eq. (4)
For convenience, let us introduce new (dimensionless) variables:
τ=σt, ζ=μz, where σ=(ˉV/ˉI)1/2, μ=a−1(ˉV/ˉK')1/2. In these variables, Eq. (4) takes the form:
φττ+βφτ−φζζ+sinφ=f0+f1cos(ωτ), (5)
where β=ˉα/(ˉIˉV)1/2, f0=M0/→V, f1=M1/ˉV,ω=Ω(ˉI/ˉV)1/2.
If the dimensionless coefficient of dissipation is small (β<<1), Eq. (5) has approximate one–soliton solution in the form of kink:
φ=4arctan[exp(γ(ζ−υ(τ)τ−ζ0] (6)
The velocity of which is defined by the equation of McLaughlin and Scott:21
dυdτ=−βυ(1−υ2)+(1−υ2)3/2+π4[f0+f1cos(ωτ)] (7)
Here γ=1/[1−(υ(τ))2]1/2, ζ0 is an arbitrary constant.
Let us determine the coordinates of the kink x by formula:
dξdτ=υ, (8)
and the total energy and the size of the kink by the formulas:
e(τ)=8√1−υ2(τ), D(τ)=1√1−υ2(τ). (9)
Differentiating (9) with respect to τ, we get two additional differential equations:
dedτ=8υ(1−υ2)3/2dυ(τ)dτ, (10)
dDdτ=υ(1−υ2)3/2dυ(τ)dτ. (11)
Eqs. (7), (8) and (10), (11) are a complete set of equations which are necessary to calculate the time dependencies of the velocity, coordinate, size and total energy of the transcription bubbles.
Here we present the results of numerical calculations of the velocity, coordinate, size and total energy of the transcription bubbles. To better understand the dynamic behavior of the bubbles, we made for three different values of the external torsion moment and for three different values of the initial velocities of the bubbles.
Bubbles dynamics under the action of constant external field
This case corresponds to f0≠0,f1=0. Dynamic characteristics of the bubbles, obtained for three different values of the torsion moment f0, are shown in Figure 3A & 3B. The bubbles trajectories in the phase plane (υ,ξ) are shown in Figure 3C. The bubble total energy and size are shown in Figures 3D & 3E. The time required for calculation of the curves in Figures 3A &3B and of the curve 1 in Figure 3(C), is equal to 2.5×103. The time required for calculation of the curves 2 and 3 in Figure 3C is equal to 5×103 and 5×104, respectively (Figure 3).
The values f01and f03 are chosen arbitrarily. The value f02is assumed to be equal to fcrit0 which in turn is determined from Eq. (7) at dυdτ=0:
fcrit0=4βυ0π√1−υ20 (12)
From Figure 3A it is seen that at f02=fcrit0(curve 2) the bubble velocity is constant and equal to the initial velocityυ0. At f03<fcrit0 (curve 3), the bubble kink velocity initially decreases and after a time period (T), equal approximately 500 (0.84ns), reaches the stationary value:
υst03=[1+(4βπf03)2]−1/2=0.010 (13)
At f01>fcrit0(curve 1), the kink velocity initially increases and after approximately the same time period T, reaches the stationary value:
υst01=[1+(4βπf01)2]−1/2=0,207 (14)
From Figure 3B it is seen that at f02=fcrit0 (line2) the bubble coordinate is a completely straight line. At f03<fcrit0 (curve 3), the coordinate initially forms a small bend, and then after a certain period of time it transforms into a straight line. Similarly, at f01>fcrit0 (curve 1), the coordinate initially forms a small bend (in opposing direction), and then it also transforms into a straight line.
From Figures 3D & 3E it is seen that at f02=fcrit0(curve 2) the bubble energy and size are constants. At f03<fcrit0 (curves with number 3), the bubble energy before reaching the stationary value is reduced from the value e01=8.04, and the bubble size is reduced from D01=0.995 till Dst,1=0.99. At f01>fcrit0 (curves with number 1), the bubble energy before reaching the stationary value is increased from the value e01=8.04 till the value est,1=8.16, and the bubble size is increased from D01=0.995 till Dst,1=1.01. Hence, by setting different values of constant external torsion field we can force the bubble to move with a certain fixed velocity. The opposite is true. If we know that bubble is moving at a given velocity, we can calculate the value of the torsion moment, which will ensure the movement of bubble at that velocity.
The above studies make it possible to estimate the value of the constant torsion moment f03, necessary for the movement of the transcription bubble at the velocity equal to the velocity of the process of transcription. To do this, we replace υst03 in formula (13) with the transcription velocity υtrand replace the constant field f03 with the value ftr0. As a result, we obtain the ratio between the transcription velocity and the sought value of the constant external field:
υtr=[1+(4βπftr0)2]−1/2 (15)
We rewrite the formula (15) in dimensional variables:
Utr=С0√1+(4ˉαπMtr0)2ˉVˉI (16)
where Utris the dimensional transcription velocity, Mtr0is the unknown constant external field, С0=√ˉK'a2/ˉI is the sound velocity in the plasmid pTTQ18. From (16) we find the sought torsion moment 0.49.10–31:
Mtr0=4ˉαπ√ˉVˉIUtrС0√1−(UtrС0)2 (17)
If we assume that the velocity of transcription is equal to 100 base pairs per second Utr=0.34.10−7m/s, then from formula (17) we find the desired estimate: Mtr0= 0.49.10–31J. Figure 4 shows the results of numerical calculations of the dynamic characteristics of the bubbles, made for three different values of the initial velocity υ0 (Figure 4).
Figure 4A shows that for any value of the initial bubble velocity υ0 the values of stationary velocity are the same. Figure 4B shows that the slopes of the bubbles trajectories reach the same value. Figure 4C & 4D shows that for different values of initial velocities the bubbles energy and size also reach the same values. Thus, we can conclude that the stationary dynamic characteristics of the bubbles do not depend on the initial velocities. On the contrary, they depend only on the value of the torsion moment.
Bubbles dynamics under the action of periodic external field
This case corresponds tof0=0,f1≠0,ω≠0. It could be realized, for example, in an experiment with a single molecule at the work of atomic force microscope in an oscillating manner.28 Figure 5A–5C show the dynamic characteristics of the bubbles, obtained for three different values of the amplitude of the external periodic field f1(Figure 5).
Analyzing the graphs of coordinate and velocity presented in Figure 5A & 5B one can notice that at the beginning of the bubbles movement there exists a short period T~500 where oscillations are establishing. After the end of the period T the bubbles continues to oscillate with the constant frequency ω. Figure 5D & 5E show the time dependence of the total energy and the size of the bubbles. It can be seen the doubling of the frequency of the stationary oscillations. Figure 6 shows the dynamic characteristics of the bubbles, obtained for three different values of the initial velocity: υ01,υ02,υ03(Figure 6).
From Figure 6 it is seen that after some short period T the velocity and coordinate of the bubbles having different initial velocities at the "start–up, begin to oscillate with constant amplitude and frequency. It can be seen also the doubling of the frequency of the oscillations of the bubbles energy and size.
Transcription bubbles have been modeled as quasi–particles – kinks, moving in the potential field of DNA. The influence of constant and periodic fields on the movement of the kinks was studied by the method McLaughlin and Scott. The time dependences of the kinks velocity, coordinate, total energy and size were obtained.
When calculating we used the following dimensionless values of the initial bubbles velocity: υ01=0.1, υ02=0.5, υ03=0.8. In the dimensional units these values correspond to: v01=187m/s, v02=935 m/s, v03=1494m/s. To model the constant torsion moment, we used the dimensionless values: f01=2.429⋅10−3, f02=1.215⋅10−3, f03=1.215⋅10−4, corresponding to the dimensional values: M01=5.37⋅10−23 J, M02=2.68⋅10−23 J, M03=2.68⋅10−24J. To simulate the amplitude and frequency of the periodic external field, we used dimensionless values: f11=2.429⋅10−3, f12=1.215⋅10,, f13=1.215⋅10−4, ω=0.02, corresponding to the dimensional values: M11=5.37⋅10−23J, M12=0.268⋅10−23J, M13=2.68⋅10−24J, Ω0=0.119⋅1011 s–1.
It was shown that in the case of constant torque, there is a small time period T~500 (0.84 ns), after which the bubble velocity reaches the stationary value. This stationary value does not depend on the initial bubble velocity and is completely determined by the value of constant torque. By setting this or that value of the constant torsion moment we can force the bubble to move at a predetermined velocity. We used this result to estimate theoretically the value of the DNA torque which can not be measured by currently existing experimental methods. We obtained that this value is equal to 0.49∙10–31 J if the velocity of transcription process is equal to 100 base pairs per second. Moreover we derived formula relating DNA torque and the velocity of transcription.
It was also shown that in the case of periodic field there is almost the same time period after which the bubble begins to oscillate with constant amplitude and frequency the value of which does not depend on the initial bubble velocity. This can lead to the destruction of the transcription bubble.
It is necessary to note, however, that all these results has been obtained in the frameworks of a rather simple model, which simulates angular oscillations of nitrous bases in one of the two polynucleotide chains, the other chain being modeled as an averaged external field. But we expect that further improvement of the model due to taking into account the transverse and longitudinal displacements of nitrous bases, the mobility of bases in the other polynucleotide chain, the helicity of the DNA structure will not change the main conclusions about the role of the constant and the periodic field in the bubble dynamics
The authors thank Dr. Dana Flavin for the valuable discussions which helped to broaden and deepen our understanding of the importance of the considered problems and its relationship to the fundamental problems of vital activity and human health.
Author declares that there is no conflict of interest.
©2018 Yakushevich, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.
2 7