Motorized chain models of the ideal chromosome.

An array of motor proteins consumes chemical energy in setting up the architectures of chromosomes. Here, we explore how the structure of ideal polymer chains is influenced by two classes of motors. The first class which we call "swimming motors" acts to propel the chromatin fiber through three-dimensional space. They represent a caricature of motors such as RNA polymerases. Previously, they have often been described by adding a persistent flow onto Brownian diffusion of the chain. The second class of motors, which we call "grappling motors" caricatures the loop extrusion processes in which segments of chromatin fibers some distance apart are brought together. We analyze these models using a self-consistent variational phonon approximation to a many-body Master equation incorporating motor activities. We show that whether the swimming motors lead to contraction or expansion depends on the susceptibility of the motors, that is, how their activity depends on the forces they must exert. Grappling motors in contrast to swimming motors lead to long-ranged correlations that resemble those first suggested for fractal globules and that are consistent with the effective interactions inferred by energy landscape analyses of Hi-C data on the interphase chromosome.

Rigidity percolation and active advection synergize in the actomyosin cortex to drive amoeboid cell motility.

Spontaneous locomotion is a common feature of most metazoan cells, generally attributed to the properties of actomyosin networks. This force-producing machinery has been studied down to the most minute molecular details, especially in lamellipodium-driven migration. Nevertheless, how actomyosin networks work inside contraction-driven amoeboid cells still lacks unifying principles. Here, using stable motile blebs from HeLa cells as a model amoeboid motile system, we imaged the dynamics of the actin cortex at the single filament level and revealed the co-existence of three distinct rheological phases. We introduce "advected percolation," a process where rigidity percolation and active advection synergize, spatially organizing the actin network's mechanical properties into a minimal and generic locomotion mechanism. Expanding from our observations on simplified systems, we speculate that this model could explain, down to the single actin filament level, how amoeboid cells, such as cancer or immune cells, can propel efficiently through complex 3D environments.

Physical modeling of nucleosome clustering in euchromatin resulting from interactions between epigenetic reader proteins.

Euchromatin is an accessible phase of genetic material containing genes that encode proteins with increased expression levels. The structure of euchromatin in vitro has been described as a 30-nm fiber formed from ordered nucleosome arrays. However, recent advances in microscopy have revealed an in vivo euchromatin architecture that is much more disordered, characterized by variable-length linker DNA and sporadic nucleosome clusters. In this work, we develop a theoretical model to elucidate factors contributing to the disordered in vivo architecture of euchromatin. We begin by developing a 1D model of nucleosome positioning that captures the interactions between bound epigenetic reader proteins to predict the distribution of DNA linker lengths between adjacent nucleosomes. We then use the predicted linker lengths to construct 3D chromatin configurations consistent with the physical properties of DNA within the nucleosome array, and we evaluate the distribution of nucleosome cluster sizes in those configurations. Our model reproduces experimental cluster-size distributions, which are dramatically influenced by the local pattern of epigenetic marks and the concentration of reader proteins. Based on our model, we attribute the disordered arrangement of euchromatin to the heterogeneous binding of reader proteins and subsequent short-range interactions between bound reader proteins on adjacent nucleosomes. By replicating experimental results with our physics-based model, we propose a mechanism for euchromatin organization in the nucleus that impacts gene regulation and the maintenance of epigenetic marks.

Compaction and Segregation of DNA in Escherichia coli.

