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

Physical models of bacterial chromosomes

The interplay between bacterial chromosome organization and functions such as transcription and replication can be studied in increasing detail using novel experimental techniques. Interpreting the resulting quantitative data, however, can be theoretically challenging. In this minireview, we discuss how connecting experimental observations to biophysical theory and modeling can give rise to new insights on bacterial chromosome organization. We consider three flavors of models of increasing complexity: simple polymer models that explore how physical constraints, such as confinement or plectoneme branching, can affect bacterial chromosome organization; bottom-up mechanistic models that connect these constraints to their underlying causes, for instance, chromosome compaction to macromolecular crowding, or supercoiling to transcription; and finally, data-driven methods for inferring interpretable and quantitative models directly from complex experimental data. Using recent examples, we discuss how biophysical models can both deepen our understanding of how bacterial chromosomes are structured and give rise to novel predictions about bacterial chromosome organization.

A Novel Fractional Brownian Dynamics Method for Simulating the Dynamics of Confined Bottle-Brush Polymers in Viscoelastic Solution

In crowded fluids, polymer segments can exhibit anomalous subdiffusion due to the viscoelasticity of the surrounding environment. Previous single-particle tracking experiments revealed that such anomalous diffusion in complex fluids (e.g., in bacterial cytoplasm) can be described by fractional Brownian motion (fBm). To investigate how the viscoelastic media affects the diffusive behaviors of polymer segments without resolving single crowders, we developed a novel fractional Brownian dynamics method to simulate the dynamics of polymers under confinement. In this work, instead of using Gaussian random numbers ("white Gaussian noise") to model the Brownian force as in the standard Brownian dynamics simulations, we introduce fractional Gaussian noise (fGn) in our homemade fractional Brownian dynamics simulation code to investigate the anomalous diffusion of polymer segments by using a simple "bottle-brush"-type polymer model. The experimental results of the velocity autocorrelation function and the exponent that characterizes the subdiffusion of the confined polymer segments can be reproduced by this simple polymer model in combination with fractional Gaussian noise (fGn), which mimics the viscoelastic media.

Scaling regimes for wormlike chains confined to cylindrical surfaces under tension

We compute the free energy of confinement [Formula: see text] for a wormlike chain (WLC), with persistence length [Formula: see text], that is confined to the surface of a cylinder of radius R under an external tension f using a mean field variational approach. For long chains, we analytically determine the behavior of the chain in a variety of regimes, which are demarcated by the interplay of [Formula: see text], the Odijk deflection length ([Formula: see text]), and the Pincus length ([Formula: see text], with [Formula: see text] being the thermal energy). The theory accurately reproduces the Odijk scaling for strongly confined chains at [Formula: see text], with [Formula: see text]. For moderate values of f, the Odijk scaling is discernible only when [Formula: see text] for strongly confined chains. Confinement does not significantly alter the scaling of the mean extension for sufficiently high tension. The theory is used to estimate unwrapping forces for DNA from nucleosomes.

Modeling the compaction of bacterial chromosomes by biomolecular crowding and the cross-linking protein H-NS

Cells orchestrate the action of various molecules toward organizing their chromosomes. Using a coarse-grained computational model, we study the compaction of bacterial chromosomes by the cross-linking protein H-NS and cellular crowders. In this work, H-NS, modeled as a mobile "binder," can bind to a chromosome-like polymer with a characteristic binding energy. The simulation results reported here clarify the relative role of biomolecular crowding and H-NS in condensing a bacterial chromosome in a quantitative manner. In particular, they shed light on the nature and degree of crowder and H-NS synergetics: while the presence of crowders enhances H-NS binding to a chromosome-like polymer, the presence of H-NS makes crowding effects more efficient, suggesting two-way synergetics in chain compaction. Also, the results show how crowding effects promote clustering of bound H-NS. For a sufficiently large concentration of H-NS, the cluster size increases with the volume fraction of crowders.

Helical insertion of polyphenylene chains into confined cylindrical slits composed of two carbon nanotubes

