Soliton Models in DNA Dynamics: An Introduction

In 1834, the Scottish engineer John Scott Russell observed a peculiar wave in a canal — one that maintained its shape and speed over a remarkable distance without dispersing. He called it a “wave of translation.” Today we call it a soliton, and it turns out this 19th-century curiosity has deep consequences for understanding the very molecule that encodes life.

The Sine-Gordon Equation

At the heart of soliton physics lies the Sine-Gordon equation:

$\frac{\partial^2 \phi}{\partial t^2} - \frac{\partial^2 \phi}{\partial x^2} + \sin(\phi) = 0$

This nonlinear PDE admits exact soliton solutions of the form:

$\phi(x,t) = 4 \arctan\!\left(\exp\!\left(\gamma(x - vt)\right)\right)$

where $\gamma = 1/\sqrt{1-v^2}$ is the relativistic Lorentz factor and $v$ is the soliton velocity. The derivative $\partial \phi / \partial x$ traces out a bell-shaped pulse — localized, stable, and traveling without dispersion.

DNA as a Coupled Pendulum Chain

The connection to DNA comes through a model introduced by Englander et al. (1980) and later refined by Peyrard and Bishop. The idea: treat each base pair in a DNA double helix as a pendulum coupled to its neighbors along the strand. When one base pair opens (rotates away from its partner), the disturbance can propagate along the chain as a soliton.

Specifically, if we let $\phi_n$ denote the angle of the $n$-th base pair from equilibrium, the equation of motion in the continuum limit becomes the Sine-Gordon equation. A soliton in this context represents a localized region of base-pair opening — sometimes called a “bubble” — that travels along the helix without spreading.

Why This Matters

This isn’t merely an elegant analogy. The opening of base pairs is a prerequisite for DNA transcription — the process by which genetic information is read and proteins are synthesized. If soliton dynamics govern how these bubbles form and propagate, then understanding solitons is, in a very real sense, understanding a mechanism of life itself.

Current research (including work ongoing at Howard’s Department of Mathematics under Prof. Roberto De Leo) investigates:

• Numerical stability of soliton propagation under realistic coupling constants
• Interaction between multiple solitons on the same strand
• Perturbative effects of thermal fluctuations on soliton trajectories

A Note on Computation

Simulating these dynamics numerically requires care. Naive finite-difference schemes can introduce artificial damping or amplification. We use symplectic integrators that preserve the Hamiltonian structure of the equations, ensuring energy is conserved over long simulation windows. Python implementations using numpy and scipy — with custom integrators — form the backbone of this computational work.

What’s Next

In a future post, I’ll walk through a concrete Python implementation of the Sine-Gordon soliton and show how varying the velocity $v$ changes the shape and energy of the pulse. The mathematics is clean, the physics is profound, and the code is surprisingly simple.

---

Umanga Baral is a Mathematics and Physics undergraduate at Howard University, working as a Research Assistant in the Department of Mathematics.