Theoretical and experimental approaches have been applied to study the polymer physics underlying the compaction of DNA in the bacterial nucleoid. Knowledge of the compaction mechanism is necessary to obtain a mechanistic understanding of the segregation process of replicating chromosome arms (replichores) during the cell cycle. The first part of this review discusses light microscope observations demonstrating that the nucleoid has a lower refractive index and thus, a lower density than the cytoplasm. A polymer physics explanation for this phenomenon was given by a theory discussed at length in this review. By assuming a phase separation between the nucleoid and the cytoplasm and by imposing equal osmotic pressure and chemical potential between the two phases, a minimal energy situation is obtained, in which soluble proteins are depleted from the nucleoid, thus explaining its lower density. This theory is compared to recent views on DNA compaction that are based on the exclusion of polyribosomes from the nucleoid or on the transcriptional activity of the cell. These new views prompt the question of whether they can still explain the lower refractive index or density of the nucleoid. In the second part of this review, we discuss the question of how DNA segregation occurs in Escherichia coli in the absence of the so-called active ParABS system, which is present in the majority of bacteria. How is the entanglement of nascent chromosome arms generated at the origin in the parental DNA network of the E. coli nucleoid prevented? Microscopic observations of the position of fluorescently-labeled genetic loci have indicated that the four nascent chromosome arms synthesized in the initial replication bubble segregate to opposite halves of the sister nucleoids. This implies that extensive intermingling of daughter strands does not occur. Based on the hypothesis that leading and lagging replichores synthesized in the replication bubble fold into microdomains that do not intermingle, a passive four-excluding-arms model for segregation is proposed. This model suggests that the key for segregation already exists in the structure of the replication bubble at the very start of DNA replication; it explains the different patterns of chromosome arms as well as the segregation distances between replicated loci, as experimentally observed.

Activity-driven chromatin organization during interphase: Compaction, segregation, and entanglement suppression.

In mammalian cells, the cohesin protein complex is believed to translocate along chromatin during interphase to form dynamic loops through a process called active loop extrusion. Chromosome conformation capture and imaging experiments have suggested that chromatin adopts a compact structure with limited interpenetration between chromosomes and between chromosomal sections. We developed a theory demonstrating that active loop extrusion causes the apparent fractal dimension of chromatin to cross-over between two and four at contour lengths on the order of 30 kilo-base pairs. The anomalously high fractal dimension [Formula: see text] is due to the inability of extruded loops to fully relax during active extrusion. Compaction on longer contour length scales extends within topologically associated domains (TADs), facilitating gene regulation by distal elements. Extrusion-induced compaction segregates TADs such that overlaps between TADs are reduced to less than 35% and increases the entanglement strand of chromatin by up to a factor of 50 to several Mega-base pairs. Furthermore, active loop extrusion couples cohesin motion to chromatin conformations formed by previously extruding cohesins and causes the mean square displacement of chromatin loci during lag times ([Formula: see text]) longer than tens of minutes to be proportional to [Formula: see text]. We validate our results with hybrid molecular dynamics-Monte Carlo simulations and show that our theory is consistent with experimental data. This work provides a theoretical basis for the compact organization of interphase chromatin, explaining the physical reason for TAD segregation and suppression of chromatin entanglements which contribute to efficient gene regulation.

Coarse-Grained Simulations on Polyethylene Crystal Network Formation and Microstructure Analysis.

Understanding and characterizing semi-crystalline models with crystalline and amorphous segments is crucial for industrial applications. A coarse-grained molecular dynamics (CGMD) simulations study probed the crystal network formation in high-density polyethylene (HDPE) from melt, and shed light on tensile properties for microstructure analysis. Modified Paul-Yoon-Smith (PYS/R) forcefield parameters are used to compute the interatomic forces among the PE chains. The isothermal crystallization at 300 K and 1 atm predicts the multi-nucleus crystal growth; moreover, the lamellar crystal stems and amorphous region are alternatively oriented. A one-dimensional density distribution along the alternative lamellar stems further confirms the ordering of the lamellar-stack orientation. Using this plastic model preparation approach, the semi-crystalline model density (ρcr) of ca. 0.913 g·cm-3 and amorphous model density (ρam) of ca. 0.856 g·cm-3 are obtained. Furthermore, the ratio of ρcr/ρam ≈ 1.06 is in good agreement with computational (≈1.096) and experimental (≈1.14) data, ensuring the reliability of the simulations. The degree of crystallinity (χc) of the model is ca. 52% at 300 K. Nevertheless, there is a gradual increase in crystallinity over the specified time, indicating the alignment of the lamellar stems during crystallization. The characteristic stress-strain curve mimicking tensile tests along the z-axis orientation exhibits a reversible sharp elastic regime, tensile strength at yield ca. 100 MPa, and a non-reversible tensile strength at break of 350%. The cavitation mechanism embraces the alignment of lamellar stems along the deformation axis. The study highlights an explanatory model of crystal network formation for the PE model using a PYS/R forcefield, and it produces a microstructure with ordered lamellar and amorphous segments with robust mechanical properties, which aids in predicting the microstructure-mechanical property relationships in plastics under applied forces.

