Forming superhelix of double stranded DNA from local deformation

Heeyuen Koh; Jae Gyung Lee; Jae Young Lee

arxiv: 2307.04597 · v19 · pith:HZK4MJOXnew · submitted 2023-07-10 · 🌊 nlin.AO

Forming superhelix of double stranded DNA from local deformation

Heeyuen Koh , Jae Young Lee , Jae Gyung Lee This is my paper

Pith reviewed 2026-05-24 07:24 UTC · model grok-4.3

classification 🌊 nlin.AO

keywords DNA superhelixgeometrical constraintskurtosisbase-pair resolutioncurvature formationDNA packagingmolecular dynamics

0 comments

The pith

Geometrical constraints on base pairs in curved DNA determine superhelix height via kurtosis rather than topology.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper derives geometrical constraints on base-pair resolution in a curved DNA strand separately from its elasticity using differential geometry. These constraints arise from the inherent helicity of base pairs and create a conditional affinity for curvature that fixes the bend-twist ratio. The central result is that kurtosis computed from an open-ended strand under these constraints sets the height of the superhelix formed around a core. This separation of geometry from elasticity matters for understanding DNA packaging because the process requires precise 1.7 turns of superhelix. Coarse-grained molecular dynamics simulations confirm the curvature formation process follows the derived constraints.

Core claim

The geometrical constraints that base pairs impose due to their inherent helicity characterize the conditional affinity for curvature formation, thereby specifying the bend-twist ratio. The result of the derivation indicates that the kurtosis derived from the geometrical constraints of an open-ended strand can be the factor that decides the height of the superhelix, rather than the topological properties of a circular chain or the binding conformation with the obstacles.

What carries the argument

geometrical constraints on base-pair-wise resolution in a curved DNA strand, derived independently via differential geometry, that impose conditional affinity for curvature and fix the bend-twist ratio

Load-bearing premise

That geometric constraints on base-pair resolution in a curved strand can be derived independently of the strand's elasticity and that an open-ended strand with a simplified core structure is representative of the packaging process.

What would settle it

A simulation or measurement in which kurtosis is held fixed while superhelix height changes with altered topology or binding conformation would show that kurtosis is not the deciding factor.

Figures

Figures reproduced from arXiv: 2307.04597 by Heeyuen Koh, Jae Gyung Lee, Jae Young Lee.

**Figure 2.** Figure 2: A-a. Arrangement between n − 1 th(blue) and n th(red) base pairs. θω is the angle for inherent helicity, which is from 32.4 ◦ . Two solid dots represent two nucleotides for each base pair, A-b. Modified twist angle, θ ′ ω with ∆Ly, B-a. Additional twist angle ∆θω for θ ′ ω = θω + ∆θω with ∆Ly and ∆Lx. Note that ∆θω is measured from the median of n th base pair. The contact point between R⃗c and the circumf… view at source ↗

**Figure 3.** Figure 3: A. Schematic figure of roll(ρ) and tilt(τ) combination along the strand for its helicity on the curved surface. Dashed line is the arclength of the strand that is marked with ∆⃗L, B. Deformation of each rotational variable in curvature unit along the contact angle. Maximum value of twist(ω) in the figure is around 10◦ . For an effective comparison, two sets of strand information are adapted from Freeman et… view at source ↗

