DNA overstretching in the presence of glyoxal: structural evidence of force-induced DNA melting

Cytosolic Interactome Protects Against Protein Unfolding in a Single Molecule Experiment

Single molecule techniques are particularly well suited for investigating the processes of protein folding and chaperone assistance. However, current assays provide only a limited perspective on the various ways in which the cellular environment can influence the folding pathway of a protein. In this study, a single molecule mechanical interrogation assay is developed and used to monitor protein unfolding and refolding within a cytosolic solution. This allows to test the cumulative topological effect of the cytoplasmic interactome on the folding process. The results reveal a stabilization against forced unfolding for partial folds, which are attributed to the protective effect of the cytoplasmic environment against unfolding and aggregation. This research opens the possibility of conducting single molecule molecular folding experiments in quasi-biological environments.

Structure and Dynamics of dsDNA in Cell-like Environments

Deoxyribonucleic acid (DNA) is a fundamental biomolecule for correct cellular functioning and regulation of biological processes. DNA's structure is dynamic and has the ability to adopt a variety of structural conformations in addition to its most widely known double-stranded DNA (dsDNA) helix structure. Stability and structural dynamics of dsDNA play an important role in molecular biology. In vivo, DNA molecules are folded in a tightly confined space, such as a cell chamber or a channel, and are highly dense in solution; their conformational properties are restricted, which affects their thermodynamics and mechanical properties. There are also many technical medical purposes for which DNA is placed in a confined space, such as gene therapy, DNA encapsulation, DNA mapping, etc. Physiological conditions and the nature of confined spaces have a significant influence on the opening or denaturation of DNA base pairs. In this review, we summarize the progress of research on the stability and dynamics of dsDNA in cell-like environments and discuss current challenges and future directions. We include studies on various thermal and mechanical properties of dsDNA in ionic solutions, molecular crowded environments, and confined spaces. By providing a better understanding of melting and unzipping of dsDNA in different environments, this review provides valuable guidelines for predicting DNA thermodynamic quantities and for designing DNA/RNA nanostructures.

Mechanical response to tension and torque of molecular chains via statistically interacting particles associated with extension, contraction, twist, and supercoiling

A methodology for the statistical mechanical analysis of polymeric chains under tension introduced previously is extended to include torque. The response of individual bonds between monomers or of entire groups of monomers to a combination of tension and torque involves, in the framework of this method of analysis, the (thermal or mechanical) activation of a specific mix of statistically interacting particles carrying quanta of extension or contraction and quanta of twist or supercoiling. The methodology, which is elucidated in applications of increasing complexity, is capable of describing the conversion between twist chirality and plectonemic chirality in quasistatic processes. The control variables are force or extension and torque or linkage (a combination of twist and writhe). The versatility of this approach is demonstrated in two applications relevant and promising for double-stranded DNA under controlled tension and torque. One application describes conformational transformations between (native) B-DNA, (underwound) S-DNA, and (overwound) P-DNA in accord with experimental data. The other application describes how the conversion between a twisted chain and a supercoiled chain accommodates variations of linkage and excess length in a buckling transition.

Constructing Free Energy Landscapes of Nucleic Acid Hairpin Unfolding

Single nucleic acid molecules form hairpins that may stabilize secondary and tertiary structures as well as perform enzymatic and other chemical functions. Considerable progress has been made in the effort to understand the contributions of various factors to the stability of a given hairpin sequence. For a given sequence, it is possible to compute both the most likely structural arrangements and their associated free energies over a range of experimental conditions. However, there are many observed hairpin irregularities for which the energies and function are not well understood. Here we examine the irregular RNA Transactivation Response (TAR) hairpin from the HIV-1 genome. Using single molecule optical tweezers, the hairpin is force unfolded, revealing the overall unfolding free energy and the character of the transition state. These measurements allow the construction of a simple energy landscape from unfolding measurements, which can be directly compared to a theoretical landscape. This method is easily adapted to other structures, including the effects of noncanonical bases and even ligand binding.

The role of the C-domain of bacteriophage T4 gene 32 protein in ssDNA binding and dsDNA helix-destabilization: Kinetic, single-molecule, and cross-linking studies

The model single-stranded DNA binding protein of bacteriophage T4, gene 32 protein (gp32) has well-established roles in DNA replication, recombination, and repair. gp32 is a single-chain polypeptide consisting of three domains. Based on thermodynamics and kinetics measurements, we have proposed that gp32 can undergo a conformational change where the acidic C-terminal domain binds internally to or near the single-stranded (ss) DNA binding surface in the core (central) domain, blocking ssDNA interaction. To test this model, we have employed a variety of experimental approaches and gp32 variants to characterize this conformational change. Utilizing stopped-flow methods, the association kinetics of wild type and truncated forms of gp32 with ssDNA were measured. When the C-domain is present, the log-log plot of k vs. [NaCl] shows a positive slope, whereas when it is absent (*I protein), there is little rate change with salt concentration, as expected for this model.A gp32 variant lacking residues 292-296 within the C-domain, ΔPR201, displays kinetic properties intermediate between gp32 and *I. The single molecule force-induced DNA helix-destabilizing activitiesas well as the single- and double-stranded DNA affinities of ΔPR201 and gp32 truncated at residue 295 also fall between full-length protein and *I. Finally, chemical cross-linking of recombinant C-domain and gp32 lacking both N- and C-terminal domains is inhibited by increasing concentrations of a short single-stranded oligonucleotide, and the salt dependence of cross-linking mirrors that expected for the model. Taken together, these results provide the first evidence in support of this model that have been obtained through structural probes.

Single-molecule mechanochemical characterization of E. coli pol III core catalytic activity