Activity-driven chromatin organization during interphase: compaction, segregation, and entanglement suppression.

In mammalian cells, the cohesin protein complex is believed to translocate along chromatin during interphase to form dynamic loops through a process called active loop extrusion. Chromosome conformation capture and imaging experiments have suggested that chromatin adopts a compact structure with limited interpenetration between chromosomes and between chromosomal sections. We developed a theory demonstrating that active loop extrusion causes the apparent fractal dimension of chromatin to cross over between two and four at contour lengths on the order of 30 kilo-base pairs (kbp). The anomalously high fractal dimension D=4 is due to the inability of extruded loops to fully relax during active extrusion. Compaction on longer contour length scales extends within topologically associated domains (TADs), facilitating gene regulation by distal elements. Extrusion-induced compaction segregates TADs such that overlaps between TADs are reduced to less than 35% and increases the entanglement strand of chromatin by up to a factor of 50 to several Mega-base pairs. Furthermore, active loop extrusion couples cohesin motion to chromatin conformations formed by previously extruding cohesins and causes the mean square displacement of chromatin loci during lag times (Δt) longer than tens of minutes to be proportional to Δt1/3. We validate our results with hybrid molecular dynamics - Monte Carlo simulations and show that our theory is consistent with experimental data. This work provides a theoretical basis for the compact organization of interphase chromatin, explaining the physical reason for TAD segregation and suppression of chromatin entanglements which contribute to efficient gene regulation.

The Challenges Facing the Current Paradigm Describing Viscoelastic Interactions in Polymer Melts.

Staudinger taught us that macromolecules were made up of covalently bonded monomer repeat units chaining up as polymer chains. This paradigm is not challenged in this paper. The main question raised in polymer physics remains: how do these long chains interact and move as a group when submitted to shear deformation at high temperature when they are viscous liquids? The current consensus is that we need to distinguish two cases: the deformation of "un-entangled chains" for macromolecules with molecular weight, M, smaller than Me, "the entanglement molecular weight", and the deformation of "entangled" chains for M > Me. The current paradigm stipulates that the properties of polymers derive from the statistical characteristics of the macromolecule itself, the designated statistical system that defines the thermodynamic state of the polymer. The current paradigm claims that the viscoelasticity of un-entangled melts is well described by the Rouse model and that the entanglement issues raised when M > Me, are well understood by the reptation model introduced by de Gennes and colleagues. Both models can be classified in the category of "chain dynamics statistics". In this paper, we examine in detail the failures and the current challenges facing the current paradigm of polymer rheology: the Rouse model for un-entangled melts, the reptation model for entangled melts, the time-temperature superposition principle, the strain-induced time dependence of viscosity, shear-refinement and sustained-orientation. The basic failure of the current paradigm and its inherent inability to fully describe the experimental reality is documented in this paper. In the discussion and conclusion sections of the paper, we suggest that a different solution to explain the viscoelasticity of polymer chains and of their "entanglement" is needed. This requires a change in paradigm to describe the dynamics of the interactions within the chains and across the chains. A brief description of our currently proposed open dissipative statistical approach, "the Grain-Field Statistics", is presented.

Toward a 3D physical model of diffusive polymer chains.

