| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Department of Chemical Engineering, University of California, Berkeley, California
Correspondence: Address reprint requests to Dr. Harvey W. Blanch, Dept. of Chemical Engineering, University of California, Berkeley, CA 94720. Tel.: 510-643-1387; Fax: 510-643-1228; E-mail: blanch{at}socrates.berkeley.edu.
| ABSTRACT |
|---|
|
|
|---|
R, scales with chain length as
R
L1.9±0.17 when analyzed by the data collapse method, and as
R
L2.17±0.17 when analyzed by multiexponential fit. | INTRODUCTION |
|---|
|
|
|---|
Epifluorescence microscopy has shown that the conformation of long DNA molecules during electrophoretic separations is highly dynamic, fluctuating between collapsed and highly extended states (Schwartz and Michael, 1989
; Smith et al., 1989
; Sunada and Blanch, 1998a
). Because DNA itself is a polymer, the description of these dynamics has often been pursued using the classical theories of polymer dynamics. Rouse (1953)
developed a bead-spring model of a polymer fluctuating around its equilibrium conformation, a model that has been used extensively since. The model assumes that each bead has the same friction coefficient (i.e., the polymer is free-draining), the restoring force of the springs is linear, and the solvent is athermal. The position of each bead with time is found from the differential equation
![]() | (1) |
the friction factor of the bead, and g is a noise term that accounts for Brownian motion. From the solution to Eq. 1, the motion of a polymer, ranging from the diffusion of the whole molecule to rearrangements of single segments, can be described. For example, the time correlation function of the end-to-end vector,
P(t) · P(0)
for a polymer with N segments, is described by a sum of independent, normal modes, each with its own relaxation time,
R:
![]() | (2) |
![]() | (3) |
-solvent, the characteristic relaxation time should scale with the square of the polymer length. However, the Rouse model neglects both hydrodynamic interactions between segments and the solvent-mediated excluded volume interaction. Zimm (1956)
s:
![]() | (4) |
Solvent quality has been explicitly introduced into the Zimm model (Doi and Edwards, 1986
), but a simple scaling argument first put forward by de Gennes (1976)
gives the same functional form for the dynamical properties. De Gennes proposed that the dynamical properties of polymers will scale similarly to their static properties, depending on the solvent quality through the scaling parameter
(which is equal to 1/2 in a
-solvent and
3/5 in a good solvent). For example, if the number of segments in a chain is changed from N to N/
(where
is a constant) then the physical properties can be held constant by simply changing the segment length from b to b
(since these properties do not depend on the local structure of the chain). By this transformation, the invariant static and dynamical properties should change from A to A
x, where x is a scaling parameter that depends on
. It can be shown that for a Rouse chain and a Zimm chain, respectively, that
![]() | (5a) |
![]() | (5b) |
The hydrodynamic features of any electrophoretic separation are going to be system-dependent. The polymer identity, concentration, and possibly polydispersity could all affect the manner in which single molecules of DNA will stretch and relax. However, whether the presence of neutral polymer at a concentration relevant to electrophoretic separations will screen these hydrodynamic interactions, and whether this screening will significantly affect the electrophoretic stretch and relaxation, remain open questions. This article addresses the magnitude and result of these hydrodynamics on DNA stretch and relaxation in a semidilute polymer solution. We use epifluorescence light microscopy to measure the end-to-end distance of DNA molecules tethered in an electric field. Data on both the steady stretching as a function of electric field and dynamic relaxation after the field is removed are presented for DNA molecules of varying length. The data are interpreted using several theories of polymer physics.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Device fabrication
We wish to examine only the effects of the neutral polymer in solution on the electrophoretic stretch and relaxation of single DNA molecules. Consequently, we need to have the DNA molecules tethered far enough away from any solid surfaces so that the hydrodynamic screening is dominated by the neutral polymer and not by the surface. As will be described in detail later, the hydrodynamic screening length at the polymer concentration used here is
20 nm. The DNA was therefore attached to a microstructure consisting of a 500-nm-high, 2-µm-wide gold wall deposited on a silicon wafer. The P-type silicon wafers were first cleaned by immersing in piranha (5:1, H2SO4:H2O2) solution for 10 min, rinsed in a tank of nanopure deionized water until resistivity reached 13 M
, and blown dry under a stream of nitrogen. A 1.5-µm-thick layer of silicon oxide was grown on the surface of the wafer by wet thermal oxidation at 1000°C. Positive photoresist was then patterned on the wafer by contact lithography to leave an array of 2-µm-wide lines of bare oxide exposed. A 7-nm layer of chromium was deposited by thermal evaporation at 5 x 106 torr, followed by a 500-nm layer of gold, without breaking the vacuum. After liftoff of the metals in a bath of acetone, the wafers were left patterned with a series of parallel 2-µm-wide lines, 500 nm in height.
The silicon oxide surfaces of the microfabricated device will generate electro-osmotic flow under our experimental conditions. Therefore, before immobilizing the DNA all surfaces were coated with a brush of polyacrylamide, which we have previously shown eliminates electroosmotic flow in microfluidic channels (Ferree and Blanch, 2003
). Before use, the microcells and glass coverslips were cleaned with piranha solution then rinsed thoroughly in deionized water. All surfaces were then coated with a self-assembled monolayer of methacryloxypropyltrimethoxysilane (MPTS, Sigma, St. Louis, MO) by immersing in a 2% solution of MPTS in 95% ethanol (0.3% acetic acid) for 10 min, rinsing with ethanol, and curing at 110°C for 5 min. A polyacrylamide brush was then polymerized from the surface by soaking overnight in a 8% solution of acrylamide (Bio-Rad Laboratories, Hercules, CA) in 0.5x TBE buffer, pH 8.0 (45 mM tris, 45 mM borate, 1 mM EDTA, all Sigma) using ammonium persulfate as a free-radical generator and TEMED (n,n,n',n'-tetramethylenediamine, Sigma) as a stabilizer. The unbound polymer was then removed by rinsing with TBE.
DNA immobilization
The chromosome of the
-phage was labeled at one end with a 12 basepair 3'-biotinylated oligonucleotide (Genemed Synthesis, South San Francisco, CA), and purified by filtering through a Microcon 50 microcentrifuge filter (Amicon, Beverly, MA). Concatamers of the molecule were made by incubating a solution of
with T4 DNA ligase (New England Biolabs, Beverly, MA) for several hours at room temperature before labeling. The biotinylated DNA was then immobilized on the gold surfaces by a technique reported in detail previously (Zimmerman and Cox, 1994
). Briefly, a self-assembled monolayer of thiolated biotin (Biotin-HPDP, Pierce-Endogen, Rockford, IL) was formed on the gold lines, to which was bound the protein neutravidin (Pierce-Endogen). Since neutravidin is a tetramer, the immobilized protein retains free biotin-binding sites that are used to bind free molecules of biotinylated DNA.
Fluorescence imaging
A 30-µm-deep channel was made by separating a glass coverslip and the silicon device with Kapton tape (3M, Minneapolis, MN) and sealing with nail polish. The channel was filled with a solution of
-molecules, labeled with the fluorophore YOYO-1 at an 8:1 basepair/dye ratio (Molecular Probes, Eugene, OR) in HEC/TBE, by capillary action. Photobleaching and photoscission of the DNA were minimized by adding 3% ß-mercaptoethanol, 15 µg/ml catalase, 50 µg/ml glucose oxidase, and 0.1% glucose (all Sigma). All experiments were performed on an upright Zeiss epifluorescence microscope (Carl Zeiss, Zurich, Switzerland) with a 100x oil immersion lens (NA = 1.3). The samples were illuminated with a xenon arc lamp using a standard FITC filter set (450 < ex < 490, em > 510), intensified and recorded to a computer hard drive with a low light vidicon camera. For steady-state measurements, images were captured at 10 fps and the extension averaged over 10 frames for each electric field. For the relaxation measurements of longer DNA (>25 µm in length) images were captured at 5 fps and 34 full relaxation curves averaged for each molecule. For the relaxation of shorter molecules of DNA (<25 µm) images were captured at 10 fps to ensure that the faster relaxation dynamics were acquired.
| RESULTS |
|---|
|
|
|---|
4.
The entanglement threshold can be found from solution viscosity measurements by a number of different techniques. Often, the log of the specific viscosity is plotted as a function of the log of the concentration and the entanglement threshold designated as the point where the data deviate from linearity (Grossman and Soane, 1991
). From the data in Fig. 1 we designated c* = 0.4 wt % as the polymer overlap concentration. As a second estimate, we also calculated the entanglement threshold from the intrinsic viscosity of the polymer in TBE. We found the intrinsic viscosity of HEC 250 from a dual Huggins-Kraemer plot to be [
] = 315 ± 20 ml/g. From geometric arguments it can be shown that c* = 2.5[
]1 = 0.79 wt % (Allcock et al., 2003
). This quantity is approximately double our estimate from Fig. 1; the difference is probably due to the sample polydispersity and the approximate nature of the geometric arguments. We therefore used a concentration of 0.7 wt % HEC 250 in TBE as the solution to study the stretch and relaxation of DNA, a concentration above the entanglement threshold, but still in the semidilute regime and a concentration relevant to electrophoretic separations (Barron and Blanch, 1995
).
|
To ensure that the HEC network did not form long-lived entanglements that would constrain DNA to a reptation tube, we performed dynamic rheology measurements on HEC 250 solutions at a concentration much higher than that used in our DNA stretching experiments. Linear viscoelastic theory predicts an entanglement plateau in a plot of the storage modulus versus frequency (Larson, 1999
). Polymers that are too short or solutions that are too dilute will instead show a steadily increasing modulus with frequency. The inset of Fig. 1 indicates that even at a concentration of 4 wt. % (
6 times more concentrated than the solutions used in the DNA stretching experiments reported here), there are no long-lived HEC entanglements.
Steady stretching in a polymer solution
In free electrophoresis DNA behaves as a free-draining polymer because the long-range hydrodynamic perturbation caused by the molecule moving through the quiescent solution is cancelled out by the electro-osmosis of the electric double layer around the DNA. However, during gel electrophoresis when a DNA molecule interacts with a gel fiber, it is temporarily held stationary. In this scenario the hydrodynamic perturbation is removed and the electro-osmotic flow around the charged DNA backbone creates a long-ranged hydrodynamic interaction along the chain. We recently reported on the steady stretching of single tethered molecules of DNA in a gel-free electrophoresis buffer (Ferree and Blanch, 2003
). It was observed that the hydrodynamic interaction results in tethered DNA stretching the same in an electric field, E, as it would stretch in a flow field with a velocity equal to µE (where µ is the free solution electrophoretic mobility). When the surrounding medium is replaced with a solution of uncharged polymers in electrophoresis buffer, much of this hydrodynamic interaction should be screened. The hydrodynamic interaction acts to lessen the force felt by a DNA molecule, so when it is screened by the polymers the molecule should stretch even further at a given electric field. We measured the extension of single DNA molecules tethered by one end to a chemically modified gold surface in a microfluidic device as a function of electric field from <1 V/cm up to 25 V/cm. Since DNA has been shown to behave as a worm-like chain (Bustamante et al., 1994
; Marko and Siggia, 1995
), which approaches its contour length with an inverse square-root dependence on force, the contour length of each molecule was found as the y intercept of a plot of extension versus E0.5. Fig. 2 is a plot of the extension of several molecules as a function of the electric field in both free solution and in 0.7% HEC. In the solution of HEC 250, the DNA extension is seen to increase above the free solution value as a result of hydrodynamic screening.
|
EL1, where x is the extension of a molecule length L in an electric field E. A plot of the fractional extension for five different chains as a function of the electric field multiplied by the contour length is shown in Fig. 3. The collapse of all of the data to one universal curve shows that the HEC solution used here does screen out nearly all of the hydrodynamic interaction along the DNA backbone.
|
25 V/cm) to a sample of tethered DNA, then suddenly switched off the field and followed the relaxation of single molecules until they reached an equilibrium size. A time lapse of the relaxation of a 51.5-µm-long DNA molecule is shown in Fig. 4. The relaxation dynamics depend on the molecular conformation, so data from three or four separate relaxation experiments were averaged. Fig. 5 is a plot of extension versus time for several molecules of different lengths in a solution of 0.7% HEC 250. By fitting a universal curve to the data in Fig. 3, we were able to estimate the contour length of each molecule from the extension before turning off the electric field.
|
|
-solvent should scale with the square of the chain length. We analyzed the relaxation data by two model independent methods: data collapse and multiexponential fit. The data collapse method has proven effective in previous studies of relaxation of DNA in free solution and relies on the different data sets having similar shapes (i.e., same functional form with a similar number of relaxation modes; Perkins et al., 1994
t, with the length of the test chain will reflect the scaling of the normal relaxation mode. In Fig. 6 a, the data for relaxation of a 37.5-µm template chain is rescaled to fit data for 16.2- and 46.8-µm test chains. A plot of log(
t) versus log(test chain length) should lie on a line with slope equal to the dynamical scaling constant as shown in the inset of Fig. 6. We repeated this technique using three template chains of different lengths: 16.2, 37.5, and 55.9 µm. We found the average of the dynamical scaling exponent for relaxation from these three calculations to be 1.93 ± 0.20.
|
R, would scale most closely with the Rouse model. Two-thirds of our data were described best by a sum of just two exponentials (plus a baseline offset), whereas the other third were split evenly between one and three exponentials. By performing a least-squares linear regression on data of log(length) versus log(
R), we found the dynamical scaling constant to be 2.17 ± 0.17 (Fig. 7).
|
| DISCUSSION |
|---|
|
|
|---|
![]() | (6) |
is the fluid viscosity and p is the pressure. Equation 6 is solved along with the continuity equation for an incompressible fluid by assuming that the prefactor of the second term is equal to 1/S2, where S is the hydrodynamic screening length. This calculation results in new equations for the friction factors of each DNA segment, as well as for the hydrodynamic perturbation around each segment. With the fluid perturbation around each segment known, the total force is equal to the bare electrophoretic force plus the hydrodynamic interaction with all other segments in the chain. The tension at each segment is then used to find a new orientation with respect to the electric field from a worm-like chain entropic elasticity formula. The calculation then begins againfinding the tension at each segment based on the new chain configurationand continues until self-consistency is found. The only unknown parameter in the theory is the screening length.
The screening length could be used as an adjustable parameter, but we would prefer to use Stigter's model to predict the behavior of DNA in polymer solutions. We can derive a value for the hydrodynamic screening length of the polymer solution using the effective medium theory developed by Freed, Edwards, and Muthukumar (Edwards and Freed, 1974
; Edwards and Muthukumar, 1984
; Freed and Edwards, 1974
). The fluid is assumed to obey the Navier-Stokes equation under the assumptions of incompressibility and creeping flow. The polymer is described by a Rouse-like bead-spring model and its motion is explicitly coupled with the motion of the fluid. At low frequencies it is found that
![]() | (7) |
-potential equal to 0.079 V (Ferree and Blanch, 2003
As described in the Introduction, the Rouse model predicts that the normal modes for relaxation of a free-draining polymer will scale with N1+2
where
is 1/2 in a
-solvent and
3/5 in a good solvent. Kantor et al. (1999)
used fluorescence microscopy to follow the relaxation of single molecules of yeast chromosomal DNA stretched by an alternating electric field in the interface between an agarose gel and a glass coverslip, a technique they termed optical contour maximization. A log-log plot of
R versus length for five different molecules resulted in a scaling exponent of 1.45. There are several aspects of the work which are questionable. First, the DNA molecules studied were very long, and therefore extremely sensitive to shearing and breakage. However, the contour length of each DNA molecule was assumed to be equal to the length of the molecule as found from pulsed field gel electrophoresis. This method of determining the contour length is very inaccurate since many molecules will break during handling and this may be the main source of the large error in their data. Second, the environment in which the DNA molecules were stretched is inhomogeneous and uncharacterized. The authors state that the molecules are trapped in a region very near to the glass where the agarose concentration is low, explaining the Zimm-like scaling of the relaxation time. However, since the solvent in these experiments is good, the scaling exponent should actually be closer to 1.8, as was found by Perkins et al. (1994)
. Because the environment around the DNA molecules is not well characterized and the technique is highly irreproducible, optical contour maximization is not a reliable method for measuring the dynamics of single DNA molecules.
We have found the relaxation scaling time in a semidilute solution of HEC 250 to be 1.93 ± 0.20 using data collapse and 2.17 ± 0.17 using multiexponential fitvalues significantly higher than found in free solution (1.7 ± 0.1; Perkins et al., 1994
). Since the electrophoresis buffer used here is a good solvent for DNA, the expected value for the scaling exponent of the relaxation time is 2.2 for Rouse dynamics. This prediction is higher than our experimental measurements, although still within the experimental uncertainty for the multiexponential fit. Since the length scale for the hydrodynamic perturbation (the persistence length,
50 nm) is of the same order of magnitude as that for the hydrodynamic screening length, hydrodynamic interactions between segments are not entirely screened and will lead to a decrease in the scaling exponent. This interaction could explain the difference between the scaling prediction and our results. It has also been proposed that the dynamical scaling constants in a good solvent may approach the theoretical value slowly as the polymer length increases (Weill and Des Cloizeaux, 1979
). Our experiments involve DNA molecules ranging from
300 to 1100 persistence lengths, which are short compared to many synthetic polymers, and may also lead to a decrease in the scaling exponent.
| CONCLUSIONS |
|---|
|
|
|---|
The dynamic relaxation of
30 separate molecules ranging in length from 300 to 1100 persistence lengths was measured. Because the HEC samples studied are too short to form long-lived reptation tubes, the DNA can be described by the bead-spring models of Rouse and Zimm. The relaxation time of the slowest mode,
R, was found to scale with length as
R
L1.93±0.20 using data collapse and
R
L2.17±0.17 using multiexponential fit. These values agree well with the Rouse model prediction.
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
This work was supported by the Department of Energy under grant FG03-94ER-14456.
Submitted on January 12, 2004; accepted for publication March 26, 2004.
| REFERENCES |
|---|
|
|
|---|
Bandrup, J., and E. H. Immergut. 1989. The Polymer Handbook. Wiley Intersciences, Ithaca, NY.
Barron, A. E., and H. W. Blanch. 1995. DNA separations by slab gel and capillary electrophoresis: theory and practice. Separ. Purific. Meth. 24:1118.
Bustamante, C., J. F. Marko, E. D. Siggia, and S. Smith. 1994. Entropic elasticity of
-phage DNA. Science. 265:15991600.
De Gennes, P. G. 1976. Scaling Concepts in Polymer Physics. Cornell University Press, Ithaca, NY .
Doi, M., and S. F. Edwards. 1986. The Theory of Polymer Dynamics. J. Birman, S. F. Edwards, C. H. Llewellyn Smith, and M. Rees, editors. Clarendon Press, Oxford, UK.
Duke, T., and J.-L. Viovy. 1994. Theory of DNA electrophoresis in physical gels and entangled polymer solutions. Phys. Rev. E. 49:24082416.[CrossRef]
Edwards, S. F., and K. F. Freed. 1974. Theory of the dynamical viscosity of polymer solutions. J. Chem. Phys. 61:11891202.[CrossRef]
Edwards, S. F., and M. Muthukumar. 1984. Brownian dynamics of polymer solutions. Macromolecules. 17:586596.[CrossRef]
Fernandes, M. X., A. Ortega, M. C. Lopez Martinez, and J. Garcia de la Torre. 2002. Calculation of hydrodynamic properties of small nucleic acids from their atomic structure. Nucl. Acids Res. 30:17821788.
Ferree, S., and H. W. Blanch. 2003. Electrokinetic stretching of DNA. Biophys. J. 85:25392546.
Freed, K. F., and S. F. Edwards. 1974. Polymer viscosity in concentrated solutions. J. Chem. Phys. 61:36263633.[CrossRef]
Grossman, P. D., and D. S. Soane. 1991. Experimental and theoretical studies of DNA separations by capillary electrophoresis in entangled polymer solutions. Biopolymers. 31:12211228.[CrossRef][Medline]
Kantor, R. M., X.-H. Guo, E. J. Huff, and D. C. Schwartz. 1999. Dynamics of DNA molecules in gel studied by fluorescence microscopy. Biochem. Biophys. Res. Commun. 258:102108.[CrossRef][Medline]
Kirkwood, J. G., and J. Riseman. 1948. The intrinsic viscosities and diffusion constants of flexible macromolecules in solution. J. Chem. Phys. 16:565573.[CrossRef]
Larson, R. G. 1999. The structure and rheology of complex fluids. K. E. Gubbins, editor. Oxford University Press, New York, NY.
Marko, J. F., and E. D. Siggia. 1995. Stretching DNA. Macromolecules. 28:87598770.[CrossRef]
Perkins, T. T., S. R. Quake, D. E. Smith, and S. Chu. 1994. Relaxation of a single DNA molecule observed by optical microscopy. Science. 264:822826.
Provencher, S. W. 1976. An eigenfunction expansion method for the analysis of exponential decay curves. J. Chem. Phys. 64:27722777.[CrossRef]
Rouse, P. E. J. 1953. A theory of the linear viscoelastic properties of dilute solutions of coiling polymers. J. Chem. Phys. 21:12721280.[CrossRef]
Sartori, A., V. Barbier, and J.-L. Viovy. 2003. Sieving mechanisms in polymeric matrices. Electrophoresis. 24:421440.[CrossRef][Medline]
Schwartz, D. C., and K. Michael. 1989. Conformational dynamics of individual DNA molecules during gel electrophoresis. Nature. 338:520522.[CrossRef][Medline]
Smith, S. B., P. K. Aldridge, and J. B. Callis. 1989. Observation of individual DNA molecules undergoing gel electrophoresis. Science. 243:203243.
Stigter, D. 2000. Influence of agarose gel on electrophoretic stretch, on trapping, and on relaxation of DNA. Macromolecules. 33:88788889.[CrossRef]
Stigter, D., and C. Bustamante. 1998. Theory for the hydrodynamic and electrophoretic stretch of tethered B-DNA. Biophys. J. 75:11971210.
Sunada, W. M., and H. W. Blanch. 1998a. Microscopy of DNA in dilute polymer solutions. Biotechnol. Prog. 14:766772.[CrossRef][Medline]
Sunada, W. M., and H. W. Blanch. 1998b. A theory for the electrophoretic separation of DNA in polymer solutions. Electrophoresis. 19:31283136.[CrossRef][Medline]
Viovy, J.-L. 2000. Electrophoresis of DNA and other polyelectrolytes: physical mechanisms. Rev. Mod. Phys. 72:813872.[CrossRef]
Weill, G., and J. Des Cloizeaux.1979. Dynamics of polymers in dilute solutions: an explanation of anomalous indices by cross-over effects. J. Phys. (Paris). 40:99105.
Zimm, B. 1956. Dynamics of polymer molecules in dilute solution: viscoelasticity, flow birefringence and dielectric loss. J. Chem. Phys. 24:269278.[CrossRef]
Zimm, B. 1998. Extension in flow of a DNA molecule tethered at one end. Macromolecules. 31:60896098.[CrossRef]
Zimmerman, R. M., and E. C. Cox. 1994. DNA stretching on functionalized gold surfaces. Nucleic Acids Res. 22:492497.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |