Kinetics of loop formation in worm-like chain polymers

Attractive crowding effect on passive and active polymer looping kinetics

Loop formation in complex environments is crucially important to many biological processes in life. In the present work, we adopt three-dimensional Langevin dynamics simulations to investigate passive and active polymer looping kinetics in crowded media featuring polymer-crowder attraction. We find polymers undergo a remarkable coil-globule-coil transition, highlighted by a marked change in the Flory scaling exponent of the gyration radius. Meanwhile, looping time as a function of the crowder's volume fraction demonstrates an apparent non-monotonic alteration. A small number of crowders induce a compact structure, which largely facilitates the looping process. While a large number of crowders heavily impede end-to-end diffusion, looping kinetics is greatly inhibited. For a self-propelled chain, we find that the attractive crowding triggers an unusual activity effect on looping kinetics. Once a globular state is formed, activity takes an effort to open the chain from the compact structure, leading to an unexpected activity-induced inhibition of looping. If the chain maintains a coil state, the dominant role of activity is to enhance diffusivity and, thus, speed up looping kinetics. The novel conformational change and looping kinetics of both passive and active polymers in the presence of attractive crowding highlight a rather distinct scenario that has no analogy in a repulsive crowding counterpart. The underlying mechanism enriches our understanding of the crucial role of attractive interactions in modulating polymer structure and dynamics.

Thermal fracture kinetics of heterogeneous semiflexible polymers

The fracture and severing of polymer chains plays a critical role in the failure of fibrous materials and the regulated turnover of intracellular filaments. Using continuum wormlike chain models, we investigate the fracture of semiflexible polymers via thermal bending fluctuations, focusing on the role of filament flexibility and dynamics. Our results highlight a previously unappreciated consequence of mechanical heterogeneity in the filament, which enhances the rate of thermal fragmentation particularly in cases where constraints hinder the movement of the chain ends. Although generally applicable to semiflexible chains with regions of different bending stiffness, the model is motivated by a specific biophysical system: the enhanced severing of actin filaments at the boundary between stiff bare regions and mechanically softened regions that are coated with cofilin regulatory proteins. The results presented here point to a potential mechanism for disassembly of polymeric materials in general and cytoskeletal actin networks in particular by the introduction of locally softened chain regions, as occurs with cofilin binding.

Heterogeneous kinetics of the loop formation of a single polymer chain in crowded and disordered media

The cytoplasmic volume of cells is occupied and crowded by a variety of macromolecules, such as proteins and cytoskeleton structures. Such diverse macromolecules make the cell cytoplasm not only structurally heterogeneous but also dynamically heterogeneous: Some macromolecules may diffuse freely inside cell cytoplasm at certain timescales while others hardly diffuse. Studies on the effects of the dynamic heterogeneity on reaction kinetics have been limited even though the effects of the crowdedness and structural heterogeneity were investigated extensively. In this study, we employ a simple model of mixtures of mobile and immobile matrix particles, tune the degree of dynamic heterogeneity by changing the fraction of immobile matrix particles, and investigate reaction kinetics in such heterogeneous media. We employ the loop formation of a single polymer chain as a model reaction and perform Langevin dynamics simulations. We find that the free-energy barrier of the loop formation is decreased as the systems become more crowded with matrix particles. But the free-energy barrier is not sensitive to the dynamic heterogeneity. As dynamic heterogeneity increases with an increase in the fraction of immobile matrix particles, however, the diffusivity of the system decreases significantly. The decrease in the diffusion (due to the dynamic heterogeneity) and the decrease in the free-energy barrier (due to the crowdedness) lead together to a complicated trend of the loop formation kinetics. As the volume fraction of immobile matrix particles reaches a critical value at the percolation transition, the reaction kinetics becomes significantly heterogeneous and the survival probability distribution of the chain loop formation becomes stretched-exponential. We also illustrate that the heterogeneous reaction rate near the percolation transition relates closely to the structures of local pores in which the polymer is located.

Loop formation and translational diffusion of intrinsically disordered proteins

The conformational flexibility and dynamics of unfolded peptide chains is of major interest in the context of protein folding and protein functioning. The rate with which amino acids at different positions along the peptide chain meet sets an upper speed limit for protein folding. By using single-molecule photo-induced energy transfer spectroscopy, we have systematically measured end-to-end and end-to-internal site contact formation rates for several intrinsically disordered protein fragments (IDPs) (11 to 41 amino acids) and have also determined their hydrodynamic radius using dual-focus fluorescence correlation spectroscopy. For interpreting the measured values, we have developed a Brownian dynamics model (based on bead-rod chain dynamics in a thermal bath including hydrodynamic interactions) which quantitatively reproduces all measured data surprisingly well while requiring only two fit parameters. The model provides a complete picture of the peptides' dynamics and allows us to translate the experimental rates and radii into molecular properties of the peptides: We find a persistence length of l_{P}=5.2±1.9Å, a hydrodynamic radius of a=3.5±0.7Å per amino acid, and that excluded volume effects play an important role in the dynamics of IDPs.

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.

Fractional Viscosity Dependence of Reaction Kinetics in Glass-Forming Liquids

The diffusion of molecules in complex systems such as glasses and cell cytoplasm is slow, heterogeneous, and sometimes nonergodic. The effects of such intriguing diffusion on the kinetics of chemical and biological reactions remain elusive. In this Letter, we report that the kinetics of the polymer loop formation reaction in a Kob-Andersen (KA) glass forming liquid is influenced significantly by the dynamic heterogeneity. The diffusion coefficient D of a KA liquid deviates from the Stokes-Einstein relation at low temperatures and D shows a fractional dependence on the solvent viscosity η_{s}, i.e., D∼η_{s}^{-ξ_{D}} with ξ_{D}=0.85. The dynamic heterogeneity of a KA liquid affects the rate constant k_{rxn} of the loop formation and leads to the identical fractional dependence of k_{rxn} on η_{s} with k_{rxn}∼η_{s}^{-ξ} and ξ=ξ_{D}, contrary to reactions in dynamically homogeneous solutions where k_{rxn}∼η_{s}^{-1}.

Interplay of Protein Binding Interactions, DNA Mechanics, and Entropy in DNA Looping Kinetics

DNA looping plays a key role in many fundamental biological processes, including gene regulation, recombination, and chromosomal organization. The looping of DNA is often mediated by proteins whose structural features and physical interactions can alter the length scale at which the looping occurs. Looping and unlooping processes are controlled by thermodynamic contributions associated with mechanical deformation of the DNA strand and entropy arising from thermal fluctuations of the conformation. To determine how these confounding effects influence DNA looping and unlooping kinetics, we present a theoretical model that incorporates the role of the protein interactions, DNA mechanics, and conformational entropy. We show that for shorter DNA strands the interaction distance affects the transition state, resulting in a complex relationship between the looped and unlooped state lifetimes and the physical properties of the looped DNA. We explore the range of behaviors that arise with varying interaction distance and DNA length. These results demonstrate how DNA deformation and entropy dictate the scaling of the looping and unlooping kinetics versus the J-factor, establishing the connection between kinetic and equilibrium behaviors. Our results show how the twist-and-bend elasticity of the DNA chain modulates the kinetics and how the influence of the interaction distance fades away at intermediate to longer chain lengths, in agreement with previous scaling predictions.

Contact Kinetics in Fractal Macromolecules

We consider the kinetics of first contact between two monomers of the same macromolecule. Relying on a fractal description of the macromolecule, we develop an analytical method to compute the mean first contact time for various molecular sizes. In our theoretical description, the non-Markovian feature of monomer motion, arising from the interactions with the other monomers, is captured by accounting for the nonequilibrium conformations of the macromolecule at the very instant of first contact. This analysis reveals a simple scaling relation for the mean first contact time between two monomers, which involves only their equilibrium distance and the spectral dimension of the macromolecule, independently of its microscopic details. Our theoretical predictions are in excellent agreement with numerical stochastic simulations.

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.

Cyclization kinetics of Gaussian semiflexible polymer chains

We consider the dynamics and the cyclization kinetics of Gaussian semiflexible chains, in which the interaction potential tends to align successive bonds. We provide asymptotic expressions for the cyclization time, for the eigenvalues and eigenfunctions, and for the mean square displacement at all time and length scales, with explicit dependence on the capture radius, on the positions of the reactive monomers in the chain, and on the finite number of beads. For the cyclization kinetics, we take into account non-Markovian effects by calculating the distribution of reactive conformations of the polymer, which are not taken into account in the classical Wilemski-Fixman theory. Comparison with numerical simulations confirms the accuracy of this non-Markovian theory.