Pol III core is the three‐subunit subassembly of the E. coli replicative DNA polymerase III holoenzyme. It contains the catalytic polymerase subunit α, the 3′ → 5′ proofreading exonuclease ε, and a subunit of unknown function, θ. We employ optical tweezers to characterize pol III core activity on a single DNA substrate. We observe polymerization at applied template forces F < 25 pN and exonucleolysis at F > 30 pN. Both polymerization and exonucleolysis occur as a series of short bursts separated by pauses. For polymerization, the initiation rate after pausing is independent of force. In contrast, the exonucleolysis initiation rate depends strongly on force. The measured force and concentration dependence of exonucleolysis initiation fits well to a two‐step reaction scheme in which pol III core binds bimolecularly to the primer‐template junction, then converts at rate k₂ into an exo‐competent conformation. Fits to the force dependence of k_init show that exo initiation requires fluctuational opening of two base pairs, in agreement with temperature‐ and mismatch‐dependent bulk biochemical assays. Taken together, our results support a model in which the pol and exo activities of pol III core are effectively independent, and in which recognition of the 3′ end of the primer by either α or ε is governed by the primer stability. Thus, binding to an unstable primer is the primary mechanism for mismatch recognition during proofreading, rather than an alternative model of duplex defect recognition.

Single-molecule studies of high-mobility group B architectural DNA bending proteins

Protein–DNA interactions can be characterized and quantified using single molecule methods such as optical tweezers, magnetic tweezers, atomic force microscopy, and fluorescence imaging. In this review, we discuss studies that characterize the binding of high-mobility group B (HMGB) architectural proteins to single DNA molecules. We show how these studies are able to extract quantitative information regarding equilibrium binding as well as non-equilibrium binding kinetics. HMGB proteins play critical but poorly understood roles in cellular function. These roles vary from the maintenance of chromatin structure and facilitation of ribosomal RNA transcription (yeast high-mobility group 1 protein) to regulatory and packaging roles (human mitochondrial transcription factor A). We describe how these HMGB proteins bind, bend, bridge, loop and compact DNA to perform these functions. We also describe how single molecule experiments observe multiple rates for dissociation of HMGB proteins from DNA, while only one rate is observed in bulk experiments. The measured single-molecule kinetics reveals a local, microscopic mechanism by which HMGB proteins alter DNA flexibility, along with a second, much slower macroscopic rate that describes the complete dissociation of the protein from DNA.

Mechanisms of small molecule-DNA interactions probed by single-molecule force spectroscopy

There is a wide range of applications for non-covalent DNA binding ligands, and optimization of such interactions requires detailed understanding of the binding mechanisms. One important class of these ligands is that of intercalators, which bind DNA by inserting aromatic moieties between adjacent DNA base pairs. Characterizing the dynamic and equilibrium aspects of DNA-intercalator complex assembly may allow optimization of DNA binding for specific functions. Single-molecule force spectroscopy studies have recently revealed new details about the molecular mechanisms governing DNA intercalation. These studies can provide the binding kinetics and affinity as well as determining the magnitude of the double helix structural deformations during the dynamic assembly of DNA–ligand complexes. These results may in turn guide the rational design of intercalators synthesized for DNA-targeted drugs, optical probes, or integrated biological self-assembly processes. Herein, we survey the progress in experimental methods as well as the corresponding analysis framework for understanding single molecule DNA binding mechanisms. We discuss briefly minor and major groove binding ligands, and then focus on intercalators, which have been probed extensively with these methods. Conventional mono-intercalators and bis-intercalators are discussed, followed by unconventional DNA intercalation. We then consider the prospects for using these methods in optimizing conventional and unconventional DNA-intercalating small molecules.

Dynamic Release of Bending Stress in Short dsDNA by Formation of a Kink and Forks

Bending with high curvature is one of the major mechanical properties of double-stranded DNA (dsDNA) that is essential for its biological functions. The emergence of a kink arising from local melting in the middle of dsDNA has been suggested as a mechanism of releasing the energy cost of bending. Herein, we report that strong bending induces two types of short dsDNA deformations, induced by two types of local melting, namely, a kink in the middle and forks at the ends, which we demonstrate using D-shaped DNA nanostructures. The two types of deformed dsDNA structures dynamically interconvert on a millisecond timescale. The transition from a fork to a kink is dominated by entropic contribution (anti-Arrhenius behavior), while the transition from a kink to a fork is dominated by enthalpic contributions. The presence of mismatches in dsDNA accelerates kink formation, and the transition from a kink to a fork is removed when the mismatch size is three base pairs.

Single aromatic residue location alters nucleic acid binding and chaperone function of FIV nucleocapsid protein

Feline immunodeficiency virus (FIV) is a retrovirus that infects domestic cats, and is an excellent animal model for human immunodeficiency virus type 1 (HIV-1) pathogenesis. The nucleocapsid (NC) protein is critical for replication in both retroviruses. FIV NC has several structural features that differ from HIV-1 NC. While both NC proteins have a single conserved aromatic residue in each of the two zinc fingers, the aromatic residue on the second finger of FIV NC is located on the opposite C-terminal side relative to its location in HIV-1 NC. In addition, whereas HIV-1 NC has a highly charged cationic N-terminal tail and a relatively short C-terminal extension, the opposite is true for FIV NC. To probe the impact of these differences on the nucleic acid (NA) binding and chaperone properties of FIV NC, we carried out ensemble and single-molecule assays with wild-type (WT) and mutant proteins. The ensemble studies show that FIV NC binding to DNA is strongly electrostatic, with a higher effective charge than that observed for HIV-1 NC. The C-terminal basic domain contributes significantly to the NA binding capability of FIV NC. In addition, the non-electrostatic component of DNA binding is much weaker for FIV NC than for HIV-1 NC. Mutation of both aromatic residues in the zinc fingers to Ala (F12A/W44A) further increases the effective charge of FIV NC and reduces its non-electrostatic binding affinity. Interestingly, switching the location of the C-terminal aromatic residue to mimic the HIV-1 NC sequence (N31W/W44A) reduces the effective charge of FIV NC and increases its non-electrostatic binding affinity to values similar to HIV-1 NC. Consistent with the results of these ensemble studies, single-molecule DNA stretching studies show that while WT FIV NC has reduced stacking capability relative to HIV-1 NC, the aromatic switch mutant recovers the ability to intercalate between the DNA bases. Our results demonstrate that altering the position of a single aromatic residue switches the binding mode of FIV NC from primarily electrostatic binding to more non-electrostatic binding, conferring upon it NA interaction properties comparable to that of HIV-1 NC.

