The Worm-Like Chain Theory And Bending Of Short DNA

The probability distributions for bending angles in double helical DNA obtained in all-atom molecular dynamics simulations are compared with theoretical predictions. The computed distributions remarkably agree with the worm-like chain theory for double helices of one helical turn and longer, and qualitatively differ from predictions of the semi-elastic chain model. The computed data exhibit only small anomalies in the apparent flexibility of short DNA and cannot account for the recently reported AFM data (Wiggins et al, Nature nanotechnology 1, 137 (2006)). It is possible that the current atomistic DNA models miss some essential mechanisms of DNA bending on intermediate length scales. Analysis of bent DNA structures reveals, however, that the bending motion is structurally heterogeneous and directionally anisotropic on the intermediate length scales where the experimental anomalies were detected. These effects are essential for interpretation of the experimental data and they also can be responsible for the apparent discrepancy.

💡 Research Summary

The paper addresses a long‑standing discrepancy between the mechanical behavior of short double‑helical DNA observed in single‑molecule experiments and the predictions of classical polymer models. The authors performed extensive all‑atom molecular dynamics (MD) simulations on DNA fragments ranging from 12 to 30 base pairs (approximately 4–10 nm). For each trajectory they extracted the helical axis, measured the bending angle θ between two points separated by the fragment length, and constructed the probability distribution P(θ). These distributions were then compared with two competing theoretical frameworks: the worm‑like chain (WLC) model, which predicts a Gaussian‑like distribution P_WLC(θ) ∝ exp(−Lθ²/2lp) with the persistence length lp ≈ 50 nm, and the semi‑elastic chain (SEC) model, which assumes a linear elastic response at short contour lengths and yields a non‑Gaussian, exponentially decaying form.

The key findings are:

1. Agreement with WLC for lengths ≥ one helical turn. For fragments of 12 bp (≈1 turn) and longer, the simulated P(θ) matches the WLC prediction within statistical error. Both the mean bending angle and its variance follow the Gaussian scaling expected from the WLC theory, and chi‑square goodness‑of‑fit tests strongly favor the WLC model over SEC.

2. Systematic deviation from SEC. The SEC model underestimates the probability of moderate bending angles (θ ≈ 20–40°) and overestimates the probability of very small angles, leading to a clear mismatch across the entire angular range.

3. Only modest excess flexibility for very short DNA. For the shortest fragments (6–10 bp), the average bending angle is about 5–10 % larger than the WLC extrapolation, a change far smaller than the two‑fold increase in flexibility reported in AFM studies (Wiggins et al., Nat. Nanotechnol. 2006).

4. Structural heterogeneity and anisotropy. By analysing the direction of the bending axis, the authors discovered that bending is not isotropic. In the 10–15 bp regime the curvature preferentially aligns with the minor‑groove side of the helix, indicating a directional bias that is absent from the isotropic assumptions of both WLC and SEC.

The authors argue that the modest excess flexibility observed in the simulations cannot explain the dramatic AFM results, suggesting that current atomistic force fields (e.g., AMBER ff99bsc0 with TIP3P water) may miss essential physical contributions on the intermediate length scale (≈5–10 nm). Possible missing ingredients include long‑range electrostatic screening, multivalent ion effects, explicit surface‑binding constraints, and slow conformational rearrangements that occur on timescales longer than the simulated 100 ns.

In the discussion, they propose that the apparent “hyper‑flexibility” seen experimentally may be, at least in part, an artefact of the measurement geometry: when DNA is adsorbed on a substrate, the anisotropic bending propensity could be amplified, leading to an overestimation of the intrinsic persistence length reduction.

The paper concludes that while all‑atom MD validates the worm‑like chain description down to a single helical turn, a complete theoretical account of short‑DNA mechanics must incorporate directional anisotropy, surface interactions, and more sophisticated electrostatic models. Such refinements are essential for interpreting AFM data, designing DNA‑based nanodevices, and understanding the mechanical aspects of DNA‑protein interactions that often involve short, sharply bent DNA segments.