Self-avoiding polymer trapped inside a cylindrical pore: Flory free energy and unexpected dynamics

Polymer architecture orchestrates the segregation and spatial organization of replicating E. coli chromosomes in slow growth

The mechanism of chromosome segregation and organization in the bacterial cell cycle of E. coli is one of the least understood aspects in its life cycle. The E. coli chromosome is often modelled as a bead spring ring polymer. We introduce cross-links in the DNA-ring polymer, resulting in the formation of loops within each replicating bacterial chromosome. We use simulations to show that the chosen polymer-topology ensures its self-organization along the cell long-axis, such that various chromosomal loci get spatially localized as seen in vivo. The localization of loci arises due to entropic repulsion between polymer loops within each daughter DNA confined in a cylinder. The cellular addresses of the loci in our model are in fair agreement with those seen in experiments as given in J. A. Cass et al., Biophys. J., 2016, 110, 2597-2609. We also show that the adoption of such modified polymer architectures by the daughter DNAs leads to an enhanced propensity of their spatial segregation. Secondly, we match other experimentally reported results, including observation of the cohesion time and the ter-transition. Additionally, the contact map generated from our simulations reproduces the macro-domain like organization as seen in the experimentally obtained Hi-C map. Lastly, we have also proposed a plausible reconciliation of the 'Train Track' and the 'Replication Factory' models which provide conflicting descriptions of the spatial organization of the replication forks. Thus, we reconcile observations from complementary experimental techniques probing bacterial chromosome organization.

Topological Friction and Relaxation Dynamics of Spatially Confined Catenated Polymers

We study catenated ring polymers confined inside channels and slits with Langevin dynamics simulations and address how the contour position and size of the interlocked or physically linked region evolve with time. We show that the catenation constraints generate a drag, or topological friction, that couples the contour motion of the interlocked regions. Notably, the coupling strength decreases as the interlocking is made tighter, but also shorter, by confinement. Though the coupling strength differs for channel and slit confinement, the data outline a single universal curve when plotted against the size of the linked region. Finally, we study how the relaxation kinetics changes after one of the rings is cut open and conclude that considering interlocked circular polymers is key for isolating the manifestations of topological friction. The results ought to be relevant for linked biomolecules in experimental or biological confining conditions.

Compression and Stretching of Single DNA Molecules under Channel Confinement

We study the compression and extension response of single dsDNA (double-stranded DNA) molecules confined in cylindrical channels by means of Monte Carlo simulations. The elastic response of micrometer-sized DNA to the external force acting through the chain ends or through the piston is markedly affected by the size of the channel. The interpretation of the force (f)-displacement (R) functions under quasi-one-dimensional confinement is facilitated by resolving the overall change of displacement ΔR into the confinement contribution ΔR_D and the force contribution ΔR_f. The external stretching of confined DNA results in a characteristic pattern of f-R functions involving their shift to the larger extensions due to the channel-induced pre-stretching ΔR_D. A smooth end-chain compression into loop-like conformations observed in moderately confined DNA can be accounted for by the relationship valid for a Gaussian chain in bulk. In narrow channels, the considerably pre-stretched DNA molecules abruptly buckle on compression by the backfolding into hairpins. On the contrary, the piston compression of DNA is characterized by a gradual reduction of the chain span S and by smooth f-S functions in the whole spatial range from the 3d near to 1d limits. The observed discrepancy between the shape of the f-R and f-S functions from two compression methods can be important for designing nanopiston experiments of compaction and knotting of single DNA in nanochannels.

Limitations of the equivalent neutral polymer assumption for theories describing nanochannel-confined DNA

The prevailing theories describing DNA confinement in a nanochannel are predicated on the assumption that wall-DNA electrostatic interactions are sufficiently short-ranged such that the problem can be mapped to an equivalent neutral polymer confined by hard walls with an appropriately reduced effective channel size. To determine when this hypothesis is valid, we leveraged a recently reported experimental data set for the fractional extension of DNA molecules in a 250-nm-wide poly(dimethyl siloxane) (PDMS) nanochannel with buffer ionic strengths between 0.075 and 48 mM. Evaluating these data in the context of the weakly correlated telegraph model of DNA confinement reveals that, at ionic strengths greater than 0.3 mM, the average fractional extension of the DNA molecules agree with theoretical predictions with a mean absolute error of 0.04. In contrast, experiments at ionic strengths below 0.3 mM produce average fractional extensions that are systematically smaller than the theoretical predictions with a larger mean absolute error of 0.15. The deviations between experiment and theory display a correlation coefficient of 0.82 with the decay length for the DNA-wall electrostatics, linking the deviations with a breakdown in approximating the DNA with an equivalent neutral polymer.

