Diffusion of macromolecules in narrow capillaries

The effects of slit-like confinement on flow-induced polymer deformation

This paper is broadly concerned with the dynamics of a polymer confined to a rectangular slit of width D and deformed by a planar elongational flow of strength γ̇. It is interested, more specifically, in the nature of the coil-stretch transition that such polymers undergo when the flow strength γ̇ is varied, and in the degree to which this transition is affected by the presence of restrictive boundaries. These issues are explored within the framework of a finitely extensible Rouse model that includes pre-averaged surface-mediated hydrodynamic interactions. Calculations of the chain's steady-state fractional extension x using this model suggest that different modes of relaxation (which are characterized by an integer p) exert different levels of control on the coil-stretch transition. In particular, the location of the transition (as identified from the graph of x versus the Weissenberg number Wi, a dimensionless parameter defined by the product of γ̇ and the time constant τ_p of a relaxation mode p) is found to vary with the choice of τ_p. In particular, when τ₁ is used in the definition of Wi, the x vs. Wi data for different D lie on a single curve, but when τ₃ is used instead (with τ₃ > τ₁) the corresponding data lie on distinct curves. These findings are in close qualitative agreement with a number of experimental results on confinement effects on DNA stretching in electric fields. Similar D-dependent trends are seen in our calculated force vs. Wi data, but force vs. x data are essentially D-independent and lie on a single curve.

Diffusive dynamics of polymer chains in an array of nanoposts

The diffusion of the head, side, and middle segments in confined polymer chains displays different dynamics in different directions.

Kirkwood diffusivity of long semiflexible chains in nanochannel confinement

We compute the axial diffusivity of asymptotically long semiflexible polymers confined in square channels. Our calculations employ the Kirkwood approximation of the mobility tensor by combining computational fluid dynamics (CFD) calculations of the hydrodynamic tensor in channel confinement with pruned-enriched Rosenbluth method (PERM) simulations of a discrete wormlike chain model. Three key results emerge from our study. First, for the classic de Gennes regime, we confirm that Brochard and de Gennes' blob theory correctly predicts the scaling of the axial diffusivity, contrary to the conclusions of previous analyses. Second, for the extended de Gennes regime, we show that a modified blob theory, which has been used to incorporate the effect of local stiffness on DNA diffusion in nanoslits, explains the deviation from the prediction of classic blob theory for diffusion in nanochannels. Third, we provide a calculation similar to the modified blob theory to explain the relative insensitivity of the diffusivity to channel size for channels between the extended de Gennes regime and the Odijk regime, which is the most relevant regime for experiments and technological applications of DNA confinement in nanochannels. Our results are not only relevant to the dynamics of confined semiflexible polymers such as DNA, but also reveal interesting analogies between confinement in channels and slits.

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.

The diffusion and relaxation of Gaussian chains in narrow rectangular slits

The confinement of a polymer to volumes whose characteristic linear dimensions are comparable to or smaller than its bulk radius of gyration R(G,bulk) can produce significant changes in its static and dynamic properties, with important implications for the understanding of single-molecule processes in biology and chemistry. In this paper, we present calculations of the effects of a narrow rectangular slit of thickness d on the scaling behavior of the diffusivity D and relaxation time τr of a Gaussian chain of polymerization index N and persistence length l0. The calculations are based on the Rouse-Zimm model of chain dynamics, with the pre-averaged hydrodynamic interaction being obtained from the solutions to Stokes equations for an incompressible fluid in a parallel plate geometry in the limit of small d. They go beyond de Gennes' purely phenomenological analysis of the problem based on blobs, which has so far been the only analytical route to the determination of chain scaling behavior for this particular geometry. The present model predicts that D ∼ dN(-1)ln (N∕d(2)) and τr ∼ N(2)d(-1)[ln (N∕d(2))](-1) in the regime of moderate confinement, where l0 ≪ d < R(G,bulk). The corresponding results for the blob model have exactly the same power law behavior, but contain no logarithmic corrections; the difference suggests that segments within a blob may actually be partially draining and not non-draining as generally assumed.

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.

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.

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.

Simulation of DNA Extension in Nanochannels

We have used a realistic model for double stranded DNA and Monte Carlo simulations to compute the extension (mean span) of a DNA molecule confined in a nanochannel over the full range of confinement in a high ionic strength buffer. The simulation data for square nanochannels resolve the apparent contradiction between prior simulation studies and the predictions from Flory theory, demonstrating the existence of two transition regimes between weak confinement (the de Gennes regime) and strong confinement (the Odijk regime). The simulation data for rectangular nanochannels support the use of the geometric mean for mapping data obtained in rectangular channels onto models developed for cylinders. The comparison of our results with experimental data illuminates the challenges in applying models for confined, neutral polymers to polyelectrolytes. Using a Flory-type approach, we also provide an improved scaling result for the relaxation time in the transition regime close to that found in experiments.