**Figure 4.** Figure 4: The angle between the contact point and the center of the base pair is measured as shown by spatiotemporal [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗

**Figure 5.** Figure 5: Total energy during wrapping process in 30 ns using oxDNA2 and new thermostat. A. c1 (pink) proves its [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗

**Figure 6.** Figure 6: A. Wrapping using the CG model with major-minor groove(oxDNA2), B. Wrapping from CG model without [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗

read the original abstract

In contrast to the sequence dependent elasticity of the DNA strand, which is revealed as nonlocal and nonlinear, the geometric constraints derived by differential geometry have not been fully elaborated in DNA strand dynamics, even though these constraints contribute independently to the quantification of related energetics. In this paper, the geometrical constraints on the base pair wise resolution in a curved DNA strand are derived separately from its elasticity, addressing the deformation characteristics during superhelix formation around a simplified core structure, which is the quintessential step in DNA packaging. The constraints derived from the given helicity of DNA strand characterize the conditional affinity for curvature formation, thereby specifying the bend-twist ratio required for superhelix formation. The result includes the conditional kurtosis, which is the deformation perpendicular to the plane defined by the curved strand, determining the height of the superhelix. Coarse-grained molecular dynamics simulation validates the description of the curvature formation process and its sequence dependent affinity.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper derives base-pair geometric constraints via differential geometry and claims their kurtosis sets superhelix height independently of elasticity or topology, but the abstract supplies no equations to check whether that separation holds.

read the letter

The central move is to treat the helicity of base pairs as imposing geometric constraints on curvature in an open-ended strand, derive a bend-twist ratio from that, and then argue that the resulting kurtosis, not linking number or binding geometry, fixes the height of the superhelix. Coarse-grained MD is used to check the curvature process. That separation of geometry from elasticity is the part worth looking at, because most DNA models fold the two together through sequence-dependent parameters. If the derivation really produces the constraints without feeding in elastic moduli, it could simplify some packaging calculations. The MD step at least gives a concrete check on whether the predicted curvature forms under the stated rules. The main weakness is that nothing in the abstract shows the independence the claim requires. Height in any real model still comes from minimizing total energy, which includes bending and twist terms; without an explicit test that varying those moduli while holding the geometric kurtosis fixed leaves height unchanged, the kurtosis cannot be shown to decide the outcome on its own. The open-ended strand also drops the topological constraints that exist in closed DNA, so the comparison to circular-chain topology is built on a setup that excludes the competing mechanism. The abstract gives no error bars, no parameter list, and no citation trail for the differential-geometry steps, which makes it impossible to judge whether the kurtosis result is new or reduces to earlier geometric models of DNA. Readers who work on coarse-grained DNA packaging or who want to test geometry-only contributions might still want the full derivation and the simulation details. The work is specific enough and the MD is present, so it clears the bar for a referee to examine the equations and the independence test. I would send it to review rather than desk-reject.

Referee Report

2 major / 1 minor

Summary. The manuscript derives geometrical constraints on base-pair resolution in a curved DNA strand independently from elasticity. It claims that the kurtosis arising from these constraints on an open-ended strand determines the height of the superhelix formed around a simplified core, rather than topological properties of circular chains or binding conformations with obstacles. Coarse-grained molecular dynamics simulations are used to validate the curvature formation process.

Significance. If the claimed independence of geometric kurtosis from elastic moduli and topology can be demonstrated explicitly, the result would provide a local geometric mechanism for selecting the characteristic 1.7-turn superhelix height in DNA packaging. This could complement existing elasticity- and topology-based models and inform understanding of nucleosome and viral packaging processes. The MD validation is a positive element, but its strength depends on the details of the geometric derivation.

major comments (2)

[Abstract / derivation] Abstract and derivation section: The central claim requires that kurtosis from base-pair geometric constraints sets superhelix height independently of elastic energy minimization. No explicit demonstration is described showing that varying bend and twist moduli (while holding the derived geometric kurtosis fixed) leaves height unchanged; without this, the separation of geometry from elasticity cannot be verified as load-bearing for the result.
[Model setup] Model setup (open-ended strand): The assertion that kurtosis replaces topological properties rests on an open-ended strand that by construction removes linking-number constraints present in real packaging. A direct comparison or control simulation with topologically closed chains is needed to establish that the height selection is not an artifact of the simplification.

minor comments (1)

[Abstract] Abstract: The phrase 'nonlocal and nonlinear' for sequence-dependent elasticity is stated without a supporting reference or brief elaboration on how the geometric constraints interact with it.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment below, clarifying the analytical basis for the claimed independence and the rationale for the model setup.

read point-by-point responses

Referee: [Abstract / derivation] Abstract and derivation section: The central claim requires that kurtosis from base-pair geometric constraints sets superhelix height independently of elastic energy minimization. No explicit demonstration is described showing that varying bend and twist moduli (while holding the derived geometric kurtosis fixed) leaves height unchanged; without this, the separation of geometry from elasticity cannot be verified as load-bearing for the result.

Authors: The derivation applies differential geometry directly to the base-pair resolution constraints imposed by DNA helicity in a curved strand. This yields the kurtosis and the associated bend-twist ratio as purely geometric quantities, without any reference to elastic moduli or energy minimization. The separation from elasticity is therefore established by construction in the analytic derivation itself; the geometric constraints operate independently of the elastic parameters. The coarse-grained MD simulations confirm that the resulting curvature formation produces the reported superhelix height under these constraints, but the independence claim rests on the geometry, not on numerical variation of moduli. revision: no
Referee: [Model setup] Model setup (open-ended strand): The assertion that kurtosis replaces topological properties rests on an open-ended strand that by construction removes linking-number constraints present in real packaging. A direct comparison or control simulation with topologically closed chains is needed to establish that the height selection is not an artifact of the simplification.

Authors: The open-ended strand is chosen precisely to isolate the local geometric mechanism arising from base-pair helicity. By removing global topological constraints such as linking number, the model demonstrates that the kurtosis derived from these local constraints is sufficient to select the superhelix height. This establishes a geometric contribution that can operate independently of topology. While simulations of closed chains would be informative for combined effects in real packaging, they are not required to substantiate the local geometric selection mechanism presented here. revision: no

Circularity Check

0 steps flagged

No circularity: geometric constraints derived independently of elasticity and topology

full rationale

The paper derives base-pair geometric constraints via differential geometry on an open-ended strand, separately from elasticity, then reports that the resulting kurtosis sets superhelix height. No equations or self-citations are shown that define the target height or bend-twist ratio into the geometric inputs by construction, nor that rename a fitted result as a prediction. The central claim rests on an explicit separation of geometry from elastic energy minimization plus external validation by coarse-grained MD simulation, satisfying the criteria for a self-contained derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

Paper relies on standard differential geometry for helicity constraints and assumes separability of geometry from elasticity; no explicit free parameters or invented entities stated in abstract.

axioms (2)

domain assumption Differential geometry applies directly to base-pair resolution in curved DNA strands
Invoked to derive constraints separately from elasticity
domain assumption Open-ended strand with simplified core is representative
Used to isolate kurtosis effect from topology or binding

pith-pipeline@v0.9.0 · 5726 in / 1269 out tokens · 30820 ms · 2026-05-24T07:24:08.181161+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

the kurtosis derived from the geometrical constraints of an open-ended strand can be the factor that decides the height of the superhelix... K = rτℓ(ω0 + Δω) sin ϕ ± 2ℓ²Ω sin 2ϕ
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

major-minor groove... θ=120° vs 180°... kurtosis accumulation... 10.3 bp period

What do these tags mean?

matches: The paper's claim is directly supported by a theorem in the formal canon.
supports: The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends: The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses: The paper appears to rely on the theorem as machinery.
contradicts: The paper's claim conflicts with a theorem or certificate in the canon.
unclear: Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.