Conformation and dynamics of single DNA molecules in parallel-plate slit microchannels

Modeling properties of chromosome territories using polymer filaments in diverse confinement geometries

Interphase chromosomes reside within distinct nuclear regions known as chromosome territories (CTs). Recent observations from Hi-C analyses, a method mapping chromosomal interactions, have revealed varied decay in contact probabilities among different chromosomes. Our study explores the relationship between this contact decay and the particular shapes of the chromosome territories they occupy. For this, we employed molecular dynamics (MD) simulations to examine how confined polymers, resembling chromosomes, behave within different confinement geometries similar to chromosome territory boundaries. Our simulations unveil so far unreported relationships between contact probabilities and end-to-end distances varying based on different confinement geometries. These findings highlight the crucial impact of chromosome territories on shaping the larger-scale properties of 3D genome organization. They emphasize the intrinsic connection between the shapes of these territories and the contact behaviors exhibited by chromosomes. Understanding these correlations is key to accurately interpret Hi-C and microscopy data, and offers vital insights into the foundational principles governing genomic organization.

Crowding and confinement act in concert to slow DNA diffusion within cell-sized droplets

Dynamics of biological macromolecules, such as DNA, in crowded and confined environments are critical to understanding cellular processes such as transcription, infection, and replication. However, the combined effects of cellular confinement and crowding on macromolecular dynamics remain poorly understood. Here, we use differential dynamic microscopy to investigate the diffusion of large DNA molecules confined in cell-sized droplets and crowded by dextran polymers. We show that confined and crowded DNA molecules exhibit universal anomalous subdiffusion with scaling that is insensitive to the degree of confinement and crowding. However, effective DNA diffusion coefficients Deff decrease up to 2 orders of magnitude as droplet size decreases—an effect that is enhanced by increased crowding. We mathematically model the coupling of crowding and confinement by combining polymer scaling theories with confinement-induced depletion effects. The generality and tunability of our system and models render them applicable to elucidating wide-ranging crowded and confined systems.

•

DNA diffusion measured in cell-sized droplets with differential dynamic microscopy

•

Combination of crowding and confinement leads to subdiffusion and slowing

•

Diffusion coefficients of DNA decrease strongly with decreasing droplet size

•

Polymer scaling theories and depletion effects predict observed dynamics

Intrinsic structure and dynamics of monolayer ring polymer melts

We present the general structural and dynamical characteristics of flexible ring polymers in narrowly confined two-dimensional (2D) melt systems using atomistic molecular dynamics simulations. The results are further analyzed via direct comparison with the 2D linear analogue as well as the three-dimensional (3D) ring and linear melt systems. It is observed that dimensional restriction in 2D confined systems results in an increase in the intrinsic chain stiffness of the ring polymer. Fundamentally, this arises from an entropic penalty on polymer chains along with a reduction in the available chain configuration states in phase space and spatial choices for individual segmental walks. This feature in combination with the intermolecular interactions between neighboring ring chains leads to an overall extended interpenetrated chain configuration for the 2D ring melt. In contrast to the generally large differences in structural and dynamical properties between ring and linear polymers in 3D melt systems, relatively similar local-to-global chain structures and dynamics are observed for the 2D ring and linear melts. This is attributed to the general structural similarity (i.e., extended double-stranded chain conformations), the less effective role of the chain ends, and the absence of complex topological constraints between chains (i.e., interchain entanglement and mutual ring threading) in the 2D confined systems compared with the corresponding 3D bulk systems.

Statistical Behaviors of Semiflexible Polymer Chains Stretched in Rectangular Tubes

We investigated the statistical behaviors of semiflexible polymer chains, which were simultaneously subjected to force stretching and rectangular tube confinement. Based on the wormlike chain model and Odijk deflection theory, we derived a new deflection length, by using which new compact formulas were obtained for the confinement free energy and force⁻confinement⁻extension relations. These newly derived formulas were justified by numerical solutions of the eigenvalue problem associated with the Fokker⁻Planck governing equation and extensive Brownian dynamics simulations based on the so-called generalized bead-rod (GBR) model. We found that, compared to classical deflection theory, these new formulas were valid for a much more extended range of the confinement size/persistence length ratio and had no adjustable fitting parameters for sufficiently long semiflexible chains in the whole deflection regime.

Structure and proton conduction in sulfonated poly(ether ether ketone) semi-permeable membranes: a multi-scale computational approach

The design of polymeric membranes for proton or ionic exchange highly depends on the fundamental understanding of the physical and molecular mechanisms that control the formation of the conduction channels. There is an inherent relation between the dynamical structure of the polymeric membrane and the electrostatic forces that drive membrane segregation and proton transport. Here, we used a multi-scale computational approach to analyze the morphology of sulfonated poly(ether ether ketone) membranes at the mesoscale. A self-consistent description of the electrostatic phenomenon was adopted, where discrete polymer chains and a continuum proton field were embedded in a continuum fluid. Brownian dynamics was used for the evolution of the suspended polymer molecules, while a convection-diffusion transport equation, including the Nernst-Planck diffusion mechanism, accounted for the dynamics of the proton concentration field. We varied the polymer concentration, the degree of sulfonation and the level of confinement to find relationships between membrane structure and proton conduction. Our results indicate that the reduced mobility of polymer chains, at concentrations above overlap, and a moderate degree of sulfonation - i.e., 30% - are essential elements for membrane segregation and proton domain connectivity. These conditions also ensure that the membrane structure is not affected by size or by potential gradients. Importantly, our analysis shows that membrane conductivity and current are linearly dependent on polymer concentration and quadratically dependent on the degree of sulfonation. We found that the optimal polymeric membrane design requires a polymer concentration above overlap and a degree of sulfonation around 50%. These conditions promote a dynamical membrane morphology with a constant density of proton channels. Our results and measurements agree with previous experimental works, thereby validating our model and observations.

Relaxation of surface-tethered polymers under moderate confinement

We study the relaxation of surface-tethered polymers in microchannels under moderate confinement (i.e. h ∼ Rg, where h is the channel height and Rg is the radius of gyration of the polymer) by experiments with fluorescence-marked DNA molecules and coupled lattice-Boltzmann/molecular dynamics simulations. The determined scaling exponent suggests that the relaxation is dominated by Zimm-dynamics with significant intra-chain hydrodynamic interactions. The relaxation of the DNA molecules is slower in shallower channels, indicating a pronounced effect of confinement on the longest relaxation time. An experimental correlation is obtained for the longest relaxation time as a function of the molecular contour length and the channel height. Good agreement between the experimental and the simulation results is found.

DNA Confined in a Nanodroplet: A Molecular Dynamics Study

As a major genetic material, the configuration and the mechanical properties of a double-stranded DNA (dsDNA) molecule in confinement are crucial for the application of nanotechnology and biological engineering. In the present paper, molecular dynamics simulation is utilized to study the configuration of dsDNA in a nanodroplet on a graphene substrate. The results show that the semiflexible dsDNA molecule changes its configuration with radius of gyration ( R_g) of a few nanometers because of the confined space, that is, the R_g of the dsDNA molecule decreases with the reduction of the nanodroplet size. In comparison, the dsDNA in the bulk usually has a persistent length of tens of nanometers. Especially, if the nanodroplet is small enough, the dsDNA molecule might form a loop structure inside. The dsDNA molecule affects the wetting properties of the graphene substrate. It is found that the graphene becomes more hydrophilic in smaller systems containing the dsDNA molecule, whereas for larger droplets, the changes of the contact angles are not significant with the presence of dsDNA. Moreover, the results indicate that for larger droplets, the line tension of the droplet containing DNA is positive and greater than that without DNA; for smaller droplets, the line tension becomes negative because the dsDNA is compressed and bent in the confinement, and has the potential to expand outwards. The worm-like chain model is used to study the bending energy of a dsDNA molecule in a droplet. The results address that the bending energy of the non-loop-structured dsDNA decreases as the droplet becomes larger, and it is larger than that of loop-structured dsDNA, as the loop structure efficiently prevents the DNA from bending in the vertical direction.

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.

Wall depletion length of a channel-confined polymer

Numerous experiments have taken advantage of DNA as a model system to test theories for a channel-confined polymer. A tacit assumption in analyzing these data is the existence of a well-defined depletion length characterizing DNA-wall interactions such that the experimental system (a polyelectrolyte in a channel with charged walls) can be mapped to the theoretical model (a neutral polymer with hard walls). We test this assumption using pruned-enriched Rosenbluth method (PERM) simulations of a DNA-like semiflexible polymer confined in a tube. The polymer-wall interactions are modeled by augmenting a hard wall interaction with an exponentially decaying, repulsive soft potential. The free energy, mean span, and variance in the mean span obtained in the presence of a soft wall potential are compared to equivalent simulations in the absence of the soft wall potential to determine the depletion length. We find that the mean span and variance about the mean span have the same depletion length for all soft potentials we tested. In contrast, the depletion length for the confinement free energy approaches that for the mean span only when depletion length no longer depends on channel size. The results have implications for the interpretation of DNA confinement experiments under low ionic strengths.

Single-molecule fluorescence resonance energy transfer in molecular biology

Single-molecule fluorescence resonance energy transfer (smFRET) is a powerful technique for studying the conformation dynamics and interactions of individual biomolecules. In this review, we describe the concept and principle of smFRET, illustrate general instrumentation and microscopy settings for experiments, and discuss the methods and algorithms for data analysis. Subsequently, we review applications of smFRET in protein conformational changes, ion channel open-close properties, receptor-ligand interactions, nucleic acid structure regulation, vesicle fusion, and force induced conformational dynamics. Finally, we discuss the main limitations of smFRET in molecular biology.

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.

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.

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.

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.

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.

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.

