Reactive conformations and non-Markovian cyclization kinetics of a Rouse polymer

Explicit analytical form for memory kernel in the generalized Langevin equation for end-to-end vector of Rouse chains

The generalized Langevin equation (GLE) provides an attractive theoretical framework for investigating the dynamics of conformational fluctuations of polymeric systems. While the memory kernel is a central function in the GLE, explicit analytical forms for this function have been challenging to obtain, even for the simple models of polymer dynamics. Here, we achieve an explicit analytical expression for the memory kernel in the GLE for the end-to-end vector of Rouse chains in the overdamped limit. Our derivation takes advantage of the finding that the dynamics of the end-to-end vector of Rouse chains with both free ends are equivalent to those of Rouse chains with one free end and the other fixed. For the latter model, we first show that the equations of motion of the Rouse modes as well as their statistical properties can be obtained under the boundary conditions where the free end is held fixed temporarily. We then analytically solve the terms associated with intrachain interactions in the GLE. By formally comparing these terms with the GLE based on the Rouse modes, we obtain an explicit expression for the memory kernel, along with analytical forms for the potential field and the random colored noise force. Our analytical memory kernel is confirmed by numerical calculations in the Laplace space and is shown to yield asymptotic behaviors that are consistent with previous studies. Finally, we utilize our analytical result to simulate the cyclization dynamics of Rouse chains and discuss the scaling of the cyclization time with chain length.

Cyclization and Relaxation Dynamics of Finite-Length Collapsed Self-Avoiding Polymers

We study the cyclization and relaxation dynamics of ideal as well as interacting polymers as a function of chain length N. For the cyclization time τ_{cyc} of ideal chains we recover the known scaling τ_{cyc}∼N^{2} for different backbone models, for a self-avoiding slightly collapsed chain we obtain from Langevin simulations and scaling theory a modified scaling τ_{cyc}∼N^{5/3}. The cyclization and relaxation dynamics of a finite-length collapsed chain scale differently; this unexpected dynamic multiscale behavior is rationalized by the crossover between swollen and collapsed chain behavior.

Redundancy principle and the role of extreme statistics in molecular and cellular biology

The paradigm of chemical activation rates in cellular biology has been shifted from the mean arrival time of a single particle to the mean of the first among many particles to arrive at a small activation site. The activation rate is set by extremely rare events, which have drastically different time scales from the mean times between activations, and depends on different structural parameters. This shift calls for reconsideration of physical processes used in deterministic and stochastic modeling of chemical reactions that are based on the traditional forward rate, especially for fast activation processes in living cells. Consequently, the biological activation time is not necessarily exponentially distributed. We review here the physical models, the mathematical analysis and the new paradigm of setting the scale to be the shortest time for activation that clarifies the role of population redundancy in selecting and accelerating transient cellular search processes. We provide examples in cellular transduction, gene activation, cell senescence activation or spermatozoa selection during fertilization, where the rate depends on numbers. We conclude that the statistics of the minimal time to activation set kinetic laws in biology, which can be very different from the ones associated to average times.

Non-Markovian closure kinetics of flexible polymers with hydrodynamic interactions

This paper presents a theoretical analysis of the closure kinetics of a polymer with hydrodynamic interactions. This analysis, which takes into account the non-Markovian dynamics of the end-to-end vector and relies on the preaveraging of the mobility tensor (Zimm dynamics), is shown to reproduce very accurately the results of numerical simulations of the complete nonlinear dynamics. It is found that Markovian treatments based on a Wilemski-Fixman approximation significantly overestimate cyclization times (up to a factor 2), showing the importance of memory effects in the dynamics. In addition, this analysis provides scaling laws of the mean first cyclization time (MFCT) with the polymer size N and capture radius b, which are identical in both Markovian and non-Markovian approaches. In particular, it is found that the scaling of the MFCT for large N is given by T ∼ N(3/2)ln(N/b(2)), which differs from the case of the Rouse dynamics where T ∼ N(2). The extension to the case of the reaction kinetics of a monomer of a Zimm polymer with an external target in a confined volume is also presented.

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.

Beyond the Markovian assumption: response–excitation probabilistic solution to random nonlinear differential equations in the long time

Uncertainty quantification for dynamical systems under non-white excitation is a difficult problem encountered across many scientific and engineering disciplines. Difficulties originate from the lack of Markovian character of system responses. The response–excitation (RE) theory, recently introduced by Sapsis & Athanassoulis (Sapsis & Athanassoulis 2008 Probabilistic Eng. Mech. 23, 289–306 ( doi:10.1016/j.probengmech.2007.12.028 )) and further studied by Venturi et al. (Venturi et al. 2012 Proc. R. Soc. A 468, 759–783 ( doi:10.1098/rspa.2011.0186 )), is a new approach, based on a simple differential constraint which is exact but non-closed. The evolution equation obtained for the RE probability density function (pdf) has the form of a generalized Liouville equation, with the excitation time frozen in the time-derivative term. In this work, the missing information of the RE differential constraint is identified and a closure scheme is developed for the long-time, stationary, limit-state of scalar nonlinear random differential equations (RDEs) under coloured excitation. The closure scheme does not alter the RE evolution equation, but collects the missing information through the solution of local statistically linearized versions of the nonlinear RDE, and interposes it into the solution scheme. Numerical results are presented for two examples, and compared with Monte Carlo simulations.

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.

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.

Gaussian semiflexible rings under angular and dihedral restrictions

Semiflexible polymer rings whose bonds obey both angular and dihedral restrictions [M. Dolgushev and A. Blumen, J. Chem. Phys. 138, 204902 (2013)], are treated under exact closure constraints. This allows us to obtain semianalytic results for their dynamics, based on sets of Langevin equations. The dihedral restrictions clearly manifest themselves in the behavior of the mean-square monomer displacement. The determination of the equilibrium ring conformations shows that the dihedral constraints influence the ring curvature, leading to compact folded structures. The method for imposing such constraints in Gaussian systems is very general and it allows to account for heterogeneous (site-dependent) restrictions. We show it by considering rings in which one site differs from the others.