DNA stretching on the wall surfaces in curved microchannels with different radii

DNA molecule conformation dynamics and stretching were made on semi-circular surfaces with different radii (500 to 5,000 μm) in microchannels measuring 200 μm × 200 μm in cross section. Five different buffer solutions - 1× Tris-acetate-EDTA (TAE), 1× Tris-borate-EDTA (TBE), 1× Tris-EDTA (TE), 1× Tris-phosphate-EDTA (TPE), and 1× Tris-buffered saline (TBS) solutions - were used with a variety of viscosity such as 40, 60, and 80 cP, with resultant 10(-4) ≤ Re ≤ 10(-3) and the corresponding 5 ≤ Wi ≤ 12. The test fluids were seeded with JOJO-1 tracer particles for flow visualization and driven through the test channels via a piezoelectric (PZT) micropump. Micro particle image velocimetry (μPIV) measuring technique was applied for the centered-plane velocity distribution measurements. It is found that the radius effect on the stretch ratio of DNA dependence is significant. The stretch ratio becomes larger as the radius becomes small due to the larger centrifugal force. Consequently, the maximum stretch was found at the center of the channel with a radius of 500 μm.

Force-induced melting of DNA--evidence for peeling and internal melting from force spectra on short synthetic duplex sequences

Overstretching of DNA occurs at about 60-70 pN when a torsionally unconstrained double-stranded DNA molecule is stretched by its ends. During the transition, the contour length increases by up to 70% without complete strand dissociation. Three mechanisms are thought to be involved: force-induced melting into single-stranded DNA where either one or both strands carry the tension, or a B-to-S transition into a longer, still base-paired conformation. We stretch sequence-designed oligonucleotides in an effort to isolate the three processes, focusing on force-induced melting. By introducing site-specific inter-strand cross-links in one or both ends of a 64 bp AT-rich duplex we could repeatedly follow the two melting processes at 5 mM and 1 M monovalent salt. We find that when one end is sealed the AT-rich sequence undergoes peeling exhibiting hysteresis at low and high salt. When both ends are sealed the AT sequence instead undergoes internal melting. Thirdly, the peeling melting is studied in a composite oligonucleotide where the same AT-rich sequence is concatenated to a GC-rich sequence known to undergo a B-to-S transition rather than melting. The construct then first melts in the AT-rich part followed at higher forces by a B-to-S transition in the GC-part, indicating that DNA overstretching modes are additive.

Transient kinetics measured with force steps discriminate between double-stranded DNA elongation and melting and define the reaction energetics

Under a tension of ∼65 pN, double-stranded DNA undergoes an overstretching transition from its basic (B-form) conformation to a 1.7 times longer conformation whose nature is only recently starting to be understood. Here we provide a structural and thermodynamic characterization of the transition by recording the length transient following force steps imposed on the λ-phage DNA with different melting degrees and temperatures (10–25°C). The shortening transient following a 20–35 pN force drop from the overstretching force shows a sequence of fast shortenings of double-stranded extended (S-form) segments and pauses owing to reannealing of melted segments. The lengthening transients following a 2–35 pN stretch to the overstretching force show the kinetics of a two-state reaction and indicate that the whole 70% extension is a B-S transition that precedes and is independent of melting. The temperature dependence of the lengthening transient shows that the entropic contribution to the B-S transition is one-third of the entropy change of thermal melting, reinforcing the evidence for a double-stranded S-form that maintains a significant fraction of the interstrand bonds. The cooperativity of the unitary elongation (22 bp) is independent of temperature, suggesting that structural factors, such as the nucleic acid sequence, control the transition.

An effective mesoscopic model of double-stranded DNA

Watson and Crick's epochal presentation of the double helix structure in 1953 has paved the way to intense exploration of DNA's vital functions in cells. Also, recent advances of single molecule techniques have made it possible to probe structures and mechanics of constrained DNA at length scales ranging from nanometers to microns. There have been a number of atomistic scale quantum chemical calculations or molecular level simulations, but they are too computationally demanding or analytically unfeasible to describe the DNA conformation and mechanics at mesoscopic levels. At micron scales, on the other hand, the wormlike chain model has been very instrumental in describing analytically the DNA mechanics but lacks certain molecular details that are essential in describing the hybridization, nano-scale confinement, and local denaturation. To fill this fundamental gap, we present a workable and predictive mesoscopic model of double-stranded DNA where the nucleotides beads constitute the basic degrees of freedom. With the inter-strand stacking given by an interaction between diagonally opposed monomers, the model explains with analytical simplicity the helix formation and produces a generalized wormlike chain model with the concomitant large bending modulus given in terms of the helical structure and stiffness. It also explains how the helical conformation undergoes overstretch transition to the ladder-like conformation at a force plateau, in agreement with the experiment.

Global force-torque phase diagram for the DNA double helix: structural transitions, triple points, and collapsed plectonemes

We present a free energy model for structural transitions of the DNA double helix driven by tensile and torsional stress. Our model is coarse grained and is based on semiflexible polymer descriptions of B-DNA, underwound L-DNA, and highly overwound P-DNA. The statistical-mechanical model of plectonemic supercoiling previously developed for B-DNA is applied to semiflexible polymer models of P- and L-DNA to obtain a model of DNA structural transitions in quantitative accord with experiment. We identify two distinct plectonemic states, one "inflated" by electrostatic repulsion and thermal fluctuations and the other "collapsed," with the two double helices inside the supercoils driven to close contact. We find that supercoiled B and L are stable only in the inflated form, while supercoiled P is always collapsed. We also predict the behavior and experimental signatures of highly underwound "Q"-DNA, the left-handed analog of P-DNA; as for P, supercoiled Q is always collapsed. Overstretched "S"-DNA and strand-separated "stress-melted" DNA are also included in our model, allowing prediction of a global phase diagram for forces up to 1000 pN and torques between ±60 pN nm, or, in terms of linking number density, from σ=-5 to +3.

Differential contribution of basic residues to HIV-1 nucleocapsid protein's nucleic acid chaperone function and retroviral replication

The human immunodeficiency virus type 1 (HIV-1) nucleocapsid (NC) protein contains 15 basic residues located throughout its 55-amino acid sequence, as well as one aromatic residue in each of its two CCHC-type zinc finger motifs. NC facilitates nucleic acid (NA) rearrangements via its chaperone activity, but the structural basis for this activity and its consequences in vivo are not completely understood. Here, we investigate the role played by basic residues in the N-terminal domain, the N-terminal zinc finger and the linker region between the two zinc fingers. We use in vitro ensemble and single-molecule DNA stretching experiments to measure the characteristics of wild-type and mutant HIV-1 NC proteins, and correlate these results with cell-based HIV-1 replication assays. All of the cationic residue mutations lead to NA interaction defects, as well as reduced HIV-1 infectivity, and these effects are most pronounced on neutralizing all five N-terminal cationic residues. HIV-1 infectivity in cells is correlated most strongly with NC’s NA annealing capabilities as well as its ability to intercalate the DNA duplex. Although NC’s aromatic residues participate directly in DNA intercalation, our findings suggest that specific basic residues enhance these interactions, resulting in optimal NA chaperone activity.

Polymerase manager protein UmuD directly regulates Escherichia coli DNA polymerase III α binding to ssDNA

Replication by Escherichia coli DNA polymerase III is disrupted on encountering DNA damage. Consequently, specialized Y-family DNA polymerases are used to bypass DNA damage. The protein UmuD is extensively involved in modulating cellular responses to DNA damage and may play a role in DNA polymerase exchange for damage tolerance. In the absence of DNA, UmuD interacts with the α subunit of DNA polymerase III at two distinct binding sites, one of which is adjacent to the single-stranded DNA-binding site of α. Here, we use single molecule DNA stretching experiments to demonstrate that UmuD specifically inhibits binding of α to ssDNA. We predict using molecular modeling that UmuD residues D91 and G92 are involved in this interaction and demonstrate that mutation of these residues disrupts the interaction. Our results suggest that competition between UmuD and ssDNA for α binding is a new mechanism for polymerase exchange.

Denaturation transition of stretched DNA

In the last two decades, single-molecule force measurements using optical and magnetic tweezers and atomic force spectroscopy have dramatically expanded our knowledge of nucleic acids and proteins. These techniques characterize the force on a biomolecule required to produce a given molecular extension. When stretching long DNA molecules, the observed force–extension relationship exhibits a characteristic plateau at approximately 65 pN where the DNA may be extended to almost twice its B-DNA length with almost no increase in force. In the present review, I describe this transition in terms of the Poland–Scheraga model and summarize recent related studies.

Revealing the competition between peeled ssDNA, melting bubbles, and S-DNA during DNA overstretching using fluorescence microscopy

Mechanical stress plays a key role in many genomic processes, such as DNA replication and transcription. The ability to predict the response of double-stranded (ds) DNA to tension is a cornerstone of understanding DNA mechanics. It is widely appreciated that torsionally relaxed dsDNA exhibits a structural transition at forces of ∼65 pN, known as overstretching, whereby the contour length of the molecule increases by ∼70%. Despite extensive investigation, the structural changes occurring in DNA during overstretching are still generating considerable debate. Three mechanisms have been proposed to account for the increase in DNA contour length during overstretching: strand unpeeling, localized base-pair breaking (yielding melting bubbles), and formation of S-DNA (strand unwinding, while base pairing is maintained). Here we show, using a combination of fluorescence microscopy and optical tweezers, that all three structures can exist, uniting the often contradictory dogmas of DNA overstretching. We visualize and distinguish strand unpeeling and melting-bubble formation using an appropriate combination of fluorescently labeled proteins, whereas remaining B-form DNA is accounted for by using specific fluorescent molecular markers. Regions of S-DNA are associated with domains where fluorescent probes do not bind. We demonstrate that the balance between the three structures of overstretched DNA is governed by both DNA topology and local DNA stability. These findings enhance our knowledge of DNA mechanics and stability, which are of fundamental importance to understanding how proteins modify the physical state of DNA.

Aromatic residue mutations reveal direct correlation between HIV-1 nucleocapsid protein's nucleic acid chaperone activity and retroviral replication

