Segregation of semiflexible macromolecules in nanochannel

Equilibrium behaviour of two cavity-confined polymers: effects of polymer width and system asymmetries

Experiments using nanofluidic devices have proven effective in characterizing the physical properties of polymers confined to small cavities. Two recent studies using such methods examined the organization and dynamics of two DNA molecules in box-like cavities with strong confinement in one direction and with square and elliptical cross sections in the lateral plane. Motivated by these experiments, we employ Monte Carlo and Brownian dynamics simulations to study the physical behaviour of two polymers confined to small cavities with shapes comparable to those used in the experiments. We quantify the effects of varying the following polymer properties and confinement dimensions on the organization and dynamics of the polymers: the polymer width, the polymer contour length ratio, the cavity cross-sectional area, and the degree of cavity elongation for cavities with rectangular and elliptical cross sections. We find that the tendency for polymers to segregate is enhanced by increasing polymer width. For sufficiently small cavities, increasing cavity elongation promotes segregation and localization of identical polymers to opposite sides of the cavity along its long axis. A free-energy barrier controls the rate of polymers swapping positions, and the observed dynamics are roughly in accord with predictions of a simple theoretical model. Increasing the contour length difference between polymers significantly affects their organization in the cavity. In the case of a large linear polymer co-trapped with a small ring polymer in an elliptical cavity, the small polymer tends to lie near the lateral confining walls, and especially at the cavity poles for highly elongated ellipses.

Knot Factories with Helical Geometry Enhance Knotting and Induce Handedness to Knots

We performed molecular dynamics simulations of DNA polymer chains confined in helical nano-channels under compression in order to explore the potential of knot-factories with helical geometry to produce knots with a preferred handedness. In our simulations, we explore mutual effect of the confinement strength and compressive forces in a range covering weak, intermediate and strong confinement together with weak and strong compressive forces. The results find that while the common metrics of polymer chain in cylindrical and helical channels are very similar, the DNA in helical channels exhibits greatly different topology in terms of chain knottedness, writhe and handedness of knots. The results show that knots with a preferred chirality in terms of average writhe can be produced by using channels with a chosen handedness.

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

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.

Equilibrium organization, conformation, and dynamics of two polymers under box-like confinement

Motivated by recent nanofluidics experiments, we use Brownian dynamics and Monte Carlo simulations to study the conformation, organization and dynamics of two polymer chains confined to a single box-like cavity. The polymers are modeled as flexible bead-spring chains, and the box has a square cross-section of side length L and a height that is small enough to compress the polymers in that dimension. For sufficiently large L, the system behaviour approaches that of an isolated polymer in a slit. However, the combined effects of crowding and confinement on the polymer organization, conformation and equilibrium dynamics become significant when where is the transverse radius of gyration for a slit geometry. In this regime, the centre-of-mass probability distribution in the transverse plane exhibits a depletion zone near the centre of the cavity (except at very small L) and a 4-fold symmetry with quasi-discrete positions. Reduction in polymer size with decreasing L arises principally from confinement rather than inter-polymer crowding. By contrast, polymer diffusion and internal motion are strongly affected by inter-polymer crowding. The two polymers tend to occupy opposite positions relative to the box centre, about which they diffuse relatively freely. Qualitatively, this static and dynamical behaviour differs significantly from that previously observed for confinement of two polymers to a narrow channel. The simulation results for a suitably chosen box width are qualitatively consistent with results from a recent experimental study of two λ-DNA chains confined to a nanofluidic cavity.

Free energy and segregation dynamics of two channel-confined polymers of different lengths

Polymers confined to a narrow channel are subject to strong entropic forces that tend to drive the molecules apart. In this study, we use Monte Carlo computer simulations to study the segregation behavior of two flexible hard-sphere polymers under confinement in a cylindrical channel. We focus on the effects of using polymers of different lengths. We measure the variation of the conformational free energy, F, with the center-of-mass separation distance, λ. The simulations reveal four different separation regimes, characterized by different scaling properties of the free energy with respect to the polymer lengths and the channel diameter, D. We propose a regime map in which the state of the system is determined by the values of the quantities N_{2}/N_{1} and λ/(N_{1}+N_{2})D^{-β}, where N_{1} and N_{2} are the polymer lengths, and where β≈0.64. The observed scaling behavior of F(λ) in each regime is in reasonable agreement with predictions using a simple theoretical model. In addition, we use MC dynamics simulations to study the segregation dynamics of initially overlapping polymers by measurement of the incremental mean first-passage time with respect to λ. For systems characterized by a wide range of λ in which a short polymer is nested within a longer one, the segregation dynamics are close to that expected for two noninteracting one-dimensional random walkers undergoing unbiased diffusion. When the free-energy gradient is large, segregation is rapid and characterized by out-of-equilibrium effects.

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.

Molecular Dynamics Simulation of Supercoiled, Knotted, and Catenated DNA Molecules, Including Modeling of Action of DNA Gyrase

A detailed protocol of molecular dynamics simulations of supercoiled DNA molecules that can be in addition knotted or catenated is described. We also describe how to model ongoing action of DNA gyrase that introduces negative supercoing into DNA molecules. The protocols provide detailed instructions about model parameters, equations of used potentials, simulation, and visualization. Implementation of the model into a frequently used molecular dynamics simulation environment, ESPResSo, is shown step by step.

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.

Manipulation of a neutral and nonpolar nanoparticle in water using a nonuniform electric field

The manipulation of nanoparticles in water is of essential importance in chemical physics, nanotechnology, medical technology, and biotechnology applications. Generally, a particle with net charges or charge polarity can be driven by an electric field. However, many practical particles only have weak and even negligible charge and polarity, which hinders the electric field to exert a force large enough to drive these nanoparticles directly. Here, we use molecular dynamics simulations to show that a neutral and nonpolar nanoparticle in liquid water can be driven directionally by an external electric field. The directed motion benefits from a nonuniform water environment produced by a nonuniform external electric field, since lower water energies exist under a higher intensity electric field. The nanoparticle spontaneously moves toward locations with a weaker electric field intensity to minimize the energy of the whole system. Considering that the distance between adjacent regions of nonuniform field intensity can reach the micrometer scale, this finding provides a new mechanism of manipulating nanoparticles from the nanoscale to the microscale.

Generation of supercoils in nicked and gapped DNA drives DNA unknotting and postreplicative decatenation

Due to the helical structure of DNA the process of DNA replication is topologically complex. Freshly replicated DNA molecules are catenated with each other and are frequently knotted. For proper functioning of DNA it is necessary to remove all of these entanglements. This is done by DNA topoisomerases that pass DNA segments through each other. However, it has been a riddle how DNA topoisomerases select the sites of their action. In highly crowded DNA in living cells random passages between contacting segments would only increase the extent of entanglement. Using molecular dynamics simulations we observed that in actively supercoiled DNA molecules the entanglements resulting from DNA knotting or catenation spontaneously approach sites of nicks and gaps in the DNA. Type I topoisomerases, that preferentially act at sites of nick and gaps, are thus naturally provided with DNA–DNA juxtapositions where a passage results in an error-free DNA unknotting or DNA decatenation.

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.

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.