Modelling of crowded polymers elucidate effects of double-strand breaks in topological domains of bacterial chromosomes

Development of a Data-Driven Integrative Model of a Bacterial Chromosome

The chromosome of archetypal bacteria E. coli is known for a complex topology with a 4.6 × 106 base pairs (bp) long sequence of nucleotides packed within a micrometer-sized cellular confinement. The inherent organization underlying this chromosome eludes general consensus due to the lack of a high-resolution picture of its conformation. Here we present our development of an integrative model of E. coli at a 500 bp resolution (https://github.com/JMLab-tifrh/ecoli_finer), which optimally combines a set of multiresolution genome-wide experimentally measured data within a framework of polymer based architecture. In particular the model is informed with an intragenome contact probability map at 5000 bp resolution derived via the Hi-C experiment and RNA-sequencing data at 500 bp resolution. Via dynamical simulations, this data-driven polymer based model generates an appropriate conformational ensemble commensurate with chromosome architectures that E. coli adopts. As a key hallmark of the E. coli chromosome the model spontaneously self-organizes into a set of nonoverlapping macrodomains and suitably locates plectonemic loops near the cell membrane. As novel extensions, it predicts a contact probability map simulated at a higher resolution than precedent experiments and can demonstrate segregation of chromosomes in a partially replicating cell. Finally, the modular nature of the model helps us devise control simulations to quantify the individual role of key features in hierarchical organization of the bacterial chromosome.

Integrative modeling of JCVI-Syn3A nucleoids with a modular approach

A lattice-based method is presented for creating 3D models of entire bacterial nucleoids integrating ultrastructural information cryoelectron tomography, genomic and proteomic data, and experimental atomic structures of biomolecules and assemblies. The method is used to generate models of the minimal genome bacterium JCVI-Syn3A, producing a series of models that test hypotheses about transcription, condensation, and overall distribution of the genome. Lattice models are also used to generate atomic models of an entire JCVI-Syn3A cell.

Cyclic-polymer grafted colloids in spherical confinement: insights for interphase chromosome organization

Interphase chromosomes are known to organize non-randomly in the micron-sized eukaryotic cell nucleus and occupy certain fraction of nuclear volume, often without mixing. Using extensive coarse-grained simulations, we model such chromosome structures as colloidal particles whose surfaces are grafted by cyclic polymers. This model system is known as Rosetta. The cyclic polymers, with varying polymerization degrees, mimic chromatin loops present in interphase chromosomes, while the rigid core models the chromocenter section of the chromosome. Our simulations show that the colloidal chromosome model provides a well-separated particle distribution without specific attraction between the chain monomers. As the polymerization degree of the grafted cyclic chains decreases while maintaining the total chromosomal length (e.g. the more potent activity of condensin-family proteins), the average chromosomal volume becomes smaller, inter-chromosomal contacts decrease, and chromocenters organize in a quasi-crystalline order reminiscent of a glassy state. This order weakens for polymer chains with a characteristic size on the order of the confinement radius. Notably, linear-polymer grafted particles also provide the same chromocenter organization scheme. However, unlike linear chains, cyclic chains result in less contact between the polymer layers of neighboring chromosome particles, demonstrating the effect of DNA breaks in altering genome-wide contacts. Our simulations show that polymer-grafted colloidal systems could help decipher 3D genome architecture along with the fractal globular and loop-extrusion models.

Equilibrium behaviour of two cavity-confined polymers: effects of polymer width and system asymmetries

Experiments using nanofluidic devices have proven effective in characterizing the physical properties of polymers confined to small cavities. Two recent studies using such methods examined the organization and dynamics of two DNA molecules in box-like cavities with strong confinement in one direction and with square and elliptical cross sections in the lateral plane. Motivated by these experiments, we employ Monte Carlo and Brownian dynamics simulations to study the physical behaviour of two polymers confined to small cavities with shapes comparable to those used in the experiments. We quantify the effects of varying the following polymer properties and confinement dimensions on the organization and dynamics of the polymers: the polymer width, the polymer contour length ratio, the cavity cross-sectional area, and the degree of cavity elongation for cavities with rectangular and elliptical cross sections. We find that the tendency for polymers to segregate is enhanced by increasing polymer width. For sufficiently small cavities, increasing cavity elongation promotes segregation and localization of identical polymers to opposite sides of the cavity along its long axis. A free-energy barrier controls the rate of polymers swapping positions, and the observed dynamics are roughly in accord with predictions of a simple theoretical model. Increasing the contour length difference between polymers significantly affects their organization in the cavity. In the case of a large linear polymer co-trapped with a small ring polymer in an elliptical cavity, the small polymer tends to lie near the lateral confining walls, and especially at the cavity poles for highly elongated ellipses.

Knot Factories with Helical Geometry Enhance Knotting and Induce Handedness to Knots

We performed molecular dynamics simulations of DNA polymer chains confined in helical nano-channels under compression in order to explore the potential of knot-factories with helical geometry to produce knots with a preferred handedness. In our simulations, we explore mutual effect of the confinement strength and compressive forces in a range covering weak, intermediate and strong confinement together with weak and strong compressive forces. The results find that while the common metrics of polymer chain in cylindrical and helical channels are very similar, the DNA in helical channels exhibits greatly different topology in terms of chain knottedness, writhe and handedness of knots. The results show that knots with a preferred chirality in terms of average writhe can be produced by using channels with a chosen handedness.

Facilitated Dissociation of Nucleoid Associated Proteins from DNA in the Bacterial Confinement

Transcription machinery depends on the temporal formation of protein-DNA complexes. Recent experiments demonstrated that not only the formation but also the lifetime of such complexes can affect the transcriptional machinery. In parallel, in vitro single-molecule studies showed that nucleoid-associated proteins (NAPs) leave the DNA rapidly as the bulk concentration of the protein increases via facilitated dissociation (FD). Nevertheless, whether such a concentration-dependent mechanism is functional in a bacterial cell, in which NAP levels and the 3d chromosomal structure are often coupled, is not clear a priori. Here, by using extensive coarse-grained molecular simulations, we model the unbinding of specific and nonspecific dimeric NAPs from a high-molecular-weight circular DNA molecule in a cylindrical structure mimicking the cellular confinement of a bacterial chromosome. Our simulations confirm that physiologically relevant peak protein levels (tens of micromolar) lead to highly compact chromosomal structures. This compaction results in rapid off rates (shorter DNA residence times) for specifically DNA-binding NAPs, such as the factor for inversion stimulation, which mostly dissociate via a segmental jump mechanism. Contrarily, for nonspecific NAPs, which are more prone to leave their binding sites via 1d sliding, the off rates decrease as the protein levels increase. The simulations with restrained chromosome models reveal that chromosome compaction is in favor of faster dissociation but only for specific proteins, and nonspecific proteins are not affected by the chromosome compaction. Overall, our results suggest that the cellular concentration level of a structural DNA-binding protein can be highly intermingled with its DNA residence time.

Equilibrium organization, conformation, and dynamics of two polymers under box-like confinement

Motivated by recent nanofluidics experiments, we use Brownian dynamics and Monte Carlo simulations to study the conformation, organization and dynamics of two polymer chains confined to a single box-like cavity. The polymers are modeled as flexible bead-spring chains, and the box has a square cross-section of side length L and a height that is small enough to compress the polymers in that dimension. For sufficiently large L, the system behaviour approaches that of an isolated polymer in a slit. However, the combined effects of crowding and confinement on the polymer organization, conformation and equilibrium dynamics become significant when where is the transverse radius of gyration for a slit geometry. In this regime, the centre-of-mass probability distribution in the transverse plane exhibits a depletion zone near the centre of the cavity (except at very small L) and a 4-fold symmetry with quasi-discrete positions. Reduction in polymer size with decreasing L arises principally from confinement rather than inter-polymer crowding. By contrast, polymer diffusion and internal motion are strongly affected by inter-polymer crowding. The two polymers tend to occupy opposite positions relative to the box centre, about which they diffuse relatively freely. Qualitatively, this static and dynamical behaviour differs significantly from that previously observed for confinement of two polymers to a narrow channel. The simulation results for a suitably chosen box width are qualitatively consistent with results from a recent experimental study of two λ-DNA chains confined to a nanofluidic cavity.

A Hi-C data-integrated model elucidates E. coli chromosome's multiscale organization at various replication stages

The chromosome of Escherichia coli is riddled with multi-faceted complexity. The emergence of chromosome conformation capture techniques are providing newer ways to explore chromosome organization. Here we combine a beads-on-a-spring polymer-based framework with recently reported Hi-C data for E. coli chromosome, in rich growth condition, to develop a comprehensive model of its chromosome at 5 kb resolution. The investigation focuses on a range of diverse chromosome architectures of E. coli at various replication states corresponding to a collection of cells, individually present in different stages of cell cycle. The Hi-C data-integrated model captures the self-organization of E. coli chromosome into multiple macrodomains within a ring-like architecture. The model demonstrates that the position of oriC is dependent on architecture and replication state of chromosomes. The distance profiles extracted from the model reconcile fluorescence microscopy and DNA-recombination assay experiments. Investigations into writhe of the chromosome model reveal that it adopts helix-like conformation with no net chirality, earlier hypothesized in experiments. A genome-wide radius of gyration map captures multiple chromosomal interaction domains and identifies the precise locations of rrn operons in the chromosome. We show that a model devoid of Hi-C encoded information would fail to recapitulate most genomic features unique to E. coli.

Free energy and segregation dynamics of two channel-confined polymers of different lengths

Polymers confined to a narrow channel are subject to strong entropic forces that tend to drive the molecules apart. In this study, we use Monte Carlo computer simulations to study the segregation behavior of two flexible hard-sphere polymers under confinement in a cylindrical channel. We focus on the effects of using polymers of different lengths. We measure the variation of the conformational free energy, F, with the center-of-mass separation distance, λ. The simulations reveal four different separation regimes, characterized by different scaling properties of the free energy with respect to the polymer lengths and the channel diameter, D. We propose a regime map in which the state of the system is determined by the values of the quantities N_{2}/N_{1} and λ/(N_{1}+N_{2})D^{-β}, where N_{1} and N_{2} are the polymer lengths, and where β≈0.64. The observed scaling behavior of F(λ) in each regime is in reasonable agreement with predictions using a simple theoretical model. In addition, we use MC dynamics simulations to study the segregation dynamics of initially overlapping polymers by measurement of the incremental mean first-passage time with respect to λ. For systems characterized by a wide range of λ in which a short polymer is nested within a longer one, the segregation dynamics are close to that expected for two noninteracting one-dimensional random walkers undergoing unbiased diffusion. When the free-energy gradient is large, segregation is rapid and characterized by out-of-equilibrium effects.

Segregation of polymers under cylindrical confinement: effects of polymer topology and crowding

Monte Carlo computer simulations are used to study the segregation behaviour of two polymers under cylindrical confinement. Using a multiple-histogram method, the conformational free energy, F, of the polymers was measured as a function of the centre-of-mass separation distance, λ. We examined the scaling of the free energy functions with the polymer length, the length and diameter of the confining cylinder, the polymer topology (i.e. linear vs. ring polymers), and the packing fraction and size of mobile crowding agents. In the absence of crowders, the observed scaling of F(λ) is similar to that predicted using a simple model employing the de Gennes blob model and the approximation that the free energy of overlapping chains in a tube is equal to that of two isolated chains each in a tube of half the cross-sectional area. Simulations were used to test the latter approximation and reveal that it yields poor quantitative predictions. This, along with generic finite-size effects, likely gives rise to the discrepancies between the predicted and measured values of scaling exponents for F(λ). For segregation in the presence of crowding agents, the free energy barrier generally decreases with increasing crowder packing fraction, thus reducing the entropic forces driving segregation. However, for fixed packing fraction, the barrier increases as the crowder/monomer size ratio decreases.

Consistent rationalization of type-2 topoisomerases' unknotting, decatenating, supercoil-relaxing actions and their scaling relation

How type-2 topoisomerases discern global topology from local properties of DNA is not known precisely but the hypothesis that the enzymes selectively pass double-helix strands at hook-like juxtapositions is promising. Building upon an investigation of unknotting and decatenating using an improved wormlike DNA model, here we focus primarily on the enzymes' action in narrowing the distribution of linking number (Lk) in supercoiled DNA. Consistent with experiments, with selective passage at a hooked juxtaposition, the simulated narrowing factor RLk diminishes with decreasing DNA circle size but approaches an asymptotic RLk ≈ 1.7-1.8 for circle size ≳3.5 kb. For the larger DNA circles, we found that (RLk - 1) ≈ 0.42log10RK ≈ 0.68log10RL and thus RK ≈ (RL)(1.6) holds for the computed RLk and knot and catenane reduction factors RK and RL attained by selective passage at different juxtaposition geometries. Remarkably, this general scaling relation is essentially identical to that observed experimentally for several type-2 topoisomerases from a variety of organisms, indicating that the different disentangling powers of the topoisomerases likely arise from variations in the hooked geometries they select. Taken together, our results suggest strongly that type-2 topoisomerases recognize not only the curvature of the G-segment but also that of the T-segment.

Friction between ring polymer brushes

Friction between ring polymer brush bilayers sliding past each other at melt densities is studied using extensive coarse-grained molecular dynamics simulations and scaling arguments, and the results are compared to the friction between bilayers of linear polymer brushes. We show that for a velocity range spanning over three decades, the frictional forces measured for ring polymer brushes are half of the corresponding friction in the case of linear brushes. In the linear-force regime, the weak inter-digitation between ring brush layers as compared to linear brushes leads also to a lower number of binary collisions between the monomers from opposing brushes. At high velocities, where the thickness of the inter-digitation between bilayers is on the order of monomer size regardless of brush topology, stretched segments of ring polymers adopt the double-stranded conformation. As a result, monomers of the double-stranded segments collide on average less with the monomers of the opposing ring brush even though a similar number of monomers occupies the inter-digitation layer for ring and linear brush bilayers. The numerical data obtained from our simulations are consistent with the proposed scaling analysis. Conformation-dependent friction reduction observed in ring brushes can have important consequences in non-equilibrium bulk systems.

Effects of physiological self-crowding of DNA on shape and biological properties of DNA molecules with various levels of supercoiling

DNA in bacterial chromosomes and bacterial plasmids is supercoiled. DNA supercoiling is essential for DNA replication and gene regulation. However, the density of supercoiling in vivo is circa twice smaller than in deproteinized DNA molecules isolated from bacteria. What are then the specific advantages of reduced supercoiling density that is maintained in vivo? Using Brownian dynamics simulations and atomic force microscopy we show here that thanks to physiological DNA–DNA crowding DNA molecules with reduced supercoiling density are still sufficiently supercoiled to stimulate interaction between cis-regulatory elements. On the other hand, weak supercoiling permits DNA molecules to modulate their overall shape in response to physiological changes in DNA crowding. This plasticity of DNA shapes may have regulatory role and be important for the postreplicative spontaneous segregation of bacterial chromosomes.

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.