Base pair opening and bubble transport in a DNA double helix induced by a protein molecule in a viscous medium

We study the nonlinear dynamics of a protein-DNA molecular system by treating DNA as a set of two coupled linear chains and protein in the form of a single linear chain sliding along the DNA at the physiological temperature in a viscous medium. The nonlinear dynamics of the above molecular system in general is governed by a perturbed nonlinear Schr"{o}dinger equation. In the non-viscous limit, the equation reduces to the completely integrable nonlinear Schr"{o}dinger (NLS) equation which admits N-soliton solutions. The soliton excitations of the DNA bases make localized base pair opening and travel along the DNA chain in the form of a bubble. This may represent the bubble generated during the transcription process when an RNA-polymerase binds to a promoter site in the DNA double helical chain. The perturbed NLS equation is solved using a perturbation theory by treating the viscous effect due to surrounding as a weak perturbation and the results show that the viscosity of the solvent in the surrounding damps out the amplitude of the soliton.

💡 Research Summary

The paper presents a theoretical study of the nonlinear dynamics of a protein–DNA system, modeling the double‑helical DNA as two coupled linear chains and the interacting protein as a single linear chain that slides along the DNA backbone. By constructing a Lagrangian that includes kinetic energy, elastic energy of each chain, and nonlinear inter‑chain potentials representing hydrogen bonding between base pairs and protein–DNA coupling, the authors derive the equations of motion. In the limit of a non‑viscous environment, the continuum approximation leads to the standard integrable nonlinear Schrödinger (NLS) equation:
i∂ψ/∂t + α∂²ψ/∂x² + β|ψ|²ψ = 0,
where ψ(x,t) is a complex field describing the collective displacement of the DNA bases, α is a dispersion coefficient, and β encodes the strength of the nonlinear interaction. Because the NLS equation is completely integrable, it admits N‑soliton solutions. The authors focus on the one‑soliton case, which corresponds to a localized excitation that opens a base‑pair “bubble” and propagates along the DNA strand. This soliton‑induced bubble is interpreted as the transient opening of the double helix that occurs when RNA polymerase binds to a promoter during transcription.

To incorporate the realistic effect of the surrounding viscous medium (e.g., cytoplasm or nucleoplasm), the authors add a weak dissipative term proportional to the solvent viscosity, resulting in a perturbed NLS equation:
i∂ψ/∂t + α∂²ψ/∂x² + β|ψ|²ψ = iεψ,
with ε ≪ 1 representing the viscous damping coefficient. Using perturbation theory, they treat ε as a small parameter and derive evolution equations for the soliton parameters (amplitude A(t), width, velocity, and phase). The analysis shows that the amplitude decays exponentially, A(t)=A₀ e⁻ᵋᵗ, while the soliton’s velocity remains essentially unchanged. Consequently, the bubble’s intensity diminishes over time, reflecting the physical damping of the excitation by the viscous environment.

Numerical simulations corroborate the analytical findings. In the non‑viscous case, the soliton propagates indefinitely without shape change, conserving its energy. When viscosity is introduced, the soliton’s peak amplitude drops, its width broadens slightly, and after a characteristic decay time the localized bubble disappears. The authors note that the decay time depends on ε, which in turn is linked to temperature, ionic strength, and the concentration of macromolecular crowding agents in the cell.

The discussion connects these results to biological transcription. In vivo, the transient opening of DNA required for RNA polymerase to access the template strand is not permanent; the surrounding medium’s viscosity imposes a finite lifetime on the opened region. The model therefore provides a physical mechanism by which cellular conditions can regulate transcriptional efficiency: higher viscosity (e.g., during stress or in crowded environments) would shorten bubble lifetimes, potentially reducing transcription rates.

In conclusion, the study demonstrates that a protein‑DNA complex can be effectively described by a perturbed NLS framework, where soliton solutions capture the essential features of base‑pair opening and bubble transport. The weak viscous perturbation accounts for realistic damping, predicting exponential attenuation of soliton amplitude. The authors suggest extensions such as incorporating three‑dimensional DNA geometry, multiple interacting proteins, and experimental validation using single‑molecule techniques. This work bridges nonlinear wave physics and molecular biology, offering a quantitative perspective on how physical forces influence genetic processes.