Worm like chain

Theoretical Considerations

The WLC model envisions an isotropic rod that is continuously flexible. This is in contrast to the freely-jointed chain model that is flexible only between discrete segments. The worm-like chain model is particularly suited for describing stiffer polymers, with successive segments displaying a sort of cooperativity: all pointing in roughly the same direction. At room temperature, the polymer adopts a conformational ensemble that is smoothly curved; at T = 0 K, the polymer adopts a rigid rod conformation.

For a polymer of length l , parametrize the path of the polymer as s ∈ ( 0 , l ) , allow t ^ ( s ) to be the unit tangent vector to the chain at s , and r → ( s ) to be the position vector along the chain. Then

It can be shown that the orientation correlation function for a worm-like chain follows an exponential decay:

where P is by definition the polymer's characteristic persistence length. A useful value is the mean square end-to-end distance of the polymer:

⟨ R 2 ⟩ = ⟨ R → ⋅ R → ⟩ = ⟨ ∫ 0 l t ^ ( s ) d s ⋅ ∫ 0 l t ^ ( s ′ ) d s ′ ⟩ = ∫ 0 l d s ∫ 0 l ⟨ t ^ ( s ) ⋅ t ^ ( s ′ ) ⟩ d s ′ = ∫ 0 l d s ∫ 0 l e − | s − s ′ | / P d s ′

⟨ R 2 ⟩ = 2 P l [ 1 − P l ( 1 − e − l / P ) ] ,

Biological Relevance

Several biologically important polymers can be effectively modeled as worm-like chains, including:

Stretching Worm-like Chain Polymers

At finite temperatures, the distance between the two ends of the polymer (end-to-end distance) will be significantly shorter than the contour length L 0 . This is caused by thermal fluctuations, which result in a coiled, random configuration of the polymer, when undisturbed. Upon stretching the polymer, the accessible spectrum of fluctuations reduces, which causes an entropic force against the external elongation. This entropic force can be estimated by considering the entropic Hamiltonian:

H = H e n t r o p i c + H e x t e r n a l = 1 2 k B T ∫ 0 L 0 P ⋅ ( ∂ 2 r → ( s ) ∂ s 2 ) 2 d s − x F .

Here, the contour length is represented by L 0 , the persistence length by P , the extension is represented by x , and external force is represented by F .

Laboratory tools such as atomic force microscopy (AFM) and optical tweezers have been used to characterize the force-dependent stretching behavior of the polymers listed above. An interpolation formula that approximates the force-extension behavior with about 15% relative error is (J. F. Marko, E. D. Siggia (1995)):

where k B is the Boltzmann constant and T is the absolute temperature. More accurate approximation for the force-extension behavior with about 1% relative error is:

Approximation for the extension-force behavior with about 1% relative error was also reported:

Extensible worm-like chain model

When extending most polymers, their elastic response cannot be neglected. As an example, for the well-studied case of stretching DNA in physiological conditions (near neutral pH, ionic strength approximately 100 mM) at room temperature, the compliance of the DNA along the contour must be accounted for. This enthalpic compliance is accounted for the material parameter K 0 , the stretch modulus. For significantly extended polymers, this yields the following Hamiltonian:

H = H e n t r o p i c + H e n t h a l p i c + H e x t e r n a l = 1 2 k B T ∫ 0 L 0 P ⋅ ( ∂ r → ( s ) ∂ s ) 2 d s + 1 2 K 0 L 0 x 2 − x F ,

with L 0 , the contour length, P , the persistence length, x the extension and F external force. This expression takes into account both the entropic term, which regards changes in the polymer conformation, and the enthalpic term, which describes the elongation of the polymer due to the external force. In the expression above, the enthalpic response is described as a linear Hookian spring. Several approximations have been put forward, dependent on the applied external force. For the low-force regime (F < about 10 pN), the following interpolation formula was derived:

F P k B T = 1 4 ( 1 − x L 0 + F K 0 ) − 2 − 1 4 + x L 0 − F K 0 .

For the higher-force regime, where the polymer is significantly extended, the following approximation is valid:

x = L 0 ( 1 − 1 2 ( k B T F P ) 1 / 2 + F K 0 ) .

A typical value for the stretch modulus of double-stranded DNA is around 1000 pN and 45 nm for the persistence length. More accurate interpolation formulas for the force-extension and extension-force behaviors are: