Weak-to-strong confinement transition of semi-flexible macromolecules in slit and in channel

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.

Density fields for branching, stiff networks in rigid confining regions

We develop a formalism to describe the equilibrium distributions for segments of confined branched networks consisting of stiff filaments. This is applicable to certain situations of cytoskeleton in cells, such as for example actin filaments with branching due to the Arp2/3 complex. We develop a grand ensemble formalism that enables the computation of segment density and polarisation profiles within the confines of the cell. This is expressed in terms of the solution to nonlinear integral equations for auxiliary functions. We find three specific classes of behaviour depending on filament length, degree of branching and the ratio of persistence length to the dimensions of the geometry. Our method allows a numerical approach for semi-flexible filaments that are networked.

Block Copolymer of Flexible and Semi-Flexible Block Confined in Nanopost Array

Coarse-grained molecular dynamics simulations of a diblock copolymer consisting of a flexible and semi-flexible block in a dense array of parallel nanoposts with a square lattice packing were performed. The mutual interactions between the two blocks of the confined diblock chain were investigated through a comparison of their size, structure, and penetration among nanoposts with the corresponding separate chains. The geometry of a nanopost array was varied at constant post separation or at constant width of the passage between nanoposts. The size of a single interstitial volume was comparable to or smaller than the size of the diblock chain. A comparison of the blocks with their separate analogous chains revealed that the mutual interactions between the blocks were shielded by the nanoposts and, thus, the blocks behaved independently. At constant passage width, competitive effects of the axial chain extension in interstitial volumes and the lateral chain expansion among interstitial volumes led to a nonmonotonic behavior of the axial span. The position of the maximum in the span plotted against the filling fraction for a diblock chain was dictated by the semi-flexible block. The semi-flexible block penetrates among the nanoposts more readily and the expansion of the whole diblock copolymer is governed by the semiflexible block. The main findings were explained using the free energy arguments when an interstitial volume was approximated by a channel geometry and a passage aperture by a slit geometry. Detail knowledge of controlled conformational behavior in a compartmentalized environment can contribute to new processes in the storage and retrieval of information.

Investigation of nematic to smectic phase transition and dynamical properties of strongly confined semiflexible polymers using Langevin dynamics

We investigated the nematic to smectic phase transition for strongly confined semiflexible polymer solutions in slit-like confinements using GPU-accelerated Langevin dynamics. We characterized the phase transitions from the nematic to smectic phases for semi-flexible polymer solutions as the polymer density increased. The dependence for the lyotropic nematic to smectic transition can be collapsed by scaling exponents between 0.2 and 0.3. The smectic C phase is found for all the cases with the polymer orientation director tilted with respect to smectic layer lateral alignment. As the chain rigidity increases, the transition density decreases for systems in which the polymer persistence length (P) to slit height (H) ratios are 1.25, 2.5, 3.75, 5 and 25. We also characterized the polymer dynamics for the isotropic-nematic-smectic transitions. The overall polymer diffusivity decreased steadily as the polymer density increased. We observed anomalous polymer diffusion along the nematic director near the isotropic-nematic transition, similar to previously reported behavior for nematic-forming ellipsoids. Polymer diffusivity decreased sharply by two orders of magnitude upon the nematic-smectic transition.

Evidence for the extended de Gennes regime of a semiflexible polymer in slit confinement

We use off-lattice, pruned-enriched Rosenbluth method (PERM) simulations to compute the confinement free energy of a real wormlike chain of effective width w and persistence length lp in a slit of height H. For slit heights much larger than the persistence length of the polymer and much smaller than the thermal blob size, the excess free energy of the confined chain is consistent with a modified version of the scaling theory for the extended de Gennes regime in a channel that reflects the blob statistics in slit confinement. Explicitly, for channel sizes [Formula: see text], the difference between the confinement free energy of the real chain and that of an ideal chain scales like w/H.

Computational approaches for investigating interfacial adhesion phenomena of polyimide on silica glass

