The shape of a flexible polymer in a cylindrical pore

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.

Nonmonotonous Translocation Time of Polymers across Pores

Polymers confined in corrugated channels, i.e., channels of varying amplitude, display multiple local maxima and minima of the diffusion coefficient upon increasing their degree of polymerization N. We propose a theoretical effective free energy for linear polymers based on a Fick-Jacobs approach. We validate the predictions against numerical data, obtaining quantitative agreement for the effective free energy, the diffusion coefficient, and the mean first passage time. Finally, we employ the effective free energy to compute the polymer lengths N_{min} at which the diffusion coefficient presents a minimum: we find a scaling expression that we rationalize with a blob model. Our results could be useful to design porous adsorbers, that separate polymers of different sizes without the action of an external flow.

Adsorption of Semiflexible Polymers in Cylindrical Tubes

Conformations of wormlike chains in cylindrical pores with attractive walls are explored for varying pore radius and strength of the attractive wall potential by molecular dynamics simulations of a coarse-grained model. Local quantities such as the fraction of monomeric units bound to the surface and the bond-orientational order parameter as well as the radial density distribution are studied, as well as the global chain extensions parallel to the cylinder axis and perpendicular to the cylinder surface. A nonmonotonic convergence of these properties to their counterparts for adsorption on a planar substrate is observed due to the conflict between pore surface curvature and chain stiffness. Also the interpretation of partially adsorbed chains in terms of trains, loops, and tails is discussed.

Cylindrical confinement of solutions containing semiflexible macromolecules: surface-induced nematic order versus phase separation

Solutions of semiflexible polymers confined in cylindrical pores with repulsive walls are studied by Molecular Dynamics simulations for a wide range of polymer concentrations. Both the case where both lengths are of the same order and the case when the persistence length by far exceeds the contour length are considered, and the enhancement of nematic order along the cylinder axis is characterized. With increasing density the character of the surface effect changes from depletion to the formation of a layered structure. For binary 50 : 50 mixtures of the two types of polymers an interplay between surface enrichment of the stiffer component and the isotropic-nematic transition is found, and a phase separated structure with cylindrical symmetry occurs, with the isotropic phase located around the cylinder axis. For melt densities the mixed nematic phase forms at the wall a layer with a screw-like structure of a tilted smectic phase. The observed behavior is tentatively interpreted in terms of the competition of the chain orientational entropy with entropy of mixing and excluded volume due to the wall.

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.

Polymers grown in cavities: Vesicles and droplets

In synthetic chemistry and biological or biomimetic systems, polymers are often grown in cavities. Polymerizations in microemulsions, biopolymers grown in cells, or in vesicles containing artificial organelles have an influence on the shape of liquid boundaries. We consider confined grand-canonical polymers to address equilibrium properties of annealed polymers. We calculate the concentration profiles established by annealed (star-) polymers inside a confining cavity. Our emphasis is on the description of pressure fields derived from the contact theorem. We further show how the pressure field exerted by a localized annealed polymer (or pair of polymers) deforms the confining vesicle/ microemulsions droplet.

Trigger factor chaperone acts as a mechanical foldase

Proteins fold under mechanical forces in a number of biological processes, ranging from muscle contraction to co-translational folding. As force hinders the folding transition, chaperones must play a role in this scenario, although their influence on protein folding under force has not been directly monitored yet. Here, we introduce single-molecule magnetic tweezers to study the folding dynamics of protein L in presence of the prototypical molecular chaperone trigger factor over the range of physiological forces (4–10 pN). Our results show that trigger factor increases prominently the probability of folding against force and accelerates the refolding kinetics. Moreover, we find that trigger factor catalyzes the folding reaction in a force-dependent manner; as the force increases, higher concentrations of trigger factor are needed to rescue folding. We propose that chaperones such as trigger factor can work as foldases under force, a mechanism which could be of relevance for several physiological processes.

Proteins fold under mechanical force during co-translational folding at the ribosome. Here, the authors use single molecule magnetic tweezers to study the influence of chaperones on protein folding and show that the ribosomal chaperone trigger factor acts as a mechanical foldase by promoting protein folding under force.

As a Nucleus Enters a Small Pore, Chromatin Stretches and Maintains Integrity, Even with DNA Breaks