The human immunodeficiency virus type 1 (HIV-1) nucleocapsid (NC) protein plays an essential role in several stages of HIV-1 replication. One important function of HIV-1 NC is to act as a nucleic acid chaperone, in which the protein facilitates nucleic acid rearrangements important for reverse transcription and recombination. NC contains only 55 amino acids, with 15 basic residues and two zinc fingers, each having a single aromatic residue (Phe16 and Trp37). Despite its simple structure, HIV-1 NC appears to have optimal chaperone activity, including the ability to strongly aggregate nucleic acids, destabilize nucleic acid secondary structure, and facilitate rapid nucleic acid annealing. Here we combine single molecule DNA stretching experiments with ensemble solution studies of protein-nucleic acid binding affinity, oligonucleotide annealing, and nucleic acid aggregation to measure the characteristics of wild-type (WT) and aromatic residue mutants of HIV-1 NC that are important for nucleic acid chaperone activity. These in vitro results are compared to in vivo HIV-1 replication for viruses containing the same mutations. This work allows us to directly relate HIV-1 NC structure with its function as a nucleic acid chaperone in vitro and in vivo. We show that replacement of either aromatic residue with another aromatic residue results in a protein that strongly resembles WT NC. In contrast, single amino acid substitutions of either Phe16Ala or Trp37Ala significantly slow down NC's DNA interaction kinetics, while retaining some helix-destabilization capability. A double Phe16Ala/Trp37Ala substitution further reduces the latter activity. Surprisingly, the ensemble nucleic acid binding, annealing, and aggregation properties are not significantly altered for any mutant except the double aromatic substitution with Ala. Thus, elimination of a single aromatic residue from either zinc finger strongly reduces NC's chaperone activity as determined by single molecule DNA stretching experiments without significantly altering its ensemble-averaged biochemical properties. Importantly, the substitution of aromatic residues with Ala progressively decreases NC's nucleic acid chaperone activity while also progressively inhibiting viral replication. Taken together, these data support the critical role of HIV-1 NC's aromatic residues, and establish a direct and statistically significant correlation between nucleic acid chaperone activity and viral replication.

Single-molecule kinetics reveal microscopic mechanism by which High-Mobility Group B proteins alter DNA flexibility

Eukaryotic High-Mobility Group B (HMGB) proteins alter DNA elasticity while facilitating transcription, replication and DNA repair. We developed a new single-molecule method to probe non-specific DNA interactions for two HMGB homologs: the human HMGB2 box A domain and yeast Nhp6Ap, along with chimeric mutants replacing neutral N-terminal residues of the HMGB2 protein with cationic sequences from Nhp6Ap. Surprisingly, HMGB proteins constrain DNA winding, and this torsional constraint is released over short timescales. These measurements reveal the microscopic dissociation rates of HMGB from DNA. Separate microscopic and macroscopic (or local and non-local) unbinding rates have been previously proposed, but never independently observed. Microscopic dissociation rates for the chimeric mutants (∼10 s⁻¹) are higher than those observed for wild-type proteins (∼0.1–1.0 s⁻¹), reflecting their reduced ability to bend DNA through short-range interactions, despite their increased DNA-binding affinity. Therefore, transient local HMGB–DNA contacts dominate the DNA-bending mechanism used by these important architectural proteins to increase DNA flexibility.

Tension induces a base-paired overstretched DNA conformation

Mixed-sequence DNA molecules undergo mechanical overstretching by approximately 70% at 60-70 pN. Since its initial discovery 15 y ago, a debate has arisen as to whether the molecule adopts a new form [Cluzel P, et al. (1996) Science 271:792-794; Smith SB, Cui Y, Bustamante C (1996) Science 271:795-799], or simply denatures under tension [van Mameren J, et al. (2009) Proc Natl Acad Sci USA 106:18231-18236]. Here, we resolve this controversy by using optical tweezers to extend small 60-64 bp single DNA duplex molecules whose base content can be designed at will. We show that when AT content is high (70%), a force-induced denaturation of the DNA helix ensues at 62 pN that is accompanied by an extension of the molecule of approximately 70%. By contrast, GC-rich sequences (60% GC) are found to undergo a reversible overstretching transition into a distinct form that is characterized by a 51% extension and that remains base-paired. For the first time, results proving the existence of a stretched basepaired form of DNA can be presented. The extension observed in the reversible transition coincides with that produced on DNA by binding of bacterial RecA and human Rad51, pointing to its possible relevance in homologous recombination.

Two distinct overstretched DNA structures revealed by single-molecule thermodynamics measurements

Double-stranded DNA is a dynamic molecule whose structure can change depending on conditions. While there is consensus in the literature about many structures DNA can have, the state of highly-stretched DNA is still not clear. Several groups have shown that DNA in the torsion-unconstrained B-form undergoes an "overstretching" transition at a stretching force of around 65 pN, which leads to approximately 1.7-fold elongation of the DNA contour length. Recent experiments have revealed that two distinct structural transitions are involved in the overstretching process: (i) a hysteretic "peeling" off one strand from its complementary strand, and (ii) a nonhysteretic transition that leads to an undetermined DNA structure. We report the first simultaneous determination of the entropy (ΔS) and enthalpy changes (ΔH) pertaining to these respective transitions. For the hysteretic peeling transition, we determined ΔS ∼ 20 cal/(K.mol) and ΔH ∼ 7 kcal/mol. In the case of the nonhysteretic transition, ΔS ∼ -3 cal/(K.mol) and ΔH ∼ 1 kcal/mol. Furthermore, the response of the transition force to salt concentration implies that the two DNA strands are spatially separated after the hysteretic peeling transition. In contrast, the corresponding response after the nonhysteretic transition indicated that the strands remained in close proximity. The selection between the two transitions depends on DNA base-pair stability, and it can be illustrated by a multidimensional phase diagram. Our results provide important insights into the thermodynamics of DNA overstretching and conformational structures of overstretched DNA that may play an important role in vivo.

A single zinc finger optimizes the DNA interactions of the nucleocapsid protein of the yeast retrotransposon Ty3

