Simulation of DNA Extension in Nanochannels

Evaluation of the Kirkwood approximation for the diffusivity of channel-confined DNA chains in the de Gennes regime

We use Brownian dynamics with hydrodynamic interactions to calculate both the Kirkwood (short-time) diffusivity and the long-time diffusivity of DNA chains from free solution down to channel confinement in the de Gennes regime. The Kirkwood diffusivity in confinement is always higher than the diffusivity obtained from the mean-squared displacement of the center-of-mass, as is the case in free solution. Moreover, the divergence of the local diffusion tensor, which is non-zero in confinement, makes a negligible contribution to the latter diffusivity in confinement. The maximum error in the Kirkwood approximation in our simulations is about 2% for experimentally relevant simulation times. The error decreases with increasing confinement, consistent with arguments from blob theory and the molecular-weight dependence of the error in free solution. In light of the typical experimental errors in measuring the properties of channel-confined DNA, our results suggest that the Kirkwood approximation is sufficiently accurate to model experimental data.

Distribution of distances between DNA barcode labels in nanochannels close to the persistence length

We obtained experimental extension data for barcoded E. coli genomic DNA molecules confined in nanochannels from 40 nm to 51 nm in width. The resulting data set consists of 1 627 779 measurements of the distance between fluorescent probes on 25 407 individual molecules. The probability density for the extension between labels is negatively skewed, and the magnitude of the skewness is relatively insensitive to the distance between labels. The two Odijk theories for DNA confinement bracket the mean extension and its variance, consistent with the scaling arguments underlying the theories. We also find that a harmonic approximation to the free energy, obtained directly from the probability density for the distance between barcode labels, leads to substantial quantitative error in the variance of the extension data. These results suggest that a theory for DNA confinement in such channels must account for the anharmonic nature of the free energy as a function of chain extension.

Using a Péclet number for the translocation of a polymer through a nanopore to tune coarse-grained simulations to experimental conditions

Coarse-grained simulations are often employed to study the translocation of DNA through a nanopore. The majority of these studies investigate the translocation process in a relatively generic sense and do not endeavor to match any particular set of experimental conditions. In this manuscript, we use the concept of a Péclet number for translocation, P(t), to compare the drift-diffusion balance in a typical experiment vs a typical simulation. We find that the standard coarse-grained approach overestimates diffusion effects by anywhere from a factor of 5 to 50 compared to experimental conditions using double stranded DNA (dsDNA). By defining a Péclet control parameter, λ, we are able to correct this and tune the simulations to replicate the experimental P(t) (for dsDNA and other scenarios). To show the effect that a particular P(t) can have on the dynamics of translocation, we perform simulations across a wide range of P(t) values for two different types of driving forces: a force applied in the pore and a pulling force applied to the end of the polymer. As P(t) brings the system from a diffusion dominated to a drift dominated regime, a variety of effects are observed including a non-monotonic dependence of the translocation time τ on P(t) and a steep rise in the probability of translocating. Comparing the two force cases illustrates the impact of the crowding effects that occur on the trans side: a non-monotonic dependence of the width of the τ distributions is obtained for the in-pore force but not for the pulling force.

Dynamic compression of single nanochannel confined DNA via a nanodozer assay

We show that a single DNA molecule confined and extended in a nanochannel can be dynamically compressed by sliding a permeable gasket at a fixed velocity relative to the stationary polymer. The gasket is realized experimentally by optically trapping a nanosphere inside a nanochannel. The trapped bead acts like a "nanodozer," directly applying compressive forces to the molecule without requirement of chemical attachment. Remarkably, these strongly nonequilibrium measurements can be quantified via a simple nonlinear convective-diffusion formalism and yield insights into the local blob statistics, allowing us to conclude that the compressed nanochannel-confined chain exhibits mean-field behavior.

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.

Knotting and Unknotting Dynamics of DNA Strands in Nanochannels