As a cell pushes or pulls its nucleus through a small constriction, the chromatin must distort and somehow maintain genomic stability despite ever-present double-strand breaks in the DNA. Here we visualize within a living cell the pore-size dependent deformation of a specific locus engineered into chromosome-1 and cleaved. An mCherry-tagged nuclease targets the submicron locus, causing DNA cleavage and recruiting repair factors such as GFP-53BP1 to a large region around the locus. Aspiration of a cell and its nucleus into a micropipette shows that chromatin aligns and stretches parallel to the pore. Extension is largest in small pores, increasing >10-fold but remaining 30-fold shorter than the DNA contour length in the locus. Brochard and de Gennes' blob model for tube geometry fits the data, with a simple modification for chromatin crowding. Continuity of the highly extended, cleaved chromatin is also maintained, consistent with folding and cross bridging of the DNA. Surprisingly, extensional integrity is unaffected by an inhibitor of the DNA repair scaffold.

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.

Electrophoretic Mobility of Polyelectrolytes within a Confining Well

We present a numerical study of polyelectrolytes electrophoresing in free solution while squeezed by an axisymmetric confinement force transverse to their net displacement. Hybrid multiparticle collision dynamics and molecular dynamics simulations with mean-field finite Debye layers show that even though the polyelectrolyte chains remain "free-draining" their electrophoretic mobility increases with confinement in nanoconfining potential wells. The primary mechanism leading to the increase in mobility above the free-solution value, despite long-range hydrodynamic screening by counterion layers, is the orientation of polymer segments within Debye layers. The observed length dependence of the electrophoretic mobility arises due to secondary effects of counterion condensation related to confinement compactification.

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.

Confined polymers in the extended de Gennes regime

We show that the problem of describing the conformations of a semiflexible polymer confined to a channel can be mapped onto an exactly solvable model in the so-called extended de Gennes regime. This regime (where the polymer is neither weakly nor strongly confined) has recently been intensively studied experimentally and by means of computer simulations. The exact solution predicts precisely how the conformational fluctuations depend upon the channel width and upon the microscopic parameters characterizing the physical properties of the polymer.

Comparison of linear and ring DNA macromolecules moderately and strongly confined in nanochannels

Understanding the mechanism of DNA extension in nanochannels is necessary for interpretation of experiments in nanofluidic channel devices that have been conducted recently with both linear and ring chains. The present article reviews the situation with linear chains and analyses the experimental results and simulations for channel-induced extension (linearization) of ring chains. Results for confined rings indicate a transition between moderate and strong confinement similar to that of linear chains. Owing to stronger self-avoidance in confined rings, the transition and chain extension is shifted relative to linear DNA. We suggest that a relationship similar to that used for the extension of linear chains may also be used for circular DNA.

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.

The entropic cost of polymer confinement

The confinement of a polymer into a small space is thermodynamically unfavorable because of the reduction in the number of conformational states. The entropic penalty affects a variety of biological processes, and it plays an important role in polymer transport properties and in microfluidic devices. We determine the entropic penalty for the confinement of elastic polymer of persistence length P in the long-chain limit. We examine three geometries: (1) parallel planes separated by a distance d (a slit); (2) a circular tube of diameter d; and (3) a sphere of diameter d. We first consider infinitely thin (ideal) chains. As d/P drops from 100 to 0.01, TΔS rises from ∼5 × 10(-4) kT to ∼30 kT per persistence length for cases (1) and (2), with the entropic penalty for case (2) being consistently about twice that for case (1). TΔS is ∼5 kT per persistence length for confinement to a sphere when d = P, about twice the value predicted by mean field theory. For all three geometries, in the limit d/P ≫ 1, the asymptotic behavior of ΔS vs d is consistent with the d(-2) behavior predicted by theory. In the limit d/P ≪ 1, the scaling of ΔS for slits and tubes is also consistent with earlier predictions (d(-2/3)). Finally, we treat volume exclusion effects, examining chains of diameter D > 0. Confinement to a narrow slit or tube (d/P ≪ 1) has the same entropic penalty as that for an ideal chain in a slit or tube with d' = d - D; in the weak confinement regime (d/P ≫ 1), the entropic penalties are significantly larger than those for infinitely thin chains. When a chain of finite diameter is forced into a sphere or other closed cavity, the entropic confinement penalty rises without limit because there are no configurations available to the chain once its volume exceeds that of the cavity.

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.

Factors governing helix formation in peptides confined to carbon nanotubes

