Strongly bent double-stranded DNA: reconciling theory and experiment

Aleksander V. Drozdetski,Abhishek Mukhopadhyay,Alexey V. Onufriev

DOI: https://doi.org/10.3389/fphy.2019.00195

2019-07-03

Abstract:The strong bending of polymers is poorly understood. We propose a general quantitative framework of polymer bending that includes both the weak and strong bending regimes on the same footing, based on a single general physical principle. As the bending deformation increases beyond a certain (polymer-specific) point, the change in the convexity properties of the effective bending energy of the polymer makes the harmonic deformation energetically unfavorable: in this strong bending regime the energy of the polymer varies linearly with the average bending angle as the system follows the convex hull of the deformation energy function. For double-stranded DNA, the effective bending deformation energy becomes non-convex for bends greater than about 2 degrees per base-pair, equivalent to the curvature of a closed circular loop of about 160 base pairs. A simple equation is derived for the polymer loop energy that covers both the weak and strong bending regimes. The theory shows quantitative agreement with recent DNA cyclization experiments on short DNA fragments, while maintaining the expected agreement with experiment in the weak bending regime. Counter-intuitively, cyclization probability (j-factor) of very short DNA loops is predicted to increase with decreasing loop length; the j-factor reaches its minimum for loops of 45 base pairs. Atomistic simulations reveal that the attractive component of the short-range Lennard-Jones interaction between the backbone atoms can explain the underlying non-convexity of the DNA effective bending energy, leading to the linear bending regime. Applicability of the theory to protein-DNA complexes, including the nucleosome, is discussed.

Biological Physics,Soft Condensed Matter,Biomolecules

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the insufficient understanding of the behavior of double - stranded DNA in a strongly bent state. Specifically, the traditional elastic rod models (such as the worm - like chain model, WLC) show significant differences from experimental results when describing the cyclization probability of short - fragment DNA (less than about 100 base pairs). These experiments show that the cyclization probability of short - fragment DNA is much higher than that predicted by the WLC model, indicating that in a strongly bent state, the elastic properties of DNA may be different from those in a weakly bent state. However, there is currently a lack of a comprehensive and consistent theoretical framework for the reasons for this difference and how to uniformly describe the behavior of DNA in different bent states.

The paper proposes a new theoretical framework, aiming to simultaneously deal with the weakly bent and strongly bent states of polymers (including DNA) by introducing a general physical principle. This framework, based on the convexity characteristics of the effective bending energy function of polymers, proposes that in a strongly bent state, the energy of the polymer has a linear relationship with its average bending angle, rather than the traditional quadratic relationship. In this way, the paper attempts to bridge the gap between theory and experiment and provide a unified model that can explain the bending behavior of DNA at different length scales.

### Main contributions of the paper
1. **Proposed a new theoretical framework**: This framework can simultaneously describe the behavior of polymers in weakly bent and strongly bent states, based on the convexity characteristics of the effective bending energy function of polymers.
2. **Explained experimentally observed phenomena**: Through the new theoretical framework, the paper successfully explained the high cyclization probability of short - fragment DNA, which is inconsistent with the results predicted by the traditional WLC model.
3. **Provided experimental verification**: Through molecular dynamics simulations and coarse - grained simulations, the correctness and applicability of the new theoretical framework were verified.

### Specific problems and solutions
- **Problem**: The traditional WLC model cannot accurately describe the cyclization probability of short - fragment DNA. Experimental results show that short - fragment DNA exhibits higher flexibility in a strongly bent state.
- **Solution**: The paper proposes a new theoretical framework. Based on the non - convexity characteristics of the effective bending energy function of polymers, this framework introduces a linear bending energy section, so it can better describe the behavior of DNA in a strongly bent state.

### Mathematical expressions
The paper proposes a new form of energy function to describe the behavior of DNA in different bent states:

\[ E(\bar{\theta}) = \begin{cases}
\frac{1}{2} k_B T L_p \bar{\theta}^2 & \text{if } \bar{\theta} \leq \theta_a \\
k_B T L_p \theta_a \left( \bar{\theta} - \frac{1}{2} \theta_a \right) & \text{if } \theta_a < \bar{\theta} < \theta_b
\end{cases} \]

where:
- \( \bar{\theta} \) is the average bending angle.
- \( L_p \) is the persistence length of DNA.
- \( \theta_a \) and \( \theta_b \) are the lower and upper limits of the non - convex region, respectively.

Through this new energy function, the paper successfully explains the experimentally observed high cyclization probability of short - fragment DNA and provides a unified theoretical framework to describe the behavior of DNA in different bent states.