The helical insertion behavior of poly(para-phenylene) (PP) chains into confined cylindrical slits constructed by two carbon nanotubes (CNTs) with different diameters is studied by molecular dynamics simulations. The contribution of system energy and each energy component to helical self-assembly is discussed to further explain the conditions, driving force and mechanism. The width and length of the slit, the diameter of the outer tube and the temperature have a great impact on the helical insertion of PP chains. Two equations are proposed to confirm the diameter and the distances between the PP helix and the inner and outer walls of the given CNTs. The helical self-assembly of PP with different numbers of chains inserted into the slits is further studied. This study has a great benefit in understanding the conformational behavior of polymers, even biological macromolecules in confinements.

Engineering Tissue-Scale Properties with Synthetic Cells: Forging One from Many

In metazoans, living cells achieve capabilities beyond individual cell functionality by assembling into multicellular tissue structures. These higher-order structures represent dynamic, heterogeneous, and responsive systems that have evolved to regenerate and coordinate their actions over large distances. Recent advances in constructing micrometer-sized vesicles, or synthetic cells, now point to a future where construction of synthetic tissue can be pursued, a boon to pressing material needs in biomedical implants, drug delivery systems, adhesives, filters, and storage devices, among others. To fully realize the potential of synthetic tissue, inspiration has been and will continue to be drawn from new molecular findings on its natural counterpart. In this review, we describe advances in introducing tissue-scale features into synthetic cell assemblies. Beyond mere complexation, synthetic cells have been fashioned with a variety of natural and engineered molecular components that serve as initial steps toward morphological control and patterning, intercellular communication, replication, and responsiveness in synthetic tissue. Particular attention has been paid to the dynamics, spatial constraints, and mechanical strengths of interactions that drive the synthesis of this next-generation material, describing how multiple synthetic cells can act as one.

Dynamics and Proton Conduction of Heterogeneously Confined Imidazole in Porous Coordination Polymers

The nanoconfinement of proton carrier molecules may contribute to the lowing of their proton dissociation energy. However, the free proton transportation does not occur as easily as in liquid due to the restricted molecular motion from surface attraction. To resolve the puzzle, herein, imidazole is confined in the channels of porous coordination polymers with tunable geometries, and their electric/structural relaxations are quantified. Imidazole confined in a square-shape channels exhibits dynamics heterogeneity of core-shell-cylinder model. The core and shell layer possess faster and slower structural dynamics, respectively, when compared to the bulk imidazole. The dimensions and geometry of the nanochannels play an important role in both the shielding of the blocking effect from attractive surfaces and the frustration filling of the internal proton carrier molecules, ultimately contributing to the fast dynamics and enhanced proton conductivity.

Topology-driven spatial organization of ring polymers under confinement

Entropic repulsion between DNA ring polymers under confinement is a key mechanism governing the spatial segregation of bacterial DNA before cell division. Here we establish that "internal" loops within a modified-ring polymer architecture enhance entropic repulsion between two overlapping polymers confined in a cylinder. Interestingly, they also induce entropy-driven spatial organization of polymer segments as seen in vivo. Here we design polymers of different architectures in our simulations by introducing a minimal number of cross-links between specific monomers along the ring-polymer contour. The cross-links are likely induced by various bridging proteins inside living cells. We investigate the segregation of two polymers with modified topologies confined in a cylinder, which initially had spatially overlapping configurations. This helps us to identify the architectures that lead to higher success rates of segregation. We also establish the mechanism that leads to localization of specific polymer segments. We use the blob model to provide a theoretical understanding of why certain architectures lead to enhanced entropic repulsive forces between the polymers. Lastly, we establish a correspondence between the organizational patterns of the chromosome of the C.crescentus bacterium and our results for a specifically designed polymer architecture. However, the principles outlined here pertaining to the organization of polymeric segments are applicable to both synthetic and biological polymers.

Confinement free energy for a polymer chain: Corrections to scaling

