Compressing a confined DNA: from nano-channel to nano-cavity

Mixing and demixing arising from compression of two semiflexible polymer chains in nanochannels

We use molecular dynamics simulation to probe the non-equilibrium physics of two nanochannel-confined semiflexible polymers in a homogeneous flow field. We find that for sufficiently stiff chains the internal organization of the two chains takes the form of interwoven folds and circular coils. This organization can lead to mixing or demixing depending on chain stiffness and flow speed. At low and intermediate flow, the two chains adopt a folded configuration, which favours mixing. At high flow, the two chains adopt a predominantly coiled configuration. The coiled configuration results in demixing when the chains are compressed from an initially demixed condition and mixing when the chains are compressed from an initially mixed condition. We find that the mixing/demixing behaviour is governed by the ratio of the number of folded segments of one chain relative to the other at low flow and by the degree of coiling in both chains at high flow. For decreasing stiffness, the chains start to aggregate locally instead of mixing smoothly at low and intermediate flow. In the limit of completely flexible chains, the two chains either completely segregate at low flow, or adopt a locally demixed configuration consisting of large aggregates of one chain relative to the other that undergo complex stochastic dynamics, diffusing, disintegrating, and reforming at intermediate flow. The transition from complete segregation to the aggregate-dominated configuration occurs when the linear intra-chain ordering breaks down.

Pressure of Linear and Ring Polymers Confined in a Cavity

Nanoscale confinement of polymers in a cavity is central to a variety of biological and nanotechnology processes. Using the discrete WLC model we simulate the compression of flexible and semiflexible polymers of linear and ring topology in a closed cavity. Simulation reveals that polymer pressure inside the cavity increases with the chain stiffness but is practically unaffected by the chain topology. For flexible polymers, the computed dependence of pressure on the cavity size and polymer concentration is consistent with the scaling behavior expected for bulk polymers in a good solvent. However, the scaling behavior of semiflexible polymers is only in partial agreement with the theory prediction, with discrepancies arising from a continuous transition between regimes in chains of moderate lengths. The computed segment density profiles endorse the propensity of semiflexible polymers to concentrate beneath the cavity surface and thus elevate the pressure. The compaction of polymers by compression into the disordered globule or growing toroidal structure is documented.

Entropic force of cone-tethered polymers interacting with a planar surface

Computer simulations are used to characterize the entropic force of one or more polymers tethered to the tip of a hard conical object that interact with a nearby hard flat surface. Pruned-enriched Rosenbluth method Monte Carlo simulations are used to calculate the variation of the conformational free energy F of a hard-sphere polymer with respect to the cone-tip-to-surface distance h from which the variation of the entropic force f≡|dF/dh| with h is determined. We consider the following cases: a single freely jointed tethered chain, a single semiflexible tethered chain, and several freely jointed chains of equal length each tethered to the cone tip. The simulation results are used to test the validity of a prediction by Maghrebi et al. [Maghrebi et al., Europhys. Lett. 96, 66002 (2011)0295-507510.1209/0295-5075/96/66002; Phys. Rev. E 86, 061801 (2012)1539-375510.1103/PhysRevE.86.061801] that f∝(γ_{∞}-γ_{0})h^{-1}, where γ_{0} and γ_{∞} are universal scaling exponents for the partition function of the tethered polymer for h=0 and h=∞, respectively. The measured functions f(h) are generally consistent with the predictions, with small quantitative discrepancies arising from the approximations employed in the theory. In the case of multiple tethered polymers, the entropic force per polymer is roughly constant, which is qualitatively inconsistent with the predictions.

Organized states arising from compression of single semiflexible polymer chains in nanochannels

We use molecular dynamics simulation to probe the nonequilibrium physics of single nanochannel-confined semiflexible polymers in a homogeneous flow field. The flow field compresses the polymer against the end of the nanochannel, simulating an experiment of a nanochannel confined chain compressed against a slit barrier. The flow-based compression gives rise to a packing of the chain against the channel end that possesses a striking organization, consisting of interweaving of folds and circular coils. For stiff chains at low flow, we find that the organization is dominated by repeated hairpin folds. For stiff chains at higher flow, we observe that circular coils arise along with the folds, with folding and coiling domains becoming interwoven at the highest flow speeds. Chain organization is retained even when the chain persistence length is on order of the channel width. We show that the global polymer organization, consisting of a number of defined folds and coiled loops, arises from the minimization of the total chain free energy.

Polymer translocation into cavities: Effects of confinement geometry, crowding, and bending rigidity on the free energy

Monte Carlo simulations are used to study the translocation of a polymer into a cavity. Modeling the polymer as a hard-sphere chain with a length up to N=601 monomers, we use a multiple-histogram method to measure the variation of the conformational free energy of the polymer with respect to the number of translocated monomers. The resulting free-energy functions are then used to obtain the confinement free energy for the translocated portion of the polymer. We characterize the confinement free energy for a flexible polymer in cavities with constant cross-sectional area A for various cavity shapes (cylindrical, rectangular, and triangular) as well as for tapered cavities with pyramidal and conical shape. The scaling of the free energy with cavity volume and translocated polymer subchain length is generally consistent with predictions from simple scaling arguments, with small deviations in the scaling exponents likely due to finite-size effects. The confinement free energy depends strongly on cavity shape anisometry and is a minimum for an isometric cavity shape with a length-to-width ratio of unity. Entropic depletion at the edges or vertices of the confining cavity are evident in the results for constant-A and pyramidal cavities. For translocation into infinitely long cones, the scaling of the free energy with taper angle is consistent with a theoretical prediction employing the blob model. We also examine the effects of polymer bending rigidity on the translocation free energy for cylindrical cavities. For isometric cavities, the observed scaling behavior is in partial agreement with theoretical predictions, with discrepancies arising from finite-size effects that prevent the emergence of well-defined scaling regimes. In addition, translocation into highly anisometric cylindrical cavities leads to a multistage folding process for stiff polymers. Finally, we examine the effects of crowding agents inside the cavity. We find that the confinement free energy increases with crowder density. At constant packing fraction the magnitude of this effect lessens with increasing crowder size for a crowder-to-monomer size ratio ≥1.