The self-knotting dynamics of DNA strands confined in nanochannels is studied with Brownian simulations. The model DNA chains are several microns long and placed inside channels that are 50-300 nm wide. This width range covers the transition between different metric scaling regimes and the concomitant drop of DNA knotting probability for channel widths below ∼75 nm. We find that knots typically originate from deep looping and backfoldings of the chain ends. Upon lowering the channel width, backfoldings become shallower and rarer and the lifetime of knots decreases while that of unknots increases. This lifetimes interplay causes the dramatic reduction of knots incidence for increasing confinement. The results can aid the design of nanochannels capable of harnessing the self-knotting dynamics to quench or relax the DNA topological state as desired.

Free energy of a long semiflexible polymer confined in a spherical cavity

The free energy and conformational properties of a wormlike chain confined inside a spherical surface are investigated. We show that in the weak-confinement limit, the wormlike chain model exactly reproduces the confinement properties of a Gaussian chain; in such a case the confinement entropy dominates the free energy; in the strong-confinement limit, the free energy is dominated by the bending energy of the chain, which is forced to wrap around the confining surface. We also present a numerical solution within the crossover region between the two limits, solving the differential equation that the probability distribution function satisfies.

Mixed confinement regimes during equilibrium confinement spectroscopy of DNA

We have used a combination of fluorescence microscopy experiments and Pruned Enriched Rosenbluth Method simulations of a discrete wormlike chain model to measure the mean extension and the variance in the mean extension of λ-DNA in 100 nm deep nanochannels with widths ranging from 100 nm to 1000 nm in discrete 100 nm steps. The mean extension is only weakly affected by the channel aspect ratio. In contrast, the fluctuations of the chain extension qualitatively differ between rectangular channels and square channels with the same cross-sectional area, owing to the "mixing" of different confinement regimes in the rectangular channels. The agreement between experiment and simulation is very good, using the extension due to intercalation as the only adjustable parameter.

Localization and stretching of polymer chains at the junction of two surfaces

We present a molecular dynamics study on the stretching of a linear polymer chain that is adsorbed at the junction of two intersecting flat surfaces of varying alignments. We observe a transition from a two-dimensional to one-dimensional (1D) structure of the adsorbed polymer when the alignment, i.e., the angle between the two surfaces that form a groove, θ, is below 135°. We show that the radius of gyration of the polymer chain Rg scales as Rg ∼ N(3/4) with the degree of polymerization N for θ = 180° (planer substrate), and the scaling changes to Rg ∼ N(1.0) for θ < 135° in good solvents. At the crossover point, θ = 135°, the exponent becomes 1.15. The 1D stretching of the polymer chain is found to be 84% of its contour length for θ ⩽ 90°. The center of mass diffusion coefficient D decreases sharply with θ. However, the diffusion coefficient scales with N as D ∼ N(-1), and is independent of θ. The relaxation time τ, for the diffusive motion, scales as τ ∼ N(2.5) for θ = 180° (planar substrate), which changes to τ ∼ N(3.0) for θ ⩽ 90°. At the crossover point, the exponent is 3.4, which is slightly higher than the 1D value of 3.0. Further, a signature of reptation-like dynamics of the polymer chain is observed at the junction for θ ⩽ 90° due to its strong 1D localization and stretching.

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.

DNA translocation through short nanofluidic channels under asymmetric pulsed electric field

Investigation of single molecule DNA dynamics in confined environments has led to important applications in DNA analysis, separation, and sequencing. Here, we studied the electrophoretic transport of DNA molecules through nanochannels shorter than the DNA contour length and calculated the associated translocation time curves. We found that the longer T4 DNA molecules required a longer time to traverse a fixed length nanochannel than shorter λ DNA molecules and that the translocation time decreased with increasing electric field which agreed with theoretical predictions. We applied this knowledge to design an asymmetric electric pulse and demonstrate the different responses of λ and T4 DNA to the pulses. We used Brownian dynamics simulations to corroborate our experimental results on DNA translocation behaviour. This work contributes to the fundamental understanding of polymer transport through nanochannels and may help in designing better separation techniques in the future.

