Conformations of DNA in Triangular Nanochannels

Crowding-induced polymer trapping in a channel

In this work, we report an intriguing phenomenon: crowding-induced polymer trapping in a channel. Using Langevin dynamics simulations and analytical calculations, we find that for a polymer confined in a channel, crowding particles can push a polymer into the channel corner through inducing an effective polymer-corner attraction due to the depletion effect. This phenomenon is referred to as polymer trapping. The occurrence of polymer trapping requires a minimum volume fraction of crowders, ϕ^{*}, which scales as ϕ^{*}∼(a_{c}/L_{p})^{1/3} for a_{c}≫a_{m} and ϕ^{*}∼(a_{c}/L_{p})^{1/3}(a_{c}/a_{m})^{1/2} for a_{c}≪a_{m}, where a_{c} is the crowder diameter, a_{m} is the monomer diameter, and L_{p} is the polymer persistence length. For DNA, ϕ^{*} is estimated to be around 0.25 for crowders with a_{c}=2nm. We find that ϕ^{*} also strongly depends on the shape of the channel cross section, and ϕ^{*} is much smaller for a triangle channel than a square channel. The polymer trapping leads to a nearly fully stretched polymer conformation along a channel corner, which may have practical applications, such as full stretching of DNA for the nanochannel-based genome mapping technology.

Stretching Wormlike Chains in Narrow Tubes of Arbitrary Cross-Sections

We considered the stretching of semiflexible polymer chains confined in narrow tubes with arbitrary cross-sections. Based on the wormlike chain model and technique of normal mode decomposition in statistical physics, we derived a compact analytical expression on the force-confinement-extension relation of the chains. This single formula was generalized to be valid for tube confinements with arbitrary cross-sections. In addition, we extended the generalized bead-rod model for Brownian dynamics simulations of confined polymer chains subjected to force stretching, so that the confinement effects to the chains applied by the tubes with arbitrary cross-sections can be quantitatively taken into account through numerical simulations. Extensive simulation examples on the wormlike chains confined in tubes of various shapes quantitatively justified the theoretically derived generalized formula on the force-confinement-extension relation of the chains.

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.

Tunable Confinement for Bridging Single-Cell Manipulation and Single-Molecule DNA Linearization

DNA linearization by nanoconfinement has offered a new avenue toward large-scale genome mapping. The ability to smoothly interface the widely different length scales from cell manipulation to DNA linearization is critical to the development of single-cell genomic mapping or sequencing technologies. Conventional nanochannel technologies for DNA analysis suffer from complex fabrication procedures, DNA stacking at the nanochannel entrance, and inefficient solution exchange. In this work, a dynamic and tunable confinement strategy is developed to manipulate and linearize genomic-length DNA molecules from a single cell. By leveraging pneumatic microvalve control and elastomeric collapse, an array of nanochannels with confining dimension down to 20 nm and length up to sub-millimeter is created and can be dynamically tuned in size. The curved edges of the microvalve form gradual transitions from microscale to nanoscale confinement, smoothly facilitating DNA entry into the nanochannels. A unified micro/nanofluidic device that integrates single-cell trapping and lysis, DNA extraction, purification, labeling, and linearization is developed based on dynamically controllable nanochannels. Mbp-long DNA molecules are extracted directly from a single cell and in situ linearized in the nanochannels. The device provides a facile and promising platform to achieve the ultimate goal of single-cell, single-genome analysis.

Spatial confinement induces hairpins in nicked circular DNA

In living cells, DNA is highly confined in space with the help of condensing agents, DNA binding proteins and high levels of supercoiling. Due to challenges associated with experimentally studying DNA under confinement, little is known about the impact of spatial confinement on the local structure of the DNA. Here, we have used well characterized slits of different sizes to collect high resolution atomic force microscopy images of confined circular DNA with the aim of assessing the impact of the spatial confinement on global and local conformational properties of DNA. Our findings, supported by numerical simulations, indicate that confinement imposes a large mechanical stress on the DNA as evidenced by a pronounced anisotropy and tangent-tangent correlation function with respect to non-constrained DNA. For the strongest confinement we observed nanometer sized hairpins and interwound structures associated with the nicked sites in the DNA sequence. Based on these findings, we propose that spatial DNA confinement in vivo can promote the formation of localized defects at mechanically weak sites that could be co-opted for biological regulatory functions.

Dynamic simulations show repeated narrowing maximizes DNA linearization in elastomeric nanochannels

This paper uses computer simulations to reveal unprecedented details about linearization of deoxyribonucleic acid (DNA) inside dynamic nanochannels that can be repeatedly widened and narrowed. We first analyze the effect of rate of channel narrowing on DNA linearization dynamics. Quick (∼0.1 s) narrowing of nanoscale channels results in rapid overstretching of the semi-flexible chain followed by a slower (∼0.1-10 s) relaxation to an equilibrium extension. Two phenomena that induce linearization during channel narrowing, namely, elongational-flow and confinement, occur simultaneously, regardless of narrowing speed. Interestingly, although elongational flow is a minimum at the mid-point of the channel and increases towards the two ends, neither the linearization dynamics nor the degree of DNA extension varies significantly with the center-of-mass of the polymer projected on the channel axis. We also noticed that there was a significant difference in time to reach the equilibrium length, as well as the degree of DNA linearization at short times, depending on the initial conformation of the biopolymer. Based on these observations, we tested a novel linearization protocol where the channels are narrowed and widened repeatedly, allowing DNA to explore multiple conformations. Repeated narrowing and widening, something uniquely enabled by the elastomeric nanochannels, significantly decrease the time to reach the equilibrium-level of stretch when performed within periods comparable to the chain relaxation time and more effectively untangle chains into more linearized biopolymers.

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.