Spatial confinement of a polymer chain results in a reduction of conformational entropy. For confinement of a flexible N-mer chain in a planar slit or cylindrical pore (confining dimension D), a blob model analysis predicts the asymptotic scaling behavior ΔF/N ∼ D^-γ with γ ≈ 1.70, where ΔF is the free energy increase due to confinement. Here, we extend this scaling analysis to include the variation of local monomer density upon confinement giving ΔF/N ∼ D^-γ(1 - h(N, D)), where the correction-to-scaling term has the form h ∼ D^y/N^Δ with exponents y = 3 - γ ≈ 1.30 and Δ = 3/γ - 1 ≈ 0.76. To test these scaling predictions, we carry out Wang-Landau simulations of confined and unconfined tangent-hard-sphere chains (bead diameter σ) in the presence of a square-well trapping potential. The fully trapped chain provides a common reference state, allowing for an absolute determination of the confinement free energy. Our simulation results for 32 ≤ N ≤ 1024 and 3 ≤ D/σ ≤ 14 are well-described by the extended scaling relation giving exponents of γ = 1.69(1), y = 1.25(2), and Δ = 0.75(6).

Confinement anisotropy drives polar organization of two DNA molecules interacting in a nanoscale cavity

There is growing appreciation for the role phase transition based phenomena play in biological systems. In particular, self-avoiding polymer chains are predicted to undergo a unique confinement dependent demixing transition as the anisotropy of the confined space is increased. This phenomenon may be relevant for understanding how interactions between multiple dsDNA molecules can induce self-organized structure in prokaryotes. While recent in vivo experiments and Monte Carlo simulations have delivered essential insights into this phenomenon and its relation to bacteria, there are fundamental questions remaining concerning how segregated polymer states arise, the role of confinement anisotropy and the nature of the dynamics in the segregated states. To address these questions, we introduce an artificial nanofluidic model to quantify the interactions of multiple dsDNA molecules in cavities with controlled anisotropy. We find that two dsDNA molecules of equal size confined in an elliptical cavity will spontaneously demix and orient along the cavity poles as cavity eccentricity is increased; the two chains will then swap pole positions with a frequency that decreases with increasing cavity eccentricity. In addition, we explore a system consisting of a large dsDNA molecule and a plasmid molecule. We find that the plasmid is excluded from the larger molecule and will exhibit a preference for the ellipse poles, giving rise to a non-uniform spatial distribution in the cavity that may help explain the non-uniform plasmid distribution observed during in vivo imaging of high-copy number plasmids in bacteria.

Channels with Helical Modulation Display Stereospecific Sensitivity for Chiral Superstructures

By means of coarse-grained molecular dynamics simulations, we explore chiral sensitivity of confining spaces modelled as helical channels to chiral superstructures represented by polymer knots. The simulations show that helical channels exhibit stereosensitivity to chiral knots localized on linear chains by effect of external pulling force and also to knots embedded on circular chains. The magnitude of the stereoselective effect is stronger for torus knots, the effect is weaker in the case of twist knots, and amphichiral knots do exhibit no chiral effects. The magnitude of the effect can be tuned by the so-far investigated radius of the helix, the pitch of the helix and the strength of the pulling force. The model is aimed to simulate and address a range of practical situations that may occur in experimental settings such as designing of nanotechnological devices for the detection of topological state of molecules, preparation of new gels with tailor made stereoselective properties, or diffusion of knotted DNA in biological conditions.

Diffusion of Short Semiflexible DNA Polymer Chains in Strong and Moderate Confinement

In many technological applications, DNA is confined within nanoenvironments that are smaller than the size of the unconfined polymer in solution. However, the dependence of the diffusion coefficient on molecular weight and characteristic confinement dimension remains poorly understood in this regime. Here, convex lens-induced confinement (CLiC) was leveraged to examine how the diffusion of short DNA fragments varied as a function of slit height by using single-molecule fluorescence tracking microscopy. The diffusion coefficient followed approximate power law behavior versus confinement height, with exponents of 0.27 ± 0.01, 0.32 ± 0.02, and 0.42 ± 0.06 for 692, 1343, and 2686 base pair chains, respectively. The weak dependence on slit height suggests that shorter semiflexible chains may adopt increasingly rodlike conformations and therefore experience weaker excluded-volume interactions as the confinement dimension is reduced. The diffusion coefficient versus molecular weight also exhibited apparent power law behavior, with exponents that varied slightly (from -0.89 to -0.85) with slit height, consistent with hydrodynamic interactions intermediate between Rouse and Zimm model predictions.

