Non-equilibrium effects of molecular motors on polymers

Excluded volume effects on tangentially driven active ring polymers

The conformational and dynamical properties of active ring polymers are studied by numerical simulations. The two-dimensionally confined polymer is modeled as a closed bead-spring chain, driven by tangential forces, put in contact with a heat bath described by the Brownian multiparticle collision dynamics. Both phantom polymers and chains comprising excluded volume interactions are considered for different bending rigidities. The size and shape are found to be dependent on persistence length, driving force, and bead mutual exclusion. The lack of excluded volume interactions is responsible for a shrinkage of active rings when increasing driving force in the flexible limit, while the presence induces a moderate swelling of chains. The internal dynamics of flexible phantom active rings shows activity-enhanced diffusive behavior at large activity values while, in the case of self-avoiding active chains, it is characterized by active ballistic motion not depending on stiffness. The long-time dynamics of active rings is marked by rotational motion whose period scales as the inverse of the applied tangential force, irrespective of persistence length and beads' self-exclusion.

Transcription-induced active forces suppress chromatin motion

The organization of interphase chromosomes in a number of species is starting to emerge thanks to advances in a variety of experimental techniques. However, much less is known about the dynamics, especially in the functional states of chromatin. Some experiments have shown that the motility of individual loci in human interphase chromosome decreases during transcription and increases upon inhibiting transcription. This is a counterintuitive finding because it is thought that the active mechanical force (F) on the order of ten piconewtons, generated by RNA polymerase II (RNAPII) that is presumably transmitted to the gene-rich region of the chromatin, would render it more open, thus enhancing the mobility. We developed a minimal active copolymer model for interphase chromosomes to investigate how F affects the dynamical properties of chromatin. The movements of the loci in the gene-rich region are suppressed in an intermediate range of F and are enhanced at small F values, which has also been observed in experiments. In the intermediate F, the bond length between consecutive loci increases, becoming commensurate with the distance at the minimum of the attractive interaction between nonbonded loci. This results in a transient disorder-to-order transition, leading to a decreased mobility during transcription. Strikingly, the F-dependent change in the locus dynamics preserves the organization of the chromosome at [Formula: see text]. Transient ordering of the loci, which is not found in the polymers with random epigenetic profiles, in the gene-rich region might be a plausible mechanism for nucleating a dynamic network involving transcription factors, RNAPII, and chromatin.

Self-Organized States from Solutions of Active Ring Polymers in Bulk and under Confinement

In the present work, we study, by means of numerical simulations, the structural and dynamical behavior of a suspension of active ring polymers in bulk and under lateral confinement. At high activity, when changing the distance between the confining planes and the polymers' density, we identify the emergence of a self-organized dynamical state, characterized by the coexistence of slowly diffusing clusters of rotating disks and faster rings moving in between them. We further assess that self-organization is robust in a range of polymer sizes, and we identify a critical value of the activity, necessary to trigger cluster formation. This system has distinctive features resembling at the same time polymers, liquid crystals, and active systems, where the interplay between activity, topology, and confinement leads to a spontaneous segregation in an initially one-component solution.

Conformation and dynamics of partially active linear polymers

We perform numerical simulations of isolated, partially active polymers, driven out-of-equilibrium by a fraction of their monomers. We show that, if the active beads are all gathered in a contiguous block, the position of the section along the chain determines the conformational and dynamical properties of the system. Notably, one can modulate the diffusion coefficient of the polymer from active-like to passive-like just by changing the position of the active block. Further, we show that a slight modification of the self-propulsion rule may give rise to an enhancement of diffusion under certain conditions, despite a decrease of the overall polymer activity. Our findings may help in the modelisation of active biophysical systems, such as filamentous bacteria or worms.

Forced and spontaneous translocation dynamics of a semiflexible active polymer in two dimensions

Polymer translocation is a fundamental topic in non-equilibrium physics and is crucially important to many biological processes in life. In the present work, we adopt two-dimensional Langevin dynamics simulations to study the forced and spontaneous translocation dynamics of an active filament. The influence of polymer stiffness on the underlying dynamics is explicitly analyzed. For the forced translocation, the results show a robust stiffness-induced inhibition, and the translocation time exhibits a dual-exponent scaling relationship with the bending modulus. Tension propagation (TP) is also examined, where we find prominent modifications in terms of both activity and stiffness. For spontaneous translocation into a pure solvent, the translocation time is almost independent of the polymer stiffness. However, when the polymer is translocated into a porous medium, an intriguing non-monotonic alteration of translocation time with increasing chain stiffness is demonstrated. The semiflexible chain is beneficial for translocation while the rigid chain is not conducive. Stiffness regulation on the diffusion dynamics of the polymer in porous media shows a consistent scenario. The interplay of activity, stiffness, and porous crowding provides a new mechanism for understanding the non-trivial translocation dynamics of an active filament in complex environments.