Inverted critical adsorption of polyelectrolytes in confinement

What are the fundamental laws for the adsorption of charged polymers onto oppositely charged surfaces, for convex, planar, and concave geometries? This question is at the heart of surface coating applications, various complex formation phenomena, as well as in the context of cellular and viral biophysics. It has been a long-standing challenge in theoretical polymer physics; for realistic systems the quantitative understanding is however often achievable only by computer simulations. In this study, we present the findings of such extensive Monte-Carlo in silico experiments for polymer-surface adsorption in confined domains. We study the inverted critical adsorption of finite-length polyelectrolytes in three fundamental geometries: planar slit, cylindrical pore, and spherical cavity. The scaling relations extracted from simulations for the critical surface charge density σc-defining the adsorption-desorption transition-are in excellent agreement with our analytical calculations based on the ground-state analysis of the Edwards equation. In particular, we confirm the magnitude and scaling of σc for the concave interfaces versus the Debye screening length 1/κ and the extent of confinement a for these three interfaces for small κa values. For large κa the critical adsorption condition approaches the known planar limit. The transition between the two regimes takes place when the radius of surface curvature or half of the slit thickness a is of the order of 1/κ. We also rationalize how σc(κ) dependence gets modified for semi-flexible versus flexible chains under external confinement. We examine the implications of the chain length for critical adsorption-the effect often hard to tackle theoretically-putting an emphasis on polymers inside attractive spherical cavities. The applications of our findings to some biological systems are discussed, for instance the adsorption of nucleic acids onto the inner surfaces of cylindrical and spherical viral capsids.

Hydrodynamics of DNA confined in nanoslits and nanochannels

Modeling the dynamics of a confined, semi exible polymer is a challenging problem, owing to the complicated interplay between the configurations of the chain, which are strongly affected by the length scale for the confinement relative to the persistence length of the chain, and the polymer-wall hydrodynamic interactions. At the same time, understanding these dynamics are crucial to the advancement of emerging genomic technologies that use confinement to stretch out DNA and "read" a genomic signature. In this mini-review, we begin by considering what is known experimentally and theoretically about the friction of a wormlike chain such as DNA confined in a slit or a channel. We then discuss how to estimate the friction coefficient of such a chain, either with dynamic simulations or via Monte Carlo sampling and the Kirk-wood pre-averaging approximation. We then review our recent work on computing the diffusivity of DNA in nanoslits and nanochannels, and conclude with some promising avenues for future work and caveats about our approach.

Stretching of DNA confined in nanochannels with charged walls

There is currently a growing interest in control of stretching of DNA inside nanoconfined regions due to the possibility to analyze and manipulate single biomolecules for applications such as DNA mapping and barcoding, which are based on stretching the DNA in a linear fashion. In the present work, we couple Finite Element Methods and Monte Carlo simulations in order to study the conformation of DNA molecules confined in nanofluidic channels with neutral and charged walls. We find that the electrostatic forces become more and more important when lowering the ionic strength of the solution. The influence of the nanochannel cross section geometry is also studied by evaluating the DNA elongation in square, rectangular, and triangular channels. We demonstrate that coupling electrostatically interacting walls with a triangular geometry is an efficient way to stretch DNA molecules at the scale of hundreds of nanometers. The paper reports experimental observations of λ-DNA molecules in poly(dimethylsiloxane) nanochannels filled with solutions of different ionic strength. The results are in good agreement with the theoretical predictions, confirming the crucial role of the electrostatic repulsion of the constraining walls on the molecule stretching.

Driven polymer transport through a periodically patterned channel

We study the driven transport of polymers in a periodically patterned channel using Langevin dynamics simulations in two dimensions. The channel walls are patterned with periodically alternating patches of attractive and non-attractive particles that act as trapping sites for the polymer. We find that the system shows rich dynamical behavior, observing giant diffusion, negative differential mobility, and several different transition mechanisms between the attractive patches. We also show that the channel can act as an efficient high-pass filter for polymers longer than a threshold length Nthr, which can be tuned by adjusting the length of the attractive patches and the driving force. Our findings suggest the possibility of fabricating polymer filtration devices based on patterned nanochannels.

Modeling the relaxation time of DNA confined in a nanochannel

Using a mapping between a Rouse dumbbell model and fine-grained Monte Carlo simulations, we have computed the relaxation time of λ-DNA in a high ionic strength buffer confined in a nanochannel. The relaxation time thus obtained agrees quantitatively with experimental data [Reisner et al., Phys. Rev. Lett. 94, 196101 (2005)] using only a single O(1) fitting parameter to account for the uncertainty in model parameters. In addition to validating our mapping, this agreement supports our previous estimates of the friction coefficient of DNA confined in a nanochannel [Tree et al., Phys. Rev. Lett. 108, 228105 (2012)], which have been difficult to validate due to the lack of direct experimental data. Furthermore, the model calculation shows that as the channel size passes below approximately 100 nm (or roughly the Kuhn length of DNA) there is a dramatic drop in the relaxation time. Inasmuch as the chain friction rises with decreasing channel size, the reduction in the relaxation time can be solely attributed to the sharp decline in the fluctuations of the chain extension. Practically, the low variance in the observed DNA extension in such small channels has important implications for genome mapping.