Modeling the Protective Role of Bacterial Lipopolysaccharides against Membrane-Rupturing Peptides

Lipopolysaccharide (LPS) is a key surface component of Gram-negative bacteria, populating the outer layer of their outer membrane. A number of experimental studies highlight its protective role against harmful molecules such as antibiotics and antimicrobial peptides (AMPs). In this work, we present a theoretical model for describing the interaction between LPS and cationic antimicrobial peptides, which combines the following two key features. The polysaccharide part is viewed as forming a polymer brush, exerting an osmotic pressure on inclusions such as antimicrobial peptides. The charged groups on LPS (those in lipid A and the two Kdo groups in the inner core) form electrostatic binding sites for cationic AMPs or cations. Using the resulting model, we offer a quantitative picture of how the brush component enhances the protective role of LPS against magainin-like peptides, in the presence of divalent cations such as Mg²⁺. The LPS brush tends to diminish the interfacial binding of the peptides, at the lipid headgroup region, by about 30%. In the presence of 5 mM of Mg²⁺, the interfacial binding does not reach a threshold value for wild-type LPS, beyond which the LPS layer is ruptured, even though it does for LPS Re (the simplest form of LPS, lacking the brush part), as long as [AMP] ≤ 20 μM, where [AMP] is the concentration of AMPs. At a low concentration of Mg²⁺ (≈1 mM), however, a smaller [AMP] value (≳2 μM) is needed to reach the threshold coverage for wild-type LPS. Our results also suggest that the interfacial binding of peptides is insensitive to their possible weak interaction with the surrounding brush chains.

Coarse-grained molecular dynamics simulations of nanoplastics interacting with a hydrophobic environment in aqueous solution

Nanoplastics (NPs) are emerging threats for marine and terrestrial ecosystems, but little is known about their fate in the environment at the molecular scale. In this work, coarse-grained molecular dynamics simulations were performed to investigate nature and strength of the interaction between NPs and hydrophobic environments. Specifically, NPs were simulated with different hydrophobic and hydrophilic polymers while carbon nanotubes (CNTs) were used to mimic surface and confinement effects of hydrophobic building blocks occurring in a soil environment. The hydrophobicity of CNTs was modified by introducing different hydrophobic and hydrophilic functional groups at their inner surfaces. The results show that hydrophobic polymers have a strong affinity to adsorb at the outer surface and to be captured inside the CNT. The accumulation within the CNT is even increased in presence of hydrophobic functional groups. This contribution is a first step towards a mechanistic understanding of a variety of processes connected to interaction of nanoscale material with environmental systems. Regarding the fate of NPs in soil, the results point to the critical role of the hydrophobicity of NPs and soil organic matter (SOM) as well as of the chemical nature of functionalized SOM cavities/voids in controlling the accumulation of NPs in soil. Moreover, the results can be related to water treatment technologies as it is shown that the hydrophobicity of CNTs and functionalization of their surfaces may play a crucial role in enhancing the adsorption capacity of CNTs with respect to organic compounds and thus their removal efficiency from wastewater.

The binding mechanisms of nanoplastics (NPs) to carbon nanotubes as hydrophobic environmental systems have been explored by coarse-grained MD simulations. The results could be closely connected to fate of NPs in soil and water treatment technologies.

Collapse transition of a heterogeneous polymer in a crowded medium

Long chain molecules can be entropically compacted in a crowded medium. We study the compaction transition of a heterogeneous polymer with ring topology by crowding effects in a free or confined space. For this, we use molecular dynamics simulations in which the effects of crowders are taken into account through effective interactions between chain segments. Our parameter choices are inspired by the Escherichia coli chromosome. The polymer consists of small and big monomers; the big monomers dispersed along the backbone are to mimic the binding of RNA polymerases. Our results show that the compaction transition is a two-step process: initial compaction induced by the association (clustering) of big monomers followed by a gradual overall compaction. They also indicate that cylindrical confinement makes the initial transition more effective; for representative parameter choices, the initial compaction accounts for about 60% reduction in the chain size. Our simulation results support the view that crowding promotes clustering of active transcription units into transcription factories.

