Cyclization of a polymer with charged reactive end groups

The role of solvent quality and chain stiffness on the end-to-end contact kinetics of semiflexible polymers

Motivated by loop closure during protein folding and DNA packing, we systemically studied the effects of the solvent quality and chain stiffness on the thermodynamics and kinetics of the end-to-end contact formation for semiflexible polymer chains with reactive ends by Langevin dynamics simulations. In thermodynamics, a rich variety of products of the end-to-end contact have been discovered, such as loop, hairpin, toroid, and rodlike bundle, the populations of which are dependent on the solvent quality and chain stiffness. In kinetics, the overall pathways to form the end-to-end contact have been identified. The change of solvent quality and chain stiffness can tune the roughness of energy landscape and modulate the kinetic partitioning of the end-to-end contact formation pathways, leading to differing kinetic behaviors. In good or poor solvents, the first end-to-end contact rate k _c decreases with increasing the strength of bending stiffness k _θ monotonically. In very poor solvents, however, the dependence of the logarithm of the first end-to-end contact rate ln k _c on k _θ exhibits erratic behavior, which stems from more rugged energy landscape due to the polymer chain getting trapped into the intermediate state composed of the rodlike bundles with two ends in separation. For semiflexible chains, with increasing chain length N, the rate k _c increases initially and then decreases: in good solvents, the rate k _c exhibits a power-law relationship to chain length N with an exponent of ∼-1.50 in the region of long chains, which is in good agreement with the value derived from the experiment in the asymptotic limit of large N; and in poor solvents, the rate k _c exhibits a significantly stronger chain length dependence than those observed in good solvents in the region of long chains due to frustration to form the end-to-end contact along a specific path, especially the scaling exponent between the rate k _c and chain length N is ∼-3.62 for the case of polymer chains with k _θ = 4 at the solvent quality ε _ij = 1, in accord with the value obtained from the experiments.

Effects of the internal friction and the solvent quality on the dynamics of a polymer chain closure

Using 3D Langevin dynamics simulations, we investigate the effects of the internal friction and the solvent quality on the dynamics of a polymer chain closure. We show that the chain closure in good solvents is a purely diffusive process. By extrapolation to zero solvent viscosity, we find that the internal friction of a chain plays a non-ignorable role in the dynamics of the chain closure. When the solvent quality changes from good to poor, the mean closure time τc decreases by about 1 order of magnitude for the chain length 20 ≤ N ≤ 100. Furthermore, τc has a minimum as a function of the solvent quality. With increasing the chain length N, the minimum of τc occurs at a better solvent. Finally, the single exponential distributions of the closure time in poor solvents suggest that the negative excluded volume of segments does not alter the nearly Poisson statistical characteristics of the process of the chain closure.

Kinetics of polymer looping with macromolecular crowding: effects of volume fraction and crowder size

The looping of polymers such as DNA is a fundamental process in the molecular biology of living cells, whose interior is characterised by a high degree of molecular crowding. We here investigate in detail the looping dynamics of flexible polymer chains in the presence of different degrees of crowding. From the analysis of the looping-unlooping rates and the looping probabilities of the chain ends we show that the presence of small crowders typically slows down the chain dynamics but larger crowders may in fact facilitate the looping. We rationalise these non-trivial and often counterintuitive effects of the crowder size on the looping kinetics in terms of an effective solution viscosity and standard excluded volume. It is shown that for small crowders the effect of an increased viscosity dominates, while for big crowders we argue that confinement effects (caging) prevail. The tradeoff between both trends can thus result in the impediment or facilitation of polymer looping, depending on the crowder size. We also examine how the crowding volume fraction, chain length, and the attraction strength of the contact groups of the polymer chain affect the looping kinetics and hairpin formation dynamics. Our results are relevant for DNA looping in the absence and presence of protein mediation, DNA hairpin formation, RNA folding, and the folding of polypeptide chains under biologically relevant high-crowding conditions.

Computation of the mean first-encounter time between the ends of a polymer chain

Using a novel theoretical approach, we study the mean first-encounter time (MFET) between the two ends of a polymer. Previous approaches used various simplifications that reduced the complexity of the problem, leading, however, to incompatible results. We construct here for the first time a general theory that allows us to compute the MFET. The method is based on estimating the mean time for a Brownian particle to reach a narrow domain in the polymer configuration space. In dimension two and three, we find that the MFET depends mainly on the first eigenvalue of the associated Fokker-Planck operator and provide precise estimates that are confirmed by Brownian simulations. Interestingly, although many time scales are involved in the encounter process, its distribution can be well approximated by a single exponential, which has several consequences for modeling chromosome dynamics in the nucleus. Another application of our result is computing the mean time for a DNA molecule to form a closed loop (when its two ends meet for the first time).

Excluded volume effects on the intrachain reaction kinetics

On the basis of the recently developed optimized Rouse-Zimm theory of chain polymers with excluded volume interactions, we calculate the long-time first-order rate constant k(1) for end-to-end cyclization of linear chain polymers. We first find that the optimized Rouse-Zimm theory provides the longest chain relaxation times tau(1) of excluded volume chains that are in excellent agreement with the available Brownian dynamics simulation results. In the free-draining limit, the cyclization rate is diffusion-controlled and k(1) is inversely proportional to tau(1), and the k(1) values calculated using the Wilemski-Fixman rate theory are in good agreement with Brownian dynamics simulation results. However, when hydrodynamic interactions are included, noticeable deviations are found. The main sources of errors are fluctuating hydrodynamic interaction and correlation hole effects as well as the non-Markovian reaction dynamic effect. The physical natures of these factors are discussed, and estimates for the magnitudes of required corrections are given. When the corrections are included, the present theory allows the prediction of accurate k(1) values for the cyclization of finite-length chains in good solvents as well as the correct scaling exponent in the long-chain limit.

Cyclization of Rouse chains at long- and short-time scales

We have investigated cyclization of a Rouse chain at long and short times by a Langevin dynamics simulation method. We measure St, the fraction of nonreacted chains, for polymerizations ranging from Z=5 to Z=800 and capture distances ranging from a=0.1b to a=8b where b is the bond length. Comparison is made with two theoretical approaches. The first is a decoupling approximation used by Wilemski and Fixman to close the relevant master equation [J. Chem. Phys. 58, 4009 (1973); 60, 866 (1974)]. The second approach is the renormalization group arguments of Friedman and O'Shaughnessy [Phys. Rev. Lett 60, 64 (1988); J. Phys. II 1, 471 (1991)]. We find that at long times St decays as a single exponential with rate k(infinity). The scaled decay rate K=k(infinity)tauR appears to approach a constant value independent of the capture distance for very large chains consistent with the predictions of both the renormalization group (RG) and Wilemski-Fixman closure approximation. We extract K*, the long chain limit of K, from the fixed point a=a* where K is independent of Z. K* is larger than both the RG and closure predictions but much closer to the RG result. More convincing evidence for the RG analysis is obtained by comparing the short-time decay of St to long-time results. The RG analysis predicts that dSdt should decay as a power law at early times and that the exponent in the power law is related to K by a simple expression with no free parameters. Our simulations find remarkable agreement with this parameter-free prediction even for relatively short chains. We discuss possible experimental consequences of our result.