Mesoscopic models for DNA stretching under force: new results and comparison to experiments

TL;DR

This paper develops a mesoscopic analytical model for DNA stretching under force, accurately fitting experimental data for ssDNA and B to S transitions, and highlights the sequence sensitivity of these structural changes.

Contribution

The authors introduce a coupled discrete worm-like chain-Ising model that captures DNA stretching behavior and provides analytical formulas for different transition regimes, connecting to previous thermal melting models.

Findings

ssDNA fits with three parameters over a nanoNewton range

monomer size of 0.2 nm is necessary for good fit

B to ssDNA transitions are not well fitted by mesoscopic models

Abstract

Single molecule experiments on B-DNA stretching have revealed one or two structural transitions, when increasing the external force. They are characterized by a sudden increase of DNA contour length and a decrease of the bending rigidity. It has been proposed that the first transition, at forces of 60--80 pN, is a transition from B to S-DNA, viewed as a stretched duplex DNA, while the second one, at stronger forces, is a strand peeling resulting in single stranded DNAs (ssDNA), similar to thermal denaturation. But due to experimental conditions these two transitions can overlap, for instance for poly(dA-dT). We derive analytical formula using a coupled discrete worm like chain-Ising model. Our model takes into account bending rigidity, discreteness of the chain, linear and non-linear (for ssDNA) bond stretching. In the limit of zero force, this model simplifies into a coupled model already developed by us for studying thermal DNA melting, establishing a connexion with previous fitting parameter values for denaturation profiles. We find that: (i) ssDNA is fitted, using an analytical formula, over a nanoNewton range with only three free parameters, the contour length, the bending modulus and the monomer size; (ii) a surprisingly good fit on this force range is possible only by choosing a monomer size of 0.2 nm, almost 4 times smaller than the ssDNA nucleobase length; (iii) mesoscopic models are not able to fit B to ssDNA (or S to ss) transitions; (iv) an analytical formula for fitting B to S transitions is derived in the strong force approximation and for long DNAs, which is in excellent agreement with exact transfer matrix calculations; (v) this formula fits perfectly well poly(dG-dC) and $\lambda$-DNA force-extension curves with consistent parameter values; (vi) a coherent picture, where S to ssDNA transitions are much more sensitive to base-pair sequence than the B to S one, emerges.