Dynamic ParB-DNA interactions initiate and maintain a partition condensate for bacterial chromosome segregation

In most bacteria, chromosome segregation is driven by the ParABS system where the CTPase protein ParB loads at the parS site to trigger the formation of a large partition complex. Here, we present in vitro studies of the partition complex for Bacillus subtilis ParB, using single-molecule fluorescence microscopy and AFM imaging to show that transient ParB-ParB bridges are essential for forming DNA condensates. Molecular Dynamics simulations confirm that condensation occurs abruptly at a critical concentration of ParB and show that multimerization is a prerequisite for forming the partition complex. Magnetic tweezer force spectroscopy on mutant ParB proteins demonstrates that CTP hydrolysis at the N-terminal domain is essential for DNA condensation. Finally, we show that transcribing RNA polymerases can steadily traverse the ParB-DNA partition complex. These findings uncover how ParB forms a stable yet dynamic partition complex for chromosome segregation that induces DNA condensation and segregation while enabling replication and transcription.

Characteristic features of self-avoiding active Brownian polymers under linear shear flow

We present Brownian dynamics simulation results of a flexible linear polymer with excluded-volume interactions under shear flow in the presence of active noise. The active noise strongly affects the polymer's conformational and dynamical properties, such as the stretching in the flow direction and compression in the gradient direction, shear-induced alignment, and shear viscosity. In the asymptotic limit of large activities and shear rates, the power-law scaling exponents of these quantities differ significantly from those of passive polymers. The chain's shear-induced stretching at a given shear rate is reduced by active noise, and it displays a non-monotonic behavior, where an initial polymer compression is followed by its stretching with increasing active force. The compression of the polymer in the gradient direction follows the relation ∼WiPe-3/4 as a function of the activity-dependent Weissenberg number WiPe, which differs from the scaling observed in passive systems ∼WiPe-1/2. The flow-induced alignment at large Péclet numbers Pe ≫ 1, where Pe is the Péclet number, and large shear rates WiPe ≫ 1 displays the scaling behavior WiPe-1/2, with an exponent differing from the passive value -1/3. Furthermore, the polymer's zero-shear viscosity displays a non-monotonic behavior, decreasing in an intermediate activity regime due to excluded-volume interactions and increasing again for large Pe. Shear thinning appears with increasing Weissenberg number with the power-laws WiPe-1/2 and WiPe-3/4 for passive and active polymers, respectively. In addition, our simulation results are compared with the results of an analytical approach, which predicts quantitatively similar behaviors for the various aforementioned physical quantities.

Structure and dynamics of chemically active ring polymers: swelling to collapse

The ring structures are common in many synthetic or natural systems and experience both local and long-range forces by chemical sensing. This work is an effort to investigate the structural and dynamical properties of a chemically active ring in an explicit solvent bath utilizing hybrid molecular dynamics (MD) and multiparticle collision dynamics (MPCD) simulation techniques. We show that by tuning the chemical properties of the ring, it can be converted from a chemo-attractant to a chemo-repellent, thereby changing the steady state to be either collapsed or swelled as compared to its passive limit. We quantify these observations by comparing the scaling laws, local structures and the dynamics of active and passive rings. Furthermore, we show the impact of varying numbers of active sites by calculating the contact probability of the collapse state that highlights diverse structures. We also analyze the dynamics of the ring by finding the relaxation time and the mean square displacement of the centre of mass. A faster relaxation with enhanced diffusion is observed for the active rings.

Polymer folding through active processes recreates features of genome organization

From proteins to chromosomes, polymers fold into specific conformations that control their biological function. Polymer folding has long been studied with equilibrium thermodynamics, yet intracellular organization and regulation involve energy-consuming, active processes. Signatures of activity have been measured in the context of chromatin motion, which shows spatial correlations and enhanced subdiffusion only in the presence of adenosine triphosphate. Moreover, chromatin motion varies with genomic coordinate, pointing toward a heterogeneous pattern of active processes along the sequence. How do such patterns of activity affect the conformation of a polymer such as chromatin? We address this question by combining analytical theory and simulations to study a polymer subjected to sequence-dependent correlated active forces. Our analysis shows that a local increase in activity (larger active forces) can cause the polymer backbone to bend and expand, while less active segments straighten out and condense. Our simulations further predict that modest activity differences can drive compartmentalization of the polymer consistent with the patterns observed in chromosome conformation capture experiments. Moreover, segments of the polymer that show correlated active (sub)diffusion attract each other through effective long-ranged harmonic interactions, whereas anticorrelations lead to effective repulsions. Thus, our theory offers nonequilibrium mechanisms for forming genomic compartments, which cannot be distinguished from affinity-based folding using structural data alone. As a first step toward exploring whether active mechanisms contribute to shaping genome conformations, we discuss a data-driven approach.