Towards a synthetic cell cycle

Recent developments in synthetic biology may bring the bottom-up generation of a synthetic cell within reach. A key feature of a living synthetic cell is a functional cell cycle, in which DNA replication and segregation as well as cell growth and division are well integrated. Here, we describe different approaches to recreate these processes in a synthetic cell, based on natural systems and/or synthetic alternatives. Although some individual machineries have recently been established, their integration and control in a synthetic cell cycle remain to be addressed. In this Perspective, we discuss potential paths towards an integrated synthetic cell cycle.

Mesoscale Simulation of Bacterial Chromosome and Cytoplasmic Nanoparticles in Confinement

In this study we investigated, using a simple polymer model of bacterial chromosome, the subdiffusive behaviors of both cytoplasmic particles and various loci in different cell wall confinements. Non-Gaussian subdiffusion of cytoplasmic particles as well as loci were obtained in our Langevin dynamic simulations, which agrees with fluorescence microscope observations. The effects of cytoplasmic particle size, locus position, confinement geometry, and density on motions of particles and loci were examined systematically. It is demonstrated that the cytoplasmic subdiffusion can largely be attributed to the mechanical properties of bacterial chromosomes rather than the viscoelasticity of cytoplasm. Due to the randomly positioned bacterial chromosome segments, the surrounding environment for both particle and loci is heterogeneous. Therefore, the exponent characterizing the subdiffusion of cytoplasmic particle/loci as well as Laplace displacement distributions of particle/loci can be reproduced by this simple model. Nevertheless, this bacterial chromosome model cannot explain the different responses of cytoplasmic particles and loci to external compression exerted on the bacterial cell wall, which suggests that the nonequilibrium activity, e.g., metabolic reactions, play an important role in cytoplasmic subdiffusion.

A Rotavirus Virus-Like Particle Confined Palladium Nanoreactor and Its Immobilization on Graphene Oxide for Catalysis

Abstract

In this work, a new viral protein cage based nanoreactor was successfully constructed via encapsulating Tween 80 stabilized palladium nanoparticles (NPs) into rotavirus capsid VP2 virus-like particles (i.e. Pd@VP2). The effects of stabilizers including CTAB, SDS, Tween 80 and PVP on controlling the particle size of Pd NPs were investigated. They were further immobilized on graphene oxide (i.e. Pd@VP2/GO) by a simple mixing method. Some characterizations including FT-IR and XPS were conducted to study adsorption mode of Pd@VP2 on GO sheets. Their catalytic performance was estimated in the reduction of 4-nitrophenol (4-NP). Results showed that Tween 80 stabilized Pd NPs with the molar ratio of Pd to Tween 80 at 1:0.1 possessed the smallest size and the best stability as well. They were encapsulated into viral protein cages (mean size 49 ± 0.26 nm) to assemble confined nanoreactors, most of which contained 1-2 Pd NPs (mean size 8.15 ± 0.26 nm). As-prepared Pd@VP2 indicated an enhanced activity (apparent reaction rate constant k _app  = (3.74 ± 0.10) × 10^-3 s^-1) for the reduction of 4-NP in comparison to non-confined Pd-Tween80 colloid (k _app  = (2.20 ± 0.06) × 10^-3 s^-1). It was logically due to confinement effects of Pd@VP2 including high dispersion of Pd NPs and high effective concentration of substrates in confined space. Pd@VP2 were further immobilized on GO surface through C-N bond. Pd@VP2/GO exhibited good reusability after recycling for four runs, confirming the strong anchoring effects of GO on Pd@VP2.

Graphic Abstract

All-or-none folding of a flexible polymer chain in cylindrical nanoconfinement