Recent studies in polymer physics have created macro-scale analogs to solute microscopic polymer chains like DNA by inducing diffusive motion on a chain of beads. These bead chains have persistence lengths of O(10) links and undergo diffusive motion under random fluctuations like vibration. We present a bead chain model within a new stochastic forcing system: an air fluidizing bed of granular media. A chain of spherical 6 mm resin beads crimped onto silk thread are buffeted randomly by the multiphase flow of grains and low density rising air "bubbles". We "thermalize" bead chains of various lengths at different fluidizing airflow rates, while X-ray imaging captures a projection of the chains' dynamics within the media. With modern 3D printing techniques, we can better represent complex polymers by geometrically varying bead connections and their relative strength, e.g., mimicking the variable stiffness between adjacent nucleotide pairs of DNA. We also develop Discrete Element Method (DEM) simulations to study the 3D motion of the bead chain, where the bead chain is represented by simulated spherical particles connected by linear and angular spring-like bonds. In experiment, we find that the velocity distributions of the beads follow exponential distributions rather than the Gaussian distributions expected from polymers in solution. Through use of the DEM simulation, we find that this difference can likely be attributed to the distributions of the forces imparted onto the chain from the fluidized bed environment. We anticipate expanding this study in the future to explore a wide range of chain composition and confinement geometry, which will provide insights into the physics of large biopolymers.

Intramolecular structural heterogeneity altered by long-range contacts in an intrinsically disordered protein.

Short-range interactions and long-range contacts drive the 3D folding of structured proteins. The proteins' structure has a direct impact on their biological function. However, nearly 40% of the eukaryotes proteome is composed of intrinsically disordered proteins (IDPs) and protein regions that fluctuate between ensembles of numerous conformations. Therefore, to understand their biological function, it is critical to depict how the structural ensemble statistics correlate to the IDPs' amino acid sequence. Here, using small-angle X-ray scattering and time-resolved Förster resonance energy transfer (trFRET), we study the intramolecular structural heterogeneity of the neurofilament low intrinsically disordered tail domain (NFLt). Using theoretical results of polymer physics, we find that the Flory scaling exponent of NFLt subsegments correlates linearly with their net charge, ranging from statistics of ideal to self-avoiding chains. Surprisingly, measuring the same segments in the context of the whole NFLt protein, we find that regardless of the peptide sequence, the segments' structural statistics are more expanded than when measured independently. Our findings show that while polymer physics can, to some level, relate the IDP's sequence to its ensemble conformations, long-range contacts between distant amino acids play a crucial role in determining intramolecular structures. This emphasizes the necessity of advanced polymer theories to fully describe IDPs ensembles with the hope that it will allow us to model their biological function.

Physics-Based Polymer Models to Probe Chromosome Structure in Single Molecules.

Human chromosomes have a complex 3D spatial organization in the cell nucleus, which comprises a hierarchy of physical interactions across genomic scales. Such an architecture serves important functional roles, as genes and their regulators have to physically interact to control gene regulation. However, the molecular mechanisms underlying the formation of those contacts remain poorly understood. Here, we describe a polymer-physics-based approach to investigate the machinery shaping genome folding and function. In silico model predictions on DNA single-molecule 3D structures are validated against independent super-resolution single-cell microscopy data, supporting a scenario whereby chromosome architecture is controlled by thermodynamics mechanisms of phase separation. Finally, as an application of our methods, the validated single-polymer conformations of the theory are used to benchmark powerful technologies to probe genome structure, such as Hi-C, SPRITE, and GAM.

Theory of chromatin organization maintained by active loop extrusion.

The active loop extrusion hypothesis proposes that chromatin threads through the cohesin protein complex into progressively larger loops until reaching specific boundary elements. We build upon this hypothesis and develop an analytical theory for active loop extrusion which predicts that loop formation probability is a nonmonotonic function of loop length and describes chromatin contact probabilities. We validate our model with Monte Carlo and hybrid Molecular Dynamics-Monte Carlo simulations and demonstrate that our theory recapitulates experimental chromatin conformation capture data. Our results support active loop extrusion as a mechanism for chromatin organization and provide an analytical description of chromatin organization that may be used to specifically modify chromatin contact probabilities.

The Bacterial Nucleoid: From Electron Microscopy to Polymer Physics-A Personal Recollection.

