Home | About | Journals | Submit | Contact Us | Français |

**|**HHS Author Manuscripts**|**PMC3775584

Formats

Article sections

Authors

Related links

Phys Chem Chem Phys. Author manuscript; available in PMC 2013 September 17.

Published in final edited form as:

Published online 2009 May 14. doi: 10.1039/b901646e

PMCID: PMC3775584

NIHMSID: NIHMS416325

The publisher's final edited version of this article is available at Phys Chem Chem Phys

See other articles in PMC that cite the published article.

As a single DNA molecule is positively supercoiled under constant tension, its extension initially increases due to a negative twist-stretch coupling^{1, 2}. The subsequent attainment of an extension maximum has previously been assumed to be indicative of the onset of a phase transition from B- to scP-DNA^{2}. Here we show that an extension maximum in fact does not coincide with the onset of a phase transition. This transition is evidenced by a direct observation of a torque plateau using an angular optical trap. Instead we find that the shape of the extension curve can be well explained with a theory by John Marko^{3} that incorporates both DNA twist-stretch coupling and bending fluctuations. This theory also provides a more accurate method of determining the value of the twist-stretch coupling modulus, which has possibly been underestimated in previous studies that did not take into consideration the bending fluctuations^{2}. Our study demonstrates the importance of torque detection in the correct identification of phase transitions as well as the contribution of the twist-stretch coupling and bending fluctuations to DNA extension.

During various cellular processes DNA molecules often experience moderate stress and torque. These perturbations may be generated by motor enzymes and provide a mechanism for the regulation of DNA replication, DNA repair, transcription, and DNA recombination^{4–8}. Therefore, understanding tensile and torsional responses of DNA is essential to understanding how mechanical perturbations may regulate cellular activities.

A single DNA molecule can be extended under force and rotated under torque, as has been investigated using optical and magnetic tweezers and micropipettes during the past two decades^{9–16}. In particular, some of these studies revealed that the tensile and torsional responses of DNA are coupled. Initial analyses suggested that DNA undertwists when extended, indicating a positive twist-stretch coupling coefficient^{3, 17}. More recent studies by Gore et al.^{1} and Lionnet et al.^{2} using magnetic tweezers, as well as our own work using angular optical trapping^{15} provide compelling evidence to the contrary: DNA overtwists when extended and, conversely, DNA extends when overtwisted. However, there remains some ambiguity in identifying the signatures of twist-stretch coupling in the extension curves when the DNA was overtwisted under constant force. While Lionnet et al.^{2} interpreted a peak in this curve to be indicative of the onset of a phase transition from B- to supercoiled P- (scP-) DNA, our earlier work indicated that the peak of the curve and the phase transition did not coincide at the particular force examined^{15}. In general, there was a lack of understanding of the nature of the extension peak location and its relation to both twist-stretch coupling and phase transitions.

The goals of this work are to differentiate between the signatures of twist-stretch coupling and those of phase transitions, and to provide an explanation for the existence of the maximum in the extension signal. To this end, we carried out DNA torsional experiments using an angular optical trap in conjunction with nanofabricated quartz cylinders so that torque, angle, force, and extension of a DNA molecule were simultaneously measured during DNA supercoiling, using previously described methods^{15, 16}. An important advantage of this approach is the direct detection of the torque signal, allowing unambiguous identification of the onset of the phase transition where torque plateaus.

The DNA template used in this study was constructed using previously described protocols^{16}. In brief, a 4218-bp piece of DNA was ligated at one end to a short (60-bp) oligonucleotide labeled with multiple digoxygenin (dig) tags and at the other end to an oligonucleotide of the same length labeled with multiple biotin tags.

The experimental configurations and procedures were similar to those described previously^{15, 16}. In brief, prior to a measurement, DNA molecules were torsionally constrained at one end to streptavidin-coated nanofabricated quartz cylinders^{15} and at the other end to an anti-dig coated coverslip. All experiments were performed in phosphate-buffered saline (157 mM Na^{+}, 4 mM K^{+}, 12 mM PO_{4}^{3−}, 140 mM Cl^{−}, pH=7.4) at 23±1 °C. The experiment began with a torsion-free DNA molecule which was held under constant tension. The DNA was first slightly undertwisted and then overwound via a steady rotation of the input laser polarization at 5 Hz, so as to explore a range of DNA supercoiling. During this time, torque, angular orientation, position, and force of the cylinder as well as the location of the coverglass were simultaneously recorded. The torque exerted on the DNA was measured from the torque exerted on the cylinder by the optical trap after subtracting the viscous drag torque of the rotating cylinder.