This manuscript provides a comprehensive study of adhesion behavior and its governing mechanisms when polyimide undergoes various modes of detachment from silica glass. Within the framework of steered molecular dynamics, we develop three different adhesion measurement techniques: pulling, peeling, and sliding. Such computational methodologies can be applied to investigate heterogeneous materials with differing interfacial adhesion modes. Here, a novel hybrid potential involving a combination of the INTERFACE force field in conjunction with ReaxFF and including Coulombic and Lennard-Jones interactions is employed to study such interfaces. The studies indicate that the pulling test requires the largest force and the shortest distance to detachment as the interfacial area is separated instantaneously, while the peeling test is observed to exhibit the largest distance for detachment because it separates via line-by-line adhesion. Two kinds of polyimides, aromatic and aliphatic type, are considered to demonstrate the rigidity dependent adhesion properties. The aromatic polyimide, which is more rigid due to the stronger charge transfer complex between chains, requires a greater force but a smaller distance at detachment than the aliphatic polyimide for all of the three methodologies.

Compressive response and helix formation of a semiflexible polymer confined in a nanochannel

Configurations of a single semiflexible polymer is studied when it is pushed into a nanochannel in the case where the polymer persistence length l_{p} is much longer than the channel diameter D:l_{p}/D≫1. Using numerical simulations, we show that the polymer undergoes a sequence of recurring structural transitions upon longitudinal compression: random deflection along the channel, a helix going around the channel wall, double-fold random deflection, double-fold helix, etc. We find that the helix transition can be understood as buckling of deflection segments, and the initial helix formation takes place at very small compression with no appreciable weak compression regime of the random deflection polymer.

Comparison of a stripe and slab confinement for ring and linear macromolecules in nanochannel

The combined effects of the channel asymmetry and the closed chain topology on the chain extension, structure factor, and the orientation correlations were studied using coarse-grained molecular dynamics simulations for moderate chain lengths. These effects are related to applications in linearization experiments with a DNA molecule in nanofluidic devices. According to the aspect ratio, the channels are classified as a stripe or slabs. The chain segments do not have any freedom to move in the direction of the narrowest stripe size, being approximately the same size as the segment size. The chains of both ring and linear topologies are extended more in a stripe than in a slab; this effect is strengthened for a ring. For a ring in a stripe, the extension-confinement strength dependence leads to effective Flory exponents even larger than 3/4, which is characteristic for a self-avoiding two-dimensional chain. While the chain extension-confinement strength dependence for both topologies conforms to the de Gennes regime in a stripe, a linear chain undergoes gradual transition to the pseudoideal regime as the slab height increases in the slab-like confinement. For a confined circle, the onset of the pseudoideal regime is shifted to larger slab heights. The structure factor confirms the absence of the pseudoideal and extended de Gennes regime in a stripe and the transition from the extended to the pseudoideal regime of a circular and linear chain upon increasing the slab heights.

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.

Modeling the relaxation of internal DNA segments during genome mapping in nanochannels

We have developed a multi-scale model describing the dynamics of internal segments of DNA in nanochannels used for genome mapping. In addition to the channel geometry, the model takes as its inputs the DNA properties in free solution (persistence length, effective width, molecular weight, and segmental hydrodynamic radius) and buffer properties (temperature and viscosity). Using pruned-enriched Rosenbluth simulations of a discrete wormlike chain model with circa 10 base pair resolution and a numerical solution for the hydrodynamic interactions in confinement, we convert these experimentally available inputs into the necessary parameters for a one-dimensional, Rouse-like model of the confined chain. The resulting coarse-grained model resolves the DNA at a length scale of approximately 6 kilobase pairs in the absence of any global hairpin folds, and is readily studied using a normal-mode analysis or Brownian dynamics simulations. The Rouse-like model successfully reproduces both the trends and order of magnitude of the relaxation time of the distance between labeled segments of DNA obtained in experiments. The model also provides insights that are not readily accessible from experiments, such as the role of the molecular weight of the DNA and location of the labeled segments that impact the statistical models used to construct genome maps from data acquired in nanochannels. The multi-scale approach used here, while focused towards a technologically relevant scenario, is readily adapted to other channel sizes and polymers.