Geometric confinement of a polymer chain results in a loss of conformational entropy. For a chain that can fold into a compact native state via a first-order-like transition, as is the case for many small proteins, confinement typically provides an entropic stabilization of the folded state, thereby shifting the location of the transition. This allows for the possibility of confinement (entropy) driven folding. Here, we investigate such confinement effects for a flexible square-well-sphere N-mer chain (monomer diameter σ) confined within a long cylindrical pore (diameter D) or a closed cylindrical box (height H = D). We carry out Wang-Landau simulations to construct the density of states, which provides access to the complete thermodynamics of the system. For a wide pore, an entropic stabilization of the folded state is observed. However, as the pore diameter approaches the size of the folded chain (D ∼ N^1/3σ), we find a destabilization effect. For pore diameters smaller than the native ground-state, the chain folds into a different, higher energy, ground state ensemble and the T vs D phase diagram displays non-monotonic behavior as the system is forced into different ground states for different ranges of D. In this regime, isothermal reduction of the confinement dimension can induce folding, unfolding, or crystallite restructuring. For the cylindrical box, we find a monotonic stabilization effect with decreasing D. Scaling laws for the confinement free energy in the athermal limit are also investigated.

A lattice kinetic Monte-Carlo method for simulating chromosomal dynamics and other (non-)equilibrium bio-assemblies

Biological assemblies in living cells such as chromosomes constitute large many-body systems that operate in a fluctuating, out-of-equilibrium environment. Since a brute-force simulation of that many degrees of freedom is currently computationally unfeasible, it is necessary to perform coarse-grained stochastic simulations. Here, we develop all tools necessary to write a lattice kinetic Monte-Carlo (LKMC) algorithm capable of performing such simulations. We discuss the validity and limits of this approach by testing the results of the simulation method in simple settings. Importantly, we illustrate how at large external forces Metropolis-Hastings kinetics violate the fluctuation-dissipation and steady-state fluctuation theorems and discuss better alternatives. Although this simulation framework is rather general, we demonstrate our approach using a DNA polymer with interacting SMC condensin loop-extruding enzymes. Specifically, we show that the scaling behavior of the loop-size distributions that we obtain in our LKMC simulations of this SMC-DNA system is consistent with that reported in other studies using Brownian dynamics simulations and analytic approaches. Moreover, we find that the irreversible dynamics of these enzymes under certain conditions result in frozen, sterically jammed polymer configurations, highlighting a potential pitfall of this approach.

Confined crowded polymers near attractive surfaces

We present results from molecular dynamics simulations of a spherically confined neutral polymer in the presence of crowding particles, studying polymer shapes and conformations as a function of the strength of the attraction to the confining wall, solvent quality, and the density of crowders. The conformations of the polymer under good solvent conditions are weakly dependent on crowder particle density, even when the polymer is strongly confined. In contrast, under poor solvent conditions, when the polymer assumes a collapsed conformation when unconfined, it can exhibit transitions to two different adsorbed phases, when either the interaction with the wall or the density of crowder particles is changed. One such transition involves a desorbed collapsed phase change to an adsorbed extended phase as the attraction of the polymer towards the confining wall is increased. Such an adsorbed extended phase can exhibit a second transition to an ordered adsorbed collapsed phase as the crowder particle density is increased. The ordered adsorbed collapsed phase of the polymer differs significantly in its structure from the desorbed collapsed phase. We revisit the earlier understanding of the adsorption of confined polymers on attractive surfaces in light of our results.

Mapping a single-molecule folding process onto a topological space

Physics of protein folding has been dominated by conceptual frameworks including the nucleation-propagation mechanism and the diffusion-collision model, and none address the topological properties of a chain during a folding process. Single-molecule interrogation of folded biomolecules has enabled real-time monitoring of folding processes at an unprecedented resolution. Despite these advances, the topology landscape has not been fully mapped for any chain. Using a novel circuit topology approach, we map the topology landscape of a model polymeric chain. Inspired by single-molecule mechanical interrogation studies, we restrained the ends of a chain and followed fold nucleation dynamics. We find that, before the nucleation, transient local entropic loops dominate. Although the nucleation length of globules is dependent on the cohesive interaction, the ultimate topological states of the collapsed polymer are largely independent of the interaction but depend on the speed of the folding process. After the nucleation, transient topological rearrangements are observed that converge to a steady-state, where the fold grows in a self-similar manner.