The effect of confinement on the stability and dynamics of peptides and proteins is relevant in the context of a number of problems in biology and biotechnology. We have examined the stability of different helix-forming sequences upon confinement to a carbon nanotube using Langevin dynamics simulations of a coarse-grained representation of the polypeptide chain. We show that the interplay of several factors that include sequence, solvent conditions, strength (lambda) of nanotube-peptide interactions, and the nanotube diameter (D) determines confinement-induced stability of helicies. In agreement with predictions based on polymer theory, the helical state is entropically stabilized for all sequences when the interaction between the peptide and the nanotube is weakly hydrophobic and D is small. However, there is a strong sequence dependence as the strength of the lambda increases. For an amphiphilic sequence, the helical stability increases with lambda, whereas for polyalanine the diagram of states is a complex function of lambda and D. In addition, decreasing the size of the "hydrophobic patch" lining the nanotube, which mimics the chemical heterogeneity of the ribosome tunnel, increases the helical stability of the polyalanine sequence. Our results provide a framework for interpreting a number of experiments involving the structure formation of peptides in the ribosome tunnel as well as transport of biopolymers through nanotubes.

Compression and stretching of a self-avoiding chain in cylindrical nanopores

Force-induced deformations of a self-avoiding chain confined inside a cylindrical cavity, with diameter D, are probed using molecular dynamics simulations, scaling analysis, and analytical calculations. We obtain and confirm a simple scaling relation -fD approximately R(-9/4) in the strong-compression regime, while for weak deformations, we find fD = -A(R/R0) + B(R/R0)(-2), where A and B are constants, f the external force, and R the chain extension (with R0 its unperturbed value). For a strong stretch, we present a universal, analytical force-extension relation. Our results can be used to analyze the behavior of biomolecules in confinement.

Shapes of semiflexible polymers in confined spaces

Biological macromolecules, living in the confines of a cell, often adopt conformations that are unlikely to occur in free space. In this paper, we investigate the effects of confinement on the shape of a semiflexible chain. Results of Monte Carlo simulations show the existence of a shape transition when the persistence length of the polymer becomes comparable to the dimensions of the box. An order parameter is introduced to quantify this behavior. A simple model is constructed to study the effect of the shape transition on the effective persistence length of the polymer.

Structure of flexible and semiflexible polyelectrolyte chains in confined spaces of slit micro/nanochannels

Understanding the behavior of a polyelectrolyte in confined spaces has direct relevance in design and manipulation of microfluidic devices, as well as transport in living organisms. In this paper, a coarse-grained model of anionic semiflexible polyelectrolyte is applied, and its structure and dynamics are fully examined with Brownian dynamics (BD) simulations both in bulk solution and under confinement between two negatively charged parallel plates. The modeling is based on the nonlinear bead-spring discretization of a continuous chain with additional long-range electrostatic, Lennard-Jones, and hydrodynamic interactions between pairs of beads. The authors also consider the steric and electrostatic interactions between the bead and the confining wall. Relevant model parameters are determined from experimental rheology data on the anionic polysaccharide xanthan reported previously. For comparison, both flexible and semiflexible models are developed accompanying zero and finite intrinsic persistence lengths, respectively. The conformational changes of the polyelectrolyte chain induced by confinements and their dependence on the screening effect of the electrolyte solution are faithfully characterized with BD simulations. Depending on the intrinsic rigidity and the medium ionic strength, the polyelectrolyte can be classified as flexible, semiflexible, or rigid. Confined flexible and semiflexible chains exhibit a nonmonotonic variation in size, as measured by the radius of gyration and end-to-end distance, with changing slit width. For the semiflexible chain, this is coupled to the variations in long-range bond vector correlation. The rigid chain, realized at low ionic strength, does not have minima in size but exhibits a sigmoidal transition. The size of confined semiflexible and rigid polyelectrolytes can be well described by the wormlike chain model once the electrostatic effects are taken into account by the persistence length measured at long length scale.

Confined space and effective interactions of multiple self-avoiding chains

We study the link between three seeming-disparate cases of self-avoiding polymers: strongly overlapping multiple chains in dilute solution, chains under spherical confinement, and the onset of semidilute solutions. Our main result is that the free energy for overlapping n chains is independent of chain length and scales as n9/4, slowly crossing over to n3, as n increases. For strongly confined polymers inside a spherical cavity, we show that rearranging the chains does not cost an additional free energy. Our results imply that, during cell cycle, global reorganization of eukaryotic chromosomes in a large cell nucleus could be readily achieved.