Experimental and theoretical studies of sequence effects on the fluctuation and melting of short DNA molecules

Understanding the melting of short DNA sequences probes DNA at the scale of the genetic code and raises questions which are very different from those posed by very long sequences, which have been extensively studied. We investigate this problem by combining experiments and theory. A new experimental method allows us to make a mapping of the opening of the guanines along the sequence as a function of temperature. The results indicate that non-local effects may be important in DNA because an AT-rich region is able to influence the opening of a base pair which is about 10 base pairs away. An earlier mesoscopic model of DNA is modified to correctly describe the time scales associated to the opening of individual base pairs well below melting, and to properly take into account the sequence. Using this model to analyze some characteristic sequences for which detailed experimental data on the melting is available [Montrichok et al. 2003 Europhys. Lett. {\bf 62} 452], we show that we have to introduce non-local effects of AT-rich regions to get acceptable results. This brings a second indication that the influence of these highly fluctuating regions of DNA on their neighborhood can extend to some distance.

💡 Research Summary

The paper tackles the melting behavior of short DNA fragments—typically 20 to 30 base pairs—in a combined experimental‑theoretical framework. While the thermodynamics of long polymers has been extensively characterized, short oligomers probe the scale of the genetic code, where individual base‑pair fluctuations become observable and biologically relevant. The authors introduce a novel chemical probing technique that exploits the oxidative sensitivity of guanine residues. When a guanine is in an “open” state (i.e., its hydrogen‑bonded pair is broken), it reacts with a specific reagent, leading to a detectable shift in electrophoretic mobility. By incrementally raising the temperature and measuring the fraction of modified guanines at each position, they obtain high‑resolution maps of base‑pair opening as a function of temperature.

These maps reveal two striking phenomena. First, AT‑rich stretches melt at lower temperatures and display large fluctuations, as expected from their weaker hydrogen‑bonding. Second, and more surprisingly, an AT‑rich region can increase the opening probability of a base pair located roughly ten nucleotides away—a clear indication of non‑local effects that are not captured by traditional, purely local melting models.

To interpret these observations, the authors revisit the Peyrard‑Bishop‑Dauxois (PBD) mesoscopic model, which treats each base pair as a particle subject to a nonlinear on‑site potential (hydrogen bonding) and a harmonic stacking interaction with its neighbors. The standard PBD framework, however, fails on two counts for short DNA: (i) it does not reproduce the experimentally observed dwell times of individual open base pairs (tens of microseconds) and (ii) it assumes that the thermodynamic parameters depend only on the local AT/GC composition, ignoring any influence that a highly fluctuating AT segment might exert on distant sites.

The authors therefore introduce two key modifications. The first is a kinetic layer based on a Markovian transition scheme that explicitly incorporates opening and reclosing rates, calibrated to match the measured lifetimes of open states. This adjustment aligns the model’s time scale with experimental data without altering the equilibrium properties. The second, more profound change is the addition of a non‑local coupling term. In this term, the energetic parameters (hydrogen‑bond depth and stacking stiffness) of an AT‑rich segment are gradually attenuated over a distance of about ten base pairs, thereby weakening neighboring GC pairs and enhancing their propensity to open. The strength and range of this coupling are determined by fitting the model to the experimental opening profiles.

When applied to several benchmark sequences for which detailed melting curves are available (e.g., the 23‑bp and 30‑bp oligomers studied by Montrichok et al., 2003), the revised model reproduces both the overall melting temperature shifts and the broadened transition regions observed experimentally. Crucially, without the non‑local term the model underestimates the opening probability of GC pairs situated beyond the immediate AT stretch, whereas inclusion of the coupling yields quantitative agreement with the measured ~10 bp influence radius.

The study concludes that the melting of short DNA is governed not only by local base‑pair energetics but also by long‑range mechanical and thermodynamic perturbations emanating from highly fluctuating AT domains. This insight has implications for biological processes such as transcription initiation, where AT‑rich promoter regions may facilitate the opening of downstream sequences, and for the design of DNA‑based nanodevices that rely on precise control of local stability. Moreover, the authors’ hybrid experimental‑theoretical approach provides a robust platform for predicting melting behavior of arbitrary short sequences, opening avenues for rational design in synthetic biology and nanotechnology.