Bacterial chromosome organization. II. Few special cross-links, cell confinement, and molecular crowders play the pivotal roles

Using a coarse-grained bead-spring model of bacterial chromosomes of Caulobacter crescentus and Escherichia coli, we show that just 33 and 38 effective cross-links in 4017 and 4642 monomer chains at special positions along the chain contour can lead to the large-scale organization of the DNA polymer, where confinement effects of the cell walls play a key role in the organization. The positions of the 33/38 cross-links along the chain contour are chosen from the Hi-C contact map of bacteria C. crescentus and E. coli. We represent 1000 base pairs as a coarse-grained monomer in our bead-spring flexible ring polymer model of the DNA polymer. Thus, 4017/4642 beads on a flexible ring polymer represent the C. crescentus/E. coli DNA polymer with 4017/4642 kilo-base pairs. Choosing suitable parameters from Paper I, we also incorporate the role of compaction of the polymer coil due to the presence of molecular crowders and the ability of the chain to release topological constraints. We validate our prediction of the organization of the bacterial chromosomes with available experimental data and also give a prediction of the approximate positions of different segments within the cell. In the absence of confinement, the minimal number of effective cross-links required to organize the DNA chains of 4017/4642 monomers was 60/82 [Agarwal et al., Europhys. Lett. 121, 18004 (2018) and Agarwal et al., J. Phys.: Condens. Matter 30, 034003 (2018)].

Relaxation of surface-tethered polymers under moderate confinement

We study the relaxation of surface-tethered polymers in microchannels under moderate confinement (i.e. h ∼ Rg, where h is the channel height and Rg is the radius of gyration of the polymer) by experiments with fluorescence-marked DNA molecules and coupled lattice-Boltzmann/molecular dynamics simulations. The determined scaling exponent suggests that the relaxation is dominated by Zimm-dynamics with significant intra-chain hydrodynamic interactions. The relaxation of the DNA molecules is slower in shallower channels, indicating a pronounced effect of confinement on the longest relaxation time. An experimental correlation is obtained for the longest relaxation time as a function of the molecular contour length and the channel height. Good agreement between the experimental and the simulation results is found.

Stretching and compression of DNA by external forces under nanochannel confinement

Mechanical deformation of dsDNA molecules inside square nanochannels is investigated using simulations based on a coarse-grained model of DNA. The combined action of confinement and weak external forces is explored in a variety of confinement regimes, including the transition zone relevant to nanofluidic experiments. The computed free energy and force profiles are markedly affected by the channel size. Effective elastic softening of confined DNA molecules relative to the bulk DNA is observed in the channels of intermediate widths. The extension of DNA from its bulk equilibrium length in nanofluidic devices is resolved into contributions from the passive extension due to confinement and from the active stretching induced by force. Potential implications of the very different energy costs computed for the two extension modes (extension by confinement takes much more free energy than stretching by force) for behavior of DNA in nanofluidic chips are indicated.

Confinement-driven organization of a histone-complexed DNA molecule in a dense array of nanoposts

The first step in the controlled storage of lengthy DNA molecules is to keep DNA molecules separated while integrated in micrometer-sized space. Herein, we present hybrid Monte Carlo simulations of a histone-complexed DNA (hcDNA) molecule confined in a dense array of nanoposts. Depending on the nanopost dimension, a single, 8.7 kilobase pair hcDNA molecule was either localized and elongated in a single inter-post space surrounded by four nanoposts or spread over several inter-post spaces through passages between two neighboring nanoposts. The conformational change of a hcDNA molecule is interpreted in terms of competitive effects of confinements in the inter-post and passage spaces. We propose that, by elaborately designing nanopost arrays, the competitive confinement effects can be adjusted such that each hcDNA molecule is localized in a single inter-post space, and thereby multiple hcDNA molecules can be physically separated from each other while stored together in the nanopost array.

