Advances in Mathematical Physics

Volume 2018, Article ID 4676281, 7 pages

https://doi.org/10.1155/2018/4676281

## Coupled Sine-Gordon Systems in DNA Dynamics

^{1}Department of Preventive Medicine, School of Medicine, Ningbo University, Ningbo 315211, China^{2}Department of Mathematics, Ningbo University, Ningbo 315211, China

Correspondence should be addressed to Chuanzhong Li; nc.ude.ubn@gnohznauhcil

Received 14 March 2018; Accepted 16 April 2018; Published 15 May 2018

Academic Editor: Zhijun Qiao

Copyright © 2018 Liya Liu and Chuanzhong Li. 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 a model called the coupled sine-Gordon equation for DNA dynamics by introducing two double helix structures. The second double helix structure is unilaterally influenced by the first one. The completely integrable coupled sine-Gordon equation admits kink-antikink solitons with increased width representing a wide base pair opening configuration in DNA. Also we propose another coupled sine-Gordon model with variable coefficients for DNA dynamics under an inhomogeneous background. We find that the inhomogeneous DNA model has many interesting localized nonrational rogue wave solutions. We can find that the appearance of the rogue waves (possibly means the genetic mutation) in the nonlinear DNA model is highly related to the inhomogeneity.

#### 1. Introduction

The sine-Gordon equation is an important integrable equation which describes many interesting phenomena including the DNA double helix molecule [1], dynamics of coupled pendulums, Josephson junction arrays [2], living cellular structures [3], and the crystalline lattices [4]. In the DNA double helix structure, the two parallel chains of nucleotides (strands) connect with dipole-dipole interaction along the strands. The two strands are bound to each other through hydrogen bonds between the complementary bases [5].

Besides in life sciences, the scalar sine-Gordon equation originates in differential geometry and has profound applications in physics which can be seen from a recent review [6]. In [7], the nonrational rogue waves can be induced by the inhomogeneity background of the sine-Gordon equation.

Similarly as the method in our previous work in [8, 9], we construct the general -sine-Gordon and -sinh-Gordon systems which contain many multicomponent sine-Gordon type and sinh-Gordon type equations in [10]. Meanwhile, we give the Bäcklund transformations of the -sine-Gordon and the nonlinear superposition formula and Lax pairs.

This paper is arranged as follows. In Section 2, we recall some basic facts about the classical sine-Gordon equation including its Lax pair and Bäcklund transformation. In Section 3, we recall the coupled sine-Gordon equation and its Lax pair. After this, we will construct the Bäcklund transformation of the inhomogeneous coupled sine-Gordon equation and, from the Bäcklund transformation, we derive many solutions including the rogue wave solutions of some special inhomogeneous coupled sine-Gordon equation.

#### 2. Sine-Gordon Equation

Under certain circumstances, the DNA can open the hydrogen bonds between conjugated pairs, exposing the unpaired bases to the action of external ligands. To model this motion each base of a strand is assumed to be coupled with the next-neighbour bases of the same strand by the elastic backbone restoring forces and with the complementary base in the opposite strand by the anharmonic potential used to model the hydrogen bond. One take as canonical variables the deflection angles and that two complementary bases form with the line passing from the attaching points of the bases to the strands.

In [11], they construct a Hamiltonian of DNA dynamics:by introducing the helical structure through twist deformation of double helix structure in analogy with the structure of helimagnet and cholesteric liquid crystal system. Here is the moment of inertia of the bases around the axes at the th site and and , respectively, represent the ferromagnetic spin-spin exchange interaction and uniaxial magnetocrystalline anisotropy with the easy axis along -direction. denotes the elastic constant associated with the twist deformation and is a constant. Having formed the Hamiltonian, the dynamics of the DNA with double helix molecule can be understood by constructing the Hamilton’s equations of motion corresponding to the Hamiltonian (1) asEquations (2a) and (2b) describe the dynamics of the DNA double helix at the discrete level when the helical nature of the molecule is represented in the form of a twist-like deformation.

Under the approximation by supposing that the difference in angular rotation of bases with respect to neighbouring bases along the two strands is small, the discrete equations of motion after redefinition of the time and the parameter reduce toin the continuum limit. After suitable rescaling of the variable , (3a) and (3b) becomewhere and we have further chosen .

By introducing two variables , , the equation can be simplified to the following well-known sine-Gordon equation:

##### 2.1. Lax Equations and Bäcklund Transformations

Let us explore the source of the matrix form of the sine-Gordon equation (5). The mathematical method of the soliton equation was derived mainly through zero curvature equation. The zero curvature equation is as follows.

We can export many soliton equations by choosing appropriate values of and . The sine-Gordon equation has the following form withHere the zero curvature equation (6) becomesFor example, we choose appropriate values of , , and aswhereThen we can get the well-known sine-Gordon equation as follows:

is a solution of (11), and, under the following Bäcklund transformation, is another solution of (11).

#### 3. Coupled Sine-Gordon Equation

In this section, we will consider the dynamics of two DNA double helix molecules by constructing the Hamilton’s equations of motion as

Equations (13a)–(14b) can possibly describe the dynamics of the interaction of two DNA double helix structures with each other at the discrete level in the form of a twist-like deformation.

Hence, under the small angle approximation, in the continuum limit, that is,the discrete equations of motion (13a)–(14b) after suitable rescaling of time and redefinition of the parameter reduce toAdding and subtracting (17a) and (17b) and after suitable rescaling of the variable , we obtainwhere , , and we have further chosen .

##### 3.1. Lax Equations of the Coupled Sine-Gordon Equation

In this section, we will consider the following coupled sine-Gordon equation: The Lax equation of the coupled sine-Gordon equation is as follows [10]:whenwe can export the Lax equation of the coupled sine-Gordon equation.

#### 4. Inhomogeneous Coupled Sine-Gordon Equation

The inhomogeneous system is modeled by the following coupled sine-Gordon equation with variable coefficients:

Here represents the inhomogeneity of the system. If , then the following transformation: , will map the inhomogeneous coupled sine-Gordon equation to the following standard coupled sine-Gordon equation:which has the following Bäcklund transformation:

There is no doubt that it is also the Bäcklund transformation of the coupled sine-Gordon equation when , .

Letting , then the first equation of the Bäcklund transformation is as follows:and the second sets of equation is as follows:Then we derive the following new solution:

From the known three solutions ,we can get the following new solution using the nonlinear superposition formula similarly as in [10]:

In this way, by the algebraic iterated operation we can get many new solutions of the coupled sine-Gordon equation and its inhomogeneous system.

When , , the solutions are exactly the solutions of the coupled sine-Gordon equation [10] whose graph can be seen in Figures 1 and 2 with a density plot of in Figure 3. From these graphs, we can find that they are all traveling wave solutions such as kink and soliton solutions with a wide base pair opening configuration in DNA.