One-Parameter Scaling Theory for DNA Extension in a Nanochannel

Experiments measuring DNA extension in nanochannels are at odds with even the most basic predictions of current scaling arguments for the conformations of confined semiflexible polymers such as DNA. We show that a theory based on a weakly self-avoiding, one-dimensional “telegraph” process collapses experimental data and simulation results onto a single master curve throughout the experimentally relevant region of parameter space and explains the mechanisms at play.

contour length L, the persistence length ℓ P quantifying the stiffness of the chain, and the effective width w that appears in the Onsager excluded volume [17]. For polyelectrolytes such as DNA, both the persistence length [18][19][20][21] and the effective polymer width [22] depend on electrostatic interactions. Recent experiments are often conducted in high ionicstrength buffers. In this case ℓ P is approximately 50 nm [23] while w is around 5 nm [24], and thus w ≪ ℓ P . This inequality emphasizes the intrinsic difficulty of describing DNA in a wide range of situations. DNA is considerably stiffer than typical synthetic polymers, yet the number of persistence lengths L/ℓ P in genomic DNA samples is large. Any theory for the conformational statistics of channel-confined DNA must account for both the local stiffness of the polymer and excluded-volume interactions. This is a formidable challenge. Matters are further complicated by the fact that most recent genome-mapping experiments are performed in nanochannels that are about 50 nm wide [6][7][8]10]. Essentially all experiments involving DNA (Fig. 1) take place in channel sizes D of the order of ℓ P and do not satisfy the criterion required for the scaling arguments of Ref. [9] to apply.
There is nomicroscopic theory for the extension of confined DNA for , despite numerous attempts [28,38,39]. Scaling arguments [13,14,16] following Refs. [9,11] yield the most useful description. If D ≫ ℓ P , they suggest that the extension X scales as X ~ D −2/3 . The problem is that the theory fails when D ≈ ℓ P , and as a result it proves to be a poor description of most recent DNA experiments in nanochannels. The earliest experiments [25], for instance, reported a much larger exponent X ~ D −0.85 , and subsequent studies [32,33,40] continue to report exponents exceeding the theoretical prediction.
We take a different approach in this Letter. We show that the DNA-confinement problem for w ≪ ℓ P and maps to the simple one-dimensional telegraph process in Fig. 2, describing the correlated walk of a particle moving with velocity v 0 along the channel axis. The velocity changes sign at rate r, creating hairpin configurations in the particle path. The process lasts from t = 0 to t = T. When the particle revisits a position it has previously visited, it incurs a penalty ε. We show that this model collapses experimental and simulation data for the extension throughout the experimentally relevant parameter range onto a universal master curve as a function of a new scaling variable α that measures the combined effects of confinement, local stiffness, and self-avoidance.
We start by considering narrow channels, D ≪ ℓ P , and later extend the arguments to channel widths up to . Since we want to compute the extension of the DNA molecule along the channel axis, it suffices to consider the projection of the three-dimensional DNA configurations x(s) to the channel axis x. Here s is the contour-length coordinate; it corresponds to time t in the telegraph process. We decompose the probability P[x(s)] of observing the projected conformation x(s) as (1) The functional P ideal [x(s)] is the probability of observing the conformation x(s) in an ensemble of ideal confined polymers, disregarding self-avoidance. The functional [x(s)] captures the effect of self-avoidance. It equals the fraction of three-dimensional polymer configurations corresponding to x(s) that contains no segments that collide with any other polymer segment.
When D < ℓ P , the local conformation statistics are determined by Odijk's theory for narrow channels [41], while the global statistics are dominated by a random sequence of direction changes (hairpins) [42]. Neglecting self-avoidance, the separation between neighboring hairpin bends is exponentially distributed [42]. On length scales much larger than the deflection length [41] λ ≡ (ℓ P D 2 ) 1/3 , the central-limit theorem assures that local alignment fluctuations are negligible [11]. These two facts imply that the ideal problem maps onto the one-dimensional telegraph process in Fig. 2. The correlation function of v x (s), the channelaxis component of the tangent vector of the ideal polymer, decays exponentially [36]: The telegraph velocity has similar correlations: Comparing Eqs. (2) and (3) we see that the contour parameter s maps to the time t in the telegraph process, whereupon the polymer-contour length L maps to the total time T in the telegraph model. The parameter a quantifies the tendency of the tangent vectors to align with the channel [31]. The parameter g is the global persistence length [42], characterizing the typical distance between hairpin turns. These parameters map to those of the telegraph process as follows: a = v 0 and g = (2r) −1 . We measured how a and g depend on the physical parameters of the full three-dimensional problem from simulations of confined ideal polymers. It turns out that it suffices to determine just two curves [Figs. 3(a) and 3(b)], since a and g/ℓ P depend on D/ℓ P only. Now consider the effect of self-avoidance. In general it is very difficult to derive an expression for [x(s)]. But for a weakly self-avoiding polymer, the problem simplifies considerably when the channel is so narrow that interactions between the polymer and the channel wall cause the three-dimensional configurations to lose correlations. We show in the Supplemental Material [36] that (4) if w ≪ ℓ P . Here ℒ(x)dx is the total amount of contour in the interval [x; x þ dx] [43]. The parameter ε penalizes overlaps. It is determined by the probability that two short polymer segments overlapping in one dimension collide in three dimensions [ The average is over the conformations of the confined ideal polymer, and y and z are the transverse channel coordinates of a short polymer segment of length ℓ. Primed coordinates belong to a second, independent segment, and v ex is the excluded volume. The excluded volume depends on the segment orientation. If ℓ ≫ w, we have v ex = 2wℓ 2 sin θ, where θ is the angle between the two segments [17]. Figure 3(c) shows how ε depends on D/ℓ P , obtained by evaluating the average in Eq. (5) from three-dimensional simulations of confined ideal polymers [36]. A single curve is sufficient to determine how ε depends on the physical parameters, because εD 2 /w is a function of D/ℓ P only (see Supplemental Material [36]).
In the telegraph model self-avoidance is incorporated in the same way. Here, ℒ has units of (time)/(position). Equation (4) then shows that ε has units of (position)/(time) 2 . Since r has units of (time) −1 , and v 0 of (position)/(time), the only dimensionless combination of ε, r, and v 0 is (6) In the limit of large T, only α can have physical significance. Invoking our mapping between telegraph model and polymer problem, we conclude that α is given by (7) This parameter measures the expected number of overlaps between the two strands of a hairpin of length g. Equation (7) has two important consequences. First, Eq. (7) allows us to generalize the mapping to all channel widths up to . To show this, consider first the ideal part. Strictly speaking, the simple picture outlined above breaks down when D ~ ℓ P because the typical hairpin length g becomes of the same order as ℓ P . But consider how α changes as D approaches ℓ P from below. For w ≪ ℓ P , the parameter α decreases below unity before g = ℓ P is reached, and for small α the precise nature of the local conformations is irrelevant. All that matters is that the ideal part is a diffusion process with exponentially decaying correlations of v x (t). Similarly, the local probability of collision is still (ε/2)ℒ 2 (x)dx, because each segment pair collides independently. The latter assumption eventually breaks down at since the transversal segment coordinates become correlated. But up to this point Eq. (4) is valid, as is Eq. (5).
Second, observables that are dimensionless in the telegraph model can only depend on α, Eq. (7), in the limit of large L. This combination α is plotted in Fig. 3(d). It turns out that αD 2 /(ℓ P w) depends only on D/ℓ P [36]. Now consider the average extension X and the variance about that average σ 2 . In the telegraph model these quantities have units of (position) and (position) 2 , and for large values of L they must be proportional to L. We therefore conclude that the data must collapse as (8) Here f X and f σ are universal scaling functions that depend only on α. We can numerically compute the form of these functions simply by simulating the telegraph model [36].
We have compared our theory to direct numerical simulations (DNS) of three-dimensional, confined, self-avoiding wormlike chains [39] using the PERM algorithm [44,45]. Figures  4(a) and 4(b) show that our theory quantitatively captures the DNS results for all values of w/ℓ P tested [36], up to w/ℓ P = 0.4. This agreement is remarkable, as the theory assumes weak self-avoidance, w ≪ ℓ P . Fig. 1 and our theory. The theory not only collapses the experimental data, but provides good quantitative agreement, in particular with the most recent experiments [10,33,34]. There is some scatter of the experimental data around the theoretical curve, but this is expected because the nanofluidic experiments are quite difficult to control.

