Theory of high-force DNA stretching and overstretching

Single-molecule experiments on single- and double-stranded DNA have sparked a renewed interest in the force versus extension of polymers. The extensible freely jointed chain (FJC) model is frequently invoked to explain the observed behavior of single-stranded DNA, but this model does not satisfactorily describe recent high-force stretching data. We instead propose a model (the discrete persistent chain) that borrows features from both the FJC and the wormlike chain, and show that it resembles the data more closely. We find that most of the high-force behavior previously attributed to stretch elasticity is really a feature of the corrected entropic elasticity; the true stretch compliance of single-stranded DNA is several times smaller than that found by previous authors. Next we elaborate our model to allow coexistence of two conformational states of DNA, each with its own stretch and bend elastic constants. Our model is computationally simple and gives an excellent fit through the entire overstretching transition of nicked, double-stranded DNA. The fit gives the first value for the bend stiffness of the overstretched state. In particular, we find the effective bend stiffness for DNA in this state to be about 12 nm k(B)T, a value quite different from either the B-form or single-stranded DNA.