Figure 1 shows representative single traces of extension and the corresponding torque as a function of number of turns added to the DNA at three different applied forces (1.9, 7.7, and 9.6 pN). The number of turns was also converted to the degree of supercoiling σ, defined as the number of turns added to dsDNA divided by the number of naturally occurring helical turns in the given dsDNA. The DNA extension is also shown as the relative extension, defined as the extension normalized to the contour length of DNA template in its B-form.

Examples of extension and torque versus turn number during DNA supercoiling. DNA molecules of 4.2 kbp in length were wound at 5 Hz under constant forces: A) 1.9 pN; B) 7.7 pN; C) 9.6 pN. Data were collected at 2 kHz and averaged with a sliding box window **...**

Some overall features of these data are summarized below. At the beginning of each trace (σ ~ −0.02) DNA was in its B-form. As the DNA was twisted, torque increased linearly while extension remained approximately constant. At lower forces (< 6 pN) this continued until the DNA buckled to form a plectoneme, evidenced by a sudden decrease in extension and a concurrent plateau in torque, whose value was force-dependent^{16}. Under this low force range, the signatures for the onset of the buckling transition may be recognized in either extension or torque. At higher forces we used (6 pN < *F* < 10 pN), as the molecule was overtwisted, DNA underwent a scP-DNA transition instead, as previously observed by Allemand et al.^{13} The onset of this transition was only evidenced by a sudden torque plateau around 40 pN nm^{14, 15}. These phase transitions may be thought of as first-order phase transitions since two separate phases may coexist and can be transformed from one to another by simply changing the twist in the DNA^{18}. Because there were no concurrent distinct features in the extension signal, the ability to monitor torque signal was essential in unambiguously locating the onset of this transition. In addition, the buckling torque showed a strong force-dependence as previously observed^{16}; whereas the scP transition torque showed little force-dependence.

Several features relevant to twist-stretch coupling were also immediately evident. First, near σ = 0, the extension curve had a small but positive slope, which implied a negative twist-stretch coupling coefficient, i.e., DNA was extended when overtwisted. This observation was in accordance with previous studies^{1, 2, 15}. Second, upon overtwisting, extension reached a maximum (indicated by an arrow in Fig. 1) at σ_{z_max}, as has also been previously observed^{2, 15}. In addition, as the force increased, σ_{z_max} increased while the magnitude in the curvature of the extension at σ_{z_max} gradually decreased. Third, the location of the extension maximum σ_{z_max} clearly did not coincide with the onset of the phase transition (indicated by a dashed line on the right) where torque began to plateau. This was the case for the range of forces we examined. Therefore the location of the maximum is not indicative of an onset of a phase transition, contrary to what has been previously reported with the magnetic tweezers experiments, where the lack of torque signal might have complicated the interpretation of the results^{2}.

To understand the nature of the extension signal for the B-form DNA, we performed a detailed analysis of both the extension and torque signals to gain insights in the twist-stretch coupling coefficient and the location of the extension maximum. We followed the analysis developed by Marko^{3}. This theory also takes into account the contribution from bending fluctuations to DNA extension, using a treatment similar to that of Moroz and Nelson^{17}. For a DNA molecule held under constant σ and force *F*, the free energy *G* includes contributions from bending, stretching, twisting, and twist-stretch coupling, and can be expressed as:

$$\frac{G}{{k}_{\mathrm{B}}T{L}_{0}}=\frac{1}{{L}_{p}}\sqrt{\frac{{L}_{p}F}{{k}_{B}T}-\frac{1}{4}{\left({C}_{0}{\omega}_{0}\sigma +\frac{gF}{{K}_{0}}\right)}^{2}}+\frac{{C}_{0}}{2}{{\omega}_{0}}^{2}{\sigma}^{2}-\frac{F}{{k}_{B}T}-\frac{{k}_{B}T}{2{K}_{0}}{\left(\frac{F}{{k}_{B}T}-g{\omega}_{0}\sigma \right)}^{2},$$

(1)

where *L*_{p} is the bending persistence length, *K*_{0} the stretch modulus, *C*_{0} the twist persistence length, *g* the twist-stretch coupling modulus (unitless), *k*_{B}*T* the thermal energy, *L*_{0} the contour length, and ω_{0} = 2π/3.57 nm^{−1} the natural twist rate. This theory predicts DNA extension *z* as a function of σ and applied force *F*:

