Revisiting blob theory for DNA diffusivity in slitlike confinement

Identifying knot types of polymer conformations by machine learning

We investigate the use of artificial neural networks (NNs) as an alternative tool to current analytical methods for recognizing knots in a given polymer conformation. The motivation is twofold. First, it is of interest to examine whether NNs are effective at learning the global and sequential properties that uniquely define a knot. Second, knot classification is an important and unsolved problem in mathematical and physical sciences, and NNs may provide insights into this problem. Motivated by these points, we generate millions of polymer conformations for five knot types: 0, 3_{1}, 4_{1}, 5_{1}, and 5_{2}, and we design various NN models for classification. Our best model achieves a five-class classification accuracy of above 99% on a polymer of 100 monomers. We find that the sequential modeling ability of recurrent NNs is crucial for this result, as it outperforms feed-forward NNs and successfully generalizes to differently sized conformations as well. We present our methods and suggest that deep learning may be used in specific applications of knot detection where some error is permissible. Hopefully, with further development, NNs can offer an alternative computational method for knot identification and facilitate knot research in mathematical and physical sciences.

Self-Avoiding Wormlike Chain Confined in a Cylindrical Tube: Scaling Behavior

Within a confining tube section, the multithreads of a strongly confined, backfolding polymer exert the excluded-volume repulsions on each other and produce physical properties that are very different from those of a confined ideal chain. The conformational properties of a such confined wormlike chain are of fundamental interest and are also practically useful in understanding the DNA confinement problems. Here, the excluded-volume effects are added to the standard wormlike-chain model by a self-consistent field theory. The numerical solutions are examined in light of their scaling properties.

Evaluation of Blob Theory for the Diffusion of DNA in Nanochannels

We have measured the diffusivity of λ-DNA molecules in approximately square nanochannels with effective sizes ranging from 117 nm to 260 nm at moderate ionic strength. The experimental results do not agree with the non-draining scaling predicted by blob theory. Rather, the data are consistent with the predictions of previous simulations of the Kirkwood diffusivity of a discrete wormlike chain model, without the need for any fitting parameters.

From bead to rod: Comparison of theories by measuring translational drag coefficients of micron-sized magnetic bead-chains in Stokes flow

Frictional drag force on an object in Stokes flow follows a linear relationship with the velocity of translation and a translational drag coefficient. This drag coefficient is related to the size, shape, and orientation of the object. For rod-like objects, analytical solutions of the drag coefficients have been proposed based on three rough approximations of the rod geometry, namely the bead model, ellipsoid model, and cylinder model. These theories all agree that translational drag coefficients of rod-like objects are functions of the rod length and aspect ratio, but differ among one another on the correction factor terms in the equations. By tracking the displacement of the particles through stationary fluids of calibrated viscosity in magnetic tweezers setup, we experimentally measured the drag coefficients of micron-sized beads and their bead-chain formations with chain length of 2 to 27. We verified our methodology with analytical solutions of dimers of two touching beads, and compared our measured drag coefficient values of rod-like objects with theoretical calculations. Our comparison reveals several analytical solutions that used more appropriate approximation and derived formulae that agree with our measurement better.

Nanoconfinement greatly speeds up the nucleation and the annealing in single-DNA collapse

Manipulating and measuring single-molecule dynamics and reactions in nanofluidics is a rapidly growing field with broad applications in developing new biotechnologies, understanding nanoconfinement effects in vivo, and exploring new phenomena in confinement. In this work, we investigate the kinetics of DNA collapse in nanoslits using single T4-DNA (165.6 kbp) and λ-DNA (48.5 kbp), with particular focus on the measurement of the nucleation and annealing times. Fixing the ethanol concentration at 35% and varying the slit height from 2000 to 31 nm, the nucleation time dramatically decreases from more than 1 hour to a few minutes or less. The increased collapsed rate results from the larger free energy experienced by coiled DNA in confinement relative to compacted DNA. Our results also shed light on other conformational transitions in confinement, such as protein folding.

Conformational Properties of a Back-Folding Wormlike Chain Confined in a Cylindrical Tube

When a semiflexible chain is confined in a narrow cylindrical tube, the formation of a polymer hairpin is a geometrical conformation that accompanies an exponentially large local free energy and, hence, is a relatively rare event. Numerical solutions of the hairpin distribution functions for persistence-length-to-tube-radius ratios over a wide range are obtained in high precision, by using the Green's function approach for the wormlike-chain model. The crossover region between the narrow and moderately narrow tubes is critically investigated in terms of the hairpin free energy, global persistence length, mean hairpin-tip distance from the tube axis, and hairpin-plane orientational properties. Accurate representations of the solutions by simple interpolation formulae are suggested.

Simulations of knotting of DNA during genome mapping

Genome mapping involves the confinement of long DNA molecules, in excess of 150 kilobase pairs, in nanochannels near the circa 50 nm persistence length of DNA. The fidelity of the map relies on the assumption that the DNA is linearized by channel confinement, which assumes the absence of knots. We have computed the probability of forming different knot types and the size of these knots for long chains (approximately 164 kilobase pairs) via pruned-enriched Rosenbluth method simulations of a discrete wormlike chain model of DNA in channel sizes ranging from 35 nm to 60 nm. Compared to prior simulations of short DNA in similar confinement, these long molecules exhibit both complex knots, with up to seven crossings, and multiple knots per chain. The knotting probability is a very strong function of channel size, ranging from 0.3% to 60%, and rationalized in the context of Odijk's theory for confined semiflexible chains. Overall, the knotting probability and knot size obtained from these equilibrium measurements are not consistent with experimental measurements of the properties of anomalously bright regions along the DNA backbone during genome mapping experiments. This result suggests that these events in experiments are either knots formed during the processing of the DNA prior to injection into the nanochannel or regions of locally high DNA concentration without a topological constraint. If so, knots during genome mapping are not an intrinsic problem for genome mapping technology.

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.

Sorting ring polymers by knot type with modulated nanochannels

In this theoretical study we discuss a novel method for sorting ring polymers according to their topological, knotted state. The proposed approach harnesses the rich dynamical behaviour of polymers confined inside spatially-modulated nanochannels. The longitudinal mobility of the rings is shown to have two key properties that are ideally suited for knot sorting. First, at fixed topology, the mobility has an intriguing oscillatory dependence on chain length. Second, the mobility ranking of different knot types is inverted upon increasing the chain length. We show that this complex interplay of channel geometry, chain length and topology can be rationalised within a simple theoretical framework based on Fick-Jacobs's diffusive theory. The results and the interpretative scheme ought to be useful for designing microfluidic devices with optimal topological sorting capabilities.

Dynamics of DNA Squeezed Inside a Nanochannel via a Sliding Gasket

We use Brownian dynamics (BD) simulation of a coarse-grained (CG) bead-spring model of DNA to study the nonequilibrim dynamics of a single DNA molecule confined inside a rectangular nanochannel being squeezed with a sliding gasket piston or "nanodozer". From our simulations we extract the nonequilibrim density profile c ( x , t ) of the squeezed molecule along the channel axis (x-coordinate) and then analyze the non-equilibrium profile using a recently introduced phenomenological Nonlinear Partial Differential Equation (NPDE) model. Since the NPDE approach also fits the experimental results well and is numerically efficient to implement, the combined BD + NPDE methods can be a powerful approach to analyze details of the confined molecular dynamics. In particular, the overall excellent agreement between the two complementary sets of data provides a strategy for carrying out large scale simulation on semi-flexible biopolymers in confinement at biologically relevant length scales.

Polymer globule with fractal properties caused by intramolecular nanostructuring and spatial constrains

By means of computer simulation, we studied macromolecules composed of N dumbbell amphiphilic monomer units with attractive pendant groups. In poor solvents, these macromolecules form spherical globules that are dense in the case of short chains (the gyration radius RG∼N(1/3)), or hollow inside and obey the RG∼N(1/2) law when the macromolecules are sufficiently long. Due to the specific intramolecular nanostructuring, the vesicle-like globules of long amphiphilic macromolecules posses some properties of fractal globules, by which they (i) could demonstrate the same scaling statistics for the entire macromolecule and for short subchains with m monomer units and (ii) possess a specific territorial structure. Within a narrow slit, the globule loses its inner cavity, takes a disk-like shape and scales as N(1/2) for much shorter macromolecules. However, the field of end-to-end distance r(m) ∼m(1/2) dependence for subchains becomes visibly smaller. The results obtained were compared with the homopolymer case.

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.

Coil-globule transition of a single semiflexible chain in slitlike confinement

Single polymer chains undergo a phase transition from coiled conformations to globular conformations as the effective attraction between monomers becomes strong enough. In this work, we investigated the coil-globule transition of a semiflexible chain confined between two parallel plates, i.e. a slit, using the lattice model and Pruned-enriched Rosenbluth method (PERM) algorithm. We find that as the slit height decreases, the critical attraction for the coil-globule transition changes non-monotonically due to the competition of the confinement free energies of the coiled and globular states. In wide (narrow) slits, the coiled state experiences more (less) confinement free energy, and hence the transition becomes easier (more difficult). In addition, we find that the transition becomes less sharp with the decreasing slit height. Here, the sharpness refers to the sensitivity of thermodynamic quantities when varying the attraction around the critical value. The relevant experiments can be performed for DNA condensation in microfluidic devices.

Experimental evidence of weak excluded volume effects for nanochannel confined DNA

We present experimental demonstration that weak excluded volume effects arise in DNA nanochannel confinement. In particular, by performing measurements of the variance in chain extension as a function of nanochannel dimension for effective channel size ranging from 305 nm to 453 nm, we show that the scaling of the variance in extension with channel size rejects the de Gennes scaling δ2X ~ D1/3 in favor of δ2X ~ D0 using uncertainty at the 95% confidence level. We also show how simulations and confinement spectroscopy can be combined to reduce molecular weight dispersity effects arising from shearing, photocleavage, and nonuniform staining of DNA.

Force-Extension for DNA in a Nanoslit: Mapping between the 3D and 2D Limits

The force-extension relation for a semiflexible polymer confined in a nanoslit is investigated. Both the effective correlation length and force-extension relation change as the chain goes from 3D (large slit heights) to 2D (tight confinement). At low forces, correlations along the polymer give an effective dimensionality. The strong force limit can be interpolated with the weak force limit for two regimes: when confinement dominates over extensile force and vice versa. These interpolations give good agreement with simulations for all slit heights and forces. We thus generalize the Marko-Siggia force-extension relation for DNA and other semiflexible biopolymers in nanoconfinement.

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.

Electrophoretic Mobility of Polyelectrolytes within a Confining Well

We present a numerical study of polyelectrolytes electrophoresing in free solution while squeezed by an axisymmetric confinement force transverse to their net displacement. Hybrid multiparticle collision dynamics and molecular dynamics simulations with mean-field finite Debye layers show that even though the polyelectrolyte chains remain "free-draining" their electrophoretic mobility increases with confinement in nanoconfining potential wells. The primary mechanism leading to the increase in mobility above the free-solution value, despite long-range hydrodynamic screening by counterion layers, is the orientation of polymer segments within Debye layers. The observed length dependence of the electrophoretic mobility arises due to secondary effects of counterion condensation related to confinement compactification.

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.

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.

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.

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.

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.

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.

Precision platform for convex lens-induced confinement microscopy

We present the conception, fabrication, and demonstration of a versatile, computer-controlled microscopy device which transforms a standard inverted fluorescence microscope into a precision single-molecule imaging station. The device uses the principle of convex lens-induced confinement [S. R. Leslie, A. P. Fields, and A. E. Cohen, Anal. Chem. 82, 6224 (2010)], which employs a tunable imaging chamber to enhance background rejection and extend diffusion-limited observation periods. Using nanopositioning stages, this device achieves repeatable and dynamic control over the geometry of the sample chamber on scales as small as the size of individual molecules, enabling regulation of their configurations and dynamics. Using microfluidics, this device enables serial insertion as well as sample recovery, facilitating temporally controlled, high-throughput measurements of multiple reagents. We report on the simulation and experimental characterization of this tunable chamber geometry, and its influence upon the diffusion and conformations of DNA molecules over extended observation periods. This new microscopy platform has the potential to capture, probe, and influence the configurations of single molecules, with dramatically improved imaging conditions in comparison to existing technologies. These capabilities are of immediate interest to a wide range of research and industry sectors in biotechnology, biophysics, materials, and chemistry.