Dimensional reduction of duplex DNA under confinement to nanofluidic slits

There has been much interest in the dimensional properties of double-stranded DNA (dsDNA) confined to nanoscale environments as a problem of fundamental importance in both biological and technological fields. This has led to a series of measurements by fluorescence microscopy of single dsDNA molecules under confinement to nanofluidic slits. Despite the efforts expended on such experiments and the corresponding theory and simulations of confined polymers, a consistent description of changes of the radius of gyration of dsDNA under strong confinement has not yet emerged. Here, we perform molecular dynamics (MD) simulations to identify relevant factors that might account for this inconsistency. Our simulations indicate a significant amplification of excluded volume interactions under confinement at the nanoscale due to the reduction of the effective dimensionality of the system. Thus, any factor influencing the excluded volume interaction of dsDNA, such as ionic strength, solution chemistry, and even fluorescent labels, can greatly influence the dsDNA size under strong confinement. These factors, which are normally less important in bulk solutions of dsDNA at moderate ionic strengths because of the relative weakness of the excluded volume interaction, must therefore be under tight control to achieve reproducible measurements of dsDNA under conditions of dimensional reduction. By simulating semi-flexible polymers over a range of parameter values relevant to the experimental systems and exploiting past theoretical treatments of the dimensional variation of swelling exponents and prefactors, we have developed a novel predictive relationship for the in-plane radius of gyration of long semi-flexible polymers under slit-like confinement. Importantly, these analytic expressions allow us to estimate the properties of dsDNA for the experimentally and biologically relevant range of contour lengths that is not currently accessible by state-of-the-art MD simulations.

Rigorous study of molecular dynamics of a single dsDNA confined in a nanochannel: Introduction of a critical mobility behaviour

The essential and effective characteristics of a double-stranded DNA (dsDNA) confined in a nanochannel is revisited by employing the rigorous full numerical approach of Molecular Dynamics (MD). The deformation of dsDNA and wall-biomolecule interaction which is critical in highly confined regime has been precisely imposed in numerical simulations. The numerical approach has been justified against available theoretical outcomes. A new and general expression for DNA electrophoretic mobility versus DNA length is extracted from numerical simulation which is out of reach of experimental methods due to practical shortcomings. The newly derived expression suggests an essential correction in the previously proposed expression for the critical case of small DNA molecules and reveals an astonishingly unbeknown trend of small DNA's mobility. Sub-molecular phenomenon of dsDNA melting under the condition of large external force is also studied. Assuming strong electric field exertion, the MD approach aptly demonstrates the elaborate melting phenomenon for dsDNA in sub-molecular scale.

Force-Extension for DNA in a Nanoslit: Mapping between the 3D and 2D Limits

The force-extension relation for a semiflexible polymer confined in a nanoslit is investigated. Both the effective correlation length and force-extension relation change as the chain goes from 3D (large slit heights) to 2D (tight confinement). At low forces, correlations along the polymer give an effective dimensionality. The strong force limit can be interpolated with the weak force limit for two regimes: when confinement dominates over extensile force and vice versa. These interpolations give good agreement with simulations for all slit heights and forces. We thus generalize the Marko-Siggia force-extension relation for DNA and other semiflexible biopolymers in nanoconfinement.

Stripe to slab confinement for the linearization of macromolecules in nanochannels