Conformation and dynamics of an active filament in crowded media

The structural and dynamical properties of active filamentous objects under macromolecular crowding have a great relevance in biology. By means of Brownian dynamics simulations, we perform a comparative study for the conformational change and diffusion dynamics of an active chain in pure solvents and in crowded media. Our result shows a robust compaction-to-swelling conformational change with the augment of the Péclet number. The presence of crowding facilitates self-trapping of monomers and, thus, reinforces the activity mediated compaction. In addition, the efficient collisions between the self-propelled monomers and crowders induce a coil-to-globulelike transition, indicated by a marked change of the Flory scaling exponent of the gyration radius. Moreover, the diffusion dynamics of the active chain in crowded solutions demonstrates activity-enhanced subdiffusion. The center of mass diffusion manifests rather new scaling relations with respect to both the chain length and Péclet number. The interplay of chain activity and medium crowding provides a new mechanism to understand the non-trivial properties of active filaments in complex environments.

Dynamics of a two-dimensional active polymer chain with a rotation-restricted active head

The dynamics of a two-dimensional active polymer composed of an active Brownian particle (ABP) at the head and a passive polymer chain is investigated using Langevin dynamics simulation. The ABP experiences a self-propulsion force f_s and a resistance torque M as the passive polymer chain is bonded to the edge of the ABP. M restricts the rotation of the ABP, and thus the dynamics of the ABP and that of the whole active polymer are influenced significantly. Due to this restriction, the persistence time τ_r, which characterizes the random rotation of the ABP, is increased significantly and changes non-monotonically with the rotational friction coefficient η_r. Our simulation results show that the effect of M on the dynamics of the active polymer can be characterized mainly by the change of τ_r. Moreover, the propulsive diffusion coefficient D_P of the whole polymer chain originated from the self-propulsion force can be described by a scaling relation D_P ∝ f_s²τ_r/N²η_t² with η_t the translational friction coefficient and N the polymer length. Our results show that the diffusion is promoted by the resistance torque M and τ_r is a key factor for the diffusion of active polymers.

Experimental verification of arcsine laws in mesoscopic nonequilibrium systems

A large number of processes in the mesoscopic world occur out of equilibrium, where the time evolution of a system becomes immensely important since it is driven principally by dissipative effects. Nonequilibrium steady states (NESS) represent a crucial category in such systems, where relaxation timescales are comparable to the operational timescales. In this study, we employ a model NESS stochastic system, which is comprised of a colloidal microparticle optically trapped in a viscous fluid, externally driven by a temporally correlated noise, and show that time-integrated observables such as the entropic current, the work done on the system or the work dissipated by it, follow the three Lévy arcsine laws [A. C. Barato et al., Phys. Rev. Lett. 121, 090601 (2018)0031-900710.1103/PhysRevLett.121.090601], in the large time limit. We discover that cumulative distributions converge faster to arcsine distributions when it is near equilibrium and the rate of entropy production is small, because in that case the entropic current has weaker temporal autocorrelation. We study this phenomenon by changing the strength of the added noise as well as by perturbing our system with a flow field produced by a microbubble at close proximity to the trapped particle. We confirm our experimental findings with theoretical simulations of the systems. Our work provides an interesting insight into the NESS statistics of the meso-regime, where stochastic fluctuations play a pivotal role.

Kick effect of enzymes causes filament compression

We investigate the influence of enzymes on the structure and dynamics of a filament by dissipative particle dynamics simulations. Enzyme exerts a kick force on the filament monomer. We pay particular attention to two factors: the magnitude of kick force and enzyme concentration. Large kick force as well as high enzyme concentration prefers a remarkable compression of the filament reminiscent of the effective depletion interaction owing to an effective increase in enzyme size and the reduction of solvent quality. Additionally, the kick effect gives rise to an increase of enzyme density from the center-of-mass of the filament to its periphery. Moreover, the increase of enzyme concentration and kick force also causes a decrease in relaxation time. Our finding is helpful to understand the role of catalytic force in chemo-mechano-biological function and the filament behavior under chemical reaction via kick-induced change of solvent quality.

Effects of alignment activity on the collapse kinetics of a flexible polymer

The dynamics of various biological filaments can be understood within the framework of active polymer models. Here we consider a bead-spring model for a flexible polymer chain in which the active interaction among the beads is introduced via an alignment rule adapted from the Vicsek model. Following quenching from the high-temperature coil phase to a low-temperature state point, we study the coarsening kinetics via molecular dynamics (MD) simulations using the Langevin thermostat. For the passive polymer case the low-temperature equilibrium state is a compact globule. The results from our MD simulations reveal that though the globular state is also the typical final state in the active case, the nonequilibrium pathways to arrive at such a state differ from the picture for the passive case due to the alignment interaction among the beads. We notice that deviations from the intermediate "pearl-necklace"-like arrangement, which is observed in the passive case, and the formation of more elongated dumbbell-like structures increase with increasing activity. Furthermore, it appears that while a small active force on the beads certainly makes the coarsening process much faster, there exists a nonmonotonic dependence of the collapse time on the strength of active interaction. We quantify these observations by comparing the scaling laws for the collapse time and growth of pearls with the passive case.

Structure and dynamics of an active polymer adsorbed on the surface of a cylinder

The structure and dynamics of an active polymer on a smooth cylindrical surface are studied by Brownian dynamics simulations. The effect of an active force on the polymer adsorption behavior and the combined effect of chain mobility, length N, rigidity κ, and cylinder radius, R, on the phase diagrams are systemically investigated. We find that complete adsorption is replaced by the irregular alternative adsorption/desorption process at a large driving force. Three typical (spiral, helix-like, and rod-like) conformations of the active polymer are observed, dependent on N, κ, and R. Dynamically, the polymer shows rotational motion in the spiral state, snake-like motion in the intermediate state, and straight translational motion without turning back in the rod-like state. In the spiral state, we find that the rotation velocity ω and the chain length follow a power-law relation ω ∼ N^-0.42, consistent with the torque-balance theory of general Archimedean spirals. And the polymer shows super-diffusive behavior along the cylinder for a long time in the helix-like and rod-like states. Our results highlight that the mobility, rigidity, and curvature of surface can be used to regulate the polymer behavior.

Spatial organization of chromosomes leads to heterogeneous chromatin motion and drives the liquid- or gel-like dynamical behavior of chromatin

Chromosome organization and dynamics are involved in regulating many fundamental processes such as gene transcription and DNA repair. Experiments unveiled that chromatin motion is highly heterogeneous inside cell nuclei, ranging from a liquid-like, mobile state to a gel-like, rigid regime. Using polymer modeling, we investigate how these different physical states and dynamical heterogeneities may emerge from the same structural mechanisms. We found that the formation of topologically associating domains (TADs) is a key driver of chromatin motion heterogeneity. In particular, we showed that the local degree of compaction of the TAD regulates the transition from a weakly compact, fluid state of chromatin to a more compact, gel state exhibiting anomalous diffusion and coherent motion. Our work provides a comprehensive study of chromosome dynamics and a unified view of chromatin motion enabling interpretation of the wide variety of dynamical behaviors observed experimentally across different biological conditions, suggesting that the "liquid" or "solid" state of chromatin are in fact two sides of the same coin.

Conformational statistics of non-equilibrium polymer loops in Rouse model with active loop extrusion

Motivated by the recent experimental observations of the DNA loop extrusion by protein motors, in this paper, we investigate the statistical properties of the growing polymer loops within the ideal chain model. The loop conformation is characterized statistically by the mean gyration radius and the pairwise contact probabilities. It turns out that a single dimensionless parameter, which is given by the ratio of the loop relaxation time over the time elapsed since the start of extrusion, controls the crossover between near-equilibrium and highly non-equilibrium asymptotics in the statistics of the extruded loop, regardless of the specific time dependence of the extrusion velocity. In addition, we show that two-sided and one-sided loop extruding motors produce the loops with almost identical properties. Our predictions are based on two rigorous semi-analytical methods accompanied by asymptotic analysis of slow and fast extrusion limits.

Nanorheology of active-passive polymer mixtures differentiates between linear and ring polymer topology

We study the motion of dispersed nanoprobes in entangled active-passive polymer mixtures. By comparing the two architectures of linear vs. unconcatenated and unknotted circular polymers, we demonstrate that novel, rich physics emerge. For both polymer architectures, nanoprobes of size smaller than the entanglement threshold of the solution move faster as activity is increased and more energy is pumped in the system. For larger nanoprobes, a surprising phenomenon occurs: while in linear solutions they move qualitatively as before, in active-passive ring solutions nanoprobes decelerate with respect to the purely passive conditions. We rationalize this effect in terms of the non-equilibrium, topology-dependent association (clustering) of nanoprobes to the cold component of the ring mixture reminiscent of the recently discovered [Weber et al., Phys. Rev. Lett., 2016, 116, 058301] phase separation in scalar active-passive mixtures. We conclude with a potential connection to the microrheology of the chromatin in the nuclei of the cells.

