Kinetics of Interior Loop Formation in Semiflexible Chains

TL;DR

This paper develops an analytical model for the kinetics of interior loop formation in semiflexible chains, such as DNA, incorporating effects of stiffness, charge interactions, and comparing predictions with simulations.

Contribution

It introduces a mean-field analytical approach to predict interior looping times in semiflexible chains, accounting for variable stiffness and electrostatic interactions.

Findings

Analytical expressions agree well with Brownian dynamics simulations.

Stiffer dangling ends increase interior looping times.

Electrostatic interactions significantly slow down loop formation at low ionic strength.

Abstract

Loop formation between monomers in the interior of semiflexible chains describes elementary events in biomolecular folding and DNA bending. We calculate analytically the interior distance distribution function for semiflexible chains using a mean-field approach. Using the potential of mean force derived from the distance distribution function we present a simple expression for the kinetics of interior looping by adopting Kramers theory. For the parameters, that are appropriate for DNA, the theoretical predictions in comparison to the case are in excellent agreement with explicit Brownian dynamics simulations of worm-like chain (WLC) model. The interior looping times ($\tau_{IC}$) can be greatly altered in cases when the stiffness of the loop differs from that of the dangling ends. If the dangling end is stiffer than the loop then $\tau_{IC}$ increases for the case of the WLC with uniform persistence length. In contrast, attachment of flexible dangling ends enhances rate of interior loop formation. The theory also shows that if the monomers are charged and interact via screened Coulomb potential then both the cyclization ($\tau_c$) and interior looping ($\tau_{IC}$) times greatly increase at low ionic concentration. Because both $\tau_c$ and $\tau_{IC}$ are determined essentially by the effective persistence length ($l_p^{(R)}$) we computed $l_p^{(R)}$ by varying the range of the repulsive interaction between the monomers. For short range interactions $l_p^{(R)}$ nearly coincides with the bare persistence length which is determined largely by the backbone chain connectivity. This finding rationalizes the efficacy of describing a number of experimental observations (response of biopolymers to force and cyclization kinetics) in biomolecules using WLC model with an effective persistence length.