Bending and Base-Stacking Interactions in Double-Stranded Semiflexible Polymer

Zhou Haijun,Ou-Yang Zhong-can

DOI: https://doi.org/10.1103/PhysRevLett.82.4560

1998-11-10

Abstract:Simple expressions for the bending and the base-stacking energy of double-stranded semiflexible biopolymers (such as DNA and actin) are derived. The distribution of the folding angle between the two strands is obtained by solving a Schrödinger equation variationally. Theoretical results based on this model on the extension versus force and extension versus degree of supercoiling relations of DNA chain are in good agreement with the experimental observations of Cluzel {\it et al.} [Science {\bf 271}, 792 (1996)], Smith {\it et al.} [{\it ibid.} {\bf 271}, 795 (1996)], and Strick {\it et al.} [{\it ibid.} {\bf 271}, 1835 (1996)].

Soft Condensed Matter,Statistical Mechanics,Biological Physics,Quantitative Biology

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the theoretical modeling of the bending and base - stacking interactions of double - stranded semi - flexible biopolymers such as DNA and actin. Specifically, the authors aim to:

1. **Derive simple expressions for bending energy**: For double - stranded semi - flexible biopolymers, derive simple expressions that describe their bending and base - stacking energies.
2. **Solve the distribution of folding angles**: Solve a Schrödinger equation by the variational method to obtain the distribution of folding angles between the two chains.
3. **Explain experimental observations**: Use the proposed model to explain the relationship between the elongation of DNA chains, external forces, and superhelicity observed in experiments by Cluzel et al., Smith et al., and Strick et al.

### Main contributions
- **Expression for bending energy**: The authors derived the expression for the total bending energy of double - stranded semi - flexible biopolymers:
\[
E_{\text{bend}}=\kappa\int\left[\left(\frac{d^{2}\mathbf{r}}{ds^{2}}\right)^{2}+R^{2}\left(\frac{d^{2}\mathbf{b}}{ds^{2}}\right)^{2}\right]ds
\]
where \(\kappa\) is the bending stiffness, \(\mathbf{r}\) is the position vector of the central axis, \(\mathbf{b}\) is the unit vector connecting the two chains, and \(R\) is the length of the base pair.

- **Distribution of folding angles**: By solving a variational Schrödinger equation, the distribution function \(\Psi(\theta, s)\) of the folding angle \(\theta\) was obtained:
\[
\frac{\partial\Psi}{\partial s}=-\frac{1}{4\ell_{p}}\left(-\frac{\partial^{2}}{\partial\theta^{2}}+4\eta\sin^{4}\theta\right)\Psi
\]
where \(\eta = (\ell_{p}/R)^{2}\) is a dimensionless constant and \(\ell_{p}\) is the bending persistence length.

- **Comparison between theory and experiment**: The authors calculated the elongation of DNA chains under different external forces and superhelicities using the model and compared it with experimental data. It was found that the theoretical results are highly consistent with the experimental observations of Cluzel et al., Smith et al., and Strick et al.

### Experimental background
In recent years, the development of experimental techniques has made it possible to directly study the mechanical properties of single macromolecules. Many biological macromolecules, such as DNA and proteins (e.g., actin), are considered semi - flexible rather than highly flexible. The traditional worm - like chain model (WLC model) can well describe the elastic behavior of DNA in the low - to - medium external force range, but it cannot explain the non - WLC elastic behavior caused by the double - stranded structure. Therefore, a new model is needed to explain these experimental phenomena.

### Conclusion
By deriving the expression for the bending energy of double - stranded semi - flexible biopolymers and solving the distribution of folding angles, the authors successfully explained the relationship between the elongation of DNA chains, external forces, and superhelicity observed in experiments. This model not only provides theoretical support but also offers a new perspective for understanding the mechanical properties of DNA and other double - stranded biopolymers.