Multivalent binding proteins can drive collapse and reswelling of chromatin in confinement

Polymer Physics Models Reveal Structural Folding Features of Single-Molecule Gene Chromatin Conformations

Here, we employ polymer physics models of chromatin to investigate the 3D folding of a 2 Mb wide genomic region encompassing the human LTN1 gene, a crucial DNA locus involved in key cellular functions. Through extensive Molecular Dynamics simulations, we reconstruct in silico the ensemble of single-molecule LTN1 3D structures, which we benchmark against recent in situ Hi-C 2.0 data. The model-derived single molecules are then used to predict structural folding features at the single-cell level, providing testable predictions for super-resolution microscopy experiments.

Connecting Chromatin Structures to Gene Regulation Using Dynamic Polymer Simulations

The transfer of regulatory information between distal loci on chromatin is thought to involve physical proximity, but key biophysical features of these contacts remain unclear. For instance, it is unknown how close and for how long two loci need to be in order to productively interact. The main challenge is that it is currently impossible to measure chromatin dynamics with high spatiotemporal resolution at scale. Polymer simulations provide an accessible and rigorous way to test biophysical models of chromatin regulation, yet there is a lack of simple and general methods for extracting the values of model parameters. Here we adapt the Nelder-Mead simplex optimization algorithm to select the best polymer model matching a given Hi-C dataset, using the MYC locus as an example. The model's biophysical parameters predict a compartmental rearrangement of the MYC locus in leukemia, which we validate with single-cell measurements. Leveraging trajectories predicted by the model, we find that loci with similar Hi-C contact frequencies can exhibit widely different contact dynamics. Interestingly, the frequency of productive interactions between loci exhibits a non-linear relationship with their Hi-C contact frequency when we enforce a specific capture radius and contact duration. These observations are consistent with recent experimental observations and suggest that the dynamic ensemble of chromatin configurations, rather than average contact matrices, is required to fully predict productive long-range chromatin interactions.

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.

Chromatin phase separated nanoregions explored by polymer cross-linker models and reconstructed from single particle trajectories

Phase separated domains (PSDs) are ubiquitous in cell biology, representing nanoregions of high molecular concentration. PSDs appear at diverse cellular domains, such as neuronal synapses but also in eukaryotic cell nucleus, limiting the access of transcription factors and thus preventing gene expression. We develop a generalized cross-linker polymer model, to study PSDs: we show that increasing the number of cross-linkers induces a polymer condensation, preventing access of diffusing molecules. To investigate how the PSDs restrict the motion of diffusing molecules, we compute the mean residence and first escaping times. Finally, we develop a method based on mean-square-displacement of single particle trajectories to reconstruct the properties of PSDs from the continuum range of anomalous exponents. We also show here that PSD generated by polymers do not induces a long-range attracting field (potential well), in contrast with nanodomains at neuronal synapses. To conclude, PSDs can result from condensed chromatin organization, where the number of cross-linkers controls molecular access.

Predicting scale-dependent chromatin polymer properties from systematic coarse-graining

Simulating chromatin is crucial for predicting genome organization and dynamics. Although coarse-grained bead-spring polymer models are commonly used to describe chromatin, the relevant bead dimensions, elastic properties, and the nature of inter-bead potentials are unknown. Using nucleosome-resolution contact probability (Micro-C) data, we systematically coarse-grain chromatin and predict quantities essential for polymer representation of chromatin. We compute size distributions of chromatin beads for different coarse-graining scales, quantify fluctuations and distributions of bond lengths between neighboring regions, and derive effective spring constant values. Unlike the prevalent notion, our findings argue that coarse-grained chromatin beads must be considered as soft particles that can overlap, and we derive an effective inter-bead soft potential and quantify an overlap parameter. We also compute angle distributions giving insights into intrinsic folding and local bendability of chromatin. While the nucleosome-linker DNA bond angle naturally emerges from our work, we show two populations of local structural states. The bead sizes, bond lengths, and bond angles show different mean behavior at Topologically Associating Domain (TAD) boundaries and TAD interiors. We integrate our findings into a coarse-grained polymer model and provide quantitative estimates of all model parameters, which can serve as a foundational basis for all future coarse-grained chromatin simulations.

Dynamics of a spherical self-propelled tracer in a polymeric medium: interplay of self-propulsion, stickiness, and crowding

We employ computer simulations to study the dynamics of a self-propelled spherical tracer particle in a viscoelastic medium, made of a long polymer chain. Here, the interplay between viscoelasticity, stickiness, and activity (self-propulsion) brings additional complexity to the tracer dynamics. Our simulations show that on increasing the stickiness of the tracer particle to the polymer beads, the dynamics of the tracer particle slows down as it gets stuck to the polymer chain and moves along with it. But with increasing self-propulsion velocity, the dynamics gets enhanced. In the case of increasing stickiness as well as activity, the non-Gaussian parameter (NGP) exhibits non-monotonic behavior, which also shows up in the re-scaled self part of the van-Hove function. Non-Gaussianity results owing to the enhanced binding events and the sticky motion of the tracer along with the chain with increasing stickiness. On the other hand, with increasing activity, initially non-Gaussianity increases as the tracer moves through the heterogeneous polymeric environment but for higher activity, the tracer escapes resulting in a negative NGP. For higher values of stickiness, the trapping time distributions of the passive tracer particle broaden and have long tails. On the other hand, for a given stickiness with increasing self-propulsion force, the trapping time distributions become narrower and have short tails. We believe that our current simulation study will be helpful in elucidating the complex motion of activity-driven probes in viscoelastic media.