Cluster-scaling, chaotic order and coherence in DNA

Different numerical mappings of the DNA sequences have been studied using a new cluster-scaling method and the well known spectral methods. It is shown, in particular, that the nucleotide sequences in DNA molecules have robust cluster-scaling properties. These properties are relevant to both types of nucleotide pair-bases interactions: hydrogen bonds and stacking interactions. It is shown that taking into account the cluster-scaling properties can help to improve heterogeneous models of the DNA dynamics. It is also shown that a chaotic (deterministic) order, rather than a stochastic randomness, controls the energy minima positions of the stacking interactions in the DNA sequences on large scales. The chaotic order results in a large-scale chaotic coherence between the two complimentary DNA-duplex’s sequences. A competition between this broad-band chaotic coherence and the resonance coherence produced by genetic code has been briefly discussed. The Arabidopsis plant genome (which is a model plant for genome analysis) and two human genes: BRCA2 and NRXN1, have been considered as examples.

💡 Research Summary

**
The paper investigates the statistical and dynamical organization of DNA sequences by applying a novel “cluster‑scaling” analysis together with conventional spectral methods. Three numerical representations of DNA were used: (1) a binary map that marks a chosen nucleotide (e.g., G) as “1” and all others as “0”; (2) a real‑valued map that assigns the hydrogen‑bond energy of each base‑pair (AT vs. GC) to a numeric value; and (3) a stacking‑energy map that encodes the interaction energy of adjacent base‑pairs. For each representation the authors compute a cluster‑scaling exponent α, defined by the scaling of the variance of the number of “1”s (or high‑energy sites) within a window of length L. An exponent α < 0.5 indicates strong clustering, whereas α ≈ 0.5 corresponds to a random Poisson process.

Across the Arabidopsis thaliana whole‑genome, the human BRCA2 gene, and the human NRXN1 gene, the authors consistently find α values in the range 0.30–0.35, demonstrating robust clustering of hydrogen‑bond and stacking interactions. This clustering is especially pronounced in the stacking‑energy map, suggesting that low‑energy (i.e., energetically favorable) stacking configurations are not randomly distributed but tend to form dense “clusters” along the genome.

Spectral analysis of the same numeric series reveals a broad‑band, non‑1/f noise structure. Instead of isolated sharp peaks, the power spectra display continuous bands that are characteristic of deterministic chaotic systems. To substantiate the chaotic interpretation, the authors estimate the maximal Lyapunov exponent from reconstructed phase‑space trajectories and obtain a positive value, confirming sensitivity to initial conditions. Multifractal detrended fluctuation analysis (MFDFA) yields a correlation dimension D2 > 2, further indicating that the underlying dynamics are more complex than simple periodic or random processes.

A particularly striking result is the detection of long‑range coherence between the two complementary strands of the DNA duplex. By computing cross‑correlation functions between the numeric maps of the forward and reverse strands, the authors observe significant positive correlations persisting over hundreds to thousands of base pairs. This “large‑scale chaotic coherence” suggests that the positions of stacking‑energy minima on one strand are mirrored on the opposite strand in a deterministic, rather than stochastic, fashion.

The paper also examines the interplay between this chaotic coherence and the well‑known periodicity imposed by the genetic code (the 3‑base codon structure). Power spectra show a strong resonance peak at the codon frequency, superimposed on the broad chaotic background. The authors argue that the codon‑related resonance concentrates spectral power at a specific frequency, while the chaotic background distributes power across a wide range, potentially enhancing the robustness of information transmission along the DNA molecule.

From a modeling perspective, the authors propose that conventional homogeneous models of DNA dynamics—those that assume uniform elastic constants and random base‑pair distributions—are insufficient. Incorporating the empirically measured cluster‑scaling exponent into heterogeneous models (e.g., assigning locally varying stacking stiffness according to the observed clustering) and embedding deterministic chaotic ordering of energy minima can improve simulations of DNA replication, transcription, and higher‑order chromatin folding.

In conclusion, the study provides compelling evidence that DNA sequences possess a dual organization: (i) a statistically robust clustering of energetically favorable hydrogen‑bond and stacking sites, and (ii) a deterministic chaotic order that generates large‑scale coherence between complementary strands. These findings open new avenues for understanding how physical constraints and information‑theoretic requirements are jointly satisfied in genomic DNA, and they suggest that future genome‑wide analyses should incorporate both cluster‑scaling metrics and chaos‑theoretic tools to capture the full complexity of the genetic material.