CC BY
22
UDC 539.12
DNA Nonlinear Dynamics With Non-classic Solitary Waves
Ma. de Lourdes Najera*, Jose A. Aguilar^, Maximo A. Aguero^
* Plantel Nezahualcoyotl ^ Facultad de Ciencias Universidad Autonoma del Estado de Mexico Instituto Literario 100 Toluca, Edo. Mex. 50000, Mexico
The nonlinear Schrodinger equation was obtained by using the method of generalized coherent states (GCS) applied to the quasi spin configuration of the DNA lattice model. In this model the DNA macromolecule is subjected to the influence of thermal phonons. By analyzing the nonlinear saturated Schrodinger equation, several non classical soliton-like solutions describing the hydrogen bonds are obtained. Among them there is a couple of compact-anticompact soliton that evolves as a healing mechanism for repairing the eventual open states of the hydrogen bonds.
Key words and phrases: nonclassical solutions, generalized coherent states, DNA, compactons.
1. Introduction
The existence of regimes where nonlinear optical propagation and DNA (or BEC) evolution models obeys the same type of equations are well known [1]. Among them the resulting for intense light pulses or spatial beams, the most fascinating is certainly the formation of solitons, which can also be considered in DNA dynamics.
The theoretical studies of DNA macromolecule dynamics in many aspects show the appearance of nonlinear traveling structures along the chain. That studies have led to many interesting types of solutions like the cnoidal, compact, peak solutions, etc.
Different approaches for variety of models of the DNA molecule have been applied. After the first proposal made by Englander [2], the model that is being deal with rotational degrees of freedom, and the first one was proposed by Yomosa [3,4]. In his pioneering work Yomosa S. presented an important DNA model which include the rotation of pair bases along the spiral model proposed by Watson and Crick. This approach takes into consideration a dynamic of a plane base rotation perpendicular to the helical axis z around the backbone structure. Further improvements were done in the papers [5-8]. Along the other type of model, the principal one in the sense of dynamics of the hydrogen bonds, is due to the Peyrard - Bishop model [9]. These two version of modeling DNA received a great improvements by many important contributions that one can find for example in the papers [5,6,10,11] and citation therein.
It is assumed that open states of the DNA is a principal requirement for arising other important features of this molecule, specifically for transcription and replication processes. Therefore, it is obvious that the opposite type of states the anti - open sates are suggested to appear inside the DNA molecule. In case of existence of these waves, this inverse process dual to the open ones, will repair the disruption or elongation of hydrogen bonds.
On the other order of things, the study of many body problems whose Hamiltonian has been explicitly written in terms of spin or quasi-spin operators, can be reduced by appropriated approximations to a quasi-classical treatment. For doing this, certain reduction procedure from the quasi - spin operator description to a classical one is
Received 27th June, 2012. We acknowledge the efforts of Secretaria de Investigación de la Universidad Autonoma del Estado de Mexico and of Secretaria de Educacion Publica (Mexico) for supporting the fundamental research under the project FE01/2012 involving several research scientific groups in the country.
needed. This procedure consists in choosing trial functions which can be used for averaging the Hamiltonian. Then it is natural to choose for this aim, coherent states since these states are the most classical and minimized the uncertainty relation [12-14].
In this contribution we study a dynamical model related to the emergence of traveling waves in the DNA molecule by using the quasi - classical approximation due to the generalized coherent states approach. This treatment was possible to implement because of the quasi spin character of the Hamiltonian of the DNA that is written in terms of the Spin operators S proposed by Takeno-Homma (TH) in [5,6]. The approach to derive analytical solutions was implemented explicitly for the case of traveling waves. In the next section we briefly expose the main features of the Takeno-Homma model and the generalized coherent state approach applied to the nonlinear lattice equation of DNA. The equation of motion derived from the application of GCS method is reported in section 3. The nonlinear soliton-like patterns are obtained by analyzing the equations in the continuum limit of the nonlinear lattice equation, that is reported in the fourth section. In the 5-th section we analyze the analytical solutions for hydrogen bond displacements. Comments and conclusion are done in the last section.
2. The Coherent States and Quasi-Spin Model of DNA
The model of the DNA can be done by constructing a Hamiltonian for description the nonlinear dynamics of the isotropically homogeneous coupled quasi-spin chain model in the nearest neighboring interactions. Additionally, since the functions of DNA are switched on under the biological temperature, it should be convenient to include thermal surrounding phonons. Thus, the total Hamiltonian that will be analyzed is written as [7,8]
H = ^ [-J(Sn ■ Sn+1 + Sn ■ Sn+1) - Gu - ai(Xn+1 - X„_i))(Sn ■ Sn)+
(1)
+ ^ + ki(Xn - Xn+1)2 2m 1
The quasi-spin operators Sn = (Sx, Sy, Sz) are related to the rotational angles and Vn as
Sn = (Sx,Sy ,SZ) = (sin cos <^n, sin sin r^n, cos $n). (2)
This representation is given for the two strands that compose the DNA. The angles and ipn(ipn) represent the angles of rotation of the bases in the nth base pair in the xz and xy plane respectively. The notation for the second strand will differ to the first one by the apostrophe in the mathematical expressions for its corresponding magnitudes. We will analyze the symmetric case when zn = zn. We denote by Xn as the displacement of the bases along the hydrogen bond at the nth site and by pn = m1X^n their corresponding momentuma. Details on how to construct the hamiltonian model for the DNA with torsion can be found in the works [5-8,15,16]. In some sense here the double strand DNA model is mapped onto a two coupled spin chain model or a spin ladder system with ferromagnetic legs.
On the other hand, there is a powerful method of the generalized coherent state. that allows to obtain a quasi-classical version of the model by averaging the quantum or quasi-quantum Hamiltonian over some states that minimizes uncertainties [17]. Here as mentioned above, we investigate the particle-like localized nonlinear excitation in the quasi quantum Hamiltonian for DNA obtained by Takeno-Homma [5,6] by making use the reduction procedure, based on the generalized coherent states (GCS) on the group SU(2)/U(1). We consider the homogeneous system that means all exchange integrals in the Hamiltonian (1) are constants and do not differ from point point to point. The quasi-spin piece of the Hamiltonian (1) is formally written in terms of spin
operators with spin value j = 1/2. Since, the spin Hamiltonian and the SU(2) coherent states are constructed on the operators of the same group, here is not necessary to carry out the bozonization procedure of the Hamiltonian. The exchange integrals in the x, y, z direction are the same so we have a quasi spin system with zero anisotropy. This type of system according to Makhankov [18] has the SU(2) symmetry. Also we can consider the easy axis of the model the axis that is the direction of helical axis of anisotropy. All these qualities permit us to use the approach of generalized coherent states based on SU(2) group.
It is known that the system of spin coherent states (GCS constructed on the SU(2)/U(1) space; which is the classical Hopf fibration S1 ^ S3 ^ S2) may be written as [12]:
№> =T(<?)№o> = 10> = (1 + H2) J>s+10> (3)
with S± = Sx + iSy, ^ = j^-Tan|a|, a, ^ are complex numbers, |0> = |j, —j> and j
defines the unitary representation of the group SU(2). Formally we have the triplet: {G,n,G/P>, where n is the continuum projector.
After averaging the Hamiltonian (1) with the aid of SU(2) GCS it is obtained obtained
Ü)
2 V (1 + I2)(1 + i^n+i!2)
+ ( CnCn+1 + £n£n+0 + (1 - 1 £n|2)(1 - 1 £n+1|2)
R _ V^ ( { (^n^n+1 +^n^n+i) + (1 - |^n|2)(1 - |^n+l|2 ) +
2\ n x l./> W, +
(1 + | £nf)(1 + |^n+112) J
1( (Y Y ^PnUn + 2jn^n + (1 - |Cn|2)(1 - |^n|2) +
- - «1(Xn+1 - Xn-1)) +
+ + k1(xn -Xn+1 )2} . (4)
2m
3. Equations of Motion
let us study the continuum version of the Hamiltonian model (4). Therefore, we consider the inter-site distance between nucleotides are quite small in comparison with the DNA length. Introducing the new variable fields Xn ^ X(z, t) with z = na and as usual making the standard procedures we have
a2
^n±i = ^(z, t) ± aipz + + •",
also
and
|^n±1|2 _ M 2 ± a(\^\)z + +...,
2!
£ - >1 dr-
The similar expansion can be straightly done for £n±1 and |£n±i|2- These last variables correspond to the complementary strand (apostrophe variables) of the double ladder of DNA. Also similar approximation is done for Xn±1. After some algebra we obtain the quasi-classical Hamiltonian
2
{t(
h = M- +
1 + W 1 + \i\2)_ _
(i - 2 aiaXz) (2(^i + ^) + (1 - \£\2)(1 - \^\2)
4 a V (1 + \£\2 )(1 + W )
+
+ —^---+ k1a(Xz)H dz + const. (5)
2 am1
The equation of motion for the unknown variables can be obtained using the standard Hamiltonian method. Thus, for finding the unknown displacements X(z, t) from
OH
the Eq.(5) we can use the standard Hamiltonian equations of motion X = —— and
op
its canonical conjugate counterpart. The equation of motion for the unknown field functions z, t) and £(z, t) can be written using the variational derivative as [19]
ifr = -(1 + M2)2
by
The similar equation of motion is directly built also for the field variable . Considering these two nonlinear equations after some algebra we can find the next nonlinear Nonlinear Schrodinger equation with saturation equation
•/ Ja /a. T ^ , (2» y \ 1 -|y|2 , ,,,
» y=+aj 1+w + v v- 4a1Xz) iw(6)
For description the hydrogen bond displacements X(z, t) one can obtain the next also nolinear equation that depend on the torsions (^ ) of the previous nonlinear equation
„1*. = 2ak№, + £ () - (7)
2 \ (1 + M2)2
4. Nonclassical Traveling Waves
We will look forward for traveling wave solutions. For doing this we change the variables a = z - vt. Applying the procedure due to Vasumathi and Daniel [7,8] i.e. making Xz = W (z, t), after some algebra it yields
X = ai 1 + M4 - 6\^\2 (8)
z 2(miv2 - 2aki) (1 + \^\2)2 ' ()
Replacing the equation (8) in to the nonlinear Schrodinger equation (6) and by consid-
[Ja
ering restriction on the exchange integrals and the change of variables z = sy — the
equation (6) can be reduced to the well known Cubic-Quintic Nonlinear Schrodinger Equation (CQNSE):
iipt + ^ss + Ki^ + «sMV - K5\*p\4^ = 0 (9)
with the following parameter relations
Ki = -2i + 3, K3 = 4i - 10/3, K5 = 4i - 34/3 (10)
with
*=i11)
m\v2 — ¿ak\
Let us analyze some particular interesting cases of the obtained equation (9)
I. Case: k5 = 0. In this particular case, the equation (9) reduces to the well known nonlinear Schrodinger equation
i pt + Pss + 7|p|2p = 0 (12)
9
with the following redefined functions and parameters ^ = p etKlt and k = —7,
8
k3 = 7 = 24*. The previous nonlinear equation is the famous nonlinear Schrodinger equation that is completely integrable. Thus there does exist a method for instance an inverse scattering or Hirota methods for finding its N-soliton solutions [20].
For obtaining some information on for the hydrogen bond displacement of the DNA macromolecule let us take some solution of the Schrodinger nonlinear equation (12). First we consider the sign of the parameter 7 to be positive. For this case as is well known [21] when 7 > 1 the solution of the NSE (12) with the trivial boundary condition that means when the field tends to zero as |íc| ^ is
/ 2V \l/2
p(s, t) = (--J Sech [V—Ds\ exp{— ivt} (13)
being an arbitrary constant. Since we are interested in obtaining traveling wave solution, we will use the well established fact that NSE is invariant under the Galilean transformations
C = s-Vt, T = t, p =p expj - ^ S- \vt\\, (14)
where Q = a/y ^, V = v/y ^ that finally gives us the traveling soliton V(s, t) = (-1/2 Sech [V-7 C] exp f - | | (V^) ' + 2^
This solution applied to DNA molecular dynamics should imply that a collective wave of bell type forming a certain open state for the angle deviation defined by |-0|2 = Ipl2 = (tan($/2)2 should be propagated along the chain due to the interaction of the DNA molecule with thermal phonons surrounding molecule. This because the collective wave is constructed by the two types of bell solutions for ^ and for £ field variables correspondingly.
II. Case. a3 = 0 In this case we have the Quintic Schrodinger equation.
i pt + pss + 72M 4P = 0. (15)
This nonlinear equation supports soliton solution in the attractive case when the coefficient 72 is positive. The solution that divides the scenario of the evolution equation in two sectors is possible to exist [22]
31/4(—w)1/4 ■ t p(s, t) = , ( w)^= e-wt,
VCosh(2V-w s)
The first sector as is well known is the sector of collapse when the amplitude of the solution infinitely grows in a finite time. The second scenario correspond to the dispersion solution when the amplitude decays and the area of its localization grows. The collapse is not physical well solution. Then, in this case it is convenient to introduce other effects perhaps dissipative ones for compensate the blow up of solutions. These can be numerous nonlinear dissipation mechanisms such as inelastic collisions of the pair of base of DNA with surrounding proteins molecules which results in loss of energy by emitting heat. This interesting issue could be studied further.
5. Hydrogen bond displacements
The mixture of both solutions ^ and — ^ will maintain a certain configuration state along the chain for angle rotations § and p. For instance, we can see a configuration made by Sech2(£) and by —Sech2( Q in such a manner that along the DNA chain should be created a wooden spinning top state, for the angle § deviation. In this case, one can observe the envelope of the edges of quasi-spin arrows. The angle that forms the arrows of quasi spins with respect to the axes Q represents the angle deviation from equilibrium position. Far from the central part the excitation zone, the angle deviations § in the middle site acquires the maximum value while outside of the center the angles vanish.
For calculating the traveling waves for displacements X along the hydrogen bonds we can use the analytical solution (13). We replace this solution into the equation (8) and integrate once with respect to the variable . After integration we have obtained a traveling solution for the displacements X( ) with an arbitrary constant of integration X0. This constant of integration can be considered by some initial condition imposed to the nonlinear equation (8). Further, the integration has been done numerically and we found the pair compacton structure represented in the Fig. (1), by virtue of satisfying the trivial boundary condition at infinity and initial conditions. Resuming, we found a solutions for the hydrogen bond displacements that can be represented by
X(0 = X(S — Vt) = djF(|p|2)dC, for — Co Co X (C) = 0, elsewhere
(16)
with F(|pi2)=1H(ilp+4|-2';f and i = f/f (17)
This solution represents a couple of compact and anti-compact fused structures. Specifically, the sector for negative values of independent variable — C0 ^ C ^ 0 corresponds to the compact piece of the solution and in the region of 0 < C ^ C0 lives the anticompacton branch of the solution. Thus, this solution is made of two different contributions of the waves. As it is well common for similar compact structures [23] we observe the discontinuity on its first derivative at points ±£0. Thus, the compacton sector ensure the formation of open states due to positive values of displacement of hydrogen bonds, while the anti-compacton sector represents the reconstruction of hydrogen bonds, see Fig.(1). The compacton and the anti-compacton parts collaborate each other in such a manner that the DNA molecule preserves its identity. We observe here a natural mechanism of healing the open states that profusely can be created along the different sector of DNA molecule.
For the second case or the quintic Schrodinger equation, the hydrogen bond displacements can be obtained analytically and we have the particular solution for example when w = —1
X(*») = 1S ( 2^ + 4^6ArcTanh
(—1 + V3) Tanh(C)
72
12Sinh(2 0 \ 73+Cosh(2 0 J
(18)
Figure 1. A couple compacton anti-compacton solution for hydrogen bond displacements along the DNA. The positive values of displacements elongate the separation of the two strands while the negative ones stretch this separation. The values of ¡3 =1 was used for numerical integration of the
corresponding equations
for -Co < c < Co-
The profile of this solution has the similar fashion of the previous solution found for the case of cubic Schrodinger equation with the Fig. (1). Again this solution represents a typical couple of compacton and anti-compacton pair solution. Although the angle deviations undergo collapse effects, the corresponding hydrogen displacement shows a regular behavior.
The energy of the compacton pair can be easily calculated by virtue of equation (5). The integration with respect to the variable C in the interval [-Co, Co] of first pair compacton anticompacton can be evaluated directly using the energy density
E = g(X2) + hXc F (M2) (19)
with g = mV , h = ai and F(j'-pj2) is defined by the relation (17). The
Ja y2Ja
compacton solution can emerge only when the velocities do not satisfy the relation mv2 — 2ak1 = 0. Otherwise, we could encounter nonphysical solutions with infinity values of energy.
6. Conclusions
We have discussed in this contribution the appearance of non classical soliton-like excitations along the DNA molecular model based on the pioneering work of [5,6]. These solutions could eventually be responsible for various fundamental processes inside the DNA macromolecule. As is well known the relation of the quasi-classical field variable of the GCS approach ^ and the angle $ being the deviation of the angle of the classical spin S see Eq. (2) from the OZ axis is the following |^|2 = tan2($/2) [12]. Thus the using of generalized coherent state approach provides us directly with the possibility to calculate the angle deviation of the "classical" spin with respect the chosen direction z. This approach allows us to derive a familiar NSE, that for very specific parameter values represents the well known completely integrable model. The angle deviation from the axis z will completely determined by the solutions of integrable NSE. These exact solutions that could serve as a good first approximation in a complex and deep theoretical and experimental analysis of the DNA molecule.
The hydrogen bond displacements, shows unexpected behavior represented by the couple of compacton anti-compacton solution. The velocities of these solutions are determined from the parameters k (10) that are restricted by the relation ft > 0 from equation (11).
The evaluation of displacements of the hydrogen bonds along the DNA can be carried out straightaway by integrating their corresponding nonlinear equations (8). We considered that for displacements that occur sufficiently far from the excitation zone, vanish. Therefore, in both direction of the axis Q we will observe the physical interesting case of vanished values of hydrogen bond displacements at long distances from the active zone. By taking in mind this physical reason and integration of the equation (8) we obtain for specific parameter values the compacton anti-compacton couple while the angle deviation evolves like a bell soliton along the DNA macromolecule.
As is well known the open states of DNA chain is considered the precursors of denat-uration of this macromolecule. As soon as the nonlinear excitations in whichever segment of the DNA molecule is being activated by external factors for example phonons in our case, nonlinear waves of coupled compacton and anti-compacton pieces could appear, describing the evolution of hydrogen bonds. The anti-compacton piece one could be interpreted as a healing wave part e that cures the bubble or open states on the DNA segments.
References
1. Serkin V. N., Hasegawa A., Belyaeva T. L. Nonautonomous Solitons in External Potentials // Phys. Rev. Lett.--Vol. 98, issue 7. — 074102.
2. Nature of the Open State in Long Polynucleotide Double Helices: Possibility of Soliton Excitations / S. W. Englander, N. R. Kallenbach, A. J. Heeger et al. // Proc. Natl. Acad. Sci (U.S.A.). — 1980. — Vol. 77 (12). — Pp. 7222-7226.
3. Yomosa S. Soliton Excitations in Deoxyribonucleic Acid (DNA) Double Helices // Phys Rev. A. — 1983. — Vol. 27. — Pp. 2120-2125.
4. Yomosa S. Solitary Excitations in Deoxyribonucleic Acid (DNA) Double Helices // Phys. Rev. A. — 1984. — Vol. 30. — Pp. 474-480.
5. Takeno S., Homma B. Topological Solitons and Modulated Structure of Bases in DNA Double Helice // Prog. Theor. Phys. — 1983. — Vol. 70, No 1. — Pp. 308311.
6. Takeno S., Homma B. A Coupled Base-Rotator Model for Structure and Dynamics of DNA // Prog. Theor. Phys. — 1984. — Vol. 72. — P. 679.
7. Daniel M., Vasumathi V. Perturbed Soliton Excitations in the DNA Double Helix // Physica D. — 2007. — Vol. 231. — Pp. 10-29.
8. Daniel M., Vasumathi V. Base-Pair Opening and Bubble Transport in a DNA Double Helix Induced by a Protein Molecule in a Viscous Medium // Physical Review E. — 2009. — Vol. 80. — P. 061904.
9. Peyrard M., Bishop A. R. Statistical Mechanics of a Nonlinear Model for DNA Denaturation // Phys. Rev. Lett. — 1989. — Vol. 62, No 23. — Pp. 2755-2758.
10. Yakushevich L. V. Nonlinear Physics of DNA. — Wiley (Chichester), 1998. — Second edit. 2004.
11. Simple Models of Non-Linear DNA Dynamics / G. Gaeta, C. Reiss, M. Peyrard, T. Dauxois // Rivista del Nuovo Cimento. — 1994. — Vol. 17, No 4. — P. 148.
12. Makhankov V. G. Soliton Phenomenology. — Kluwer Academic Publishers, 1990.
13. Perelomov A. Generalized Coherent States and Their Applications. — Berlin: Springer, 1986.
14. Klauder J. R., Skagerstam B. S. Coherent States. Applications in Physics and Mathematical Physics. — World Scient, 2010.
15. Takeno S., Homma S. Kinks and Breathers Associated with Collective Sugar Puckering in DNA // Prog. Theor. Phys. — 1987. — Vol. 77, No 3.
16. Takeno S. Compacton-Like Modes in Model DNA Systems and Their Bearing on Biological Functioning // Phys. Lett. — 2005. — Vol. A 339, No 3-5. — Pp. 352360.
17. Kosevich A. M., Ivanov B. A., Kovalyov A. S. Magnetic Solitons // Phys. Rep. — 1990. — Vol. 194, No 3-4. — Pp. 117-238.
18. Makhankov V. G., Makhankov A. V. Spin Coherent States // Phys. Stat. Sol. B. — 1988. — Vol. 145. — Pp. 669-678.
19. Fadeev L., Takhtajan L. Hamiltonian Methods in the Theory of Solitons. — Springer, 1987.
20. Sulem C., Sulem P. L. The Nonlinear Schrodinger Equation: Self-focusing and Wave Collapse. — New York: Springer-Verlag, 1999.
21. Debnath L. Nonlinear Partial Differential Equations. — Birkhauser, 1997.
22. Alfimov G. L., Konotop V. V., Pacciani P. Stationary Localized Modes in the Quintic Nonlinear Schrodinger Equation with a Periodic Potential // Arxiv:nlin/0605035v2. — 2006.
23. Rossenau P., Hyman J. M. Compactons: Solitons with Finite Wavelength // Phys. Rev. Lett. — 1993. — Vol. 70. — Pp. 564-567.
УДК 539.12
Нелинейная динамика неклассических уединённых волн
Мария Лурдес Нахера*, Хосе А. Агиллар^, Максимо А. Агуэро^
* Институт 'растениеводства ^ Факультет наук Государственный автономный университет г. Мехико, Мексика
Методом обобщённых когерентных состояний (ОКС) в применении к решеточной модели ДНК получено нелинейное уравнение Шрёдингера. В этой модели макромолекула ДНК подвержена влиянию тепловых фононов. Путем анализа уравнения Шрёдингера с насыщением получено несколько солитоно-подобных решений, описывающих водородные связи. Среди них имеется пара «компактный и антикомпактный солитоны», которая эволюционирует в соответствии с механизмом восстановления открытых состояний с водородными связями.
Ключевые слова: неклассические решения, обобщённые когерентные состояния, ДНК, компактоны.