Activity-Induced Collapse and Arrest of Active Polymer Rings

We investigate, using numerical simulations, the conformations of isolated active ring polymers. We find that their behavior depends crucially on their size: Short rings (N≲100) swell, whereas longer rings (N≳200) collapse, at sufficiently high activity. By investigating the nonequilibrium process leading to the steady state, we find a universal route driving both outcomes; we highlight the central role of steric interactions, at variance with linear chains, and of topology conservation. We further show that the collapsed rings are arrested by looking at different observables, all underlining the presence of an extremely long timescales at the steady state, associated with the internal dynamics of the collapsed section. Finally, we found that in some circumstances the collapsed state spins about its axis.

A semiflexible polymer in a gliding assay: reentrant transition, role of turnover and activity

We consider a model of an extensible semiflexible filament moving in two dimensions on a motility assay of motor proteins represented explicitly as active harmonic linkers. Their heads bind stochastically to polymer segments within a capture radius, and extend along the filament in a directed fashion before detaching. Both the extension and detachment rates are load-dependent and generate an active drive on the filament. The filament undergoes a first order phase transition from the open chain to spiral conformation and shows a reentrant behavior in both the active extension and the turnover, defined as the ratio of attachment-detachment rates. Associated with the phase transition, the size and shape of the polymer change non-monotonically, and the relevant autocorrelation functions display a double-exponential decay. The corresponding correlation times show a maximum signifying the dominance of spirals. The orientational dynamics captures the rotation of spirals, and its correlation time decays with activity as a power law.

Dynamics as a cause for the nanoscale organization of the genome

Chromatin ‘blobs’ were recently identified by live super-resolution imaging of labeled nucleosomes as pervasive but fleeting structural entities. However, the mechanisms leading to the formation of these blobs and their functional implications are unknown. We explore here whether causal relationships exist between parameters that characterize the chromatin blob dynamics and structure, by adapting a framework for spatio-temporal Granger-causality inference. Our analysis reveals that chromatin dynamics is a key determinant for both blob area and local density. Such causality, however, could be demonstrated only in 10–20% of the nucleus, suggesting that chromatin dynamics and structure at the nanometer scale are dominated by stochasticity. We show that the theory of active semiflexible polymers can be invoked to provide potential mechanisms leading to the organization of chromatin into blobs. Our results represent a first step toward elucidating the mechanisms that govern the dynamic and stochastic organization of chromatin in the cell nucleus.

Non-Equilibrium Living Polymers

Systems of "living" polymers are ubiquitous in industry and are traditionally realised using surfactants. Here I first review the theoretical state-of-the-art of living polymers and then discuss non-equilibrium extensions that may be realised with advanced synthetic chemistry or DNA functionalised by proteins. These systems are not only interesting in order to realise novel "living" soft matter but can also shed insight into how genomes are (topologically) regulated in vivo.

How a local active force modifies the structural properties of polymers

We study the dynamics of a polymer, described as a variant of a Rouse chain, driven by an active terminal monomer (head). The local active force induces a transition from a globule-like to an elongated state, as revealed by the study of the end-to-end distance, the variance of which is analytically predicted under suitable approximations. The change in the relaxation times of the Rouse-modes produced by the local self-propulsion is consistent with the transition from globule to elongated conformations. Moreover, also the bond-bond spatial correlation for the chain head are affected by the self-propulsion and a gradient of over-stretched bonds along the chain is observed. We compare our numerical results both with the phenomenological stiff-polymer theory and several analytical predictions in the Rouse-chain approximation.

Conformation and dynamics of a self-avoiding active flexible polymer

We investigate conformations and dynamics of a polymer considering its monomers to be active Brownian particles. This active polymer shows very intriguing physical behavior which is absent in an active Rouse chain. The chain initially shrinks with active force, which starts swelling on further increase in force. The shrinkage followed by swelling is attributed purely to excluded-volume interactions among the monomers. In the swelling regime, the chain shows a crossover from the self-avoiding behavior to the Rouse behavior with scaling exponent ν_{a}≈1/2 for end-to-end distance. The nonmonotonicity in the structure is analyzed through various physical quantities; specifically, radial distribution function of monomers, scattering time, as well as various energy calculations. The chain relaxes faster than the Rouse chain in the intermediate force regime, with a crossover in variation of relaxation time at large active force as given by a power law τ_{r}∼Pe^{-4/3} (Pe is Péclet number).