Effects of molecular crowding and confinement on the spatial organization of a biopolymer

A chain molecule can be entropically collapsed in a crowded medium in a free or confined space. Here, we present a unified view of how molecular crowding collapses a flexible polymer in three distinct spaces: free, cylindrical, and (two-dimensional) slit-like. Despite their seeming disparities, a few general features characterize all these cases, even though the ϕ_c-dependence of chain compaction differs between the two cases: a > a_c and a < a_c, where ϕ_c is the volume fraction of crowders, a is the monomer size, and a_c is the crowder size. For a > a_c (applicable to a coarse-grained model of bacterial chromosomes), chain size depends on the ratio aϕ_c/a_c, and "full" compaction occurs universally at aϕ_c/a_c ≈ 1; for a_c > a (relevant for protein folding), it is controlled by ϕ_c alone and crowding has a modest effect on chain size in a cellular environment (ϕ_c ≈ 0.3). Also for a typical parameter range of biological relevance, molecular crowding can be viewed as effectively reducing the solvent quality, independent of confinement.

Semiflexible macromolecules in quasi-one-dimensional confinement: Discrete versus continuous bond angles

The conformations of semiflexible polymers in two dimensions confined in a strip of width D are studied by computer simulations, investigating two different models for the mechanism by which chain stiffness is realized. One model (studied by molecular dynamics) is a bead-spring model in the continuum, where stiffness is controlled by a bond angle potential allowing for arbitrary bond angles. The other model (studied by Monte Carlo) is a self-avoiding walk chain on the square lattice, where only discrete bond angles (0° and ±90°) are possible, and the bond angle potential then controls the density of kinks along the chain contour. The first model is a crude description of DNA-like biopolymers, while the second model (roughly) describes synthetic polymers like alkane chains. It is first demonstrated that in the bulk the crossover from rods to self-avoiding walks for both models is very similar, when one studies average chain linear dimensions, transverse fluctuations, etc., despite their differences in local conformations. However, in quasi-one-dimensional confinement two significant differences between both models occur: (i) The persistence length (extracted from the average cosine of the bond angle) gets renormalized for the lattice model when D gets less than the bulk persistence length, while in the continuum model it stays unchanged. (ii) The monomer density near the repulsive walls for semiflexible polymers is compatible with a power law predicted for the Kratky-Porod model in the case of the bead-spring model, while for the lattice case it tends to a nonzero constant across the strip. However, for the density of chain ends, such a constant behavior seems to occur for both models, unlike the power law observed for flexible polymers. In the regime where the bulk persistence length ℓp is comparable to D, hairpin conformations are detected, and the chain linear dimensions are discussed in terms of a crossover from the Daoud/De Gennes "string of blobs"-picture to the flexible rod picture when D decreases and/or the chain stiffness increases. Introducing a suitable further coarse-graining of the chain contours of the continuum model, direct estimates for the deflection length and its distribution could be obtained.

Confinement and partitioning of a single polymer chain in a dense array of nanoposts

We present a Brownian dynamics simulation study on the confinement and partitioning of a single, flexible polymer chain in a dense array of nanoposts with different sizes and separations, especially, when the volume of an interstitial space formed among four nanoposts is less than the volume of the polymer chain. As the interstitial volume decreases by either increasing the nanopost diameter or decreasing the separation between nanoposts, the chain conformation becomes elongated in the direction parallel to the nanoposts. Interestingly, however, the degree of chain elongation varies in a non-monotonic fashion as the interstitial volume decreases while keeping the passage width between two nanoposts constant at a small value. We calculate the free energy of chain partitioning over several interstitial spaces from the partitioning probability, and find that the non-monotonic dependence of the chain elongation results from an interplay between the confinement-driven chain elongation along the direction parallel to the nanoposts and the chain spreading perpendicular to the nanoposts by partitioning chain segments over several interstitial spaces. These results present the possibility of utilizing a dense array of nanoposts as a template to control polymer conformations.

Polymers under confinement: single polymers, how they interact, and as model chromosomes