Reverse transcription in retroviruses and retrotransposons requires nucleic acid chaperones, which drive the rearrangement of nucleic acid conformation. The nucleic acid chaperone properties of the human immunodeficiency virus type-1 (HIV-1) nucleocapsid (NC) protein have been extensively studied, and nucleic acid aggregation, duplex destabilization and rapid binding kinetics have been identified as major components of its activity. However, the properties of other nucleic acid chaperone proteins, such as retrotransposon Ty3 NC, a likely ancestor of HIV-1 NC, are not well understood. In addition, it is unclear whether a single zinc finger is sufficient to optimize the properties characteristic of HIV-1 NC. We used single-molecule DNA stretching as a method for detailed characterization of Ty3 NC chaperone activity. We found that wild type Ty3 NC aggregates single- and double-stranded DNA, weakly stabilizes dsDNA, and exhibits rapid binding kinetics. Single-molecule studies in the presence of Ty3 NC mutants show that the N-terminal basic residues and the unique zinc finger at the C-terminus are required for optimum chaperone activity in this system. While the single zinc finger is capable of optimizing Ty3 NC's DNA interaction kinetics, two zinc fingers may be necessary in order to facilitate the DNA destabilization exhibited by HIV-1 NC.

Paired mutations abolish and restore the balanced annealing and melting activities of ORF1p that are required for LINE-1 retrotransposition

Retrotransposition amplifies LINE-1 (L1) to high copy number in mammalian genomes. The L1 protein encoded by ORF1 (ORF1p) is required for retrotransposition. This dependence on ORF1p was investigated by mutating three highly conserved residues, R238, R284 and Y318 to alanine, thereby inactivating retrotransposition. R284A and Y318A were rescued by further substituting the alanine with the appropriate conservative amino acid, e.g. lysine or phenylalanine, respectively, whereas R238K remained inactive. Quantification of the steady-state levels of L1 RNA and ORF1p failed to discriminate active from inactive variants, indicating loss of L1 retrotransposition resulted from loss of function rather than reduced expression. The two biochemical properties known for ORF1p are high-affinity RNA binding and nucleic acid chaperone activity. Only R238A/K exhibited significantly reduced RNA affinities. The nucleic acid chaperone activities of the remaining paired mutants were assessed by single-molecule DNA stretching and found to mirror retrotransposition activity. To further examine ORF1p chaperone function, their energetic barriers to DNA annealing and melting were derived from kinetic work. When plotted against each other, the ratio of these two activities distinguished functional from non-functional ORF1p variants. These findings enhance our understanding of the requirements for ORF1p in LINE-1 retrotransposition and, more generally, nucleic acid chaperone function.

Quantitative High-Resolution Sensing of DNA Hybridization Using Magnetic Tweezers with Evanescent Illumination

We applied the combined approach of evanescent nanometry and force spectroscopy using magnetic tweezers to quantify the degree of hybridization of a single synthetic single-stranded DNA oligomer to a resolution approaching a single-base. In this setup, the 200 nucleotide long DNA was covalently attached to the surface of an optically transparent solid support at one end and to the surface of a superparamagnetic fluorescent microsphere (force probe) at the other end. The force was applied to the probes using an electromagnet. The end-to-end molecular distance (i.e. out-of-image-plane position of the force probe) was determined from the intensity of the probe fluorescent image observed with total-internal reflectance microscopy. An equation of state for single stranded DNA molecules under tension (extensible freely jointed chain) was used to derive the penetration depth of the evanescent field and to calibrate the magnetic properties of the force probes. The parameters of the magnetic response of the force probes obtained from the equation of state remained constant when changing the penetration depth, indicating a robust calibration procedure. The results of such a calibration were also confirmed using independently measured probe-surface distances for probes mounted onto cantilevers of an atomic force microscope. Upon hybridization of the complementary 50 nucleotide-long oligomer to the surface-bound 200-mer, the changes in the force-distance curves were consistent with the quantitative conversion of 25% of the original single-stranded DNA to its double-stranded form, which was modeled as an elastic rod. The method presented here for quantifying the hybridization state of the single DNA molecules has potential for determining the degree of hybridization of individual molecules in a single molecule array with high accuracy.

Measuring DNA-protein binding affinity on a single molecule using optical tweezers

DNA-protein interactions may be observed on single molecules with a variety of techniques. However, quantifying the binding affinity is difficult and often requires many (∼100) individual events to characterize the interaction. We use a single λ DNA molecule that provides a lattice of binding sites for proteins. Extending and relaxing the tethered molecule reversibly melts DNA, allowing it to be converted between double-stranded (ds) and single-stranded (ss) forms. By monitoring changes in the properties of the DNA as a function of added protein concentration and fitting to binding models, the DNA-protein interaction may be characterized and quantified. As an example, the high mobility group protein HMGB1(box A  +  B) is observed to stabilize dsDNA. Measuring the strength of this effect allows us to determine the equilibrium association constant for HMGB1(box A  +  B) binding to dsDNA.

Transition dynamics and selection of the distinct S-DNA and strand unpeeling modes of double helix overstretching

Recent studies have revealed two distinct pathways for the DNA overstretching transition near 65 pN: ‘unpeeling’ of one strand from the other, and a transition from B-DNA to an elongated double-stranded ‘S-DNA’ form. However, basic questions concerning the dynamics of these transitions, relative stability of the two competing overstretched states, and effects of nicks and free DNA ends on overstretching, remain open. In this study we report that: (i) stepwise extension changes caused by sequence-defined barriers occur during the strand-unpeeling transition, whereas rapid, sequence-independent extension fluctuations occur during the B to S transition; (ii) the secondary transition that often occurs following the overstretching transition is strand-unpeeling, during which the extension increases by 0.01–0.02 nm per base pair of S-DNA converted to single-stranded DNA at forces between 75 and 110 pN; (iii) even in the presence of nicks or free ends, S-DNA can be stable under physiological solution conditions; (iv) distribution of small GC-rich islands in a large DNA plays a key role in determining the transition pathways; and (v) in the absence of nicks or free ends, torsion-unconstrained DNA undergoes the overstretching transition via creation of S-DNA. Our study provides a new, high-resolution understanding of the competition between unpeeling and formation of S-DNA.

