Mobile breathers in a Nonlinear model for DNA breathing

71 Eclética Química Journal, 42, 2017 ISSN: 1678-4618 DOI: 10.26850/1678-4618eqj.v42.1.2017.p71-75 ABSTRACT: Objectives. Analyze the DNA dynamics in PeyrardBishop-Dauxois model (PBD) with different control parameters using its energy center of the mobile “breather”. Materials and methods. We used the Peyrard-Bishop-Dauxois mathematical model and the MATLAB software for studying the DNA dynamic using Morse potential, Symmetric Morse and the “hump” potential for simulating the interactions which arise the pile up. Results. It has been observed that the analytical and computational methods allow to detect the influence of the potentials of the PBD model in the behavior of the energy center in the presence of a couple of base A(adenine) or T(thymine) using the control of parameter α=-0.30 and velocity of mobile breather: v0=0.1. In the case of Morse potential, the center of energy respect to the mobile breather undergoes a change in its trajectory and produce a DNA breathing. Conclusions. Analytical and computational approaches can be used for obtaining differences respect to the DNA dynamics using different control parameters: velocity of BM and inhomogeneity. The potential “hump” may decrease the reflective effect with the indicated parameters to the effect on the energy center to the mobile breather. Mobile breathers in a Nonlinear model for DNA breathing


Introduction
DNA is an important basis of molecular biology 1 .The dynamic in DNA depends on the physical parameters that modify the behavior of the base pairs and show large amplitude fluctuations known as the "breathing" of DNA 2 .We have some references about this dynamics in the model of Peyrard and Bishop [3][4][5][6][7][8][9][10][11][12][13][14] .This model exhibits spatially localized oscillations denoted as "discrete breathers" 15 .Flach, A. Gorbach 16 give the theory of the discrete breathers (DB) and applications.
Forinash, Peyrard and Malomed 17 proved the influence of an impurity in the movement of discrete breathers in Klein-Gordon (KG) chains.The method used for calculating discrete breathers in (KG) lattices consists on approximating the dynamical equations by a Nonlinear Schrodinger Equation (NLS) and moving the breather by an envelope solution 17 .
Cuevas, Palmero, Romero and Archilla have found that the moving breather can cross the impurity, and can be reflected by it, or can be trapped 18,19 .
The problem is to determine the parameters of the model to generate a trapped breather (TB) which is very important for the DNA transcription.The results obtained in our simulation verify the existence of breather trapped with the conditions described by Cuevas, Palmero and Archilla 18 .The trapped breather (TB) depends on the impurity and the initial velocity of mobile breather.
In addition, we have also observed that the PBD model with "hump" Morse potential and inappropriate initial velocity have different vibration.We conjecture that "hump" Morse potential with appropriate initial velocity of mobile breather and impurity can generate a trapped breather.The results of the simulations can be applied to study of moving breather collisions in the PBD model (for a review, see 20 ).

The PBD model with modified Morse potentials and the inhomogeneity for the Morse potential
At a first stage, the dynamics of this PBD model is analyzed by the equations: . u velocity  ; u = stretching of the hydrogen bonds N= Number of the pairs of basis. ( The associate equations for (4) are a system of equations (n=1, 2... N) being the derivative of the potential) We consider potentials of the form: V is the potential and D n represents a discontinuity in the bases.In case of a discontinuity on the site only n 0 , D n is the form and the parameter α is the non-homogeneity (magnitude of the impurity) and considered in the interval [-1,∞ ) and D 0 a constant.
In the case of the Potential of Morse are: D n is the depth of potential and a n represents the amplitude of the potential of Morse of the site n.
We consider the potentials:

Numerical solutions of the equations (5)
Using the approximation for the oscillator n and and substituting in (5) one has (2 ) 0 which depend on the parameter K, and V' n k is the k th Fourier coefficient for the periodic function V'(u n (t)).
In Figure 1   the maximum amplitude of the vibrations is 0.6.In biological terms, this change can be the effect of a disease or pharmaceutical control treatment.

STABILITY: Harmonic bifurcation
We analyze the stability of the breather solution.
Let us introduce a function ũ n (t) = u n (t) +  n (t), where u n (t) is the periodic breather solution shown in Figure 1 and ε n (t) is a perturbation: ũ n (t) must satisfy the system (equations 5) and expanding around u n (t) to first order (linearization), we obtain the following system of equations for ε n (t) We can associate a monodromy matrix for this equation with Floquet multipliers 7 .

Remark 2
The solution is stable if the modules of Floquet multipliers are one.The especial instability ("harmonic bifurcation") in our case happens when a pair of Floquet multipliers merges at λ=1 and splits off circle onto the positive real axis, in Figure 2.

Breather mobile and breather trapped
Once the system of equations ( 5) is worked out by Runge Kutta method and has seen its stability till the parameter K=0.13, one construct in breather mobile considering the initial conditions of the position and average speed of each position "n" respect to the harmonic oscillation corresponding to the DNA and one get construct its corresponding center of energy of the breather mobile 13 (see Figure 3).The center of energy of the breather mobile is given by 4 where the density energy is the form For the case of Morse potential the breather trapped is obtained with the additional condition: α, ε (-0.49, -0.02) and the initial velocity of BM: v o =0.1 18 .
It is very important the initial velocity of the BM for the displacement a long of sites of DNA and can be produced of DNA breathing.
This couple of parameters (α, v 0 ) is transcendental for DNA breathing.

Test of simulation for the evolution of the center of energy with modified Morse potential and Morse potential
The dynamic respect to the DNA with N=251 a couple of bases is watched more precisely analyzing its center of energy for two types of potential considering an impurity in the position n 0 = 26 and with the control parameter α = -0.30,D n =0.5 (n≠ n 0) .The interaction respect to the center of energy for two types of potential between the model of Morse and the model "hump" Morse is showed in the  Remark 4 In the reference 13 we have the complete study of the Figure 4 considering the sequence of Insulin with N=251 pairs of basis.In Figure 4 we have the trapped case for the interaction of the center of energy when we consider the Morse potential and the position of the inhomogeneity.
Remark 5 In Figure 4 we do not have the reflected case for the interaction of the center of energy when we consider the "hump" Morse potential.

Discussion
In the Remark 3: The maximum amplitude of the breather mobile has 0.6 angstrom for the symmetric Morse potential.In this case is interesting to simulate interactions of the center of energy with inhomogeneity.We have obtained an appropriate center of energy for the initial velocity of 0.1 12 .
In the Remark 4: We have shown the center of energy is trapped with the control in the inhomogeneity n 0 =26 for the Morse potential.This case is favorable for the DNA transcription because we need the breather trapping 18 .
In the Remark 5: We have shown the center of energy is not trapped for the "hump" Morse potential.This case is not favorable for the DNA transcription because we need the breather trapping 18 .In this case is necessary to introduce the control parameter appropriate of inhomogeneity for obtaining a breather trapping using the center of energy 18 .Another control parameter is the initial velocity of the breather mobile as showed in reference 16 .For Figure 4, it is firstly necessary an appropriate initial velocity for the simulations of the center of energy in interactions with inhomogeneity (for a review, see 12 and Figure 5).Our results are similar as obtained the reference 18 with the control of parameters for the BT (Breather trapped): 1. Initial velocity for the BM with the parameter λ =0.1 2. Inhomogeneity α =-0.30 for the Morse potential in the PBD model.
The proposed method can be applied to verify the breather rebounds when the parameters α, ε (-0.54, -0.49) 18 and using another initial velocity for mobile breather.
Even though, a reflection is observed after a short period of time in figure 4. In fact, many breathers trapped collide and facilitate the transcript process (for a review, see 20 ).
The theory of control optimum can be used for explore the effective breather trapped.
Finally the breather trapped act as a precursor of the transcription bubble of DNA 21 .These BTs can interact with another breathers and collect energy 22 .

Conclusions
This work has proposed state influence in the dynamic of DNA using the potentials, initial velocity of the breather mobile and the inhomogeneity as parameters of control.Next we have obtained different amplitudes of the breather solution depending of the potential.Finally, it is very important to emphasize the importance of the analytical and computational methods in the study of the components of DNA.
Next, the interaction of the breather mobile with the inhomogeneity of the DNA has been proposed for performing changes in the trajectory of the center of energy for the Morse potential with implications in the transcription process of DNA.
is depicted the example for the potential (equation 9) the solution of equations (5) denominated breather.

Figure 1 .Remark 1
Figure 1.The solution of equations (5) for the case of potential (equation 9) with K=0.3 and wb=0.8.Remark 1 The dynamics of the DNA is a set of coupled oscillators, and the vibrational motion is equivalent to equations (5) which depend on the several control parameters: potential modified of Morse and constant K of coupling.Parametric changes in the potentials imply changes in the amplitudes of the vibrations.For example, for the potential (equation 9),

Figure 3 .
Figure 3. Breather mobile for the Symmetric Morse potential.Remark 3 The maximum amplitude of the breather mobile has 0.6 angstrom for the symmetric Morse potential.

Figure 4 .
Figure 4. Evolution of the center of energy considering two types of Morse potentials with the parameters: N=251, K=0.13 and the position of the inhomogeneity n 0 =26 for the Morse potential.

Figure 5 .
Figure 5. Evolution of the center of energy considering of symmetric Morse potential.with the velocities: 0.01, 0.1 and 1.