Unzipping DNA with optical tweezers: high sequence sensitivity and force flips

Force measurements are performed on single DNA molecules with an optical trapping interferometer that combines subpiconewton force resolution and millisecond time resolution. A molecular construction is prepared for mechanically unzipping several thousand-basepair DNA sequences in an in vitro configuration. The force signals corresponding to opening and closing the double helix at low velocity are studied experimentally and are compared to calculations assuming thermal equilibrium. We address the effect of the stiffness on the basepair sensitivity and consider fluctuations in the force signal. With respect to earlier work performed with soft microneedles, we obtain a very significant increase in basepair sensitivity: presently, sequence features appearing at a scale of 10 basepairs are observed. When measured with the optical trap the unzipping force exhibits characteristic flips between different values at specific positions that are determined by the base sequence. This behavior is attributed to bistabilities in the position of the opening fork; the force flips directly reflect transitions between different states involved in the time-averaging of the molecular system.     

Probing protein-DNA interactions by unzipping a single DNA double helix

We present unzipping force analysis of protein association (UFAPA) as a novel and versatile method for detection of the position and dynamic nature of protein-DNA interactions. A single DNA double helix was unzipped in the presence of DNA-binding proteins using a feedback-enhanced optical trap. When the unzipping fork in a DNA reached a bound protein molecule we observed a dramatic increase in the tension in the DNA, followed by a sudden tension reduction. Analysis of the unzipping force throughout an unbinding "event" revealed information about the spatial location and dynamic nature of the protein-DNA complex. The capacity of UFAPA to spatially locate protein-DNA interactions is demonstrated by noncatalytic restriction mapping on a 4-kb DNA with three restriction enzymes (BsoBI, XhoI, and EcoRI). A restriction map for a given restriction enzyme was generated with an accuracy of approximately 25 bp. UFAPA also allows direct determination of the site-specific equilibrium association constant (K(A)) for a DNA-binding protein. This capability is demonstrated by measuring the cation concentration dependence of K(A) for EcoRI binding. The measured values are in good agreement with previous measurements of K(A) over an intermediate range of cation concentration. These results demonstrate the potential utility of UFAPA for future studies of site-specific protein-DNA interactions.     

Free-Energy Landscape and Characteristic Forces for the Initiation of DNA Unzipping

DNA unzipping, the separation of its double helix into single strands, is crucial in modulating a host of genetic processes. Although the large-scale separation of double-stranded DNA has been studied with a variety of theoretical and experimental techniques, the minute details of the very first steps of unzipping are still unclear. Here, we use atomistic molecular-dynamics simulations, coarse-grained simulations, and a statistical-mechanical model to study the initiation of DNA unzipping by an external force. Calculation of the potential of mean force profiles for the initial separation of the first few terminal basepairs in a DNA oligomer revealed that forces ranging between 130 and 230 pN are needed to disrupt the first basepair, and these values are an order of magnitude larger than those needed to disrupt basepairs in partially unzipped DNA. The force peak has an echo of ∼50 pN at the distance that unzips the second basepair. We show that the high peak needed to initiate unzipping derives from a free-energy basin that is distinct from the basins of subsequent basepairs because of entropic contributions, and we highlight the microscopic origin of the peak. To our knowledge, our results suggest a new window of exploration for single-molecule experiments.     

Nanopore unzipping of individual DNA hairpin molecules

We have used the nanometer scale alpha-Hemolysin pore to study the unzipping kinetics of individual DNA hairpins under constant force or constant loading rate. Using a dynamic voltage control method, the entry rate of polynucleotides into the pore and the voltage pattern applied to induce hairpin unzipping are independently set. Thus, hundreds of unzipping events can be tested in a short period of time (few minutes), independently of the unzipping voltage amplitude. Because our method does not entail the physical coupling of the molecules under test to a force transducer, very high throughput can be achieved. We used our method to study DNA unzipping kinetics at small forces, which have not been accessed before. We find that in this regime the static unzipping times decrease exponentially with voltage with a characteristic slope that is independent of the duplex region sequence, and that the intercept depends strongly on the duplex region energy. We also present the first nanopore dynamic force measurements (time varying force). Our results are in agreement with the approximately logV dependence at high V (where V is the loading rate) observed by other methods. The extension of these measurements to lower loading rates reveals a much weaker dependence on V.     