A three-state model with loop entropy for the overstretching transition of DNA

We introduce a three-state model for a single DNA chain under tension that distinguishes among B-DNA, S-DNA, and M (molten or denatured) segments and at the same time correctly accounts for the entropy of molten loops, characterized by the exponent c in the asymptotic expression S approximately -c ln n for the entropy of a loop of length n. Force extension curves are derived exactly by employing a generalized Poland-Scheraga approach and then compared to experimental data. Simultaneous fitting to force-extension data at room temperature and to the denaturation phase transition at zero force is possible and allows us to establish a global phase diagram in the force-temperature plane. Under a stretching force, the effects of the stacking energy (entering as a domain-wall energy between paired and unpaired bases) and the loop entropy are separated. Therefore, we can estimate the loop exponent c independently from the precise value of the stacking energy. The fitted value for c is small, suggesting that nicks dominate the experimental force extension traces of natural DNA.

An in silico design for a DNA nanomechanical switch

We have calculated how the charge transfer properties of DNA change in response to the application of an external stretching force. Since charge transfer occurs through the DNA π-stack, any disruption to this stacking causes dramatic changes in the transport properties of the biomolecule, as our calculations demonstrate. We therefore propose that the mechanical response of DNA to an applied stretching force might be used in the design of a nanomechanical switch.

Torsional sensing of small-molecule binding using magnetic tweezers

DNA-binding small molecules are widespread in the cell and heavily used in biological applications. Here, we use magnetic tweezers, which control the force and torque applied to single DNAs, to study three small molecules: ethidium bromide (EtBr), a well-known intercalator; netropsin, a minor-groove binding anti-microbial drug; and topotecan, a clinically used anti-tumor drug. In the low-force limit in which biologically relevant torques can be accessed (<10 pN), we show that ethidium intercalation lengthens DNA ∼1.5-fold and decreases the persistence length, from which we extract binding constants. Using our control of supercoiling, we measure the decrease in DNA twist per intercalation to be 27.3 ± 1° and demonstrate that ethidium binding delays the accumulation of torsional stress in DNA, likely via direct reduction of the torsional modulus and torque-dependent binding. Furthermore, we observe that EtBr stabilizes the DNA duplex in regimes where bare DNA undergoes structural transitions. In contrast, minor groove binding by netropsin affects neither the contour nor persistence length significantly, yet increases the twist per base of DNA. Finally, we show that topotecan binding has consequences similar to those of EtBr, providing evidence for an intercalative binding mode. These insights into the torsional consequences of ligand binding can help elucidate the effects of small-molecule drugs in the cellular environment.

The nature of the force-induced conformation transition of dsDNA studied by using single molecule force spectroscopy

Single-stranded DNA binding proteins (SSB) interact with single-stranded DNA (ssDNA) specifically. Taking advantage of this character, we have employed Bacillus subtilis SSB protein to investigate the nature of force-induced conformation transition of double-stranded DNA (dsDNA) by using AFM-based single molecule force spectroscopy (SMFS) technique. Our results show that, when a dsDNA is stretched beyond its contour length, the dsDNA is partially melted, producing some ssDNA segments which can be captured by SSB proteins. We have also systematically investigated the effects of stretching length, waiting time, and salt concentration on the conformation transition of dsDNA and SSB-ssDNA interactions, respectively. Furthermore, the effect of proflavine, a DNA intercalator, on the SSB-DNA interactions has been investigated, and the results indicate that the proflavine-saturated dsDNA can be stabilized to the extent that the dsDNA will no longer melt into ssDNA under the mechanical force even up to 150 pN, and no SSB-DNA interactions are detectable.

Biophysical characterization of DNA binding from single molecule force measurements

Single molecule force spectroscopy is a powerful method that uses the mechanical properties of DNA to explore DNA interactions. Here we describe how DNA stretching experiments quantitatively characterize the DNA binding of small molecules and proteins. Small molecules exhibit diverse DNA binding modes, including binding into the major and minor grooves and intercalation between base pairs of double-stranded DNA (dsDNA). Histones bind and package dsDNA, while other nuclear proteins such as high mobility group proteins bind to the backbone and bend dsDNA. Single-stranded DNA (ssDNA) binding proteins slide along dsDNA to locate and stabilize ssDNA during replication. Other proteins exhibit binding to both dsDNA and ssDNA. Nucleic acid chaperone proteins can switch rapidly between dsDNA and ssDNA binding modes, while DNA polymerases bind both forms of DNA with high affinity at distinct binding sites at the replication fork. Single molecule force measurements quantitatively characterize these DNA binding mechanisms, elucidating small molecule interactions and protein function.

Different pulling modes in DNA overstretching: a theoretical analysis

We study the thermally driven denaturation of a double-stranded polymer in the presence of a stretching force via Monte-Carlo simulations. When one strand only is stretched, the denaturation transition is first order, while when both strands are stretched, melting is second order. By revisiting the Poland-Scheraga model for DNA melting, we show that at room temperature, the most likely scenario is that DNA melts as it overstretches. Our results are in general agreement with the most recent experiments and suggest how varying temperature and stretching mode may help settle the question whether S-DNA exists or not.

Force spectroscopy and fluorescence microscopy of dsDNA-YOYO-1 complexes: implications for the structure of dsDNA in the overstretching region

