We consider the nonlinear dynamics in a double-chain model of DNA which consists of two long elastic homogeneous strands connected with each other by an elastic membrane. By using the method of dynamical systems, the bounded traveling wave solutions such as bell-shaped solitary waves and periodic waves for the coupled nonlinear dynamical equations of DNA model are obtained and simulated numerically. For the same wave speed, bell-shaped solitary waves of different heights are found to coexist.
1. Introduction
In 1953, Waston and Crick discovered the structure of deoxyribonucleic acid (DNA) double helix [1], which opened the area of molecular biology. The significance of the double helix model not only means the proven structure of DNA, but, more importantly, also helps to reveal the replication mechanism of DNA. To further study the characteristics of DNA, nonlinear science has been used to deal with DNA system, because its properties can be investigated by nonlinear models that combine the methods of physics with biological tools [2]. The soliton theory especially which pertains to nonlinear science has been widely applied in the study of DNA [3–6].
To study DNA structure from the view of nonlinear science, it is necessary to look for the right nonlinear mathematical models. A number of researchers have tried to establish mathematical models to describe DNA system. At the beginning, it is difficult to use a specific mathematical model to simulate DNA system due to its complex structure and the presence of various movements [7]. So some simplified models were used to study only some internal movement in DNA [8–12]. Further, two movements which made the main contribution to the DNA denaturation process were taken into account together in one model. For example, a two-dimensional discrete model was used to describe the two movements in DNA [13]. Another new double-chain model of DNA that consists of two long elastic homogeneous strands was given to describe transverse movements along the hydrogen bond and longitudinal movements along the backbone direction [14, 15].
In [14, 15], the nonlinear dynamical equations were derived and nonlinear dynamics of DNA were studied. The exact solutions of the general nonlinear dynamic system in the double-chain model of DNA were obtained and discussed numerically and analytically under some special approximate conditions. The nonlinear dynamical equations for the double-chain model were given as follows:
(1)utt-c12uxx=λ1u+γ1uv+μ1u3+β1uv3,vtt-c22vxx=λ2v+γ2u2+μ2u2v+β2v3+c0,
where u denotes the difference of the longitudinal displacements of the bottom and top strands, that is, the displacements of the bases from their equilibrium positions along the direction of the phosphodiester bridge that connects the two bases of the same strands; v represents the difference of the transverse displacements of the bottom and top strands, that is, the displacements of the bases from their equilibrium positions along the direction of the hydrogen bond that connects the two bases of the base pair. The constants c1, c2, λ1, λ2, γ1, γ2, μ1, μ2, β1, β2, and c0 in (1) are written as
(2)c1=±Yρ,c2=±Fρ,λ1=-2μρσh(h-l0),λ2=-2μρσ,γ1=2γ2=22μl0ρσh2,μ1=μ2=-2μl0ρσh3,β1=β2=4μl0ρσh3,c0=2μ(h-l0)ρσh,
where ρ, σ, Y, F, μ, l0, and h are, respectively, the mass density, the area of transverse cross-section, the Young modulus, the tension density of the strand, the rigidity of the elastic membrane, the distance of the two strands, and the height of the membrane in the equilibrium position.
In [16], introducing the transformation v=au+b (a and b are constants) in (1) yields
(3)utt-c12uxx=u3(μ1+β1a2)+u2(2β1ab+aγ1)+u(λ1+bγ1+β1b2),utt-c22uxx=u3(μ2+β2a2)+u2(γ2a+μ2ba+3β2ab)+u(λ2+3β2b2)+λ2ba+β2b3a+c0a.
Comparing (3) and using the physical quantities (2), it follows that b=h/2 and F=Y; (3) can be written as
(4)utt-c12uxx=Au3+Bu2+Cu+D,
where A=((-2α/h3)+(4a2α/h3)), B=62aα/h2, C=((-2α/l0)+(6α/h)), and D=0 with α=μl0/ρσ.
Reference [16] studied the traveling wave solutions of (4) and obtained one particular family of solitary kink-type solutions for different values of Riccati parameter of (4). Unfortunately the results are not complete, because (4) admits tanh-type, periodic-type, and bell-type solutions; [16] only presented the tanh-type solutions.
In this paper, we use the bifurcation method of dynamical systems [17–20] to study bell-type and periodic-type solutions of (6), which can improve the results in [16]. We look for solutions of (4) in the form of u(x,t)=φ(ξ) with ξ=(x/c1)-2t. Substituting u(x,t)=φ((x/c1)-2t) into (4), we have
(5)φξξ=Aφ3+Bφ2+Cφ,
which is equivalent to the following two-dimensional system:
(6)dφdξ=y,dydξ=Aφ3+Bφ2+Cφ.
Obviously, system (6) has the first integral
(7)H(φ,y)=14Aφ4+13Bφ3+12Cφ2-12y2=M,
where M is the Hamiltonian constant.
In fact, system (6) is a planar dynamical system whose phase orbits defined by the vector field of (6) determine traveling wave solutions, so we first need to investigate some different phase portraits of system (6) in different parameter regions. Then, by using the theory of dynamical system, we will identify some different wave profiles and compute exact representations for solitary wave solutions and periodic wave solutions.
This paper is organized as follows. In Section 2, we draw bifurcation phase portraits of system (6), where explicit parametric conditions will be derived. In Section 3, we give exact representations of solitary wave and periodic wave solutions of (6) in explicit form. A short conclusion will be given in Section 4.
2. Properties of Singular Points and Bifurcation Phase Portraits
Firstly, we compute the equilibrium points of system (6) under different parametric conditions. It is clear that neither the parameters B nor C takes value zero from the forms of parameters A, B, and C; then A equals zero if a=±2/2. So we give equilibrium points in different parametric cases as follows.
When A=0, there are two equilibrium points (0,0) and (-C/B,0).
When A≠0 and δ1>0, there are three equilibrium points (0,0), (φ-,0), and (φ+,0) of (6) on the φ-axis, where δ1=B2-4AC, φ-=(-B-δ1)/2A, and φ+=(-B+δ1)/2A.
Define f0(φ)=Aφ3+Bφ2+Cφ; then the type of equilibrium point (φ,0) can be decided by λ±=f0′(φ), which are eigenvalues of linearized system of (6) at point (φ,0).
Secondly, we compute the intersection points of phase portraits of (6) on the φ-axis when the orbits pass through the equilibrium point (0,0), namely, when the Hamiltonian constant M=h(0,0)=0. Let f(φ)=(1/2)Aφ4+(2/3)Bφ3+Cφ2, δ2=(4/9)B2-2AC, φ0=-3A/2B, φ1=((-2/3)B-δ2)/A, and φ2=((-2/3)B+δ2)/A. When A=0, there are a double zero point 0 and a zero point φ0 of f(φ) on the φ-axis. When A≠0 and δ2>0, there are a double zero point 0 and two zero points φ1 and φ2 of f(φ) on the φ-axis. In order to draw bifurcation phase portraits of system (6), we need to discuss the sign of parameters A, B, and C. The intersection points of phase portraits of (6) on the φ-axis are given as follows.
If A=0, there are two intersection points (0,0) and (φ0,0) of (6) on the φ-axis (see Figure 1).
If A>0, there are three intersection points (0,0), (φ1,0), and (φ2,0) of (6) on the φ-axis and φ1<φ2 (see Figure 2).
If A<0, there are three intersection points (0,0), (φ1,0), and (φ2,0) of (6) on the φ-axis and φ2<φ1 (see Figure 3).
The bifurcation phase portraits of system (6) with A=0.
B>0, C>0
B>0, C<0
B<0, C>0
B<0, C<0
The bifurcation phase portraits of system (6) with A>0.
B>0, C>0
B>0, C<0
B<0, C>0
B<0, C<0
The bifurcation phase portraits of system (6) with A<0.
B>0, C>0
B>0, C<0
B<0, C>0
B<0, C<0
Finally, from the above analysis about system (6), we give twelve different phase portraits of (6) which are shown in Figures 1–3. There exist homoclinic orbits and periodic orbits in some phase portraits of Figures 1–3, which correspond to bell-type and periodic-type waves and are of special interest to us.
3. Exact Explicit Representations of Solitary Wave Solutions
The case A=0, B>0, and C>0 (see Figure 1(a)).
Corresponding to the homoclinic orbit of (6) passing through points (φ0,0) and (0,0) on the φ-y plane, there exists a smooth solitary wave solution. Now, (7) becomes
(8)y2=23Bφ3+Cφ2(φ0⩽φ⩽0).
Choosing φ(0)=φ0, substituting (8) into the first equation of (6), and integrating along the homoclinic orbit, we have
(9)∫φφ01s(2/3)Bs+Cds=|ξ|.
Solving (9) yields
(10)u1=-3CB(cosh(C|ξ|)+1).
The case A=0, B<0, and C>0 (see Figure 1(c)).
Similarly, in Figure 1(c) the homoclinic orbit has the following expression:
(11)y2=23Bφ3+Cφ2(0⩽φ⩽φ0).
Thus the corresponding smooth solitary wave solution is obtained as (10) and its profile is shown in Figure 4(b).
The case A>0, B>0, and C>0 (see Figure 2(a)).
Smooth solitary wave and periodic wave solutions.
In this case, on the φ-y plane the homoclinic orbit of (6) passing through points (φ2,0) and (0,0) has the following expression:
(12)y2=12Aφ4+23Bφ3+Cφ2(φ2⩽φ⩽0).
Choosing φ(0)=φ2, substituting (12) into the first equation of (6), and integrating along the homoclinic orbit, we have
(13)∫φφ21s(1/2)As2+(2/3)Bs+Cds=|ξ|.
Solving (13), it follows that
(14)u2=4Cα1α12e-C|ξ|+((4/9)B2-2AC)eC|ξ|-(4/3)Bα1,
where α1=(2(-2B2+9AC+B4B2-18AC))/(-6B+34B2-18AC). The profiles of solutions u1 and u2 both have figure caption Figure 4(a).
The case A>0, B<0, and C>0 (see Figure 2(c)).
In this case, on the φ-y plane the homoclinic orbit of (6) passing through points (φ1,0) and (0,0) has the following expression:
(15)y2=12Aφ4+23Bφ3+Cφ2(0⩽φ⩽φ1).
Choosing φ(0)=φ1, substituting (15) into the first equation of (6), and integrating along the homoclinic orbit, we have
(16)∫φφ11s(1/2)As2+(2/3)Bs+Cds=|ξ|.
Solving (16), we get
(17)u3=4Cα2α22eC|ξ|+((4/9)B2-2AC)e-C|ξ|-(4/3)Bα2,
where α2=(2(2B2-9AC+B4B2-18AC))/(6B+34B2-18AC). The profile of u3 is shown in Figure 4(d).
The case A<0, B>0, and C>0 and A<0, B<0, and C>0 (see Figures 3(a) and 3(c)).
In this case, there exist two homoclinic orbits and their expressions are (12) and (15), respectively. Similarly we get the corresponding solitary wave solutions as u2 and u3, which coexist for the same speed. The plane graphs of the coexisting solitary wave solutions are shown in Figure 4(c), when A<0, B>0, and C>0, and Figure 4(d), when A<0, B<0, and C>0.
The case A<0, B>0, and C<0 (see Figure 3(b)).
In this case, on the φ-y plane the orbit of (6) which passes through points (φ1,0) and (φ2,0) is a periodic orbit; it has the following expression:
(18)y2=12Aφ4+23Bφ3+Cφ2,(0<φ2⩽φ⩽φ1).
Choosing φ(0)=φ1, substituting (18) into the first equation of (6), and integrating along the periodic orbit, we have
(19)∫φφ11s(1/2)As2+(2/3)Bs+Cds=|ξ|.
Solving (19) we obtain the following periodic wave solution:
(20)u4=6C4B2-18ACcos(-Cξ)-2B.
The profile of u4 is shown in Figure 4(e).
The case A<0, B<0, and C<0 (see Figure 3(d)).
Similarly, corresponding to the periodic orbit of (6) which passes through points (φ1,0) and (φ2,0), there exists a periodic wave solution. Now, the expression of the periodic orbit is
(21)y2=12Aφ4+23Bφ3+Cφ2,(φ2⩽φ⩽φ1<0).
Choosing φ(0)=φ2, substituting (21) into the first equation of (6), and integrating along the periodic orbit, we have
(22)∫φ2φ1s(1/2)As2+(2/3)Bs+Cds=|ξ|.
Solving (22), it follows that
(23)u5=6C-4B2-18ACcos(-Cξ)-2B.
The profile of u5 is shown in Figure 4(f).
4. Conclusion
In this paper, we employ the method of dynamical systems to study the nonlinear dynamical equations of DNA model. The bifurcation phase portraits of the DNA under some parametric conditions are drawn, and explicit exact expressions of bell-shaped solitary waves and periodic waves are obtained via some special homoclinic and periodic orbits; their planar graphs are simulated. These solutions are new and different from those in [14–16]; some existing results are improved. Furthermore, in general, the larger the wave speed, the higher the wave crest. So for the same wave speed, it is interesting that the bell-shaped solitary waves of different heights are found to coexist.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (no. 11226303) and Guangdong Province (no. 2013KJCX0189). The authors thank the editors for their hard working and also gratefully acknowledge helpful comments and suggestions by reviewers.
WatsonJ. D.CrickF. H. C.Molecular structure of nucleic acids: a structure for deoxyribose nucleic acid195317143567377382-s2.0-003849754210.1038/171737a0PeyrardM.LópezS. C.JamesG.Modelling DNA at the mesoscale: a challenge for nonlinear science?200821T9110.1088/0951-7715/21/6/T02ZBL1162.92308LipniackiT.Chemically driven traveling waves in DNA199960725310.1103/PhysRevE.60.7253YomosaS.Solitary excitations in deoxyribonuclei acid (DNA) double helices19843014744802-s2.0-2654445199010.1103/PhysRevA.30.474ZhangC. T.Soliton excitations in deoxyribonucleic acid (DNA) double helices198735288689110.1103/PhysRevA.35.886MR870770ForinashK.Nonlinear dynamics in a double-chain model of DNA1991431073410.1103/PhysRevB.43.10743YakushevichL. V.2004Berlin, GermanyJohn Wiley & SonsEnglanderS. W.KallenbanchN. R.HeegerA. J.KrumhanslJ. A.LitwinS.Nature of the open state in long polynucleotide double helices: possibility of soliton excitations1980777222722610.1073/pnas.77.12.7222YomosaS.Soliton excitations in deoxyribonucleic acid (DNA) double helices19832742120212510.1103/PhysRevA.27.2120MR700064TakenoS.HommaS.A coupled base-rotator model for structure and dynamics of DNA—local fluctuations in helical twist angles and topological solitons198472467969310.1143/PTP.72.679MR772719ZBL1074.92501PeyrardM.BishopA. R.Statistical mechanics of a nonlinear model for DNA denaturation19896223275527582-s2.0-424318731910.1103/PhysRevLett.62.2755DauxoisT.PeyrardM.BishopA. R.Entropy-driven DNA denaturation1993471R44R472-s2.0-1714440497210.1103/PhysRevE.47.R44MutoV.LomdahlP. S.ChristiansenP. L.Two-dimensional discrete model for DNA dynamics: longitudinal wave propagation and denaturation19904212745274582-s2.0-000134896010.1103/PhysRevA.42.7452KongD. X.LouS. Y.ZengJ.Nonlinear dynamics in a new double Chain-model of DNA2001366737742QianX.-M.LouS.-Y.Exact solutions of nonlinear dynamics equation in a new double-chain model of DNA2003394501505MR2001933AlkaW.GoyalA.Nagaraja KumarC.Nonlinear dynamics of DNA—riccati generalized solitary wave solutions201137534804832-s2.0-7865088074610.1016/j.physleta.2010.11.017ZBL1241.35205ChowS. N.HaleJ. K.1981Berlin, GermanySpringerMR660633LiJ. B.LiuZ. R.Smooth and non-smooth traveling waves in a nonlinearly dispersive equation200025415610.1016/S0307-904X(00)00031-7ZBL0985.37072LiJ. B.LiuZ. R.Traveling wave solutions for a class of nonlinear dispersive equations20022339710.1142/S0252959902000365ZBL1011.35014LiuZ. R.OuyangZ. Y.A note on solitary waves for modified forms of Camassa-Holm and Degasperis-Procesi equations20073664-537738110.1016/j.physleta.2007.01.074MR2570695ZBL1203.35234