Probing the statistics of sequence-dependent DNA conformations in solution using SAXS

SAXS studies of four 60 base-pair DNA duplexes with sequences closely related to part of the GAGE6 (G-antigen 6) promoter have been performed to study the role of DNA conformations in solution and their potential relationship to DNA-protein binding. We show that the SAXS data can be analysed using a simple polymer model which nevertheless quantitatively describes the average persistence length and torsional rigidity of the DNA double helix to determine the statistical distribution of local conformations of the DNA in solution to a high accuracy. Although the SAXS data is averaged over time and all spatial orientations of the molecules, for sequences which have some asymmetry in the data we show that the conformations can be oriented with respect to the sequence. This allows specific features detected by the analysis to be precisely related to the DNA sequence, opening up new opportunities for SAXS to investigate the properties of DNA in solution. The biological implications of these results are discussed.

💡 Research Summary

This study investigates how DNA sequence influences the three‑dimensional conformation of short duplexes in solution and how such conformational variability might affect protein binding. Four 60‑base‑pair double‑stranded DNA fragments derived from the GAGE6 promoter were synthesized: the native sequence and three variants in which AT‑rich tracts were progressively replaced by GC‑rich tracts (GAGE6‑1, GAGE6‑2, GAGE6‑3). All samples were measured by size‑exclusion‑chromatography coupled to small‑angle X‑ray scattering (SEC‑SAXS) at physiological salt conditions (150 mM KCl, 5 mM MgCl₂, 20 mM HEPES, pH 7.4).

The authors argue that analysis in real space, i.e., the pair‑distance distribution function P(r), is more sensitive to local bends and twists than the traditional reciprocal‑space intensity I(q). By performing indirect Fourier transforms with the AutoGNOM routine, they obtained robust P(r) curves for each construct, carefully checking regularisation effects and confirming that the main features are stable across different parameter sets.

A simple polymer model was then employed to extract two global mechanical parameters: the persistence length (≈ 50 nm) and the torsional rigidity (≈ 2.5 × 10⁻¹⁹ J·rad⁻²). Monte‑Carlo sampling generated thousands of possible 3‑D conformations consistent with these parameters; each conformation’s P(r) was calculated and compared to the experimental P(r) to identify the most probable ensemble. This approach yields a quantitative distribution of local bending angles and twist deviations along the 60‑bp chain. Notably, the native AT‑rich regions display an increased probability of bends of 3–5° centered around base pairs 30–35, whereas the GC‑substituted variants show reduced bending and a modest increase in overall stiffness.

A particularly innovative aspect is the detection of asymmetry in the P(r) curves for sequences with non‑uniform AT‑tract placement. Because SAXS averages over all orientations, such asymmetry implies that the ensemble retains a preferred orientation relative to the sequence. The authors exploit this to map the direction of the dominant bend to a specific segment of the DNA, effectively “orienting” the molecule in silico despite the rotational averaging inherent to SAXS.

Biologically, the GAGE6 fragment is known to bind the multifunctional nuclear protein SFPQ (also called PSF). SFPQ contains an arginine‑glycine‑rich (RGG) domain that preferentially interacts with exposed bases. The authors propose that the AT‑rich tracts, by locally softening the helix and promoting transient openings or kinks, create a more accessible binding platform for SFPQ. Their SAXS‑derived conformational maps support this hypothesis, linking sequence‑dependent flexibility to protein‑DNA recognition.

In summary, the paper demonstrates that (1) SEC‑SAXS combined with a minimal polymer model can accurately quantify both persistence length and torsional rigidity of short DNA duplexes; (2) subtle, sequence‑dependent variations in local curvature and twist can be resolved at the level of a few base pairs; (3) asymmetries in the real‑space P(r) allow orientation of the DNA relative to its sequence; and (4) these structural nuances provide a plausible mechanistic basis for the selective binding of SFPQ to the GAGE6 promoter. The methodology is relatively straightforward, yet powerful enough to be applied to a broad range of DNA‑protein interaction studies where solution‑phase dynamics are critical.