$$\frac{z}{{L}_{0}}=-\frac{\partial (G/{L}_{0})}{\partial F}{|}_{\sigma}=1+\frac{F}{{K}_{0}}-\frac{{k}_{\mathrm{B}}Tg{\omega}_{0}\sigma}{{K}_{0}}-\sqrt{\frac{{k}_{\mathrm{B}}T}{4{L}_{\mathrm{p}}F}}\xb7\frac{1-\frac{g{C}_{0}\sigma {\omega}_{0}+{g}^{2}F/{K}_{0}}{2{L}_{\mathrm{p}}{K}_{0}}{k}_{\mathrm{B}}T}{\sqrt{1-\frac{{k}_{\mathrm{B}}T}{4{L}_{\mathrm{p}}F}{\left({C}_{0}{\omega}_{0}\sigma -\frac{gF}{{K}_{0}}\right)}^{2}}}.$$

(2)

We used Equation (1) to obtain a prediction of torque τ as a function of σ and applied force *F*:

$$\tau =\frac{1}{{\omega}_{0}}\frac{\partial (G/{L}_{0})}{\partial \sigma}{|}_{F}={k}_{\mathrm{B}}T{\omega}_{0}\sigma \left({C}_{0}-\frac{{k}_{\mathrm{B}}T{g}^{2}}{{K}_{0}}\right)+\frac{{k}_{\mathrm{B}}TgF}{{K}_{0}}-\frac{\frac{{k}_{\mathrm{B}}T{C}_{0}}{4{L}_{\mathrm{p}}}\left({C}_{0}{\omega}_{0}\sigma +\frac{gF}{{K}_{0}}\right)}{\sqrt{\frac{{L}_{\mathrm{p}}F}{{k}_{\mathrm{B}}T}-\frac{1}{4}{\left({C}_{0}{\omega}_{0}\sigma +\frac{gF}{{K}_{0}}\right)}^{2}}}.$$

(3)

Thus, Equations (2) and (3) form a complete set of relations that fully describe force, extension, torque, and twist for a B-form DNA.

Below, we will focus on the analysis of the extension signal. Equation (2) contains several terms and its last term is due to the contribution of DNA bending fluctuations in the presence of the twist-stretch coupling. We found that Equation (2) predicts the existence of a maximum in the extension curve, reflecting an interplay between twist-stretch coupling and bending fluctuations. In the absence of any twist-stretch coupling, the maximum is centrally located, i.e., if *g* = 0, σ_{z_max} = 0. Twist-stretch coupling shifts the location of the maximum away from the center (Supplementary Materials). If *g* > 0, then σ_{z_max} < 0; and if *g* < 0, then σ_{z_max} > 0. The absolute value of σ_{z_max} increases with an increase in force.

Notice that if bending fluctuations are neglected, Equation (2) is then simplified to:

$$\frac{z}{{L}_{0}}=1+\frac{F}{{K}_{0}}-\frac{{k}_{\mathrm{B}}Tg{\omega}_{0}\sigma}{{K}_{0}}.$$

(4)

In contrast to Equation (2), Equation (4) predicts a linear relation between extension and the degree of supercoiling: if *g* > 0, then slope < 0 ; and if *g* < 0, then slope > 0 It does not predict a maximum in the extension curve.

We performed a fit of Equation (2) to our extension data for the B-form DNA with *g* as the only fit parameter. Other parameters were taken from previous measurements obtained under similar experimental conditions: *K*_{0} = 1200 pN^{9, 11}, *C*_{0} = 100 nm^{14, 16}, and *L*_{p} = 43 nm^{9, 11, 16}. Our data were well fit by Equation (2) and some examples are shown in Figure 1 (red). In particular, the experimental observations that σ_{z_max} > 0 and increases with increasing force were indicative of a negative value for the twist-stretch coupling modulus *g*. Figure 2A shows *g* values obtained from the fits under various forces. Over the force range examined, *g* was essentially independent of the force: *g* = −21±1 (mean ± sem, N = 41). The *g* value is sensitive to the value used for *C*. For example, a ±10% uncertainty in *C* will result in a ±15% uncertainty in *g*.

A) Measurement of the twist-stretch coupling modulus *g*. For each force, the twist-stretch coupling modulus was determined using Equation (2) from 7–11 traces of data. The mean of the modulus *g* is shown as the solid horizontal line. For comparison, **...**

If Equation (4) is used to obtain a *g* value instead, the magnitude of *g* is underestimated by ~ 20% over the range of forces examined (dashed line in Figure 2A). This difference does not vanish even with an increase in force, since the geometric coupling between bending and writhe fluctuations leads to a decrease in extension^{3}. Interestingly, Lionnet et al.^{2} used Equation (4) to fit their extension data and obtained *g* = −16 ± 7 (mean ± sd, N > 36), whose magnitude is ~ 20% lower than our measured value. On the other hand, Gore et al.^{1} made a measurement of twist angle as a function of extension change when the DNA was not torsionally constrained and obtained *g* = −22±5 (mean ± sem, N = 4), more in accord with our measured value. In the analysis by Gore et al.^{1}, the bending fluctuations were also neglected resulting in a similar linear approximation. Thus one might expect that the magnitude of *g* was similarly underestimated. Careful analysis using Equations (2) and (4) reveals that this effect was almost completely canceled by another linear approximation made to convert force to extension change. Taking together all these results (summarized in Fig. 2B), the twist-stretch coupling modulus *g* should be about − 21.