How confinement or a physical constraint modifies polymer chains is not only a classical problem in polymer physics but also relevant in a variety of contexts such as single-molecule manipulations, nanofabrication in narrow pores, and modelling of chromosome organization. Here, we review recent progress in our understanding of polymers in a confined (and crowded) space. To this end, we highlight converging views of these systems from computational, experimental, and theoretical approaches, and then clarify what remains to be clarified. In particular, we focus on exploring how cylindrical confinement reshapes individual chains and induces segregation forces between them - by pointing to the relationships between intra-chain organization and chain segregation. In the presence of crowders, chain molecules can be entropically phase-separated into a condensed state. We include a kernel of discussions on the nature of chain compaction by crowders, especially in a confined space. Finally, we discuss the relevance of confined polymers for the nucleoid, an intracellular space in which the bacterial chromosome is tightly packed, in part by cytoplasmic crowders.

A polymer in a crowded and confined space: effects of crowder size and poly-dispersity

DNA compaction in a bacterial cell is in part carried out by entropic (depletion) forces induced by "free" proteins or crowding particles in the cytoplasm. Indeed, recent in vitro experiments highlight these effects by showing that they alone can condense the E. coli chromosome to its in vivo size. Using molecular dynamics simulations and a theoretical approach, we study how a flexible chain molecule can be compacted by crowding particles with variable sizes in a (cell-like) cylindrical space. Our results show that with smaller crowding agents the compaction occurs at a lower volume fraction but at a larger concentration such that doubling their size is equivalent to increasing their concentration fourfold. Similarly, the effect of polydispersity can be correctly mimicked by adjusting the size of crowders in a homogeneous system. Under different conditions, however, crowding particles can induce chain adsorption onto the cylinder wall, stretching the chain, which would otherwise remain condensed.

Accuracy of the blob model for single flexible polymers inside nanoslits that are a few monomer sizes wide

The de Gennes' blob model is extensively used in different problems of polymer physics. This model is theoretically applicable when the number of monomers inside each blob is large enough. For confined flexible polymers, this requires the confining geometry to be much larger than the monomer size. In this paper, the opposite limit of polymer in nanoslits with one to several monomers width is studied, using molecular dynamics simulations. Extension of the polymer inside nanoslits, confinement force on the plates, and the effective spring constant of the confined polymer are investigated. Despite the theoretical limitations of the blob model, the simulation results are explained with the blob model very well. The agreement is observed for the static properties and the dynamic spring constant of the polymer. A theoretical description of the conditions under which the dynamic spring constant of the polymer is independent of the small number of monomers inside blobs is given. Our results on the limit of applicability of the blob model can be useful in the design of nanotechnology devices.

Localization and stretching of polymer chains at the junction of two surfaces

We present a molecular dynamics study on the stretching of a linear polymer chain that is adsorbed at the junction of two intersecting flat surfaces of varying alignments. We observe a transition from a two-dimensional to one-dimensional (1D) structure of the adsorbed polymer when the alignment, i.e., the angle between the two surfaces that form a groove, θ, is below 135°. We show that the radius of gyration of the polymer chain Rg scales as Rg ∼ N(3/4) with the degree of polymerization N for θ = 180° (planer substrate), and the scaling changes to Rg ∼ N(1.0) for θ < 135° in good solvents. At the crossover point, θ = 135°, the exponent becomes 1.15. The 1D stretching of the polymer chain is found to be 84% of its contour length for θ ⩽ 90°. The center of mass diffusion coefficient D decreases sharply with θ. However, the diffusion coefficient scales with N as D ∼ N(-1), and is independent of θ. The relaxation time τ, for the diffusive motion, scales as τ ∼ N(2.5) for θ = 180° (planar substrate), which changes to τ ∼ N(3.0) for θ ⩽ 90°. At the crossover point, the exponent is 3.4, which is slightly higher than the 1D value of 3.0. Further, a signature of reptation-like dynamics of the polymer chain is observed at the junction for θ ⩽ 90° due to its strong 1D localization and stretching.

The multifork Escherichia coli chromosome is a self-duplicating and self-segregating thermodynamic ring polymer

