Discrete-breather-assisted charge transport along DNA-like molecular wires

Bubble Relaxation Dynamics in Homopolymer DNA Sequences

Understanding the inherent timescales of large bubbles in DNA is critical to a thorough comprehension of its physicochemical characteristics, as well as their potential role on helix opening and biological function. In this work, we employ the coarse-grained Peyrard-Bishop-Dauxois model of DNA to study relaxation dynamics of large bubbles in homopolymer DNA, using simulations up to the microsecond time scale. By studying energy autocorrelation functions of relatively large bubbles inserted into thermalised DNA molecules, we extract characteristic relaxation times from the equilibration process for both adenine-thymine (AT) and guanine-cytosine (GC) homopolymers. Bubbles of different amplitudes and widths are investigated through extensive statistics and appropriate fittings of their relaxation. Characteristic relaxation times increase with bubble amplitude and width. We show that, within the model, relaxation times are two orders of magnitude longer in GC sequences than in AT sequences. Overall, our results confirm that large bubbles leave a lasting impact on the molecule's dynamics, for times between 0.5-500 ns depending on the homopolymer type and bubble shape, thus clearly affecting long-time evolutions of the molecule.

Discrete breathers in a triangular β-Fermi-Pasta-Ulam-Tsingou lattice

A practical approach to the search for (quasi-) discrete breathers (DBs) in a triangular β-FPUT lattice (after Fermi, Pasta, Ulam, and Tsingou) is proposed. DBs are obtained by superimposing localizing functions on delocalized nonlinear vibrational modes (DNVMs) having frequencies above the phonon spectrum of the lattice. Zero-dimensional and one-dimensional DBs are obtained. The former ones are localized in both spatial dimensions, and the latter ones are only in one dimension. Among the one-dimensional DBs, two families are considered: the first is based on the DNVMs of a triangular lattice, and the second is based on the DNVMs of a chain. We speculate that our systematic approach on the triangular β-FPUT lattice reveals all possible types of spatially localized oscillations with frequencies bifurcating from the upper edge of the phonon band as all DNVMs with frequencies above the phonon band are analyzed.

Distributions of bubble lifetimes and bubble lengths in DNA

We investigate the distribution of bubble lifetimes and bubble lengths in DNA at physiological temperature, by performing extensive molecular dynamics simulations with the Peyrard-Bishop-Dauxois (PBD) model, as well as an extended version (ePBD) having a sequence-dependent stacking interaction, emphasizing the effect of the sequences' guanine-cytosine (GC)/adenine-thymine (AT) content on these distributions. For both models we find that base pair-dependent (GC vs AT) thresholds for considering complementary nucleotides to be separated are able to reproduce the observed dependence of the melting temperature on the GC content of the DNA sequence. Using these thresholds for base pair openings, we obtain bubble lifetime distributions for bubbles of lengths up to ten base pairs as the GC content of the sequences is varied, which are accurately fitted with stretched exponential functions. We find that for both models the average bubble lifetime decreases with increasing either the bubble length or the GC content. In addition, the obtained bubble length distributions are also fitted by appropriate stretched exponential functions and our results show that short bubbles have similar likelihoods for any GC content, but longer ones are substantially more likely to occur in AT-rich sequences. We also show that the ePBD model permits more, longer-lived, bubbles than the PBD system.

A review on nonlinear DNA physics

The study and the investigation of structural and dynamical properties of complex systems have attracted considerable interest among scientists in general and physicists and biologists in particular. The present review paper represents a broad overview of the research performed over the nonlinear dynamics of DNA, devoted to some different aspects of DNA physics and including analytical, quantum and computational tools to understand nonlinear DNA physics. We review in detail the semi-discrete approximation within helicoidal Peyrard-Bishop model and show that localized modulated solitary waves, usually called breathers, can emerge and move along the DNA. Since living processes occur at submolecular level, we then discuss a quantum treatment to address the problem of how charge and energy are transported on DNA and how they may play an important role for the functioning of living cells. While this problem has attracted the attention of researchers for a long time, it is still poorly understood how charge and energy transport can occur at distances comparable to the size of macromolecules. Here, we review a theory based on the mechanism of 'self-trapping' of electrons due to their interaction with mechanical (thermal) oscillation of the DNA structure. We also describe recent computational models that have been developed to capture nonlinear mechanics of DNA in vitro and in vivo, possibly under topological constraints. Finally, we provide some conjectures on potential future directions for this field.