The extension maximum, however, can only be explained using Equation (2) and not Equation (4). Figure 3 plots σ_{z_max} obtained as a function of force. For comparison, we also plotted the critical σ* _{c}* values for phase transitions to plectoneme (scB-) or scP-DNA. As this figure indicates, when DNA is positively supercoiled under moderate forces (2 pN <

Twist-stretch coupling should also alter the torque signal. We found that consideration of the twist-stretch coupling as in Equation (3) will lower the expected torque value at most by ~ 1 pN nm (Supplementary Materials). This is below the uncertainty of our experimental determination of torque.

The analysis described in this work may not be valid at forces significantly higher than those used in the current work. In the absence of torsional constraints, Equation (3) predicts a monotonic increase in overtwisting angle with an increase in force. However, Gore et al.^{1} found that the twist angle starts to decrease with force above ~ 30 pN.

It is also worth mentioning that we do not consider sequence-dependent effects on twist-stretch coupling or phase transitions here. Prior simulation work ^{2, 19} suggests that DNA sequence may modulate the twist-stretch coupling. Future experiments with more refined measurements may help verify this prediction.

We thank members of the Wang lab for critical reading of the manuscript, S. Forth for help with the angular trap instrument and DNA template preparations, and T. Lionnet, V. Croquette, and D. Bensimon for helpful comments. M.D.W. wishes to acknowledge support from NSF grant (MCB-0820293), NIH grant (R01 GM059849), and the Cornell Nanobiotechnology Center.

1. Gore J, Bryant Z, Nollmann M, Le MU, Cozzarelli NR, Bustamante C. Nature. 2006;442:836–839. [PubMed]

2. Lionnet T, Joubaud S, Lavery R, Bensimon D, Croquette V. Phys Rev Lett. 2006;96:178102. [PubMed]

3. Marko JF. Physical Review E. 1998;57:2134–2149.

4. Liu LF, Wang JC. Proc Natl Acad Sci U S A. 1987;84:7024–7027. [PubMed]

5. Travers AA, Thompson JM. Philos Transact A Math Phys Eng Sci. 2004;362:1265–1279. [PubMed]

6. Koster DA, Croquette V, Dekker C, Shuman S, Dekker NH. Nature. 2005;434:671–674. [PubMed]

7. Cozzarelli NR, Cost GJ, Nollmann M, Viard T, Stray JE. Nat Rev Mol Cell Biol. 2006;7:580–588. [PubMed]

8. Garcia HG, Grayson P, Han L, Inamdar M, Kondev J, Nelson PC, Phillips R, Widom J, Wiggins PA. Biopolymers. 2007;85:115–130. [PMC free article] [PubMed]

9. Smith SB, Cui Y, Bustamante C. Science. 1996;271:795–799. [PubMed]

10. Strick TR, Allemand JF, Bensimon D, Bensimon A, Croquette V. Science. 1996;271:1835–1837. [PubMed]

11. Wang MD, Yin H, Landick R, Gelles J, Block SM. Biophys J. 1997;72:1335–1346. [PubMed]

12. Strick TR, Allemand JF, Bensimon D, Croquette V. Biophys J. 1998;74:2016–2028. [PubMed]

13. Allemand JF, Bensimon D, Lavery R, Croquette V. Proc Natl Acad Sci U S A. 1998;95:14152–14157. [PubMed]

14. Bryant Z, Stone MD, Gore J, Smith SB, Cozzarelli NR, Bustamante C. Nature. 2003;424:338–341. [PubMed]

15. Deufel C, Forth S, Simmons CR, Dejgosha S, Wang MD. Nat Methods. 2007;4:223–225. [PubMed]

16. Forth S, Deufel C, Sheinin MY, Daniels B, Sethna JP, Wang MD. Phys Rev Lett. 2008;100:148301. [PMC free article] [PubMed]

17. Moroz JD, Nelson P. Proc Natl Acad Sci U S A. 1997;94:14418–14422. [PubMed]

18. Marko JF, Siggia ED. Physical Review E. 1995;52:2912–2938. [PubMed]

19. Lionnet T, Lankas F. Biophys J. 2007;92:L30–L32. [PubMed]

20. Marko JF. Phys Rev E Stat Nonlin Soft Matter Phys. 2007;76:021926. [PubMed]

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |