Energy localization in the Peyrard-Bishop DNA model

High-order nonlinear excitations in the Joyeux-Buyukdagli model of DNA

By means of the semidiscrete multiple-scale method, we study the existence and properties of high-order envelope solitons and discrete breathers in a homogeneous DNA chain model that is based on pairing enthalpies and site-dependent finite stacking. We obtain the analytical solutions for an envelope soliton, and find that at the Brillouin zone center, discrete breather solutions can appear below the bottom of the phonon band. The behavior of two solitons in collisions and the stability of discrete breathers are confirmed by numerical simulations of the exact equations of the system.

Compact-envelope bright solitary wave in a DNA double strand

We study the nonlinear dynamics of a homogeneous DNA chain that is based on site-dependent finite stacking and pairing enthalpies. We introduce an extended nonlinear Schrödinger equation describing the dynamics of modulated wave in DNA model. We obtain envelope bright solitary waves with compact support as a solution. Analytical criteria of existence of this solution are derived. The stability of bright compactons is confirmed by numerical simulations of the exact equations of the lattice. The impact of the finite stacking energy is investigated, and we show that some of these compact bright solitary waves are robust, while others decompose very quickly depending on the finite stacking parameters.

Multifractal analysis of thermal denaturation based on the Peyrard-Bishop-Dauxois model

The theory of DNA dynamics is exceedingly complex and not easily explained. In the past two decades, by adapting methods of statistical physics, the dynamics of DNA in contact with a thermal bath is widely studied. In this paper, the thermal denaturation of DNA in the framework of the Peyrard-Bishop-Dauxois (PBD) model through the Rényi dimension is investigated. As a result, the Rényi dimension spectrum of the melting transition process reveals the multifractal nature of the dynamics of the Peyrard-Bishop-Dauxois model. Also, it can be concluded that the Rényi dimension (D(q)) at negative values of q is the characteristic signature of pre-melting and thermal denaturation of DNA. Furthermore, this approach is in excellent agreement with previous experimental studies.

Localization and delocalization of energy in a Peyrard-Bishop chain

This paper studies energy localization conditions in lattices of the type proposed by Peyrard and Bishop. Homogeneous and inhomogeneous lattices are analyzed and the role of interfaces in the latter is emphasized. Simulations allowed us to identify critical energy values for the existence of localization. After a certain energy value, it is possible to observe the loss of energy localization along the chain.

Thermal equivalence of DNA duplexes for probe design

We present the theory of thermal equivalence in the framework of the Peyrard-Bishop model and some of its anharmonic variants. The thermal equivalence gives rise to a melting index τ which maps closely the experimental DNA melting temperatures for short DNA sequences. We show that the efficient calculation of the melting index can be used to analyse the parameters of the Peyrard-Bishop model and propose an improved set of Morse potential parameters. With this new set we are able to calculate some of the experimental melting temperatures to ± 1.2 °C. We review some of the concepts of sequencing probe design and show how to use the melting index to explore the possibilities of gene coverage by tuning the model parameters.

A model coupling vibrational and rotational motion for the DNA molecule

We investigate a mechanical model for the DNA molecule using an extension of the Peyrard and Bishop model. In the present model, there are two chains of oscillators linked by a Morse potential, which represent the hydrogen bonds. The rotational and vibrational motions of each base pair are considered and the coupling for these motions are introduced by a nonlinear combination of them in the Morse potential. In this context, thermodynamics and structural properties are discussed.

The relevance of nonlinear stacking interactions in simple models of double-stranded DNA

Single molecule DNA experiments provide interesting data that allow a better understanding of the mechanical interactions between the strands and the nucleotides of this molecule. In some sense, these experiments complement the classical ones about DNA thermal denaturation. It is well known that the original Peyrard-Bishop (PB) model by means of a harmonic stacking potential and a nonlinear substrate potential has been able to predict the existence of a critical temperature of full denaturation of the molecule. In the present paper, driven by the findings of single molecule experiments, we substitute the original harmonic intra-strand stacking potential with a Duffing type potential. By elementary and analytical arguments, we show that with this choice it is possible to obtain a sharp transition in the classical domain wall solution of the PB model and the compactification of the classical solitary wave solutions of other models for the dynamics of DNA. We discuss why these solutions may improve our knowledge of the DNA dynamics in several directions.

Composite model for DNA torsion dynamics

DNA torsion dynamics is essential in the transcription process; a simple model for it, in reasonable agreement with experimental observations, has been proposed by Yakushevich (Y) and developed by several authors; in this, the nucleotides (the DNA subunits made of a sugar-phosphate group and the attached nitrogen base) are described by a single degree of freedom. In this paper we propose and investigate, both analytically and numerically, a "composite" version of the Y model, in which the sugar-phosphate group and the base are described by separate degrees of freedom. The model proposed here contains as a particular case the Y model and shares with it many features and results, but represents an improvement from both the conceptual and the phenomenological point of view. It provides a more realistic description of DNA and possibly a justification for the use of models which consider the DNA chain as uniform. It shows that the existence of solitons is a generic feature of the underlying nonlinear dynamics and is to a large extent independent of the detailed modeling of DNA. The model we consider supports solitonic solutions, qualitatively and quantitatively very similar to the Y solitons, in a fully realistic range of all the physical parameters characterizing the DNA.

Energy localization in the phi4 oscillator chain

We study energy localization in a finite one-dimensional phi(4) oscillator chain with initial energy in a single oscillator of the chain. We numerically calculate the effective number of degrees of freedom sharing the energy on the lattice as a function of time. We find that for energies smaller than a critical value, energy equipartition among the oscillators is reached in a relatively short time. On the other hand, above the critical energy, a decreasing number of particles sharing the energy is observed. We give an estimate of the effective number of degrees of freedom as a function of the energy. Our results suggest that localization is due to the appearance, above threshold, of a breather-like structure. Analytic arguments are given, based on the averaging theory and the analysis of a discrete nonlinear Schrödinger equation approximating the dynamics, to support and explain the numerical results.