Compressing Θ-Chain in Slit Geometry

When compressed in a slit of width D, a Θ-chain that displays the scaling of size R₀ (diameter) with respect to the number of monomers N, R₀ ∼ aN^1/2, expands in the lateral direction as R_|| ∼ aN^ν(a/D)^2ν-1. Provided that the Θ condition is strictly maintained throughout the compression, the well-known scaling exponent of Θ-chain in two dimensions, ν = 4/7, is anticipated in a perfect confinement. However, numerics shows that upon increasing compression from R₀/D < 1 to R₀/D ≫ 1, ν gradually deviates from ν = 1/2 and plateaus at ν = 3/4, the exponent associated with the self-avoiding walk in two dimensions. Using both theoretical considerations and numerics, we argue that it is highly nontrivial to maintain the Θ condition under confinement because of two major effects. First, as the dimension is reduced from three to two dimensions, the contributions of higher order virial terms, which can be ignored in three dimensions at large N, become significant, making the perturbative expansion used in Flory-type approach inherently problematic. Second and more importantly, the geometrical confinement, which is regarded as an applied external field, alters the second virial coefficient (B₂) changes from B₂ = 0 (Θ condition) in free space to B₂ > 0 (good-solvent condition) in confinement. Our study provides practical insight into how confinement affects the conformation of a single polymer chain.

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)].

Confinement induces helical organization of chromosome-like polymers

Helical organization is commonly observed for a variety of biopolymers. Here we study the helical organization of two types of biopolymers, i.e., DNA-like semiflexible and bottle-brush polymers, in a cell-like confined space. A bottle-brush polymer consists of a backbone and side chains emanating from the backbone, resembling a supercoiled bacterial chromosome. Using computer simulations, we calculate 'writhe' distributions of confined biopolymers for a wide range of parameters. Our effort clarifies the conditions under which biopolymers are helically organized. While helical organization is not easily realized for DNA-like biomolecules, cylindrical confinement can induce spiral patterns in a bottle brush, similarly to what was observed with bacterial chromosomes. They also suggest that ring-shape bottle brushes have a stronger tendency for helical organization. We discuss how our results can be used to interpret chromosome experiments. For instance, they suggest that experimental resolution has unexpected consequences on writhe measurements (e.g., narrowing of the writhe distribution and kinetic separation of opposite helical states).

Revealing the critical role of template functional group ordering in the template-directed crystallization of pyrazinamide

The nucleation of γ form pyrazinamide can be directed by the ordering and specific orientation of the template functional groups.

Comparative study of the crowding-induced collapse effect in hard-sphere, flexible polymer and rod-like polymer systems

A systematic Langevin simulation is performed to study the crowding-induced collapse effect on a probed chain in three typical systems: hard sphere (HS), flexible polymer and rod-like polymer.

Probing the organization and dynamics of two DNA chains trapped in a nanofluidic cavity

Here we present a pneumatically-actuated nanofluidic platform that has the capability of dynamically controlling the confinement environment of macromolecules in solution. Using a principle familiar from classic devices based on soft-lithography, the system uses pneumatic pressure to deflect a thin nitride lid into a nanoslit, confining molecules in an array of cavities embedded in the slit. We use this system to quantify the interactions of multiple confined DNA chains, a key problem in polymer physics with important implications for nanofluidic device performance and DNA partitioning/organization in bacteria and the eukaryotes. In particular, we focus on the problem of two-chain confinement, using differential staining of the chains to independently assess the chain conformation, determine the degree of partitioning/mixing in the cavities and assess coupled diffusion of the chain center-of-mass positions. We find that confinement of more than one chain in the cavity can have a drastic impact on the polymer dynamics and conformation.

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.

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.