Transport of DNA in hydrophobic microchannels: a dissipative particle dynamics simulation

In this work, we numerically study a new means of manipulating single DNA chains in microchannels. The method is based on the effect of finite slip at hydrophobic walls on the hydrodynamics and, consequently, on the dynamics of the DNA in microchannels. We use dissipative particle dynamics to study DNA transport as a function of chain length and the Reynolds number in two dimensional parallel plate channels. We show how an asymmetric velocity profile in a channel with hydrophobic and hydrophilic walls can be used to manipulate the location of the DNA molecules. Using this effect, we propose a simple arrangement of hydrophobic and hydrophilic strips which can be exploited to separate long and short DNA chains.

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.

Simulations of DNA stretching by flow field in microchannels with complex geometry

Recently, we have reported the experimental results of DNA stretching by flow field in three microchannels (C. H. Lee and C. C. Hsieh, Biomicrofluidics 7(1), 014109 (2013)) designed specifically for the purpose of preconditioning DNA conformation for easier stretching. The experimental results do not only demonstrate the superiority of the new devices but also provides detailed observation of DNA behavior in complex flow field that was not available before. In this study, we use Brownian dynamics-finite element method (BD-FEM) to simulate DNA behavior in these microchannels, and compare the results against the experiments. Although the hydrodynamic interaction (HI) between DNA segments and between DNA and the device boundaries was not included in the simulations, the simulation results are in fairly good agreement with the experimental data from either the aspect of the single molecule behavior or from the aspect of ensemble averaged properties. The discrepancy between the simulation and the experimental results can be explained by the neglect of HI effect in the simulations. Considering the huge savings on the computational cost from neglecting HI, we conclude that BD-FEM can be used as an efficient and economic designing tool for developing new microfluidic device for DNA manipulation.

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.

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.

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.

Simulation of single DNA molecule stretching and immobilization in a de-wetting two-phase flow over micropillar-patterned surface

We investigate single DNA stretching dynamics in a de-wetting flow over micropillars using Brownian dynamics simulation. The Brownian dynamics simulation is coupled with transient flow field computation through a numerical particle tracking algorithm. The droplet formation on the top of the micropillar during the de-wetting process creates a flow pattern that allows DNA to stretch across the micropillars. It is found that DNA nanowire forms if DNA molecules could extend across the stagnation point inside the connecting water filament before its breakup. It also shows that DNA locates closer to the top wall of the micropillar has higher chance to enter the flow pattern of droplet formation and thus has higher chance to be stretched across the micropillars. Our simulation tool has the potential to become a design tool for DNA manipulation in complex biomicrofluidic devices.

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.

Critical adsorption of a flexible polymer confined between two parallel interacting surfaces

Critical adsorption of a lattice self-avoiding bond fluctuation polymer chain confined between two parallel impenetrable surfaces is studied using the Monte Carlo method. The dependence of the mean contact number <M> on the temperature T and on the chain length N is simulated for a polymer-surface interaction E=-1. A critical adsorption of the polymer is found at T(c)=1.65 for large surface separation distance D>N(ν)b, whereas no critical adsorption is observed for small distance D<N(ν)b, where ν ≈ 0.58 is the Flory exponent and b is the mean bond length. The critical adsorption point T(c)=1.65 is the same as that of a grafted polymer. Normal diffusion is observed for the confined polymer; however, the diffusion rate is dependent on the temperature and surface separation distance.

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.

Stretching DNA by electric field and flow field in microfluidic devices: An experimental validation to the devices designed with computer simulations

We examined the performance of three microfluidic devices for stretching DNA. The first device is a microchannel with a contraction, and the remaining two are the modifications to the first. The modified designs were made with the help of computer simulations [C. C. Hsieh and T. H. Lin, Biomicrofluidics 5(4), 044106 (2011) and C. C. Hsieh, T. H. Lin, and C. D. Huang, Biomicrofluidics 6, 044105 (2012)] and they were optimized for operating with electric field. In our experiments, we first used DC electric field to stretch DNA. However, the experimental results were not even in qualitative agreement with our simulations. More detailed investigation revealed that DNA molecules adopt a globular conformation in high DC field and therefore become more difficult to stretch. Owing to the similarity between flow field and electric field, we turned to use flow field to stretch DNA with the same devices. The evolution patterns of DNA conformation in flow field were found qualitatively the same as our prediction based on electric field. We analyzed the maximum values, the evolution and the distributions of DNA extension at different Deborah number in each device. We found that the shear and the hydrodynamic interaction have significant influence on the performance of the devices.

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.