Kinetics of Interior Loop Formation in Semiflexible Chains

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.