Interplay between chain stiffness and excluded volume of semiflexible polymers confined in nanochannels

The properties of channel-confined semiflexible polymers are determined by a complicated interplay of chain stiffness and excluded volume effects. Using Pruned-Enriched Rosenbluth Method (PERM) simulations, we study the equilibrium properties of channel-confined polymers by systematically controlling chain stiffness and excluded volume. Our calculations of chain extension and confinement free energy for freely jointed chains with and without excluded volume show excellent agreement with theoretical predictions. For ideal wormlike chains, the extension is seen to crossover from Odijk behavior in strong confinement to zero-stretching, bulk-like behavior in weak confinement. In contrast, for self-avoiding wormlike chains, we always observe that the linear scaling of the extension with the contour length is valid in the long-chain limit irrespective of the regime of confinement, owing to the coexistence of stiffness and excluded volume effects. We further propose that the long-chain limit for the extension corresponds to chain lengths wherein the projection of the end-to-end distance along the axis of the channel is nearly equal to the mean span parallel to the axis. For DNA in nanochannels, this limit was identified using PERM simulations out to molecular weights of more than 1 megabase pairs; the molecular weight of λ-DNA is found to exhibit nearly asymptotic fractional extension for channels sizes used commonly in experiments.

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.

Is DNA a Good Model Polymer?

The details surrounding the cross-over from wormlike-specific to universal polymeric behavior has been the subject of debate and confusion even for the simple case of a dilute, unconfined wormlike chain. We have directly computed the polymer size, form factor, free energy and Kirkwood diffusivity for unconfined wormlike chains as a function of molecular weight, focusing on persistence lengths and effective widths that represent single-stranded and double-stranded DNA in a high ionic strength buffer. To do so, we use a chain-growth Monte Carlo algorithm, the Pruned-Enriched Rosenbluth Method (PERM), which allows us to estimate equilibrium and near-equilibrium dynamic properties of wormlike chains over an extremely large range of contour lengths. From our calculations, we find that very large DNA chains (≈ 1,000,000 base pairs depending on the choice of size metric) are required to reach flexible, swollen non-draining coils. Furthermore, our results indicate that the commonly used model polymer λ-DNA (48,500 base pairs) does not exhibit "ideal" scaling, but exists in the middle of the transition to long-chain behavior. We subsequently conclude that typical DNA used in experiments are too short to serve as an accurate model of long-chain, universal polymer behavior.

Presentation of large DNA molecules for analysis as nanoconfined dumbbells

The analysis of very large DNA molecules intrinsically supports long-range, phased sequence information, but requires new approaches for their effective presentation as part of any genome analysis platform. Using a multi-pronged approach that marshaled molecular confinement, ionic environment, and DNA elastic properties-but tressed by molecular simulations-we have developed an efficient and scalable approach for presentation of large DNA molecules within nanoscale slits. Our approach relies on the formation of DNA dumbbells, where large segments of the molecules remain outside the nanoslits used to confine them. The low ionic environment, synergizing other features of our approach, enables DNA molecules to adopt a fully stretched conformation, comparable to the contour length, thereby facilitating analysis by optical microscopy. Accordingly, a molecular model is proposed to describe the conformation and dynamics of the DNA molecules within the nanoslits; a Langevin description of the polymer dynamics is adopted in which hydrodynamic effects are included through a Green's function formalism. Our simulations reveal that a delicate balance between electrostatic and hydrodynamic interactions is responsible for the observed molecular conformations. We demonstrate and further confirm that the "Odijk regime" does indeed start when the confinement dimensions size are of the same order of magnitude as the persistence length of the molecule. We also summarize current theories concerning dumbbell dynamics.

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.

