Persistence length of short homopolymeric single-stranded DNA sequences in polyvalent cations

We used simulations of short single stranded DNA (ssDNA) homopolymers, based on the sequence dependent Three Interaction Site (TIS) model, to calculate the persistence length ( l _p ) in polyvalent cations. The TIS model accounts for stacking interactions and electrostatic interactions are treated using the Coulomb potential. We find that l _p for dT ₃₀ (T is thymine) and dA ₃₀ (A is adenine) is quantitatively fit using ( is the bare persistence length, λ is a dimensionless constant, and κ ⁻¹ is the inverse Debye length) in the divalent cations Mg ²⁺ and Ca ²⁺ . The dependence of l _p on κ is surprising because it was derived for long flexible polyelectrolytes in which the charges interact via the Debye-Hückel potential. The values are 0.4 nm and 1.1 nm for polyT and polyA, respectively. Strikingly, l _p is almost independent of the tetravalent spermine concentration. There is no clear theoretical explanation although simulations suggest that the number of spermine molecules that bind to the ssDNA saturates at a small value. A qualitative picture, based on the restrictions of access to the phosphate groups due to volume exclusion of the anisotropic structure of Spm ⁴⁺ , rationalizes the simulation results. The predicted dependence of l _p in spermine awaits experimental test.