We investigated the recently suggested advantageous analysis of chain linearization experiments with macromolecules confined in a stripe-like channel (Huang and Battacharya, EPL, 2014, 106, 18004) using Monte Carlo simulations. The enhanced chain extension in a stripe, which is due to the significant excluded volume interactions between the monomers in two dimensions, weakens considerably on transition to an experimentally feasible slit-like channel. Based on the chain extension-confinement strength dependence and the structure factor behavior for a chain in a stripe, we infer the excluded volume regime (de Gennes regime) typical for two-dimensional systems. On widening of the stripe in a direction perpendicular to the stripe plane, i.e. on the transition to the slab geometry, the advantageous chain extension decreases and a Gaussian regime is observed for not very long semiflexible chains. The evidence for pseudo-ideality in confined chains is based on four indicators: the extension curves, variation of the extension with the persistence length P, estimated limits for the regimes in the investigated systems, and the structure factor behavior. The slab behavior can be observed when the two-dimensional stripe (originally of a one-monomer thickness) reaches a reduced thickness D larger than approximately D/P ≈ 0.2 in the third dimension. This maximum height of a slab at which the advantage of a stripe is retained is very low and has implications for DNA linearization experiments.

Nanoslit Confined DNA at Low Ionic Strengths

We present nanoslit confined DNA conformations at very low ionic strengths and a theory to explain most measurements for single DNA molecule size under strong nanoslit confinement. Very low ionic strength conditions not only increase the DNA persistence length dramatically, but also cause DNA molecules to swell to the extent that the effective diameter of DNA becomes larger than the nanoslit height. By accounting for these effects, our results and theory provide a reasonable clue for a current controversy regarding the dependence of the DNA conformation on slit height (h), persistence length (p), and effective diameter (w).

Probing transient protein-mediated DNA linkages using nanoconfinement

We present an analytic technique for probing protein-catalyzed transient DNA loops that is based on nanofluidic channels. In these nanochannels, DNA is forced in a linear configuration that makes loops appear as folds whose size can easily be quantified. Using this technique, we study the interaction between T4 DNA ligase and DNA. We find that T4 DNA ligase binding changes the physical characteristics of the DNApolymer, in particular persistence length and effective width. We find that the rate of DNA fold unrolling is significantly reduced when T4 DNA ligase and ATP are applied to bare DNA. Together with evidence of T4 DNA ligase bridging two different segments of DNA based on AFM imaging, we thus conclude that ligase can transiently stabilize folded DNA configurations by coordinating genetically distant DNA stretches.

Extension of DNA in a nanochannel as a rod-to-coil transition

DNA confinement in nanochannels is emerging as an important tool for genomics and an excellent platform for testing the theories of confined wormlike polymers. Using cutting-edge, large scale Monte Carlo simulations of asymptotically long wormlike chains, we show that, in analogy to the rod-to-coil transition for free wormlike polymers, there exists a universal, Gauss-de Gennes regime that connects the classic Odijk and de Gennes regimes of channel-confined chains. For DNA in a nanochannel, this Gauss-de Gennes regime spans practically the entire experimentally relevant range of channel sizes, including the nanochannels used in an incipient genome mapping technology.

Entropic depletion of DNA in triangular nanochannels

Using Monte Carlo simulations of a touching-bead model of double-stranded DNA, we show that DNA extension is enhanced in isosceles triangular nanochannels (relative to a circular nanochannel of the same effective size) due to entropic depletion in the channel corners. The extent of the enhanced extension depends non-monotonically on both the accessible area of the nanochannel and the apex angle of the triangle. We also develop a metric to quantify the extent of entropic depletion, thereby collapsing the extension data for circular, square, and various triangular nanochannels onto a single master curve for channel sizes in the transition between the Odijk and de Gennes regimes.

Conformation Model of Back-Folding and Looping of a Single DNA Molecule Confined Inside a Nanochannel

Currently, no theory is available to describe the conformation of DNA confined in a channel when the nanochannel diameter is around the persistence length. Back-folded hairpins in the undulating wormlike chain conformation result in the formation of loops, which reduces the stretch of the molecule in the longitudinal direction of the channel. A cooperativity model is used to quantify the frequency and size of the loop domains. The predictions agree with results from the Monte Carlo simulation.