At all but the slowest growth rates, Escherichia coli cell cycles overlap, and its nucleoid is segregated to daughter cells as a forked DNA circle with replication ongoing-a state fundamentally different from eukaryotes. We have solved the chromosome organization, structural dynamics, and segregation of this constantly replicating chromosome. It is locally condensed to form a branched donut, compressed so that the least replicated DNA spans the cell center and the newest DNA extends toward the cell poles. Three narrow zones at the cell center and quarters contain both the replication forks and nascent DNA and serve to segregate the duplicated chromosomal information as it flows outward. The overall pattern is smoothly self-replicating, except when the duplicated terminus region is released from the septum and recoils to the center of a sister nucleoid. In circular cross-section of the cell, the left and right arms of the chromosome form separate, parallel structures that lie in each cell half along the radial cell axis. In contrast, replication forks and origin and terminus regions are found mostly at the center of the cross section, balanced by the parallel chromosome arms. The structure is consistent with the model in which the nucleoid is a constrained ring polymer that develops by spontaneous thermodynamics. The ring polymer pattern extrapolates to higher growth rates and also provides a structural basis for the form of the chromosome during very slow growth.

DNA confinement in nanochannels: physics and biological applications

DNA is the central storage molecule of genetic information in the cell, and reading that information is a central problem in biology. While sequencing technology has made enormous advances over the past decade, there is growing interest in platforms that can readout genetic information directly from long single DNA molecules, with the ultimate goal of single-cell, single-genome analysis. Such a capability would obviate the need for ensemble averaging over heterogeneous cellular populations and eliminate uncertainties introduced by cloning and molecular amplification steps (thus enabling direct assessment of the genome in its native state). In this review, we will discuss how the information contained in genomic-length single DNA molecules can be accessed via physical confinement in nanochannels. Due to self-avoidance interactions, DNA molecules will stretch out when confined in nanochannels, creating a linear unscrolling of the genome along the channel for analysis. We will first review the fundamental physics of DNA nanochannel confinement--including the effect of varying ionic strength--and then discuss recent applications of these systems to genomic mapping. Apart from the intense biological interest in extracting linear sequence information from elongated DNA molecules, from a physics view these systems are fascinating as they enable probing of single-molecule conformation in environments with dimensions that intersect key physical length-scales in the 1 nm to 100 µm range.

Fluctuation modes of nanoconfined DNA

We report an experimental investigation of the magnitude of length and density fluctuations in DNA that has been stretched in nanofluidic channels. We find that the experimental data can be described using a one-dimensional overdamped oscillator chain with nonzero equilibrium spring length and that a chain of discrete oscillators yields a better description than a continuous chain. We speculate that the scale of these discrete oscillators coincides with the scale at which the finite extensibility of the polymer manifests itself. We discuss how the measurement process influences the apparent measured dynamic properties, and outline requirements for the recovery of true physical quantities.

Simulation of DNA Extension in Nanochannels

We have used a realistic model for double stranded DNA and Monte Carlo simulations to compute the extension (mean span) of a DNA molecule confined in a nanochannel over the full range of confinement in a high ionic strength buffer. The simulation data for square nanochannels resolve the apparent contradiction between prior simulation studies and the predictions from Flory theory, demonstrating the existence of two transition regimes between weak confinement (the de Gennes regime) and strong confinement (the Odijk regime). The simulation data for rectangular nanochannels support the use of the geometric mean for mapping data obtained in rectangular channels onto models developed for cylinders. The comparison of our results with experimental data illuminates the challenges in applying models for confined, neutral polymers to polyelectrolytes. Using a Flory-type approach, we also provide an improved scaling result for the relaxation time in the transition regime close to that found in experiments.

Entropically driven motion of polymers in nonuniform nanochannels

In nanofluidic devices, nonuniform confinement induces an entropic force that automatically drives biopolymers toward less-confined regions to gain entropy. To understand this phenomenon, we first analyze the diffusion of an entropy-driven particle system. The derived Fokker-Planck equation reveals an effective driving force as the negative gradient of the free energy. The derivation also shows that both the diffusion constant and drag coefficient are location dependent on an arbitrary free-energy landscape. As an application, DNA motion and deformation in nonuniform channels are investigated. Typical solutions reveal large gradients of stress on the polymer where the channel width changes rapidly. Migration of DNA in several nonuniform channels is discussed.