In the 1960s, electron microscopy did not provide a clear answer regarding the compact or dispersed organization of the bacterial nucleoid. This was due to the necessary preparation steps of fixation and dehydration (for embedding) and freezing (for freeze-fracturing). Nevertheless, it was possible to measure the lengths of nucleoids in thin sections of slow-growing Escherichia coli cells, showing their gradual increase along with cell elongation. Later, through application of the so-called agar filtration method for electron microscopy, we were able to perform accurate measurements of cell size and shape. The introduction of confocal and fluorescence light microscopy enabled measurements of size and position of the bacterial nucleoid in living cells, inducing the concepts of "nucleoid occlusion" for localizing cell division and of "transertion" for the final step of nucleoid segregation. The question of why the DNA does not spread throughout the cytoplasm was approached by applying polymer-physical concepts of interactions between DNA and proteins. This gave a mechanistic insight in the depletion of proteins from the nucleoid, in accordance with its low refractive index observed by phase-contrast microscopy. Although in most bacterial species, the widely conserved proteins of the ParABS-system play a role in directing the segregation of newly replicated DNA strands, the basis for the separation and opposing movement of the chromosome arms was proposed to lie in preventing intermingling of nascent daughter strands already in the early replication bubble. E. coli, lacking the ParABS system, may be suitable for investigating this basic mechanism of DNA strand separation and segregation.

Structural changes in chromosomes driven by multiple condensin motors during mitosis.

We create a computational framework that utilizes loop extrusion (LE) by multiple condensin I/II motors to predict changes in chromosome organization during mitosis. The theory accurately reproduces the experimental contact probability profiles for the mitotic chromosomes in HeLa and DT40 cells. The LE rate is smaller at the start of mitosis and increases as the cells approach metaphase. Condensin II-mediated mean loop size is about six times larger than loops because of condensin I. The loops, which overlap each other, are stapled to a central dynamically changing helical scaffold formed by the motors during the LE process. A polymer physics-based data-driven method that uses the Hi-C contact map as the only input shows that the helix is characterized as random helix perversions (RHPs) in which the handedness changes randomly along the scaffold. The theoretical predictions, which are testable using imaging experiments, do not contain any parameters.

Unveiling the Machinery behind Chromosome Folding by Polymer Physics Modeling.

Understanding the mechanisms underlying the complex 3D architecture of mammalian genomes poses, at a more fundamental level, the problem of how two or multiple genomic sites can establish physical contacts in the nucleus of the cells. Beyond stochastic and fleeting encounters related to the polymeric nature of chromatin, experiments have revealed specific, privileged patterns of interactions that suggest the existence of basic organizing principles of folding. In this review, we focus on two major and recently proposed physical processes of chromatin organization: loop-extrusion and polymer phase-separation, both supported by increasing experimental evidence. We discuss their implementation into polymer physics models, which we test against available single-cell super-resolution imaging data, showing that both mechanisms can cooperate to shape chromatin structure at the single-molecule level. Next, by exploiting the comprehension of the underlying molecular mechanisms, we illustrate how such polymer models can be used as powerful tools to make predictions in silico that can complement experiments in understanding genome folding. To this aim, we focus on recent key applications, such as the prediction of chromatin structure rearrangements upon disease-associated mutations and the identification of the putative chromatin organizing factors that orchestrate the specificity of DNA regulatory contacts genome-wide.

Consistencies and contradictions in different polymer models of chromatin architecture.

Genetic information is stored in very long DNA molecules, which are folded to form chromatin, a similarly long polymer fibre that is ultimately organised into chromosomes. The organisation of chromatin is fundamental to many cellular functions, from the expression of the genetic information to cell division. As a long polymer, chromatin is very flexible and may adopt a myriad of shapes. Globally, the polymer physics governing chromatin dynamics is very well understood. But chromatin is not uniform and regions of it, with chemical modifications and bound effectors, form domains and compartments through mechanisms not yet clear. Polymer models have been successfully used to investigate these mechanisms to explain cytological observations and build hypothesis for experimental validation. Many different approaches to conceptualise chromatin in polymer models can be envisioned and each reflects different aspects. Here, we compare recent approaches that aim at reproducing prominent features of interphase chromatin organisation: the compartmentalisation into eu- and heterochromatin compartments, the formation of a nucleolus, chromatin loops and the rosette and Rabl conformations of interphase chromosomes. We highlight commonalities and contradictions that point to a modulation of the mechanisms involved to fine degree. Consolidating models will require the inclusion of yet hidden or neglected parameters.

