(ProQuest: ... denotes non-US-ASCII text omitted.)
Zheng-yong Ouyang 1 and Shan Zheng 2
Academic Editor:Yonghui Xia
1, Department of Mathematics, Foshan University, Foshan Guangdong 528000, China
2, Guangzhou Maritime College, Guangzhou Guangdong 510725, China
Received 6 January 2014; Revised 24 February 2014; Accepted 10 March 2014; 3 April 2014
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.
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: [figure omitted; refer to PDF] 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 c 1 , c 2 , λ 1 , λ 2 , γ 1 , γ 2 , μ 1 , μ 2 , β 1 , β 2 , and c 0 in (1) are written as [figure omitted; refer to PDF] where ρ , σ , Y , F , μ , l 0 , 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 = a u + b ( a and b are constants) in (1) yields [figure omitted; refer to PDF] Comparing (3) and using the physical quantities (2), it follows that b = h / 2 and F = Y ; (3) can be written as [figure omitted; refer to PDF] where A = ( ( - 2 α / h 3 ) + ( 4 a 2 α / h 3 ) ) , B = 6 2 a α / h 2 , C = ( ( - 2 α / l 0 ) + ( 6 α / h ) ) , and D = 0 with α = μ l 0 / ρ σ .
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 ) = [straight phi] ( ξ ) with ξ = ( x / c 1 ) - 2 t . Substituting u ( x , t ) = [straight phi] ( ( x / c 1 ) - 2 t ) into (4), we have [figure omitted; refer to PDF] which is equivalent to the following two-dimensional system: [figure omitted; refer to PDF] Obviously, system (6) has the first integral [figure omitted; refer to PDF] 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.
(1) When A = 0 , there are two equilibrium points ( 0,0 ) and ( - C / B , 0 ) .
(2) When A ...0; 0 and δ 1 > 0 , there are three equilibrium points ( 0,0 ) , ( [straight phi] - , 0 ) , and ( [straight phi] + , 0 ) of (6) on the [straight phi] -axis, where δ 1 = B 2 - 4 A C , [straight phi] - = ( - B - δ 1 ) / 2 A , and [straight phi] + = ( - B + δ 1 ) / 2 A .
Define f 0 ( [straight phi] ) = A [straight phi] 3 + B [straight phi] 2 + C [straight phi] ; then the type of equilibrium point ( [straight phi] , 0 ) can be decided by λ ± = f 0 [variant prime] ( [straight phi] ) , which are eigenvalues of linearized system of (6) at point ( [straight phi] , 0 ) .
Secondly, we compute the intersection points of phase portraits of (6) on the [straight phi] -axis when the orbits pass through the equilibrium point ( 0,0 ) , namely, when the Hamiltonian constant M = h ( 0,0 ) = 0 . Let f ( [straight phi] ) = ( 1 / 2 ) A [straight phi] 4 + ( 2 / 3 ) B [straight phi] 3 + C [straight phi] 2 , δ 2 = ( 4 / 9 ) B 2 - 2 A C , [straight phi] 0 = - 3 A / 2 B , [straight phi] 1 = ( ( - 2 / 3 ) B - δ 2 ) / A , and [straight phi] 2 = ( ( - 2 / 3 ) B + δ 2 ) / A . When A = 0 , there are a double zero point 0 and a zero point [straight phi] 0 of f ( [straight phi] ) on the [straight phi] -axis. When A ...0; 0 and δ 2 > 0 , there are a double zero point 0 and two zero points [straight phi] 1 and [straight phi] 2 of f ( [straight phi] ) on the [straight phi] -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 [straight phi] -axis are given as follows.
(1) If A = 0 , there are two intersection points ( 0,0 ) and ( [straight phi] 0 , 0 ) of (6) on the [straight phi] -axis (see Figure 1).
(2) If A > 0 , there are three intersection points ( 0,0 ) , ( [straight phi] 1 , 0 ) , and ( [straight phi] 2 , 0 ) of (6) on the [straight phi] -axis and [straight phi] 1 < [straight phi] 2 (see Figure 2).
(3) If A < 0 , there are three intersection points ( 0,0 ) , ( [straight phi] 1 , 0 ) , and ( [straight phi] 2 , 0 ) of (6) on the [straight phi] -axis and [straight phi] 2 < [straight phi] 1 (see Figure 3).
The bifurcation phase portraits of system (6) with A = 0 .
(a) B > 0 , C > 0
[figure omitted; refer to PDF]
(b) B > 0 , C < 0
[figure omitted; refer to PDF]
(c) B < 0 , C > 0
[figure omitted; refer to PDF]
(d) B < 0 , C < 0
[figure omitted; refer to PDF]
The bifurcation phase portraits of system (6) with A > 0 .
(a) B > 0 , C > 0
[figure omitted; refer to PDF]
(b) B > 0 , C < 0
[figure omitted; refer to PDF]
(c) B < 0 , C > 0
[figure omitted; refer to PDF]
(d) B < 0 , C < 0
[figure omitted; refer to PDF]
The bifurcation phase portraits of system (6) with A < 0 .
(a) B > 0 , C > 0
[figure omitted; refer to PDF]
(b) B > 0 , C < 0
[figure omitted; refer to PDF]
(c) B < 0 , C > 0
[figure omitted; refer to PDF]
(d) B < 0 , C < 0
[figure omitted; refer to PDF]
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
(1) The case A = 0 , B > 0 , and C > 0 (see Figure 1(a)).
Corresponding to the homoclinic orbit of (6) passing through points ( [straight phi] 0 , 0 ) and ( 0,0 ) on the [straight phi] - y plane, there exists a smooth solitary wave solution. Now, (7) becomes [figure omitted; refer to PDF] Choosing [straight phi] ( 0 ) = [straight phi] 0 , substituting (8) into the first equation of (6), and integrating along the homoclinic orbit, we have [figure omitted; refer to PDF] Solving (9) yields [figure omitted; refer to PDF]
(2) 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: [figure omitted; refer to PDF] Thus the corresponding smooth solitary wave solution is obtained as (10) and its profile is shown in Figure 4(b).
(3) The case A > 0 , B > 0 , and C > 0 (see Figure 2(a)).
Smooth solitary wave and periodic wave solutions.
(a) [figure omitted; refer to PDF]
(b) [figure omitted; refer to PDF]
(c) [figure omitted; refer to PDF]
(d) [figure omitted; refer to PDF]
(e) [figure omitted; refer to PDF]
(f) [figure omitted; refer to PDF]
In this case, on the [straight phi] - y plane the homoclinic orbit of (6) passing through points ( [straight phi] 2 , 0 ) and ( 0,0 ) has the following expression: [figure omitted; refer to PDF] Choosing [straight phi] ( 0 ) = [straight phi] 2 , substituting (12) into the first equation of (6), and integrating along the homoclinic orbit, we have [figure omitted; refer to PDF] Solving (13), it follows that [figure omitted; refer to PDF] where α 1 = ( 2 ( - 2 B 2 + 9 A C + B 4 B 2 - 18 A C ) ) / ( - 6 B + 3 4 B 2 - 18 A C ) . The profiles of solutions u 1 and u 2 both have figure caption Figure 4(a).
(4) The case A > 0 , B < 0 , and C > 0 (see Figure 2(c)).
In this case, on the [straight phi] - y plane the homoclinic orbit of (6) passing through points ( [straight phi] 1 , 0 ) and ( 0,0 ) has the following expression: [figure omitted; refer to PDF] Choosing [straight phi] ( 0 ) = [straight phi] 1 , substituting (15) into the first equation of (6), and integrating along the homoclinic orbit, we have [figure omitted; refer to PDF] Solving (16), we get [figure omitted; refer to PDF] where α 2 = ( 2 ( 2 B 2 - 9 A C + B 4 B 2 - 18 A C ) ) / ( 6 B + 3 4 B 2 - 18 A C ) . The profile of u 3 is shown in Figure 4(d).
(5) 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 u 2 and u 3 , 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 .
(6) The case A < 0 , B > 0 , and C < 0 (see Figure 3(b)).
In this case, on the [straight phi] - y plane the orbit of (6) which passes through points ( [straight phi] 1 , 0 ) and ( [straight phi] 2 , 0 ) is a periodic orbit; it has the following expression: [figure omitted; refer to PDF] Choosing [straight phi] ( 0 ) = [straight phi] 1 , substituting (18) into the first equation of (6), and integrating along the periodic orbit, we have [figure omitted; refer to PDF] Solving (19) we obtain the following periodic wave solution: [figure omitted; refer to PDF] The profile of u 4 is shown in Figure 4(e).
(7) 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 ( [straight phi] 1 , 0 ) and ( [straight phi] 2 , 0 ) , there exists a periodic wave solution. Now, the expression of the periodic orbit is [figure omitted; refer to PDF] Choosing [straight phi] ( 0 ) = [straight phi] 2 , substituting (21) into the first equation of (6), and integrating along the periodic orbit, we have [figure omitted; refer to PDF] Solving (22), it follows that [figure omitted; refer to PDF] The profile of u 5 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.
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.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
[1] J. D. Watson, F. H. C. Crick, "Molecular structure of nucleic acids: a structure for deoxyribose nucleic acid," Nature , vol. 171, no. 4356, pp. 737-738, 1953.
[2] M. Peyrard, S. C. López, G. James, "Modelling DNA at the mesoscale: a challenge for nonlinear science?," Nonlinearity , vol. 21, pp. T91, 2008.
[3] T. Lipniacki, "Chemically driven traveling waves in DNA," Physical Review E , vol. 60, pp. 7253, 1999.
[4] S. Yomosa, "Solitary excitations in deoxyribonuclei acid (DNA) double helices," Physical Review A , vol. 30, no. 1, pp. 474-480, 1984.
[5] C. T. Zhang, "Soliton excitations in deoxyribonucleic acid (DNA) double helices," Physical Review A , vol. 35, no. 2, pp. 886-891, 1987.
[6] K. Forinash, "Nonlinear dynamics in a double-chain model of DNA," Physical Review B , vol. 43, pp. 10734, 1991.
[7] L. V. Yakushevich Nonlinear Physics of DNA , John Wiley & Sons, Berlin, Germany, 2004.
[8] S. W. Englander, N. R. Kallenbanch, A. J. Heeger, J. A. Krumhansl, S. Litwin, "Nature of the open state in long polynucleotide double helices: possibility of soliton excitations," Proceedings of the National Academy of Sciences of the United States of America , vol. 77, pp. 7222-7226, 1980.
[9] S. Yomosa, "Soliton excitations in deoxyribonucleic acid (DNA) double helices," Physical Review A , vol. 27, no. 4, pp. 2120-2125, 1983.
[10] S. Takeno, S. Homma, "A coupled base-rotator model for structure and dynamics of DNA--local fluctuations in helical twist angles and topological solitons," Progress of Theoretical Physics , vol. 72, no. 4, pp. 679-693, 1984.
[11] M. Peyrard, A. R. Bishop, "Statistical mechanics of a nonlinear model for DNA denaturation," Physical Review Letters , vol. 62, no. 23, pp. 2755-2758, 1989.
[12] T. Dauxois, M. Peyrard, A. R. Bishop, "Entropy-driven DNA denaturation," Physical Review E , vol. 47, no. 1, pp. R44-R47, 1993.
[13] V. Muto, P. S. Lomdahl, P. L. Christiansen, "Two-dimensional discrete model for DNA dynamics: longitudinal wave propagation and denaturation," Physical Review A , vol. 42, no. 12, pp. 7452-7458, 1990.
[14] D. X. Kong, S. Y. Lou, J. Zeng, "Nonlinear dynamics in a new double Chain-model of DNA," Communications in Theoretical Physics , vol. 36, no. 6, pp. 737-742, 2001.
[15] X.-M. Qian, S.-Y. Lou, "Exact solutions of nonlinear dynamics equation in a new double-chain model of DNA," Communications in Theoretical Physics , vol. 39, no. 4, pp. 501-505, 2003.
[16] W. Alka, A. Goyal, C. Nagaraja Kumar, "Nonlinear dynamics of DNA--riccati generalized solitary wave solutions," Physics Letters A: General, Atomic and Solid State Physics , vol. 375, no. 3, pp. 480-483, 2011.
[17] S. N. Chow, J. K. Hale Methods of Bifurcation Theory , Springer, Berlin, Germany, 1981.
[18] J. B. Li, Z. R. Liu, "Smooth and non-smooth traveling waves in a nonlinearly dispersive equation," Applied Mathematical Modelling , vol. 25, pp. 41-56, 2000.
[19] J. B. Li, Z. R. Liu, "Traveling wave solutions for a class of nonlinear dispersive equations," Chinese Annals of Mathematics , vol. 23, pp. 397, 2002.
[20] Z. R. Liu, Z. Y. Ouyang, "A note on solitary waves for modified forms of Camassa-Holm and Degasperis-Procesi equations," Physics Letters A , vol. 366, no. 4-5, pp. 377-381, 2007.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2014 Zheng-yong Ouyang and Shan Zheng. Zheng-yong Ouyang 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
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.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer