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

Electrical Transport and Dynamics of Confined DNA through Highly Conductive 2D Graphene Nanochannels

Confining DNA in nanochannels is an important approach to studying its structure and transportation dynamics. Graphene nanochannels are particularly attractive for studying DNA confinement due to their atomic flatness, precise height control, and excellent mechanical strength. Here, using femtosecond laser etching and wetting transfer, we fabricate graphene nanochannels down to less than 4.3 nm in height, with the length-to-height ratios up to 103. These channels exhibit high stability, low noise, and self-cleaning ability during the long-term ionic current recording. We report a clear linear relationship between DNA length and the residence time in the channel and further utilize this relationship to differentiate DNA fragments based on their lengths, ranging widely from 200 bps to 48.5 kbps. The graphene nanochannel presented here provides a potential platform for label-free analyses and reveals fundamental insights into the conformational dynamics of DNA and proteins in confined space.

Compression and Stretching of Single DNA Molecules under Channel Confinement

We study the compression and extension response of single dsDNA (double-stranded DNA) molecules confined in cylindrical channels by means of Monte Carlo simulations. The elastic response of micrometer-sized DNA to the external force acting through the chain ends or through the piston is markedly affected by the size of the channel. The interpretation of the force (f)-displacement (R) functions under quasi-one-dimensional confinement is facilitated by resolving the overall change of displacement ΔR into the confinement contribution ΔR_D and the force contribution ΔR_f. The external stretching of confined DNA results in a characteristic pattern of f-R functions involving their shift to the larger extensions due to the channel-induced pre-stretching ΔR_D. A smooth end-chain compression into loop-like conformations observed in moderately confined DNA can be accounted for by the relationship valid for a Gaussian chain in bulk. In narrow channels, the considerably pre-stretched DNA molecules abruptly buckle on compression by the backfolding into hairpins. On the contrary, the piston compression of DNA is characterized by a gradual reduction of the chain span S and by smooth f-S functions in the whole spatial range from the 3d near to 1d limits. The observed discrepancy between the shape of the f-R and f-S functions from two compression methods can be important for designing nanopiston experiments of compaction and knotting of single DNA in nanochannels.

Stretching Wormlike Chains in Narrow Tubes of Arbitrary Cross-Sections

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

Force-displacement relations at compression of dsDNA macromolecules

The elasticity of dsDNA molecules is investigated by Monte Carlo simulations based on a coarse-grained model of DNA. The force-displacement (f-r) curves are computed under the constraints of the constant force (Gibbs) or the constant length (Helmholtz) ensemble. Particular attention was paid to the compressional (negative) and weak tensile forces. It was confirmed that simulations using the vector Gibbs ensemble fail to represent the compression behavior of polymers. Simulations using the scalar Gibbs protocol resulted in a qualitatively correct compressional response of DNA provided that the quadratic averages of displacements were employed. Furthermore, a well-known shortcoming of the popular Marko-Siggia relation for DNA elasticity at weak tensile forces is elucidated. Conversely, the function f-r from the simulation at the constant length constraint, as well as the new closed-form expressions, provides a realistic depiction of the DNA elasticity over the wide range of negative and positive forces. Merely a qualitative resemblance of the compression functions f-r predicted by the employed approaches supports the notion that the elastic response of DNA molecules may be greatly affected by the specifics of the experimental setups and the kind of averaging of the measured variable.

Correlation anisotropy and stiffness of DNA molecules confined in nanochannels

The anisotropy of orientational correlations in DNA molecules confined in cylindrical channels is explored by Monte Carlo simulations using a coarse-grained model of double-stranded (ds) DNA. We find that the correlation function ⟨C(s)⟩_⊥ in the transverse (confined) dimension exhibits a region of negative values in the whole range of channel sizes. Such a clear-cut sign of the opposite orientation of chain segments represents a microscopic validation of the Odijk deflection mechanism in narrow channels. At moderate-to-weak confinement, the negative ⟨C(s)⟩_⊥ correlations imply a preference of DNA segments for transverse looping. The inclination for looping can explain a reduction of stiffness as well as the enhanced knotting of confined DNA relative to that detected earlier in bulk at some channel sizes. Furthermore, it is shown that the orientational persistence length P_or fails to convey the apparent stiffness of DNA molecules in channels. Instead, correlation lengths P_∥ and P_⊥ in the axial and transverse directions, respectively, encompass the channel-induced modifications of DNA stiffness.

Stretching and compression of DNA by external forces under nanochannel confinement

Mechanical deformation of dsDNA molecules inside square nanochannels is investigated using simulations based on a coarse-grained model of DNA. The combined action of confinement and weak external forces is explored in a variety of confinement regimes, including the transition zone relevant to nanofluidic experiments. The computed free energy and force profiles are markedly affected by the channel size. Effective elastic softening of confined DNA molecules relative to the bulk DNA is observed in the channels of intermediate widths. The extension of DNA from its bulk equilibrium length in nanofluidic devices is resolved into contributions from the passive extension due to confinement and from the active stretching induced by force. Potential implications of the very different energy costs computed for the two extension modes (extension by confinement takes much more free energy than stretching by force) for behavior of DNA in nanofluidic chips are indicated.

Sorting polymers by size via an array of viscous posts

DNA fragments can be sorted according to size by forcing them through an array of nanoposts. Whereas previous studies have explored solid nanoposts, this work examines nanoposts constructed out of viscous inclusions. Langevin dynamics simulations are used to study the dynamics of polymers driven through arrays of these viscous nanoposts for a range of post viscosities. The results are compared to the solid post case. Increasing post viscosity causes a decrease in the mobility of polymers traversing the array. In the limit of high post viscosity, the mobility becomes lower than in the solid post arrays, rather than converging to it. Analysis of the distributions of event times also shows that the viscous case is fundamentally different from the solid post case. The decrease in mobility in the viscous case arises from slowing down the polymer as it interacts with or even moves through the nanoposts, whereas the solid post case exhibits wrapping and unwrapping dynamics, yielding escape-like statistics. This work suggests that it may be possible to use viscous inclusions within nanofluidic and microfluidic devices to sort biomolecules with high resolution.

Stretching of Tethered DNA in Nanoslits

External forces and confinement are two fundamental and complementary approaches for biopolymer stretching. By employing micro- and nanofluidics, we study the force-extension dynamics by simultaneously applying external forces and confinement to single-DNA molecules. In particular, we apply external electric fields to stretch single DNA molecules that are attached to microspheres anchored at a nanoslit entrance. Using this method, we measure the force-extension relation of tethered DNA and describe this relation with modified wormlike chain models. This allowed experimental validations of several theoretical predictions, including the increase in the global persistence length of confined DNA with increasing degree of confinement and the "confined Pincus" regime in slit confinement.

How to determine local stretching and tension in a flow-stretched DNA molecule

We determine the nonuniform stretching of and tension in a mega base pairs-long fragment of deoxyribonucleic acid (DNA) that is flow stretched in a nanofluidic chip. We use no markers, do not know the contour length of the DNA, and do not have the full DNA molecule inside our field of view. Instead, we analyze the transverse thermal motion of the DNA. Tension at the center of the DNA adds up to 16 pN, giving almost fully stretched DNA. This method was devised for optical mapping of DNA, specifically, DNA denaturation patterns. It may be useful also for other studies, e.g., DNA-protein interactions, specifically, their tension dependence. Generally, wherever long strands of DNA-e.g., native DNA extracted from human cells or bacteria-must be stretched with ease for inspection, this method applies.

Effect of excluded volume on the force-extension of wormlike chains in slit confinement

We use pruned-enriched Rosenbluth method simulations to develop a quantitative phase diagram for the stretching of a real wormlike chain confined in a slit. Our simulations confirm the existence of a "confined Pincus" regime in slit confinement, analogous to the Pincus regime in free solution, where excluded volume effects are sensible. The lower bound for the confined Pincus regime in the force-molecular weight plane, as well as the scaling of the extension with force and slit size, agree with an existing scaling theory for this regime. The upper bound of the confined Pincus regime depends on the strength of the confinement. For strong confinement, the confined Pincus regime ends when the contour length in the Pincus blob is too short to have intrablob excluded volume. As a result, the chain statistics become ideal and the confined Pincus regime at low forces is connected directly to ideal chain stretching at large forces. In contrast, for weak confinement, the confined Pincus regime ends when the Pincus blobs no longer fit inside the slit, even though there is sufficient contour length to have excluded volume inside the Pincus blob. As a result, weak confinement leads to a free-solution Pincus regime intervening between the confined Pincus regime for weak forces and ideal chain stretching at strong forces. Our results highlight shortcomings in existing models for the stretching of wormlike chains in slits.