Figures 4(c) and 4(d) show the comparison between the results of the experiments summarized in
In the limit α ≪ 1 our theory allows us to map the problem to an uncorrelated weakly selfavoiding diffusion process [43,47]. This follows from the fact that the correlation function in the telegraph model, Eq. (3), decays to zero before the next collision occurs, for α ≪ 1. As a result, the ideal random walk is simply diffusive, with diffusion constant . This has two consequences.
First, for α ≪ 1 observables depend on v 0 and r only through the combination . Since the extension is linear in L for large L, we deduce that X/L can only depend on ε and in this limit. Since has units of (position) 2 /(time) while ε has units of (position)/(time) 2 in the telegraph model, we see that the only possible combination is X/L ∝ (ε ) 1/3 . This gives X/ (La)=f X (α)∝α 1/3 , explaining the power law in Figs. 4(a) and 4(c). For the variance we conclude that σ 2 /L ∝ , independent of ε [ Figs. 4(b) and 4(d)]. Alternatively, we can deduce these scalings by a mean-field argument, indicating that fluctuations of ℒ are negligible when α ≪ 1. Assuming that ℒ ~ T/X, we find for the extension distribution , whereupon d log P/dX = 0 yields X/(v 0 T) ~ α 1/3 .
For the variance we obtain .
Second, we can use the exact mathematical results derived in Ref.
Now consider the limit of large α. The extension X tends to La in this limit [11], since the frequency of hairpins tends to zero. The variance decays as σ 2 ~ α −3 , as Fig. 4 shows. To deduce this power law, we estimate the variance of the strongly extended polymer as (number of hairpins)× (hairpin extension) 2 . To determine the number of hairpins, note that the expected number of collisions for a hairpin of contour length h is ~αh/g. The resulting hairpin is therefore likely to survive the collision check only if h is of the order g/α or smaller. This requires a second switch of direction within the length g/α. This occurs with probability (g/α)/g = α −1 , so that the number of hairpins is (L/g)α −1 . To obtain the hairpin extension we multiply its contour length h ~ g/α by its alignment a, so that the typical hairpin extension becomes ~ga/α. Therefore, σ 2 ∝ (L/g)α −1 (ga/α) 2 = Lga 2 α −3 for large α.
For very large values of α, the theory fails [ Fig. 4(b)] because hairpins are so rare that alignment fluctuations (not included in the telegraph model) dominate the variance [11]. This correction is taken into account simply by adding the variance in the extreme Odijk limit [48]: Here d σ is a universal constant. By fitting the solid line in Fig. 4(b) for α > 10, we find d σ ≈ 0.09. We observe excellent agreement between this refined theory and the simulation data for ℓ P /w = 12. For the stiffer polymers (ℓ P /w = 36), still longer contour lengths are required to reach the large-L limit and to reduce the statistical error.
Finally, we show that our theory contains scaling laws derived earlier as particular asymptotic limits. In very narrow channels, a ≈ 1 and 〈sin θ〉 ≈ (D/ℓ P ) 1/3 . Using these approximations in Eqs. (5) and (6) gives (11) where C ≈ 1.95 is a constant [36]. The parameter ξ appears in Odijk's scaling theory [11] and the extension scales as X ~ ξ 1/3 [11] in this special limit. In wide channels, for , we have that and g ≈ ℓ P . Using diffusion approximations for the distribution of the polymer in the channel [36] gives (12) This is the result of Ref. [15], implying the same scaling X ~ D −2/3 that Odijk's scaling arguments [11] predict in this asymptotic limit. At first glance it is perhaps surprising that the two distinct limits Eqs. (11) and (12) are described by the same random-walk process. After all, the three-dimensional polymer conformations are entirely different in the two regimes, described by invoking deflection segments, hairpins, and blobs. Our universal theory, by contrast, rests on the fact that the macroscopic statistics of a weakly interacting random walk do not depend on the microscopic details of the process [49].
We can also conclude that the DNA experiments shown in Fig. 1 cannot obey the scalings X ~ D −2/3 or X ~ ξ 1/3 because the experiments do not satisfy the strong inequalities D ≫ ℓ P or D ≪ ℓ P (see Fig. S-5 in the Supplemental Material [36]), and therefore do not reach the asymptotic limits required for these power laws to emerge. Our theory shows, and Fig. 4(a) confirms, that X ~ α 1/3 for small values of α. But the parameter α depends upon D and ℓ P in an intricate way via Eq. (7), in general not in a power-law fashion.
In summary, we have shown that DNA confinement in nanochannels can be modeled by a telegraph process, collapsing all of the data in terms of a scaling variable α. Our theory brings to the fore universal properties of confined polymers in a good solvent in a way that is obscured by the prevailing scaling theories [9,11,12,14,15,41]. The distinction between deflection segments, hairpins, and blobs, which leads to the need to define separate regimes, is not necessary. Rather, the statistics of the confined polymer chain for adopt a universal behavior at sufficiently long length scales, independent of the details of the microscopic model.

Supplementary Material
Refer to Web version on PubMed Central for supplementary material.

FIG. 2.
Illustration of the telegraph process along the channel axis (x axis). The walk is one dimensional, but for clarity it is expanded vertically, to show the changes in direction that create hairpin configurations of the confined DNA molecule.