Unzipping DNA from the condensed globule state-effects of unraveling

We have studied theoretically the unzipping of a double-stranded DNA from a condensed globule state by an external force. At constant force, we found that the double-stranded DNA unzips an at critical force Fc and the number of unzipped monomers M goes as M approximately (Fc - F)-3, for both the homogeneous and heterogeneous double-stranded DNA sequence. This is different from the case of unzipping from an extended coil state in which the number of unzipped monomers M goes as M approximately (Fc - F)-chi, where the exponent chi is either 1 or 2 depending on whether the double-stranded DNA sequence is homogeneous or heterogeneous, respectively. In the case of unzipping at constant extension, we found that for a double-stranded DNA with a very large number N of base pairs, the force remains almost constant as a function of the extension, before the unraveling transition, at which the force drops abruptly to zero. Right at the unraveling transition, the number of base pairs remaining in the condensed globule state is still very large and goes as N(3/4), in agreement with theoretical predictions of the unraveling transition of polymers stretched by an external force.     

Mechanical stability of single DNA molecules

Using a modified atomic force microscope (AFM), individual double-stranded (ds) DNA molecules attached to an AFM tip and a gold surface were overstretched, and the mechanical stability of the DNA double helix was investigated. In lambda-phage DNA the previously reported B-S transition at 65 piconewtons (pN) is followed by a second conformational transition, during which the DNA double helix melts into two single strands. Unlike the B-S transition, the melting transition exhibits a pronounced force-loading-rate dependence and a marked hysteresis, characteristic of a nonequilibrium conformational transition. The kinetics of force-induced melting of the double helix, its reannealing kinetics, as well as the influence of ionic strength, temperature, and DNA sequence on the mechanical stability of the double helix were investigated. As expected, the DNA double helix is considerably destabilized under low salt buffer conditions (相似文献   

Sequence-dependent mechanics of single DNA molecules

Atomic force microscope-based single-molecule force spectroscopy was employed to measure sequence-dependent mechanical properties of DNA by stretching individual DNA double strands attached between a gold surface and an AFM tip. We discovered that in lambda-phage DNA the previously reported B-S transition, where 'S' represents an overstretched conformation, at 65 pN is followed by a nonequilibrium melting transition at 150 pN. During this transition the DNA is split into single strands that fully recombine upon relaxation. The sequence dependence was investigated in comparative studies with poly(dG-dC) and poly(dA-dT) DNA. Both the B-S and the melting transition occur at significantly lower forces in poly(dA-dT) compared to poly(dG-dC). We made use of the melting transition to prepare single poly(dG-dC) and poly(dA-dT) DNA strands that upon relaxation reannealed into hairpins as a result of their self-complementary sequence. The unzipping of these hairpins directly revealed the base pair-unbinding forces for G-C to be 20 +/- 3 pN and for A-T to be 9 +/- 3 pN.     

Equilibrium, kinetic, and footprinting studies of the Tus-Ter protein-DNA interaction.

Arrest of DNA replication in the terminus region of the Escherichia coli chromosome is mediated by protein-DNA complexes composed of the Tus protein and 23 base pair sequences generically called Ter sites. We have characterized the in vitro binding of purified Tus protein to a 37-base pair oligodeoxyribonucleotide containing the TerB sequence. The measured equilibrium binding constant (KD) for the chromosomal TerB site in KG buffer (50 mM Tris-Cl, 150 mM potassium glutamate, 25 degrees C, pH 7.5, 0.1 mM dithiothreitol, 0.1 mM EDTA, and 100 micrograms/ml bovine serum albumin) was 3.4 x 10(-13) M. Kinetic measurements in the same buffer revealed that the Tus-TerB complex was very stable, with a half-life of 550 min, a dissociation rate constant of 2.1 x 10(-5) s-1, and an association rate constant of 1.4 x 10(8) M-1 s-1. Similar measurements of Tus protein binding to the TerR2 site of the plasmid R6K showed an affinity 30-fold lower than the Tus-TerB interaction. This difference was due primarily to a more rapid dissociation of the Tus-TerR2 complex. Using standard chemical modification techniques, we also examined the DNA-protein contacts of the Tus-TerB interaction. Extensive contacts between the Tus protein and the TerB sequence were observed in the highly conserved 11 base-pair "core" sequence common to all identified Ter sites. In addition, protein-DNA contact sites were observed in the region of the Ter site where DNA replication is arrested. Projection of the footprinting data onto B-form DNA indicated that the majority of the alkylation interference and hydroxyl radical-protected sites were arranged on one face of the DNA helix. We also observed dimethyl sulfate protection of 2 guanine residues on the opposite side of the helix, suggesting that part of the Tus protein extends around the double helix. The distribution of contacts along the TerB sequence was consistent with the functional polarity of the Tus-Ter complex and suggested possible mechanisms for the impediment of protein translocation along DNA.     

Computer graphics program to reveal the dependence of the gross three-dimensional structure of the B-DNA double helix on primary structure.

Programs are presented to plot the gross three-dimensional structure of the DNA double helix with the base sequence as input information. The rules that determine the overall structure of the double helix are those that predict the dependence of local helix parameters (specifically, helix twist angle and relative basepair roll angle) on sequence. For this purpose, the user can select either the Calladine-Dickerson parameters or the Tung-Harvey parameters. These programs can be used as tools to investigate the variation of DNA tertiary structure with sequence, which may play an important role in the sequence-specific recognition of DNA by proteins.     

Long lifetime of hydrogen-bonded DNA basepairs by force spectroscopy

Electron-tunneling data suggest that a noncovalently-bonded complex of three molecules, two recognition molecules that present hydrogen-bond donor and acceptor sites via a carboxamide group, and a DNA base, remains bound for seconds. This is surprising, given that imino-proton exchange rates show that basepairs in a DNA double helix open on millisecond timescales. The long lifetime of the three-molecule complex was confirmed using force spectroscopy, but measurements on DNA basepairs are required to establish a comparison with the proton-exchange data. Here, we report on a dynamic force spectroscopy study of complexes between the bases adenine and thymine (A-T, two-hydrogen bonds) and 2-aminoadenine and thymine (2AA-T, three-hydrogen bonds). Bases were tethered to an AFM probe and mica substrate via long, covalently linked polymer tethers. Data for bond-survival probability versus force and the rupture-force distributions were well fitted by the Bell model. The resulting lifetime of the complexes at zero pulling force was ~2 s for two-hydrogen bonds (A-T) and ~4 s for three-hydrogen bonds (2AA-T). Thus, DNA basepairs in an AFM pulling experiment remain bonded for long times, even without the stabilizing influence of base-stacking in a double helix. This result suggests that the pathways for opening, and perhaps the open states themselves, are very different in the AFM and proton-exchange measurements.     

Microscopic Mechanics of Hairpin DNA Translocation through Synthetic Nanopores

Nanoscale pores have proved useful as a means to assay DNA and are actively being developed as the basis of genome sequencing methods. Hairpin DNA (hpDNA), having both double-helical and overhanging coil portions, can be trapped in a nanopore, giving ample time to execute a sequence measurement. In this article, we provide a detailed account of hpDNA interaction with a synthetic nanopore obtained through extensive all-atom molecular dynamics simulations. For synthetic pores with minimum diameters from 1.3 to 2.2 nm, we find that hpDNA can translocate by three modes: unzipping of the double helix and—in two distinct orientations—stretching/distortion of the double helix. Furthermore, each of these modes can be selected by an appropriate choice of the pore size and voltage applied transverse to the membrane. We demonstrate that the presence of hpDNA can dramatically alter the distribution of ions within the pore, substantially affecting the ionic current through it. In experiments and simulations, the ionic current relative to that in the absence of DNA can drop below 10% and rise beyond 200%. Simulations associate the former with the double helix occupying the constriction and the latter with accumulation of DNA that has passed through the constriction.     

Protein Unfolding by Biological Unfoldases: Insights from Modeling

The molecular determinants of the high efficiency of biological machines like unfoldases (e.g., the proteasome) are not well understood. We propose a model to study protein translocation into the chamber of biological unfoldases represented as a funnel. It is argued that translocation is a much faster way of unfolding a protein than end-to-end stretching, especially in a low-force regime, because it allows for a conformational freedom while concentrating local tension on consecutive regions of a protein chain and preventing refolding. This results in a serial unfolding of the protein structures dominated by unzipping. Thus, pulling against the unfoldase pore is an efficient catalyst of the unfolding reaction. We also show that the presence of the funnel makes the tension along the backbone of the substrate protein nonuniform even when the protein gets unfolded. Hence, the stalling force measured by single-molecule force spectroscopy techniques may be smaller than the traction force of the unfoldase motor.     

Protein Unfolding by Biological Unfoldases: Insights from Modeling

The molecular determinants of the high efficiency of biological machines like unfoldases (e.g., the proteasome) are not well understood. We propose a model to study protein translocation into the chamber of biological unfoldases represented as a funnel. It is argued that translocation is a much faster way of unfolding a protein than end-to-end stretching, especially in a low-force regime, because it allows for a conformational freedom while concentrating local tension on consecutive regions of a protein chain and preventing refolding. This results in a serial unfolding of the protein structures dominated by unzipping. Thus, pulling against the unfoldase pore is an efficient catalyst of the unfolding reaction. We also show that the presence of the funnel makes the tension along the backbone of the substrate protein nonuniform even when the protein gets unfolded. Hence, the stalling force measured by single-molecule force spectroscopy techniques may be smaller than the traction force of the unfoldase motor.     

The double helix is dehydrated: evidence from the hydrolysis of acridinium ester-labeled probes

A highly chemiluminescent reporter molecule, acridinium ester (AE), was tethered to single-stranded oligonucleotide probes and hybridized to complementary as well as mismatched target sequences. When tethered to single-stranded probes, AE was readily hydrolyzed by water or hydroxide ion. In contrast, when hybridized to a complementary target, hydrolysis of the AE probe was markedly inhibited. Mismatches near AE eliminated the ability of the double helix to strongly inhibit AE hydrolysis. To establish the molecular basis for these remarkable hydrolysis properties of AE-labeled probes, the binding and hydrolysis mechanisms of AE-labeled probes were examined. When tethered to single- or double-stranded nucleic acids, hydrolysis of AE was found to proceed by generalized base catalysis in which a base abstracts a proton from water and the resulting hydroxide ion then hydrolyzes AE. Analysis of the hydrolysis rates of AE bound to DNA revealed that AE binds the minor groove of DNA and that its hydrolysis is inhibited by low water activity within the minor groove of the helix. Depending upon the sequence of the DNA, the water activity of the minor groove was estimated to be at least 2-4-fold lower than bulk solution. Hydrolysis measurements of AE tethered to RNA as well as RNA/DNA hybrids argued that the grooves of these double helices are also dehydrated relative to bulk solution. Remarkably, mismatched bases, regardless of their structure or sequence context, enhanced hydrolysis of AE by inducing hydration of the double helix that spread approximately five base pairs on either side of the mismatch.     

DNA structural variations in the E. coli tyrT promoter

X-ray studies have established that the structure of a right-handed, Watson-Crick double helix can change from place to place along its length as a function of base sequence. The base pairs transmit deformations out to the phosphate backbone, where they can then be recognized by proteins and other DNA-binding reagents. Here we have examined at single-bond resolution the interactions of three commonly used nucleases (DNAase I, DNAase II, and copper-phenanthroline) with a DNA of natural origin, the 160 bp tyrT promoter. All three of these reagents seem sensitive to DNA backbone geometry rather than base sequence per se. Their sequence-dependent patterns of cleavage provide evidence for structural polymorphism of several sorts: global variation in helix groove width, global variation in radial asymmetry, and local variation in phosphate accessibility. These findings explain how sequence zones of a certain base composition, or purine-pyrimidine asymmetry, can influence the recognition of DNA by protein molecules.     

Signatures of DNA flexibility, interactions and sequence-related structural variations in classical X-ray diffraction patterns

The theory of X-ray diffraction from ideal, rigid helices allowed Watson and Crick to unravel the DNA structure, thereby elucidating functions encoded in it. Yet, as we know now, the DNA double helix is neither ideal nor rigid. Its structure varies with the base pair sequence. Its flexibility leads to thermal fluctuations and allows molecules to adapt their structure to optimize their intermolecular interactions. In addition to the double helix symmetry revealed by Watson and Crick, classical X-ray diffraction patterns of DNA contain information about the flexibility, interactions and sequence-related variations encoded within the helical structure. To extract this information, we have developed a new diffraction theory that accounts for these effects. We show how double helix non-ideality and fluctuations broaden the diffraction peaks. Meridional intensity profiles of the peaks at the first three helical layer lines reveal information about structural adaptation and intermolecular interactions. The meridional width of the fifth layer line peaks is inversely proportional to the helical coherence length that characterizes sequence-related and thermal variations in the double helix structure. Analysis of measured fiber diffraction patterns based on this theory yields important parameters that control DNA structure, packing and function.     

Torsional tension in the DNA double helix measured with trimethylpsoralen in living E. coli cells: analogous measurements in insect and human cells

The rate of covalent photobinding of trimethylpsoralen to DNA is greater when the DNA is wound with negative superhelical tension than when it is relaxed. In vitro the rate of photobinding is directly proportional to the negative superhelical density of the DNA. Thus measurement of the rate of photobinding provides an assay for probing in vivo unrestrained tension in the winding of the DNA double helix. This approach has been applied to measure torsional tension in DNA as it is packaged in living E. coli. Drosophila and HeLa cells. A method is described for measuring the rate of photobinding to intracellular DNA and rRNA, and for using the latter measurement as an internal control of the rate of me3-psoralen photobinding in vivo. This permits more accurate and reproducible measurement of changes in the DNA-psoralen photobinding reaction. The me3-psoralen probe interacts with intracellular bacterial DNA as expected for a purified DNA duplex wound with superhelical density sigma = -0.05 +/- 0.01. This superhelical tension is relaxed in cells when multiple single-strand breaks are introduced into the chromosomal DNA by gamma-irradiation. Similar relaxation occurs when cells are treated with the DNA gyrase inhibitor coumermycin. The results suggest that the DNA double helix is wound with torsional tension in vivo and that DNA supercoils which are equilibrated with this tension are not completely restrained in nucleosome-like structures. Torsional tension in the DNA of eucaryotic cells is not detectable in analogous measurements of the packaged DNA of HeLa and Drosophila cells. The simplest interpretation of this finding is that, within the limits of detection, all superhelical turns in the DNA are restrained in nucleosomes or nucleosome-like structures in these eucaryotic cells.     

Equal Sensitivity of the Two Strands of øX174 Replicative DNA to Breakage by Ionizing Radiation

The supertwisted, double-stranded, replicative intermediate of øX174 DNA (RFI) has been used to determine whether one of the two strands of the double helix is uniquely sensitive to induction of single-strand breaks by ionizing radiation. This could result from a particularly sensitive base sequence or a transfer of energy to a specific location of the DNA molecule. The results indicate that both strands of the double helix are equally broken, even though their base compositions are significantly different. If there are “hot spots” in the strands, then they are present in equal amounts in each strand.