Revisiting blob theory for DNA diffusivity in slitlike confinement

Blob theory has been widely applied to describe polymer conformations and dynamics in nanoconfinement. In slit confinement, blob theory predicts a scaling exponent of 2/3 for polymer diffusivity as a function of slit height, yet a large body of experimental studies using DNA produce a scaling exponent significantly less than 2/3. In this work, we develop a theory that predicts that this discrepancy occurs because the segment correlation function for a semiflexible chain such as DNA does not follow the Flory exponent for length scales smaller than the persistence length. We show that these short length scale effects contribute significantly to the scaling for the DNA diffusivity, but do not appreciably affect the scalings for static properties. Our theory is fully supported by Monte Carlo simulations, quantitative agreement with DNA experiments, and the results reconcile this outstanding problem for confined polymers.

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.

Electro-entropic excluded volume effects on DNA looping and relaxation in nanochannels

We investigate the fluctuation-relaxation dynamics of entropically restricted DNA molecules in square nanochannels ranging from 0.09 to 19.9 times the persistence length. In nanochannels smaller than the persistence length, the chain relaxation time is found to have cubic dependence on the channel size. It is found that the effective polymer width significantly alter the chain conformation and relaxation time in strong confinement. For thinner chains, looped chain configurations are found in channels with height comparable to the persistence length, with very slow relaxation compared to un-looped chains. Larger effective chain widths inhibit the formation of hairpin loops.

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.

Electrokinetically-driven transport of DNA through focused ion beam milled nanofluidic channels

The electrophoretically driven transport of double-stranded λ-phage DNA through focused ion beam (FIB) milled nanochannels is described. Nanochannels were fabricated having critical dimensions (width and depth) corresponding to 0.5×, 1×, and 2× the DNA persistence length, or 25 nm, 50 nm, and 100 nm, respectively. The threshold field strength required to drive transport, the threading mobility, and the transport mobility were measured as a function of nanochannel size. As the nanochannel dimensions decreased, the entropic barrier to translocation increased and transport became more constrained. Equilibrium models of confinement provide a framework in which to understand the observed trends, although the dynamic nature of the experiments resulted in significant deviations from theory. It was also demonstrated that the use of dynamic wall coatings for the purpose of electroosmotic flow suppression can have a significant impact on transport dynamics that may obfuscate entropic contributions. The nonintermittent DNA transport through the FIB milled nanochannels demonstrates that they are well suited for use in nanofluidic devices. We expect that an understanding of the dynamic transport properties reported here will facilitate the incorporation of FIB-milled nanochannels in devices for single molecule and ensemble analyses.

Orientational correlations in confined DNA

We study how the orientational correlations of DNA confined to nanochannels depend on the channel diameter D by means of Monte Carlo simulations and a mean-field theory. This theory describes DNA conformations in the experimentally relevant regime where the Flory-de Gennes theory does not apply. We show how local correlations determine the dependence of the end-to-end distance of the DNA molecule upon D. Tapered nanochannels provide the necessary resolution in D to study experimentally how the extension of confined DNA molecules depends upon D. Our experimental and theoretical results are in qualitative agreement.

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.

DNA conformation in nanochannels: Monte Carlo simulation studies using a primitive DNA model

