Propagation of twist solitons in real DNA chains

We report on numerical investigations concerning the propagation of solitons in a real DNA chain (the Human Adenovirus 2) using a realistic model of DNA torsional dynamics; this takes fully into account the inhomogeneities in the real chain. We find that twist solitons propagate for considerable distances (2-10 times their diameters) before stopping due to phonon emission. Our results show that twist solitons may exist in real DNA chains; and on a more general level that solitonic propagation can take place in highly inhomogeneous media.

💡 Research Summary

The paper investigates the propagation of twist solitons—non‑linear torsional excitations—in a realistic model of DNA that fully incorporates the sequence‑dependent inhomogeneities of an actual genome, specifically the Human Adenovirus 2 (HAdV‑2) DNA. The authors begin by constructing an “inhomogeneous chain” model in which each base pair is assigned distinct mechanical parameters (torsional stiffness, mass, coupling constants, and electrostatic charge) derived from its AT or GC composition. This contrasts with earlier theoretical work that typically assumes a homogeneous lattice or uses averaged parameters, thereby neglecting the complex variability present in real DNA.

Using the sine‑Gordon (SG) type potential to describe the torsional interaction between neighboring base pairs, the authors embed a classic SG soliton as the initial condition. The soliton’s width is set to roughly five to ten base pairs (≈1 nm), and its initial velocity is varied across a realistic range (0.5–2 nm ps⁻¹). The equations of motion for the discrete chain are integrated with a fourth‑order Runge‑Kutta scheme, employing a small time step (Δt≈0.01 ps) and absorbing boundary layers to prevent artificial reflections. Energy conservation is monitored throughout, and Fourier analysis is applied to detect phonon (linear wave) emission.

The simulations reveal two principal findings. First, despite the pronounced heterogeneity of the chain, the soliton can travel distances equal to 2–10 times its own diameter (≈10–50 base pairs) before its amplitude decays appreciably. This indicates that the non‑linear restoring force of the soliton dominates over the dispersive effects introduced by the variable parameters, allowing relatively long‑range transport of torsional stress. Second, the ultimate stopping distance is governed by the gradual transfer of soliton energy to phonon modes. As the soliton moves, it excites small‑amplitude vibrations in the surrounding lattice; these phonons accumulate and eventually erode the soliton’s coherent structure, leading to its dissipation. The rate of this energy leakage is sequence‑dependent: GC‑rich regions, which possess higher torsional stiffness, cause faster attenuation, whereas AT‑rich stretches permit longer propagation.

A systematic parameter sweep identifies a “stability window” for soliton propagation. When the ratio of torsional stiffness (K) to the coupling constant (γ) lies between approximately 1.2 and 1.8, solitons maintain their shape and travel at least five soliton widths before significant decay. Outside this window, the initial soliton either destabilizes immediately or its propagation distance collapses to less than one or two widths. The authors also quantify the phonon spectrum at the point where the soliton halts, observing a pronounced low‑frequency peak that confirms the energy transfer mechanism.

These results have several important implications. They demonstrate that non‑linear torsional excitations can exist and move over biologically relevant distances in real, highly heterogeneous DNA, challenging the conventional view that DNA torsional stress dissipates solely by linear diffusion. This suggests a possible mechanism for the rapid transmission of supercoiling or torque generated during transcription, replication, or chromatin remodeling, where a soliton could convey torsional information across several nucleotides without immediate loss. Moreover, the sensitivity of soliton dynamics to local base composition implies that specific genomic regions (e.g., promoters, replication origins) might be tuned to either facilitate or impede soliton propagation, potentially influencing the regulation of genetic processes.

The paper concludes by proposing experimental avenues to validate the computational predictions. High‑speed atomic force microscopy (AFM) or optical tweezers could be employed to impose controlled torsional twists on single DNA molecules and monitor the ensuing dynamics. By measuring the velocity, attenuation length, and associated phonon emission, researchers could directly observe soliton‑like behavior. Coupling such measurements with fluorescence‑based torque sensors would further enable real‑time mapping of torsional stress propagation along the DNA contour.

In summary, the study provides a rigorous numerical demonstration that twist solitons can travel appreciable distances in a realistic DNA chain, despite the presence of strong sequence‑dependent inhomogeneities. It bridges non‑linear physics and molecular biology, offering a fresh perspective on how mechanical signals might be transmitted within the genome and opening new experimental directions to explore soliton‑mediated DNA dynamics.