When individual dsDNA molecules are stretched beyond their B-form contour length, they reveal a structural transition in which the molecule extends 1.7 times its contour length. The nature of this transition is still a subject of debate. In the first model, the DNA helix unwinds and combined with the tilting of the base pairs (which remain intact), results in a stretched form of DNA (also known as S-DNA). In the second model the base pairs break resulting effectively in two single-strands, which is referred to as force-induced melting. Here a combination of optical tweezers force spectroscopy with fluorescence microscopy was used to study the structure of dsDNA in the overstretching regime. When dsDNA was stretched in the presence of 10 nM YOYO-1 an initial increase in total fluorescence intensity of the dye–DNA complex was observed and at an extension where the dsDNA started to overstretch the fluorescence intensity leveled off and ultimately decreased when stretched further into the overstretching region. Simultaneous force spectroscopy and fluorescence polarization microscopy revealed that the orientation of dye molecules did not change significantly in the overstretching region (78.0°± 3.2°). These results presented here clearly suggest that, the structure of overstretched dsDNA can be explained accurately by force induced melting.

C-terminal domain modulates the nucleic acid chaperone activity of human T-cell leukemia virus type 1 nucleocapsid protein via an electrostatic mechanism

Retroviral nucleocapsid (NC) proteins are molecular chaperones that facilitate nucleic acid (NA) remodeling events critical in viral replication processes such as reverse transcription. Surprisingly, the NC protein from human T-cell leukemia virus type 1 (HTLV-1) is an extremely poor NA chaperone. Using bulk and single molecule methods, we find that removal of the anionic C-terminal domain (CTD) of HTLV-1 NC results in a protein with chaperone properties comparable with that of other retroviral NCs. Increasing the ionic strength of the solution also improves the chaperone activity of full-length HTLV-1 NC. To determine how the CTD negatively modulates the chaperone activity of HTLV-1 NC, we quantified the thermodynamics and kinetics of wild-type and mutant HTLV-1 NC/NA interactions. The wild-type protein exhibits very slow dissociation kinetics, and removal of the CTD or mutations that eliminate acidic residues dramatically increase the protein/DNA interaction kinetics. Taken together, these results suggest that the anionic CTD interacts with the cationic N-terminal domain intramolecularly when HTLV-1 NC is not bound to nucleic acids, and similar interactions occur between neighboring molecules when NC is NA-bound. The intramolecular N-terminal domain-CTD attraction slows down the association of the HTLV-1 NC with NA, whereas the intermolecular interaction leads to multimerization of HTLV-1 NC on the NA. The latter inhibits both NA/NC aggregation and rapid protein dissociation from single-stranded DNA. These features make HTLV-1 NC a poor NA chaperone, despite its robust duplex destabilizing capability.

Stretching and twisting of the DNA duplexes in coarse-grained dynamical models

Three coarse-grained molecular dynamics models of the double-stranded DNA are proposed and compared in the context of single molecule mechanical manipulation such as twisting and various schemes of stretching-unzipping, shearing, two-strand stretching and stretching of only one strand. The models differ in the number of effective beads (between two and five) representing each nucleotide. They all show similar behaviour, but the bigger the resolution, the more details in the force patterns. The models incorporate the effective Lennard-Jones potentials in the couplings between two strands and harmonic potentials to describe the structure of a single strand. The force patterns are shown to depend on the sequence studied. In particular, both shearing and unzipping for an all-AT sequence lead to lower forces than for an all-CG sequence. The unzipping patterns and the corresponding scenario diagrams for the contact rupture events are found to reflect the sequential information if the temperature is moderate and initial transients are discarded. The derived torque-force phase diagram is found to be qualitatively consistent with experiments and all-atom simulations.

Overstretching of a 30 bp DNA duplex studied with steered molecular dynamics simulation: effects of structural defects on structure and force-extension relation

Single-molecule experiments on polymeric DNA show that the molecule can be overstretched at nearly constant force by about 70% beyond its relaxed contour length. In this publication we use steered molecular dynamics (MD) simulation to study the effect of structural defects on force-extension curves and structures at high elongation in a 30 base pair duplex pulled by its torsionally unconstrained 5' -5' ends. The defect-free duplex shows a plateau in the force-extension curve at 120 pN in which large segments with inclined and paired bases ("S-DNA") near both ends of the duplex coexist with a central B-type segment separated from the former by small denaturation bubbles. In the presence of a base mismatch or a nick, force-extension curves are very similar to the ones of the defect-free duplex. For the duplex with a base mismatch, S-type segments with highly inclined base pairs are not observed; rather, the overstretched duplex consists of B-type segments separated by denaturation bubbles. The nicked duplex evolves, via a two-step transition, into a two-domain structure characterized by a large S-type segment coexisting with several short S-type segments which are separated by short denaturation bubbles. Our results suggest that in the presence of nicks the force-extension curve of highly elongated duplex DNA might reflect locally highly inhomogeneous stretching.

Unraveling the structure of DNA during overstretching by using multicolor, single-molecule fluorescence imaging

Single-molecule manipulation studies have revealed that double-stranded DNA undergoes a structural transition when subjected to tension. At forces that depend on the attachment geometry of the DNA (65 pN or 110 pN), it elongates approximately 1.7-fold and its elastic properties change dramatically. The nature of this overstretched DNA has been under debate. In one model, the DNA cooperatively unwinds, while base pairing remains intact. In a competing model, the hydrogen bonds between base pairs break and two single DNA strands are formed, comparable to thermal DNA melting. Here, we resolve the structural basis of DNA overstretching using a combination of fluorescence microscopy, optical tweezers, and microfluidics. In DNA molecules undergoing the transition, we visualize double- and single-stranded segments using specific fluorescent labels. Our data directly demonstrate that overstretching comprises a gradual conversion from double-stranded to single-stranded DNA, irrespective of the attachment geometry. We found that these conversions favorably initiate from nicks or free DNA ends. These discontinuities in the phosphodiester backbone serve as energetically favorable nucleation points for melting. When both DNA strands are intact and no nicks or free ends are present, the overstretching force increases from 65 to 110 pN and melting initiates throughout the molecule, comparable to thermal melting. These results provide unique insights in the thermodynamics of DNA and DNA-protein interactions.