Self-avoiding polymer trapped inside a cylindrical pore: Flory free energy and unexpected dynamics

We study the elastic and dynamic behavior of a self-avoiding chain confined inside a cylindrical pore using a Flory-type approach and molecular-dynamics simulations. In the Hookean regime, we find that the effective spring constant of the chain is given by keff approximately N(-1)D(-gamma), where N is the number of monomers and D the diameter of the pore. While the Flory approach reproduces the earlier scaling result gamma=1/3, our simulations confirm a more recent numerical result gamma approximately 0.9 for the computationally accessible regimes. In the absence of hydrodynamic interactions, the relaxation dynamics of a stretched-and-released chain in this regime is characterized by a global relaxation time tauR approximately N2Dgamma with the same exponent gamma for keff. We also discuss how chain relaxation under confinement is influenced by hydrodynamic interactions. In the presence (or absence) of the hydrodynamic interaction, the finite-size effect observed in keff is shown to persist in chain relaxation, resulting in tauR markedly different from previous results.

Diffusion of PEG confined between lamellae of negatively magnetically aligned bicelles: pulsed field gradient 1H NMR measurements

The diffusion of various molecular weight poly(ethyleneglycol)s (PEG) confined between the lamellae of magnetically aligned bicelles has been measured using stimulated echo (STE) pulsed field gradient (PFG) 1H nuclear magnetic resonance (NMR) spectroscopy. Bicelles were formulated to contain dimyristoylphosphatidylcholine (DMPC), dimyristoylphosphatidylglycerol (DMPG), and dihexanoylphosphatidylcholine (DHPC) in the proportion DMPG/DMPC = 0.05 and q = (DMPC + DMPG)/DHPC = 4.5. PEG diffusion within the interlamellar spaces between such bicelles was found to be unrestricted over diffusion distances of tens of microns. Two confinement regimes could be differentiated according to the dependence of the reduced PEG diffusivity D/D0, where D0 is the unconfined PEG diffusion coefficient, on the relative confinement Rh/H, where Rh is the unperturbed hydration radius of the particular PEG and H approximately 60 A is the separation between apposing lamellae of the magnetically aligned bicelles. In the regime Rh/H < 0.4, the reduced PEG diffusivity was altered only in proportion to the viscosity increase associated with the bicelle dispersion relative to bulk solution. In the regime Rh/H > 0.4, the reduced PEG diffusivity scaled as (Rh/H)-2/3, in agreement with scaling theories for confined polymers.

Lattice-Boltzmann simulations of the dynamics of polymer solutions in periodic and confined geometries

A numerical method to simulate the dynamics of polymer solutions in confined geometries has been implemented and tested. The method combines a fluctuating lattice-Boltzmann model of the solvent [Ladd, Phys. Rev. Lett. 70, 1339 (1993)] with a point-particle model of the polymer chains. A friction term couples the monomers to the fluid [Ahlrichs and Dunweg, J. Chem. Phys. 111, 8225 (1999)], providing both the hydrodynamic interactions between the monomers and the correlated random forces. The coupled equations for particles and fluid are solved on an inertial time scale, which proves to be surprisingly simple and efficient, avoiding the costly linear algebra associated with Brownian dynamics. Complex confined geometries can be represented by a straightforward mapping of the boundary surfaces onto a regular three-dimensional grid. The hydrodynamic interactions between monomers are shown to compare well with solutions of the Stokes equations down to distances of the order of the grid spacing. Numerical results are presented for the radius of gyration, end-to-end distance, and diffusion coefficient of an isolated polymer chain, ranging from 16 to 1024 monomers in length. The simulations are in excellent agreement with renormalization group calculations for an excluded volume chain. We show that hydrodynamic interactions in large polymers can be systematically coarse-grained to substantially reduce the computational cost of the simulation. Finally, we examine the effects of confinement and flow on the polymer distribution and diffusion constant in a narrow channel. Our results support the qualitative conclusions of recent Brownian dynamics simulations of confined polymers [Jendrejack et al., J. Chem. Phys. 119, 1165 (2003) and Jendrejack et al., J. Chem. Phys. 120, 2513 (2004)].