Biophysical characterization of DNA binding from single molecule force measurements

Alba shapes the archaeal genome using a delicate balance of bridging and stiffening the DNA

Architectural proteins have an important role in shaping the genome and act as global regulators of gene expression. How these proteins jointly modulate genome plasticity is largely unknown. In archaea, one of the most abundant proteins, Alba, is considered to have a key role in organizing the genome. Here we characterize the multimodal architectural properties and interplay of the Alba1 and Alba2 proteins using single-molecule imaging and manipulation techniques. We demonstrate that the two paralogues can bridge and rigidify DNA and that the interplay between the two proteins influences the balance between these effects. Our data yield a structural model that explains the multimodal behaviour of Alba proteins and its impact on genome folding.

How the genome is physically organized is less understood in archaea than in eubacteria or eukaryotes. Laurens et al. measure DNA binding by the Sulfolobus solfataricus proteins Alba1 and Alba2 using single-molecule techniques and conclude that the presence of Alba2 leads to more bridging between DNA.

Strong bending of the DNA double helix

During the past decade, the issue of strong bending of the double helix has attracted a lot of attention. Here, we overview the major experimental and theoretical developments in the field sorting out reliably established facts from speculations and unsubstantiated claims. Theoretical analysis shows that sharp bends or kinks have to facilitate strong bending of the double helix. It remains to be determined what is the critical curvature of DNA that prompts the appearance of the kinks. Different experimental and computational approaches to the problem are analyzed. We conclude that there is no reliable evidence that any anomalous behavior of the double helix happens when DNA fragments in the range of 100 bp are circularized without torsional stress. The anomaly starts at the fragment length of about 70 bp when sharp bends or kinks emerge in essentially every molecule. Experimental data and theoretical analysis suggest that kinks may represent openings of isolated base pairs, which had been experimentally detected in linear DNA molecules. The calculation suggests that although the probability of these openings in unstressed DNA is close to 10⁻⁵, it increases sharply in small DNA circles reaching 1 open bp per circle of 70 bp.

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.

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.

Optical tweezers to study viruses

A virus is a complex molecular machine that propagates by channeling its genetic information from cell to cell. Unlike macroscopic engines, it operates in a nanoscopic world under continuous thermal agitation. Viruses have developed efficient passive and active strategies to pack and release nucleic acids. Some aspects of the dynamic behavior of viruses and their substrates can be studied using structural and biochemical techniques. Recently, physical techniques have been applied to dynamic studies of viruses in which their intrinsic mechanical activity can be measured directly. Optical tweezers are a technology that can be used to measure the force, torque and strain produced by molecular motors, as a function of time and at the single-molecule level. Thanks to this technique, some bacteriophages are now known to be powerful nanomachines; they exert force in the piconewton range and their motors work in a highly coordinated fashion for packaging the viral nucleic acid genome. Nucleic acids, whose elasticity and condensation behavior are inherently coupled to the viral packaging mechanisms, are also amenable to examination with optical tweezers. In this chapter, we provide a comprehensive analysis of this laser-based tool, its combination with imaging methods and its application to the study of viruses and viral molecules.

Quantitative assessment of the interplay between DNA elasticity and cooperative binding of ligands

Binding of ligands to DNA can be studied by measuring the change of the persistence length of the complex formed, in single-molecule assays. We propose a methodology for persistence length data analysis based on a quenched disorder statistical model and describing the binding isotherm by a Hill-type equation. We obtain an expression for the effective persistence length as a function of the total ligand concentration, which we apply to our data of the DNA-cationic β-cyclodextrin and to the DNA-HU protein data available in the literature, determining the values of the local persistence lengths, the dissociation constant, and the degree of cooperativity for each set of data. In both cases the persistence length behaves nonmonotonically as a function of ligand concentration and based on the results obtained we discuss some physical aspects of the interplay between DNA elasticity and cooperative binding of ligands.

Range of interaction between DNA-bending proteins is controlled by the second-longest correlation length for bending fluctuations

When a DNA molecule is stretched, the zero-force correlation length for its bending fluctuations-the persistence length A-bifurcates into two different correlation lengths-the shorter "longitudinal" correlation length ξ_{∥}(f) and the longer "transverse" correlation length ξ_{⊥}(f). In the high-force limit, ξ_{∥}(f)=ξ_{⊥}(f)/2=sqrt[k_{B}TA/f]/2. When DNA-bending proteins bind to the DNA molecule, there is an effective interaction between the protein-generated bends mediated by DNA elasticity and bending fluctuations. Surprisingly, the range of this interaction is not the longest correlation length associated with transverse fluctuations of the tangent vector along the polymer, but instead is the second longest longitudinal correlation length ξ_{∥}(f,μ). The effect arises from the protein-bend contribution to the Hamiltonian having an axial rotational symmetry which eliminates its coupling to the transverse fluctuations.

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.

