Marginally compact fractal trees with semiflexibility

Effect of loops on the mean-square displacement of Rouse-model chromatin

Chromatin polymer dynamics are commonly described using the classical Rouse model. The subsequent discovery, however, of intermediate-scale chromatin organization known as topologically associating domains (TADs) in experimental Hi-C contact maps for chromosomes across the tree of life, together with the success of loop extrusion factor (LEF) model in explaining TAD formation, motivates efforts to understand the effect of loops and loop extrusion on chromatin dynamics. This paper seeks to fulfill this need by combining LEF-model simulations with extended Rouse-model polymer simulations to investigate the dynamics of chromatin with loops and dynamic loop extrusion. We show that loops significantly suppress the averaged mean-square displacement (MSD) of a gene locus, consistent with recent experiments that track fluorescently labeled chromatin loci. We also find that loops reduce the MSD's stretching exponent from the classical Rouse-model value of 1/2 to a loop-density-dependent value in the 0.45-0.40 range. Remarkably, stretching exponent values in this range have also been observed in recent experiments [Weber et al., Phys. Rev. Lett. 104, 238102 (2010)0031-900710.1103/PhysRevLett.104.238102; Bailey et al., Mol. Biol. Cell 34, ar78 (2023)1059-152410.1091/mbc.E23-04-0119]. We also show that the dynamics of loop extrusion itself negligibly affects chromatin mobility. By studying static "rosette" loop configurations, we also demonstrate that chromatin MSDs and stretching exponents depend on the location of the locus in question relative to the position of the loops and on the local friction environment.

Influence of Branching on the Configurational and Dynamical Properties of Entangled Polymer Melts

We probe the influence of branching on the configurational, packing, and density correlation function properties of polymer melts of linear and star polymers, with emphasis on molecular masses larger than the entanglement molecular mass of linear chains. In particular, we calculate the conformational properties of these polymers, such as the hydrodynamic radius R h , packing length p, pair correlation function g ( r ) , and polymer center of mass self-diffusion coefficient, D, with the use of coarse-grained molecular dynamics simulations. Our simulation results reproduce the phenomenology of simulated linear and branched polymers, and we attempt to understand our observations based on a combination of hydrodynamic and thermodynamic modeling. We introduce a model of "entanglement" phenomenon in high molecular mass polymers that assumes polymers can viewed in a coarse-grained sense as "soft" particles and, correspondingly, we model the emergence of heterogeneous dynamics in polymeric glass-forming liquids to occur in a fashion similar to glass-forming liquids in which the molecules have soft repulsive interactions. Based on this novel perspective of polymer melt dynamics, we propose a functional form for D that can describe our simulation results for both star and linear polymers, covering both the unentangled to entangled polymer melt regimes.

Many-body contacts in fractal polymer chains and fractional Brownian trajectories

We calculate the probabilities that a trajectory of a fractional Brownian motion with arbitrary fractal dimension d_{f} visits the same spot n≥3 times, at given moments t_{1},...,t_{n}, and obtain a determinant expression for these probabilities in terms of a displacement-displacement covariance matrix. Except for the standard Brownian trajectories with d_{f}=2, the resulting many-body contact probabilities cannot be factorized into a product of single-loop contributions. Within a Gaussian network model of a self-interacting polymer chain, which we suggested recently [K. Polovnikov et al., Soft Matter 14, 6561 (2018)1744-683X10.1039/C8SM00785C], the probabilities we calculate here can be interpreted as probabilities of multibody contacts in a fractal polymer conformation with the same fractal dimension d_{f}. This Gaussian approach, which implies a mapping from fractional Brownian motion trajectories to polymer conformations, can be used as a semiquantitative model of polymer chains in topologically stabilized conformations, e.g., in melts of unconcatenated rings or in the chromatin fiber, which is the material medium containing genetic information. The model presented here can be used, therefore, as a benchmark for interpretation of the data of many-body contacts in genomes, which we expect to be available soon in, e.g., Hi-C experiments.

Influence of knot complexity on glass-formation in low molecular mass ring polymer melts

We perform molecular dynamics simulations on a coarse-grained polymer melt to study the dynamics of glass-formation in ring polymer melts of variable knot complexity. After generating melts of non-concatenated polymeric rings having a range of minimum crossing number values, mc, we compute the coherent intermediate scattering function, the segmental α-relaxation time, fragility, and the glass transition temperature as a function of mc. Variation of knot complexity is found to have a pronounced effect on the dynamics of polymer melts since both molecular rigidity and packing are altered, primary physical factors governing glass-formation in polymeric materials.

Effective Hamiltonian of topologically stabilized polymer states

Topologically stabilized polymer conformations in melts of nonconcatenated polymer rings and crumpled globules are considered to be a good candidate for the description of the spatial structure of mitotic chromosomes. Despite significant efforts, the microscopic Hamiltonian capable of describing such systems still remains unknown. We describe a polymer conformation by a Gaussian network - a system with a Hamiltonian quadratic in all coordinates - and show that by tuning interaction constants, one can obtain equilibrium conformations with any fractal dimension between 2 (an ideal polymer chain) and 3 (a crumpled globule). Monomer-to-monomer distances in topologically stabilized states, according to available numerical data, fit very well the Gaussian distribution, giving an additional argument in support of the quadratic Hamiltonian model. Mathematically, the polymer conformations are mapped onto the trajectories of a subdiffusive fractal Brownian particle. Moreover, we explicitly show that the quadratic Hamiltonian with a hierarchical set of coupling constants provides the microscopic background for the description of the path integral of the fractional Brownian motion with an algebraically decaying kernel.

A comparative study of thermodynamic, conformational, and structural properties of bottlebrush with star and ring polymer melts

Thermodynamic, conformational, and structural properties of bottlebrush polymer melts are investigated with molecular dynamics simulations and compared to linear, regular star, and unknotted ring polymer melts to gauge the influence of molecular topology on polymer melt properties. We focus on the variation of the backbone chain length, the grafting density along the backbone, and the length of the side chains at different temperatures above the melt glass transition temperature. Based on these comparisons, we find that the segmental density, isothermal compressibility, and isobaric thermal expansion of bottlebrush melts are quantitatively similar to unknotted ring polymer melts and star polymer melts having a moderate number ( f = 5 to 6) of arms. These similarities extend to the mass scaling of the chain radius of gyration. Our results together indicate that the configurational properties of bottlebrush polymers in their melt state are more similar to randomly branched polymers than linear polymer chains. We also find that the average shape of bottlebrush polymers having short backbone chains with respect to the side chain length is also rather similar to the unknotted ring and moderately branched star polymers in their melt state. As a general trend, the molecular shape of bottlebrush polymers becomes more spherically symmetric when the length of the side chains has a commensurate length as the backbone chain. Finally, we calculate the partial static structure factor of the backbone segments and we find the emergence of a peak at the length scales that characterizes the average distance between the backbone chains. This peak is absent when we calculate the full static structure factor. We characterize the scaling of this peak with parameters characterizing the bottlebrush molecular architecture to aid in the experimental characterization of these molecules by neutron scattering.