Topological Considerations in Biomolecular Condensation.

Biomolecular condensation and phase separation are increasingly understood to play crucial roles in cellular compartmentalization and spatiotemporal regulation of cell machinery implicated in function and pathology. A key aspect of current research is to gain insight into the underlying physical mechanisms of these processes. Accordingly, concepts of soft matter and polymer physics, the thermodynamics of mixing, and material science have been utilized for understanding condensation mechanisms of multivalent macromolecules resulting in viscoelastic mesoscopic supramolecular assemblies. Here, we focus on two topological concepts that have recently been providing key mechanistic understanding in the field. First, we will discuss how percolation provides a network-topology-related framework that offers an interesting paradigm to understand the complex networking of dense 'connected' condensate structures and, therefore, their phase behavior. Second, we will discuss the idea of entanglement as another topological concept that has deep roots in polymer physics and important implications for biomolecular condensates. We will first review some historical developments and fundamentals of these concepts, then we will discuss current advancements and recent examples. Our discussion ends with a few open questions and the challenges to address them, hinting at unveiling fresh possibilities for the modification of existing knowledge as well as the development of new concepts relevant to condensate science.

Nucleic Acid Structure Prediction Including Pseudoknots Through Direct Enumeration of States: A User's Guide to the LandscapeFold Algorithm.

Here we detail the LandscapeFold secondary structure prediction algorithm and how it is used. The algorithm was previously described and tested in (Kimchi O et al., Biophys J 117(3):520-532, 2019), though it was not named there. The algorithm directly enumerates all possible secondary structures into which up to two RNA or single-stranded DNA sequences can fold. It uses a polymer physics model to estimate the configurational entropy of structures including complex pseudoknots. We detail each of these steps and ways in which the user can adjust the algorithm as desired. The code is available on the GitHub repository https://github.com/ofer-kimchi/LandscapeFold .

Molecular Dynamics Simulation of Entangled Melts at High Rates: Identifying Entanglement Lockup Mechanism Leading to True Strain Hardening.

In the present work, molecular dynamics simulations are carried out based on the bead-spring model to indicate how the entanglement lockup manifests in the late stage of fast Rouse-Weissnberg number (WiR >>1) uniaxial melt stretching of entangled polymer melts. At high strains, distinct features show up to reveal the emergence of an increasingly tightened entanglement network. Chain tension can build up, peaking at the middle of the chain, to a level for chain scission, through accumulated interchain interactions, as if there is a tug-of-war ongoing for each load-bearing chain. Thanks to the interchain uncrossability, network junctions form by the pairing of two or more hairpins. It is hypothesized that the interchain entanglement at junctions can lockup through prevailing twist-like interchain couplings as long as WiR > 9. In this limit, a significant fraction of chains act like cyclic chains to form a network held by interchain uncrossability, and appreciable chain tension emerges.

How enzymatic activity is involved in chromatin organization.

Spatial organization of chromatin plays a critical role in genome regulation. Previously, various types of affinity mediators and enzymes have been attributed to regulate spatial organization of chromatin from a thermodynamics perspective. However, at the mechanistic level, enzymes act in their unique ways and perturb the chromatin. Here, we construct a polymer physics model following the mechanistic scheme of Topoisomerase-II, an enzyme resolving topological constraints of chromatin, and investigate how it affects interphase chromatin organization. Our computer simulations demonstrate Topoisomerase-II's ability to phase separate chromatin into eu- and heterochromatic regions with a characteristic wall-like organization of the euchromatic regions. We realized that the ability of the euchromatic regions to cross each other due to enzymatic activity of Topoisomerase-II induces this phase separation. This realization is based on the physical fact that partial absence of self-avoiding interaction can induce phase separation of a system into its self-avoiding and non-self-avoiding parts, which we reveal using a mean-field argument. Furthermore, motivated from recent experimental observations, we extend our model to a bidisperse setting and show that the characteristic features of the enzymatic activity-driven phase separation survive there. The existence of these robust characteristic features, even under the non-localized action of the enzyme, highlights the critical role of enzymatic activity in chromatin organization.