Single-polymer dynamics under constraints: scaling theory and computer experiment

The relaxation, diffusion and translocation dynamics of single linear polymer chains in confinement is briefly reviewed with emphasis on the comparison between theoretical scaling predictions and observations from experiment or, most frequently, from computer simulations. Besides cylindrical, spherical and slit-like constraints, related problems such as the chain dynamics in a random medium and the translocation dynamics through a nanopore are also considered. Another particular kind of confinement is imposed by polymer adsorption on attractive surfaces or selective interfaces--a short overview of single-chain dynamics is also contained in this survey. While both theory and numerical experiments consider predominantly coarse-grained models of self-avoiding linear chain molecules with typically Rouse dynamics, we also note some recent studies which examine the impact of hydrodynamic interactions on polymer dynamics in confinement. In all of the aforementioned cases we focus mainly on the consequences of imposed geometric restrictions on single-chain dynamics and try to check our degree of understanding by assessing the agreement between theoretical predictions and observations.

Transition between two regimes describing internal fluctuation of DNA in a nanochannel

We measure the thermal fluctuation of the internal segments of a piece of DNA confined in a nanochannel about 50-100 nm wide. This local thermodynamic property is key to accurate measurement of distances in genomic analysis. For DNA in ~100 nm channels, we observe a critical length scale ~10 m for the mean extension of internal segments, below which the de Gennes' theory describes the fluctuations with no fitting parameters, and above which the fluctuation data falls into Odijk's deflection theory regime. By analyzing the probability distributions of the extensions of the internal segments, we infer that folded structures of length 150-250 nm, separated by ~10 m exist in the confined DNA during the transition between the two regimes. For ~50 nm channels we find that the fluctuation is significantly reduced since the Odijk regime appears earlier. This is critical for genomic analysis. We further propose a more detailed theory based on small fluctuations and incorporating the effects of confinement to explicitly calculate the statistical properties of the internal fluctuations. Our theory is applicable to polymers with heterogeneous mechanical properties confined in non-uniform channels. We show that existing theories for the end-to-end extension/fluctuation of polymers can be used to study the internal fluctuations only when the contour length of the polymer is many times larger than its persistence length. Finally, our results suggest that introducing nicks in the DNA will not change its fluctuation behavior when the nick density is below 1 nick per kbp DNA.

Size and shape of excluded volume polymers confined between parallel plates

A number of recent experiments have provided detailed observations of the configurations of long DNA strands under nano-to-micrometer-sized confinement. We therefore revisit the problem of an excluded volume polymer chain confined between two parallel plates with varying plate separation. We show that the nonmonotonic behavior of the overall size of the chain as a function of plate separation, seen in computer simulations and reproduced by earlier theories, can already be predicted on the basis of scaling arguments. However, the behavior of the size in a plane parallel to the plates, a quantity observed in recent experiments, is predicted to be monotonic, in contrast to the experimental findings. We analyze this problem in depth with a mean-field approach that maps the confined polymer onto an anisotropic Gaussian chain, which allows the size of the polymer to be determined separately in the confined and unconfined directions. The theory allows the analytical construction of a smooth crossover between the small-plate-separation de Gennes regime and the large-plate-separation Flory regime. The results show good agreement with molecular dynamics simulations in the presence of a Langevin heat bath and confirm the scaling predictions.

Overlapping two self-avoiding polymers in a closed cylindrical pore: Implications for chromosome segregation in a bacterial cell

We study the spatial organization and segregation of two self-avoiding polymers trapped inside a closed cylindrical pore. Using molecular-dynamics simulations, we show how confinement shapes the chains, especially their mutual (entropic) force, chain miscibility, and segregation dynamics. Under strong confinement, the chains are shown to repel more strongly and thus segregate better if they are shorter and the confining space is more asymmetric, in contrast to the spherically confined case, where nonlinear chain topology is required for chain partitioning in equilibrium. When applied to bacterial chromosomes, our results imply that chromosome miscibility depends on how they are compacted and structured inside the cell (by proteins and supercoiling). Finally, longitudinal confinement is shown to have nontrivial effects on segregation dynamics by randomizing and thus slowing down the segregation process, which would otherwise be assisted with entropic forces.