Mesoscopic models for DNA stretching under force: New results and comparison with experiments

Single-molecule experiments on double-stranded B-DNA stretching have revealed one or two structural transitions, when increasing the external force. They are characterized by a sudden increase of DNA contour length and a decrease of the bending rigidity. The nature and the critical forces of these transitions depend on DNA base sequence, loading rate, salt conditions and temperature. It has been proposed that the first transition, at forces of 60-80 pN, is a transition from B to S-DNA, viewed as a stretched duplex DNA, while the second one, at stronger forces, is a strand peeling resulting in single-stranded DNAs (ssDNA), similar to thermal denaturation. But due to experimental conditions these two transitions can overlap, for instance for poly(dA-dT). In an attempt to propose a coherent picture compatible with this variety of experimental observations, we derive an analytical formula using a coupled discrete worm-like chain-Ising model. Our model takes into account bending rigidity, discreteness of the chain, linear and non-linear (for ssDNA) bond stretching. In the limit of zero force, this model simplifies into a coupled model already developed by us for studying thermal DNA melting, establishing a connection with previous fitting parameter values for denaturation profiles. Our results are summarized as follows: i) ssDNA is fitted, using an analytical formula, over a nano-Newton range with only three free parameters, the contour length, the bending modulus and the monomer size; ii) a surprisingly good fit on this force range is possible only by choosing a monomer size of 0.2 nm, almost 4 times smaller than the ssDNA nucleobase length; iii) mesoscopic models are not able to fit B to ssDNA (or S to ss) transitions; iv) an analytical formula for fitting B to S transitions is derived in the strong force approximation and for long DNAs, which is in excellent agreement with exact transfer matrix calculations; v) this formula fits perfectly well poly(dG-dC) and λ-DNA force-extension curves with consistent parameter values; vi) a coherent picture, where S to ssDNA transitions are much more sensitive to base-pair sequence than the B to S one, emerges. This relatively simple model might allow one to further study quantitatively the influence of salt concentration and base-pairing interactions on DNA force-induced transitions.

Determining the DNA stability parameters for the breathing dynamics of heterogeneous DNA by stochastic optimization

We suggest that the thermodynamic stability parameters (nearest neighbor stacking and hydrogen bonding free energies) of double-stranded DNA molecules can be inferred reliably from time series of the size fluctuations (breathing) of local denaturation zones (bubbles). On the basis of the reconstructed bubble size distribution, this is achieved through stochastic optimization of the free energies in terms of simulated annealing. In particular, it is shown that even noisy time series allow the identification of the stability parameters at remarkable accuracy. This method will be useful to obtain the DNA stacking and hydrogen bonding free energies from single bubble breathing assays rather than equilibrium data.

Force spectroscopy reveals the DNA structural dynamics that govern the slow binding of Actinomycin D

Actinomycin D (ActD) is a small molecule with strong antibiotic and anticancer activity. However, its biologically relevant DNA-binding mechanism has never been resolved, with some studies suggesting that the primary binding mode is intercalation, and others suggesting that single-stranded DNA binding is most important. To resolve this controversy, we develop a method to quantify ActD’s equilibrium and kinetic DNA-binding properties as a function of stretching force applied to a single DNA molecule. We find that destabilization of double stranded DNA (dsDNA) by force exponentially facilitates the extremely slow ActD-dsDNA on and off rates, with a much stronger effect on association, resulting in overall enhancement of equilibrium ActD binding. While we find the preferred ActD–DNA-binding mode to be to two DNA strands, major duplex deformations appear to be a pre-requisite for ActD binding. These results provide quantitative support for a model in which the biologically active mode of ActD binding is to pre-melted dsDNA, as found in transcription bubbles. DNA in transcriptionally hyperactive cancer cells will therefore likely efficiently and rapidly bind low ActD concentrations (∼10 nM), essentially locking ActD within dsDNA due to its slow dissociation, blocking RNA synthesis and leading to cell death.

EMSA and single-molecule force spectroscopy study of interactions between Bacillus subtilis single-stranded DNA-binding protein and single-stranded DNA

In this article, interactions between Bacillus subtilis single-stranded DNA binding proteins (BsSSB) and single-stranded DNA (ssDNA) were systematically studied. The effect of different molar ratios between BsSSB and ssDNA on their binding modes was first investigated by electrophoretic mobility shift assays (EMSAs). It is found that a high molar ratio of BsSSB to ssDNA can produce BsSSB-ssDNA complexes formed in the mode of two proteins binding one 65-nt (nucleotide) ssDNA whereas a low molar ratio facilitates the formation of BsSSB-ssDNA complexes in the mode of one protein binding one 65-nt ssDNA. Furthermore, two binding modes are in dynamic equilibrium. The unbinding force of BsSSB-ssDNA complexes was measured quantitatively in solutions with different salt concentrations by using AFM-based single-molecule force spectroscopy (SMFS). Our results show that the unbinding force is about 10 pN higher at high salt concentration (0.5 M NaCl) than at low salt concentration (0.1 M NaCl) and the lifetime of BsSSB-ssDNA complexes at high salt concentration is twice as long as that at low salt concentration. These results indicate that more tightly packed BsSSB-ssDNA complexes can form at high salt (0.5 M NaCl) concentration. In addition, the results of EMSA show that ssDNA, which is bound to BsSSB, can dissociate from BsSSB in the presence of the cDNA strand, indicating the dynamic nature of BsSSB-ssDNA interactions.

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.