A ring-polymer model shows how macromolecular crowding controls chromosome-arm organization in Escherichia coli

Macromolecular crowding influences various cellular processes such as macromolecular association and transcription, and is a key determinant of chromosome organization in bacteria. The entropy of crowders favors compaction of long chain molecules such as chromosomes. To what extent is the circular bacterial chromosome, often viewed as consisting of “two arms”, organized entropically by crowding? Using computer simulations, we examine how a ring polymer is organized in a crowded and cylindrically-confined space, as a coarse-grained bacterial chromosome. Our results suggest that in a wide parameter range of biological relevance crowding is essential for separating the two arms in the way observed with Escherichia coli chromosomes at fast-growth rates, in addition to maintaining the chromosome in an organized collapsed state. Under different conditions, however, the ring polymer is centrally condensed or adsorbed onto the cylindrical wall with the two arms laterally collapsed onto each other. We discuss the relevance of our results to chromosome-membrane interactions.

How are molecular crowding and the spatial organization of a biopolymer interrelated

In a crowded cellular interior, dissolved biomolecules or crowders exert excluded volume effects on other biomolecules, which in turn control various processes including protein aggregation and chromosome organization. As a result of these effects, a long chain molecule can be phase-separated into a condensed state, redistributing the surrounding crowders. Using computer simulations and a theoretical approach, we study the interrelationship between molecular crowding and chain organization. In a parameter space of biological relevance, the distributions of monomers and crowders follow a simple relationship: the sum of their volume fractions rescaled by their size remains constant. Beyond a physical picture of molecular crowding it offers, this finding explains a few key features of what has been known about chromosome organization in an E. coli cell.

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.

Polymer globule with fractal properties caused by intramolecular nanostructuring and spatial constrains

By means of computer simulation, we studied macromolecules composed of N dumbbell amphiphilic monomer units with attractive pendant groups. In poor solvents, these macromolecules form spherical globules that are dense in the case of short chains (the gyration radius RG∼N(1/3)), or hollow inside and obey the RG∼N(1/2) law when the macromolecules are sufficiently long. Due to the specific intramolecular nanostructuring, the vesicle-like globules of long amphiphilic macromolecules posses some properties of fractal globules, by which they (i) could demonstrate the same scaling statistics for the entire macromolecule and for short subchains with m monomer units and (ii) possess a specific territorial structure. Within a narrow slit, the globule loses its inner cavity, takes a disk-like shape and scales as N(1/2) for much shorter macromolecules. However, the field of end-to-end distance r(m) ∼m(1/2) dependence for subchains becomes visibly smaller. The results obtained were compared with the homopolymer case.

The polymer physics of single DNA confined in nanochannels

In recent years, applications and experimental studies of DNA in nanochannels have stimulated the investigation of the polymer physics of DNA in confinement. Recent advances in the physics of confined polymers, using DNA as a model polymer, have moved beyond the classic Odijk theory for the strong confinement, and the classic blob theory for the weak confinement. In this review, we present the current understanding of the behaviors of confined polymers while briefly reviewing classic theories. Three aspects of confined DNA are presented: static, dynamic, and topological properties. The relevant simulation methods are also summarized. In addition, comparisons of confined DNA with DNA under tension and DNA in semidilute solution are made to emphasize universal behaviors. Finally, an outlook of the possible future research for confined DNA is given.

Finite-size corrections for confined polymers in the extended de Gennes regime

Theoretical results for the extension of a polymer confined to a channel are usually derived in the limit of infinite contour length. But experimental studies and simulations of DNA molecules confined to nanochannels are not necessarily in this asymptotic limit. We calculate the statistics of the span and the end-to-end distance of a semiflexible polymer of finite length in the extended de Gennes regime, exploiting the fact that the problem can be mapped to a one-dimensional weakly self-avoiding random walk. The results thus obtained compare favorably with pruned-enriched Rosenbluth method (PERM) simulations of a three-dimensional discrete wormlike chain model of DNA confined in a nanochannel. We discuss the implications for experimental studies of linear λ-DNA confined to nanochannels at the high ionic strengths used in many experiments.