We have performed canonical ensemble Monte Carlo simulations of a primitive DNA model to study the conformation of 2.56 ~ 21.8 μm long DNA molecules confined in nanochannels at various ionic concentrations with the comparison of our previous experimental findings. In the model, the DNA molecule is represented as a chain of charged hard spheres connected by fixed bond length and the nanochannels as planar hard walls. System potentials consist of explicit electrostatic potential along with short-ranged hard-sphere and angle potentials. Our primitive model system provides valuable insight into the DNA conformation, which cannot be easily obtained from experiments or theories. First, the visualization and statistical analysis of DNA molecules in various channel dimensions and ionic strengths verified the formation of locally coiled structures such as backfolding or hairpin and their significance even in highly stretched states. Although the folding events mostly occur within the region of ~0.5 μm from both chain ends, significant portion of the events still take place in the middle region. Second, our study also showed that two controlling factors such as channel dimension and ionic strength widely used in stretching DNA molecules have different influence on the local DNA structure. Ionic strength changes local correlation between neighboring monomers by controlling the strength of electrostatic interaction (and thus the persistence length of DNA), which leads to more coiled local conformation. On the other hand, channel dimension controls the overall stretch by applying the geometric constraint to the non-local DNA conformation instead of directly affecting local correlation. Third, the molecular weight dependence of DNA stretch was observed especially in low stretch regime, which is mainly due to the fact that low stretch modes observed in short DNA molecules are not readily accessible to much longer DNA molecules, resulting in the increase in the stretch of longer DNA molecules.

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.

Effect of Nanoslit Confinement on the Knotting Probability of Circular DNA

Monte Carlo simulations are used to study the knotting probability of circular DNA confined in a slit. We systematically vary the slit height, the width, and the contour length of the DNA molecule. We find that the trend in knotting probability with respect to slit height can be monotonic or nonmonotonic, depending on the width and contour length. The nonmonotonic trend is caused by two competing factors: the increase of the effective persistence length and the increase of segment density by slit confinement. These factors are antagonistic, in the sense that the increase in effective persistence length disfavors knot formation, whereas the increase in segment density favors the knotting probability. Our simulation results bring to light the importance of both chain length and width for slit-confined circular DNA and can be used to guide future experiments which aim to produce populations of knotted DNA through cyclization or catalyzed double-strand passage reactions in confinement.

Mobility of a semiflexible chain confined in a nanochannel

The classic results of de Gennes and Odijk describe the mobility of a semiflexible chain confined in a nanochannel only in the limits of very weak and very strong confinement, respectively. Using Monte Carlo sampling of the Kirkwood diffusivity with full hydrodynamic interactions, we show that the mobility of a semiflexible chain exhibits a broad plateau as a function of extension before transitioning to an Odijk regime, and that the width of the plateau depends on the anisotropy of the monomers. For the particular case of DNA in a high ionic strength buffer, which has highly anisotropic monomers, we predict that this Rouse-like behavior will be observed over most of the measurable chain extensions seen in experiments.

Resolution limit for DNA barcodes in the Odijk regime

We develop an approximation for the probability of optically resolving two fluorescent labels on the backbone of a DNA molecule confined in a nanochannel in the Odijk regime as a function of the fluorescence wavelength, channel size, and the properties of the DNA (persistence length and effective width). The theoretical predictions agree well with equivalent data produced by Monte Carlo simulations of a touching wormlike bead model of DNA in a high ionic strength buffer. Although the theory is only strictly valid in the limit where the effective width of the nanochannel is small compared with the persistence length of the DNA, simulations indicate that the theoretical predictions are reasonably accurate for channel widths up to two-thirds of the persistence length. Our results quantify the conjecture that DNA barcoding has kilobase pair resolution-provided the nanochannel lies in the Odijk regime.

Intermediate filaments in small configuration spaces

Intermediate filaments play a key role in cell mechanics. Apart from their great importance from a biomedical point of view, they also act as a very suitable micrometer-sized model system for semiflexible polymers. We perform a statistical analysis of the thermal fluctuations of individual filaments confined in microchannels. The small channel width and the resulting deflections at the walls give rise to a reduction of the configuration space by about 2 orders of magnitude. This circumstance enables us to precisely measure the intrinsic persistence length of vimentin intermediate filaments and to show that they behave as ideal wormlike chains; we observe that small fluctuations in perpendicular planes decouple. Furthermore, the inclusion of results for confined actin filaments demonstrates that the Odijk confinement regime is valid over at least 1 order of magnitude in persistence length.