B-DNA to zip-DNA: simulating a DNA transition to a novel structure with enhanced charge-transport characteristics

The forced extension of a DNA segment is studied in a series of steered molecular dynamics simulations, employing a broad range of pulling forces. Throughout the entire force range, the formation of a zipper-like (zip-) DNA structure is observed. In that structure, first predicted by Lohikoski et al., the bases of the DNA strands interdigitate with each other and form a single-base aromatic stack. Similar motifs, albeit only a few base pairs in extent, have been observed in experimental crystal structures. Analysis of the dynamics of structural changes in pulled DNA shows that S-form DNA, thought to be adopted by DNA under applied force, serves as an intermediate between B-DNA and zip-DNA. Therefore, the phase transition plateau observed in force-extension curves of DNA is suggested to reflect the B-DNA to zip-DNA structural transition. Electronic structure analysis of purine bases in zip-DNA indicates a several-fold to order of magnitude increase in the π-π electronic coupling among nearest-neighbor nucleobases, compared to B-DNA. We further observe that zip-DNA does not require base pair complementarity between DNA strands, and we predict that the increased electronic coupling in zip-DNA will result in a much higher rate of charge transfer through an all-purine zip-DNA compared to B-DNA of equal length.

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.

Combining optical trapping, fluorescence microscopy and micro-fluidics for single molecule studies of DNA-protein interactions

Complexity and heterogeneity are common denominators of the many molecular events taking place inside the cell. Single-molecule techniques are important tools to quantify the actions of biomolecules. Heterogeneous interactions between multiple proteins, however, are difficult to study with these technologies. One solution is to integrate optical trapping with micro-fluidics and single-molecule fluorescence microscopy. This combination opens the possibility to study heterogeneous/complex protein interactions with unprecedented levels of precision and control. It is particularly powerful for the study of DNA-protein interactions as it allows manipulating the DNA while at the same time, individual proteins binding to it can be visualized. In this work, we aim to illustrate several published and unpublished key results employing the combination of fluorescence microscopy and optical tweezers. Examples are recent studies of the structural properties of DNA and DNA-protein complexes, the molecular mechanisms of nucleo-protein filament assembly on DNA and the motion of DNA-bound proteins. In addition, we present new results demonstrating that single, fluorescently labeled proteins bound to individual, optically trapped DNA molecules can already be tracked with localization accuracy in the sub-10 nm range at tensions above 1 pN. These experiments by us and others demonstrate the enormous potential of this combination of single-molecule techniques for the investigation of complex DNA-protein interactions.

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.

Single-molecule stretching studies of RNA chaperones

RNA chaperone proteins play significant roles in diverse biological contexts. The most widely studied RNA chaperones are the retroviral nucleocapsid proteins (NC), also referred to as nucleic acid (NA) chaperones. Surprisingly, the biophysical properties of the NC proteins vary significantly for different viruses, and it appears that HIV-1 NC has optimal NA chaperone activity. In this review we discuss the physical nature of the NA chaperone activity of NC. We conclude that the optimal NA chaperone must saturate NA binding, leading to strong NA aggregation and slight destabilization of all NA duplexes. Finally, rapid kinetics of the chaperone protein interaction with NA is another primary component of its NA chaperone activity. We discuss these characteristics of HIV-1 NC and compare them with those of other NA binding proteins and ligands that exhibit only some characteristics of NA chaperone activity, as studied by single molecule DNA stretching.

Temperature dependence of DNA persistence length

We have determined the temperature dependence of DNA persistence length, a, using two different methods. The first approach was based on measuring the j-factors of short DNA fragments at various temperatures. Fitting the measured j-factors by the theoretical equation allowed us to obtain the values of a for temperatures between 5°C and 42°C. The second approach was based on measuring the equilibrium distribution of the linking number between the strands of circular DNA at different temperatures. The major contribution into the distribution variance comes from the fluctuations of DNA writhe in the nicked circular molecules which are specified by the value of a. The computation-based analysis of the measured variances was used to obtain the values of a for temperatures up to 60°C. We found a good agreement between the results obtained by these two methods. Our data show that DNA persistence length strongly depends on temperature and accounting for this dependence is important in quantitative comparison between experimental results obtained at different temperatures.