RecA filament sliding on DNA facilitates homology search

Open accessCopyright infoDownload PDFRelated contenteLife Lens

RecA filament sliding on DNA facilitates homology search

Affiliation details

This article has been corrected
DOI: http://dx.doi.org/10.7554/eLife.00067Published December 13, 2012 Cite as eLife 2012;1:e00067

Abstract

During homologous recombination, RecA forms a helical filament on a single stranded (ss) DNA that searches for a homologous double stranded (ds) DNA and catalyzes the exchange of complementary base pairs to form a new heteroduplex. Using single molecule fluorescence imaging tools with high spatiotemporal resolution we characterized the encounter complex between the RecA filament and dsDNA. We present evidence in support of the ‘sliding model’ wherein a RecA filament diffuses along a dsDNA track. We further show that homology can be detected during sliding. Sliding occurs with a diffusion coefficient of approximately 8000 bp2/s allowing the filament to sample several hundred base pairs before dissociation. Modeling suggests that sliding can accelerate homology search by as much as 200 fold. Homology recognition can occur for as few as 6 nt of complementary basepairs with the recognition efficiency increasing for higher complementarity. Our data represents the first example of a DNA bound multi-protein complex which can slide along another DNA to facilitate target search.

DOI: http://dx.doi.org/10.7554/eLife.00067.001

eLife digest

The DNA molecules in cells are continuously bombarded with radiation, chemicals and other agents, and it is important for cells to repair the damage caused by these before the process of cell division begins. Most DNA molecules consist of two single strands of DNA that are held together by hydrogen bonds in the familiar double-helix structure. Of the various types of damage that DNA molecules are prone to, double-strand breaks are among the most dangerous because they can lead to cancer if they are not repaired.

DNA molecules use four bases—adenine, cytosine, guanine, and thymine—to store genetic information. In single-stranded DNA these bases are attached to a backbone made of alternating sugar and phosphate groups. A crucial feature of double-stranded DNA is that the sequences of bases in the two strands are complementary to each other—adenine is always paired with thymine, and cytosine is always paired with guanine. However, the hydrogen bonds that hold the pairs of bases together are quite weak, which means that the two strands of the double helix can be pulled apart quite easily. The ease with which these bonds can be formed and broken is crucial for many genetic processes.

One way to repair a double strand break is to replace the damaged stretch of DNA with an undamaged stretch from another DNA molecule. This process of swapping DNA molecules, which is called strand exchange, is catalyzed by a protein that is able to interact with two DNA molecules at the same time. An important first step within this process is identifying the stretch of DNA that can be used to repair the break.

Ragunathan et al. now report evidence from experiments on Escherichia coli that support a model in which the protein catalyst (RecA in the case of E. coli) combines with a single strand of DNA to form a filamentous DNA–protein complex (RecA filament) that can then slide along a double-stranded DNA molecule to search for a complementary sequence of base pairs. High-resolution fluorescent imaging reveals that the RecA filament is able to sample several hundred base pairs before the filament dissociates from the DNA and rebinds at a different location. The sliding was largely driven by electrostatic interactions between the RecA filament and the double-stranded DNA, and the filament was capable of identifying matching sequences that contained as few as six matching bases.

Ragunathan et al. estimate that sliding is about two orders of magnitude faster at finding matching sequences compared to mechanisms that do not involve sliding, such as models that rely solely on chance encounters between DNA molecules and the RecA filament. By showing that a DNA–protein complex can slide along another DNA molecule to search for a target, these results could lead to new insights into other systems in which it is necessary for protein-nucleic acid complexes to locate a particular sequence of bases.

DOI: http://dx.doi.org/10.7554/eLife.00067.002

Main text

Introduction

The ubiquitous presence of DNA damaging agents poses a constant threat to genome integrity and protein machineries are required to repair DNA damage. Homologous recombination is one of the pathways involved in double strand break repair (Cox et al., 2000). An important step in homologous recombination is the reciprocal exchange of basepairs between complementary DNA molecules during a reaction called strand exchange, which is catalyzed by RecA in E. coli. Homologs of RecA, Rad51 and Dmc1 in eukaryotes and RadA in archaea, carry out similar functions during DNA repair underscoring the central role of proteins catalyzing strand exchange across all forms of life (Bianco et al., 1998).

The strand exchange reaction involves three steps: (1) pre-synapsis, (2) synapsis, and (3) heteroduplex extension via branch migration. Pre-synapsis involves the assembly of RecA monomers on single stranded (ss) DNA in the presence of ATP with a stoichiometry of 3 nt per monomer (Di Capua et al., 1982; Dombroski et al., 1983). RecA forms a filament which stretches the ssDNA to a length of 1.5 times the length of B-form DNA (Stasiak et al., 1981; Dunn et al., 1982; Stasiak and Egelman, 1986). During synapsis, the RecA filament finds a homologous double stranded (ds) DNA and catalyzes the exchange of complementary base pairs to form a new heteroduplex product. Subsequently, branch migration mediates the extension of the heteroduplex product (Cox and Lehman, 1981).

Homology search is the first step in synapsis and is arguably the most mysterious aspect of the strand exchange reaction (Barzel and Kupiec, 2008). A RecA filament has to rapidly and accurately find a homologous sequence in the presence of a vast excess of non-homologous dsDNA. The multivalent nature of the RecA filament enables it to sustain contact with a long dsDNA while allowing for frequent dissociation events that permits rapid sampling of different segments of the DNA, thus facilitating homology search (Forget and Kowalczykowski, 2012). Experiments involving the mechanical manipulation of the incoming dsDNA upon binding to the RecA filament identified structural intermediates which facilitate homology recognition and provided insight into mechanisms by which RecA can discriminate between homologous and non-homologous DNA sequences (Danilowicz et al., 2012; De Vlaminck et al., 2012; Peacock-Villada et al., 2012). However, random 3D collision between a RecA filament and the dsDNA even when aided by intersegmental transfer is unlikely to result in precise alignment of matching sequences. For this reason, RecA filament sliding on dsDNA, even over a short distance range of ten to hundreds of bp would significantly accelerate the search for homology. Although many proteins are now known to diffuse along the DNA via 1D sliding (Blainey et al., 2006; Wang et al., 2006; Gorman et al., 2007; Bonnet et al., 2008; Liu et al., 2008; Roy et al., 2009; Zhou et al., 2011) RecA filament diffusion on DNA has not been observed before and in fact an early study ruled out a role for long range (approximately several kb) 1D sliding during homology search (Adzuma, 1998). Here, we utilized the high spatio-temporal resolution of single molecule FRET (Ha et al., 1996) to examine the possibility of RecA filament sliding during homology search.

Using two and three color FRET measurements, we show that RecA filament slides along dsDNA, primarily mediated by electrostatic interactions. Furthermore, by using sequences with short stretches of homology, we could monitor repeated events of homology recognition and disengagement in real time without full dissociation of RecA filament from dsDNA. RecA-mediated homology search is the first example of a DNA bound multi protein complex which can slide on another DNA during target search process. It may also serve as a canonical example for other proteins (e.g., telomerase and Argonaute) which are bound to nucleic acid sequences that act as ‘guide’ strands conferring target site specificity.

Results

Dynamic interactions between RecA filament and non-homologous dsDNA

For short (<80 bp) homologous dsDNA substrates, homology search is completed rapidly (within 30 ms) once the dsDNA encounters a RecA filament (Ragunathan et al., 2011). Therefore, we used a non-homologous dsDNA to avoid stable product formation and monitor RecA filament in the act of homology search.

We immobilized a partial dsDNA with a 5′ 39 nt ssDNA tail on a passivated quartz surface via biotin–neutravidin interaction (Figure 1A). The DNA is labeled with a FRET acceptor (Cy5) at the ssDNA/dsDNA junction. We first formed a stable RecA filament on the DNA by using ATPγS as the cofactor. Then, we added a solution containing non-homologous dsDNA (LdsDNA = 39 bp) and ATPγS while simultaneously removing free RecA from solution, allowing us to observe the interaction solely between the incoming non-homologous dsDNA (free of RecA) and a single isolated RecA filament. RecA filaments formed under these conditions are stable and can carry out the strand exchange reaction with a homologous dsDNA (Ragunathan et al., 2011).

Docking of non-homologous dsDNA to the RecA filament is detected as an abrupt appearance of fluorescence signal from the background level. After docking, we observed large and rapid fluctuations in FRET detected as anti-correlated changes of donor and acceptor intensities. The FRET fluctuations are indicative of extensive distance changes between the donor on the dsDNA and the acceptor on the RecA/ssDNA filament. Single molecule time traces show multiple dsDNA binding and dissociation events to the same filament because the observed interactions are transient without forming a stable product (Figure 1B). The lifetime of the binding events is exponentially distributed with an average lifetime of 3.5 s (Figure 1C). Filaments formed with ATP displayed similar FRET fluctuations upon docking of non-homologous dsDNA and exhibit comparable dissociation times (Figure 1—figure supplement 1A–C). Hence, the fluctuations observed here do not require ATP hydrolysis. Using different lengths for the ssDNA tail (Lfilament = 50 or 99 nt) did not significantly change the lifetime of the encounter complex (Figure 1—figure supplement 2A). In contrast, increasing dsDNA length increased the lifetime presumably due to a larger number of contacts between the dsDNA and the RecA filament (Figure 1—figure supplement 2B).

We determined the lifetime of the transient encounter complex as a function of magnesium or sodium concentrations. Rapid FRET fluctuations persisted under all the solution conditions tested while the lifetime of the complex decreased with increasing magnesium or sodium concentrations suggesting that the interaction is electrostatic in nature and can be weakened by increased screening (Figure 1D–F). Indeed, eliminating negative charge from the RecA C-terminus (Lusetti et al., 2003) by deletion of acidic residues resulted in an increased lifetime of dsDNA bound to the RecA filament (Figure 1—figure supplement 3).

FRET fluctuations are due to sliding of RecA filament

We considered two possible explanations for the large FRET fluctuations observed upon docking of dsDNA to the RecA filament. The first involves limited unwinding of the dsDNA (Bianchi et al., 1985) by the RecA filament resulting in separation of the duplex ends. Repeated melting and annealing transitions at the labeled end of the dsDNA could in principle result in FRET fluctuations. We tested this possibility using a dsDNA (LdsDNA = 39 bp) labeled at one duplex end (donor and acceptor fluorophores on the two opposing strands) so that local melting would cause a FRET decrease and vice versa. When the dsDNA docks to an unlabeled RecA filament immobilized on the surface, we observed stable high FRET (Figure 2—figure supplement 1) suggesting that separation of the duplex ends is not the source of FRET fluctuations described in Figure 1.

The second explanation involves 1D sliding (or diffusion) of the RecA filament along dsDNA. In order to test the sliding model, we examined whether changes in the length of the RecA filament would affect the time scale of FRET fluctuations since such a change would modulate the encounter frequency between the donor and the acceptor fluorophores. In contrast, conformational changes either within the dsDNA or the protein would presumably exhibit the same time scale of FRET fluctuations independent of RecA filament length.

Shorter RecA filaments (Lfilament = 21 nt) exhibited more rapid FRET changes and of smaller amplitudes (Figure 2A) compared to longer filaments (Lfilament = 99 nt) which showed larger and slower changes in FRET (Figure 2B). To quantify the time scale of FRET fluctuations, we determined cross-correlation of the donor and acceptor fluorescence intensities for four different filament lengths (Lfilament = 21, 39, 69 and 99 nt) and a single length of non-homologous dsDNA (LdsDNA = 39 bp). The average cross correlation time increased with increasing filament lengths (Figure 2C), supporting the 1D sliding model. Furthermore, the histogram of FRET efficiencies (E) shifted towards lower FRET values for increasing filament lengths (Figure 2D). This observation also supports the 1D sliding model because the two fluorophores would spend a smaller fraction of time in close proximity if the filament is longer. To estimate the diffusion coefficient of the 1D sliding process, we performed Monte Carlo simulations of dsDNA diffusing along a RecA filament (‘Experimental procedures’ and Figure 2—figure supplement 2). Given its large persistence length of about 800 nm (Hegner et al., 1999), we treated the RecA filament as a rigid rod. We then simulated time traces of donor and acceptor intensities of dsDNA bound to RecA filament of various lengths using different pre-assigned diffusion coefficients (Dslide) and then calculated the average cross correlation times for each case. By comparing the simulation results with the experimental data we estimated the diffusion constant Dslide for 1D sliding of dsDNA relative to the RecA filament to be approximately 0.9 × 10−3 μm2/s or 7700 bp2/s. Additionally the time scale of sliding (as measured by the cross correlation time) does not change as a function of sodium or magnesium ion concentrations in solution while the sliding time scale increases twofold when ATPγS is replaced by ATP (Figure 2—figure supplement 3A–B). However, the use of ATP as a co-factor, especially in experiments involving short ssDNA substrates is complicated by the dissociation of RecA monomers from the ssDNA ends leading to filament instabilities which could potentially disrupt sliding along the RecA filament track.

Three color experiments support the sliding model

To further test the 1D sliding model, we designed a three-color FRET assay. Here, the immobilized RecA filament is labeled with two different acceptor fluorophores, Cy5 at the ssDNA/dsDNA junction and Cy7 in the middle of the ssDNA embedded within the RecA filament (Figure 3A). The two acceptors are separated by 33 nt, which upon RecA binding results in a large separation and negligible FRET between them (Joo et al., 2006). The sliding model predicts anticorrelated changes between the two FRET efficiencies, one between Cy3 and Cy5 (ECy3-Cy5) and the other between Cy3 and Cy7 (ECy3-Cy7) because when Cy3 on the dsDNA approaches Cy5, it should move away from Cy7 and vice versa. We observed large and rapid fluctuations in FRET from Cy3 to Cy5 (ECy3-Cy5) and FRET from Cy3 to Cy7 (ECy3-Cy7) (Figure 3B, bottom panel). Consistent with the sliding model, the time trace of ECy3-Cy5 and ECy3-Cy7 exhibits anticorrelation between the two FRET efficiencies. The timescale of fluctuations determined from the cross-correlation of the two FRET efficiencies (Txcorr-3color = 0.04 s) (Figure 3C) is similar to that measured using the two color assay (Figure 2C). The scatter plot of ECy3-Cy5 vs ECy3-Cy7 and calculation of Pearson's correlation coefficient for the two FRET efficiencies (rpearson = −0.6) provide further support for a negative correlation between ECy3-Cy5 and ECy3-Cy7 (Figure 3D). Hence, cumulatively the two and three color FRET results support the sliding model.

Figure 3.
Download figureOpen in new tabDownload powerpointFigure 3. Three color FRET observations support RecA filament sliding.

(A) A schematic of the single molecule three color FRET assay to measure RecA filament sliding. ssDNA (Lfilament = 99 nt) labeled with two acceptor fluorophores (Cy5-red and Cy7-black) with a separation of 33 nt between the fluorophores, was immobilized on the surface. Upon docking of non-homologous donor (Cy3) labeled dsDNA to the pre-formed RecA filament formation, sliding predicts anticorrelated emissions between the two acceptors. (B) Single molecule time traces of Cy3 (green), Cy5 (red) and Cy7 (black) intensities (top panel). Corresponding FRET time traces of FRET between Cy3 and Cy5 (ECy3-Cy5-blue) and FRET between Cy3 and Cy7 (ECy3-Cy7-grey). (C) Normalized cross correlation plot of ECy3-Cy5 and ECy3-Cy7 averaged over 30 molecules and a single exponential fit of the data is overlaid (black). (D). Scatter plot of ECy3-Cy5 and ECy3-Cy7 for 30 molecules showing unique high FRET regions along both axes.

DOI: http://dx.doi.org/10.7554/eLife.00067.011

Homology detection during sliding

Can the sliding of RecA filament along dsDNA be a physiologically relevant activity, that is, can the RecA filament recognize a homologous sequence during sliding? To answer this question, we embedded two repeats of an identical sequence at positions HS1 and HS2 within an otherwise non-homologous ssDNA (Figure 4A). If the target dsDNA contains base pairs which are homologous to the short repeat sequence, it may be possible to observe back and forth sliding events between the two homology sites (HS1 and HS2) without full dissociation of the dsDNA from the RecA filament. Such an observation would indicate that RecA filament sliding allows for homology recognition and base pairing. In all cases, HS1 and HS2 are complementary to a sequence in close proximity to the donor labeled end of the target dsDNA and the location and spacing between the two homology sites were chosen to be within a FRET sensitive regime (approximately 20–80 Å).

We confirmed that dsDNA docking to a RecA filament formed on a ssDNA homopolymer sequence (poly-T sequence, Lfilament = 50 nt) preserves the large and rapid FRET fluctuations (Figure 4—figure supplement 1 and Figure 2B). We then introduced two identical sequences of length, Lh = 5 nt, at HS1 and HS2 positions which are complementary to a 5 bp sequence within the target dsDNA (LdsDNA = 39 bp). Binding of the dsDNA to the RecA filament formed on the ssDNA with the two 5 nt repeats displayed rapid FRET fluctuations across a broad range of FRET values (Figure 4B,E) similar to that observed using a poly-T DNA sequence (Figure 4—figure supplement 1). In contrast, when we increased the length of homology at HS1 and HS2 by a single nucleotide to Lh = 6 nt, single molecule time traces displayed transitions between discrete FRET states (Figure 4C). The same trend was observed upon further increasing the homology length to Lh = 7 nt (Figure 4D). The resulting FRET histograms for 6 and 7 nt homology lengths showed distinct FRET peaks (Figure 4F,G).

In order to confirm that the observed FRET peaks arise from base pairing and recognition at specific homology sites, we analyzed different ssDNA sequences containing only one of the two 6 nt homology sites (HS1 or HS2) and obtained a distinct FRET peak at either approximately 0.9 or 0.5, respectively (Figure 4H,I). Thus, we assigned the highest FRET state (approximately 0.9) to homology recognition and base pairing at HS1 and the mid FRET state (approximately 0.5) to homology recognition and base pairing at HS2. The lowest FRET state, E ∼ 0.1, is likely to correspond to dsDNA sliding outside the boundary of HS1 and HS2 in a FRET insensitive regime. Given that this location lacks stable base pairing interactions, we referred to the low FRET state as the non-homologous state (NH). Although the gap between HS1 and HS2 is 5 nt, it is noteworthy that the total distance traversed by dsDNA to exhibit complete basepairing is (Lh + 5) nt. Thus, for the smallest Lh of 6 nt, this translates to a distance of approximately 50 Ǻ (11 nt × 3.4 Ǻ × 1.5). Given the rapid movements of the dsDNA between adjacent homology sites and the fact that both the dsDNA and RecA filament are stiff at the length scales of our experiments, we cannot think of mechanisms other than sliding to explain the observed transitions. We used a statistical approach based on Hidden Markov Model (HMM) analysis to make unbiased assignments of the various FRET states present within each single molecule time trace. The FRET transitions obtained from HMM analysis were then plotted in the form of a transition density plot (TDP) which is a 2D histogram reflecting the frequency of transitions between the various FRET states. For Lh = 6 nt, the TDP obtained from 236 molecules exhibiting 13,020 transitions displays three distinct FRET states, E ∼ 0.1, 0.5 and 0.9 (Figure 4J). Furthermore, the HMM/TDP analysis gave the transition rates and the number of transitions between the different FRET states. The rate of departing HS1 or HS2 was threefold lower than that of leaving NH, showing that even 6 bp homology can significantly stabilize joint molecule formation (Figure 4L). Similarly, we fitted the data for Lh = 7 nt using HMM analysis (191 molecules, 3819 transitions) and extracted the transition rates and number of transitions from the resulting TDP (Figure 4K). As expected, we observed a slower rate of leaving HS1 and HS2 for Lh = 7 nt compared to Lh = 6 nt (Figure 4L). Also, the dissociation rate of dsDNA bound to HS1 was higher than that of dsDNA bound to HS2 possibly because HS1 is located at the filament end where the terminal RecA monomer at the ssDNA–dsDNA junction binds with altered kinetics (Joo et al., 2006) as opposed to HS2 which is embedded within the RecA filament and benefits from complementary basepairing with the homology site and stabilizing interactions with neighboring RecA monomers.

Interestingly, transitions did not occur solely between neighboring sites (NH<->HS2 and HS1<->HS2 type transitions) but also involved transitions of the type NH->HS1 indicating that the RecA filament might overlook a region of homology during the sliding process. Homology recognition and base pairing of dsDNA at either HS1 or HS2 during sliding is a stochastic process. For a dsDNA initially bound to a non-homologous site (NH), the efficiency with which the first homology site the dsDNA encounters, HS2, is recognized before the more distant homology site, HS1, is recognized is given by,ErecognitionNH->HS2=NNH->HS2(NNH->HS2+NNH->HS1),where NNH->HS1 and NNH->HS2 denote the number of transitions from NH->HS1 and NH->HS2 respectively. For 6 nt homology (Lh = 6 nt) the efficiency of homology site recognition, ErecognitionNH->HS2 is approximately 55%. This efficiency increased for 7 nt homology (ErecognitionNH->HS2 for Lh = 7 nt is approximately 77%). Cumulatively, our data suggest that the efficiency of homology recognition by the RecA filament increases with the number of homologous nucleotides available for basepairing. We could not perform a similar analysis for Lh = 8 nt since dsDNA bound to HS1 and HS2 positions in this case exhibited stable FRET states and transitions between them were rare (Figure 4—figure supplement 2).

Discussion

Our study establishes RecA filament sliding as a possible mechanism to assist homology search during strand exchange reaction. Until now, there has been no example of a multi-protein complex bound to DNA such as the RecA filament which is capable of sliding along a second DNA strand to locate a matching sequence of bases. We estimated the diffusion constant, Dslide, for RecA filament sliding to be approximately 0.9 × 10−3 μm2/s or 7700 bp2/s. Because the lifetime of heterologous synapses ranges from 0.5 to 10 s (Figure 1—figure supplement 2B) (Mani et al., 2009; Forget and Kowalczykowski, 2012), we estimate that dsDNA would diffuse over a length of 60–300 bp during each encounter with a RecA filament prior to its dissociation.

Using the estimated diffusion coefficient and sliding distance, we calculated the effective rate enhancement of target search due to sliding which is given by the ratio of the rate of association of RecA filament with the target site in the presence of sliding to the rate of association with the target site in the absence of sliding (D = 0) (Hammar et al., 2012). In this hypothetical model, RecA filaments can bind non-specifically to DNA but do not slide (Figure 4—figure supplement 3). Though, we do not know the exact microscopic rate constants which determine the binding probability of the RecA filament to the target site (pbind), we can estimate that even for a recognition probability as low as approximately 1%, sliding leads to a rate enhancement in homology search by approximately 200-fold. If no sliding is allowed and homology must be recognized via fortuitous 3D contact that aligns the dsDNA and filament in perfect registry, homology recognition would be about two orders of magnitude slower. Therefore, 1D sliding combined with intersegmental transfer proposed by Forget and Kowalcyzkowski (2012) can dramatically accelerate homology search.

We cannot rule out the possibility that changes in angle between the dsDNA and the RecA filament (off-axis motions) might affect the observed FRET changes. Incorporating such effects of off-axis motions into the model would increase the time scale of linear sliding but the effect will be small because control experiments such as the DNA length dependence (Figure 2C) showed that on-axis motion is the dominant source of FRET fluctuations. Although our experimental scheme cannot probe for small and dynamic off-axis motions between the RecA filament and dsDNA, we can at least rule out the extreme possibility of ‘dsDNA flipping’ in which case the dsDNA may undergo a 180° change in orientation while remaining bound to the RecA filament (Figure 4—figure supplement 4). Therefore, a dsDNA that binds to a RecA filament in the wrong orientation cannot establish the correct orientation required for homology recognition and basepairing without full dissociation and rebinding.

We also demonstrated that homology recognition and base pairing processes can occur during sliding. Our results showed that as few as 6 nt of complementary base pairs are sufficient to act as the unit of homology recognition for RecA mediated homology search. It is noteworthy that recognition of a homology site is a stochastic process and the efficiency of homology recognition by the RecA filament depends on the length of homology which is encountered at a particular site. Furthermore, we show that adding even a single nucleotide to the minimum recognition unit of 6 nt improves recognition efficiency by approximately 1.5-fold.

Proteins mediate target search reactions by a combination of 1D sliding and 3D diffusion based processes (Gorman and Greene, 2008). While intersegmental transfer serves to bring non-contiguous segments of dsDNA close to each other, 1D sliding during homology search could help RecA reorganize the initial synaptic complex in cis without the need for a RecA filament to fully dissociate and rebind to a nearby homology site located within a short distance (up to few hundred base pairs). Given the ability of the RecA filament to slide between nearby homology sites (Figure 4), sliding could serve as a mechanism to rapidly scan neighboring sequences for the existence of an optimal seed sequence from which base pair propagation and heteroduplex extension reactions can proceed. An earlier report demonstrated sliding of Rad51 oligomers (not a Rad51 filament bound to ssDNA) without any clear functional role attributed to the observed sliding activity (Graneli et al., 2006). Because of significant structural and functional similarities between RecA and Rad51 (Conway et al., 2004), it is likely that sliding might also play a role in homology search mediated by Rad51.

Materials and methods

DNA preparation

DNA oligos used in our measurements were purchased from Integrated DNA Technologies (IDT, Coralville, IA). The oligos were suspended in T50 buffer (10 mM Tris-Cl, 50 mM NaCl, pH 8.0). The DNA sequences used in our measurements are specified in Supplementary file 1. Double strand DNA was prepared by mixing complementary DNA molecules and heating to 90°C followed by slow cooling to room temperature over a period of 2 hr. dsDNA was purified from free ssDNA using a 12% native PAGE gel to ensure the absence of free ssDNA in our dsDNA preparations. The partial duplex DNA molecules for immobilization were prepared by annealing with the biotin_DNA sequence. All DNA molecules were labeled at the terminal ends (5′ or 3′) with Cy3 or Cy5 (labeling performed by IDT) as specified in the experimental scheme. For internally labeled DNA oligos, Cy3 N-hydroxysuccinimido (NHS) ester and Cy5 NHS ester (GE Healthcare) were internally labeled to a dT of ssDNA modified via a C6 amino linker (IDT, Coralville, IA).

Single molecule RecA filament and DNA binding assay

The quartz (Finkenbeiner, Waltham, MA) surface is passivated with Polyethylene glycol (m-PEG-5000; Laysan Bio Inc.) and 1–2% biotinylated PEG (biotin-PEG-5000; Laysan Bio Inc., Arab, AL). The coating of the quartz imaging surface with PEG (Roy et al., 2008) eliminated effects due to non-specific binding of proteins. Acceptor labeled reference ssDNA molecules were immobilized on the passivated surface by means of a biotin–neutravidin interaction. After washing away excess of acceptor molecules, the reference ssDNA was incubated with 1 μM RecA (Epicenter biotechnologies, Madison, WI) and 1 mM ATPγS (EMD Calbiochem, Billerica, MA) in an incubation buffer containing 25 mM Tris Acetate pH 7.5, 100 mM Sodium Acetate and 1 mM Magnesium Acetate. In some cases, 1 mM ATP was used instead of ATPγS. RecA ΔC17 protein was a generous gift from Dr. M. M. Cox (University of Wisconsin). After incubation for 15 min to ensure complete filament formation on ssDNA (Lfilament) molecules, the buffer in the chamber was exchanged with a solution of non-homologous dsDNA (1 nM) and 1 mM ATPγS in a strand exchange buffer (25 mM Tris Acetate pH 7.5, 100 mM Sodium Acetate, 10 mM or 1 mM Magnesium Acetate) supplemented with an oxygen scavenging system (1 mg/mL glucose oxidase, 0.8% glucose, 0.04 mg/mL catalase and 3 mM Trolox). Imaging was initiated as soon as the buffer exchange was complete. All measurements were carried out at room temperature (23 ± 1°C).

Single molecule data acquisition

Excitation of the donor, Cy3, was carried out using a Nd:YAG laser (532 nm, 75 mW; Crystalaser, Reno, NV) by means of prism type total internal reflection microscopy (Roy et al., 2008). After filtering the scattered excitation light using a 550 nm long pass filter, fluorescence emission from the donor and the acceptor was refocused onto an EMCCD camera (Andor, UK). The Cy3 and Cy5 emissions were split into two channels using a 630 nm dichroic mirror. In the case of three color measurements involving Cy3, Cy5 and Cy7, an additional dichroic mirror (730 nm) was used to split the emission across three emission channels. The time resolution for all single molecule strand exchange experiments was 30 ms unless otherwise specified. The data acquisition was carried out using home built software written in Visual C++. The movies obtained with the CCD were analyzed first using IDL and the intensities of the fluorophores and the time traces was visualized using customized MATLAB programs (Joo and Ha, 2008; Roy et al., 2008).

Single molecule data analysis

The dwell time analysis was carried out by a home-written MATLAB program. The background intensity in the donor and acceptor channel was subtracted followed by leakage subtraction of the donor signal to the acceptor channel. Details regarding the acquisition and analysis are based on previously published methods. After visually inspecting the acquired data, we manually selected the relevant time periods for analysis and used Origin 8.0 to plot the data. Cross-correlation analysis was performed as previously described (Kim et al., 2002) by calculating cross-correlation functions for donor and acceptor time traces for a given molecule. Fitting the averaged cross-correlation functions to a single exponential function allows to estimate the average cross correlation time for each measurement. Three color data analysis was performed as previously described (Hohng et al., 2004; Roy et al., 2009). Following leakage correction for Cy3, Cy5 and Cy7 intensities, we additionally corrected Cy7 intensity by utilizing a gamma factor (to correct for differences in optics and detection efficiency).γ=ΔICy5ΔICy7ΔICy5 and ΔICy7 represent the ratio of the change in intensity in Cy5 and Cy7 channels which is calculated by following the change in intensity in the two channels following Cy7 photobleaching. ECy3->Cy5, the FRET efficiency corresponding to Cy3–Cy5 interaction and ECy3->Cy7 corresponding to the Cy3-Cy7 interaction was calculated asED-A,i=IA,iΣIA,i+ID,where ID and IA are the donor and acceptor intensities respectively. The suffix, i, denotes the two acceptors Cy5 and Cy7 utilized in this measurement. We carefully designed the three color measurement to ensure that the two acceptors Cy5 and Cy7 spaced far apart from each other. such that ECy5->Cy7 ∼ 0 and the presence of a stable RecA filament (by using ATPγS) results in the spacing between Cy5 and Cy7 remaining unchanged during the course of the measurement. Software for acquiring and analyzing single molecule FRET data is freely available for download from https://physics.illinois.edu/cplc/software.

Monte Carlo simulation for calibration

Monte Carlo (MC) simulation of the free diffusion along dsDNA was written in MATLAB. Since the persistence length of RecA-ssDNA was reported to be approximately 784 nm (approximately 2300 bp) (Hegner et al., 1999), the RecA-DNA filaments in our experiments can be simply treated as a rod in one dimension. Equal probabilities of moving forward or backward were generated by the program. The step size of the MC simulation was set to be small enough (0.1 or 0.01 bp, in our case) to mimic continuous sliding of the RecA filament along dsDNA. The distance between the donor and acceptor fluorophores and their corresponding FRET value was calculated after each step. Since RecA binding to ssDNA results in an average rise between successive nucleotides of d = 5.1 Å/bp, L could be expressed byL=r2+(d×|xxo|)2,where x and xo respectively represent the position of non-homologous dsDNA and the static acceptor along the dsDNA. Total intensity of each data bin, 30 ms, is set to be 500 a.u., matching our typical data intensity in TIRF experiments. Gaussian white noise with a peak intensity of 40 a.u. is applied. Finally, same cross correlation (CC) analysis as described before was applied to the MC simulated traces for calibration purpose (τ, dwell time of CC analysis; D, diffusion coefficient).

References

Acknowledgements

We dedicate this paper to the memory of our brilliant colleague, Robert MClegg (1945-2012), a pioneer in the application of fluorescence and FRET approaches to the life sciences, a great mentor and a remarkable human being. We thank K Lee, J Yoo and C Joo for discussions and experimental help. We are grateful to A Jain for help with data acquisition. We thank A Jain and R Zhou for careful reading of the manuscript. We also thank A Jain, J Fei, H Koh, Y Ishitsuka, S Doganay and X Shi for advice and comments related to the project.

Decision letter

Xiaowei Zhuang, Reviewing editor, Harvard University, United States

eLife posts the editorial decision letter and author response on a selection of the published articles (subject to the approval of the authors). An edited version of the letter sent to the authors after peer review is shown, indicating the substantive concerns or comments; minor concerns are not usually shown. Reviewers have the opportunity to discuss the decision before the letter is sent (see review process). Similarly, the author response typically shows only responses to the major concerns raised by the reviewers.

Thank you for choosing to send your work entitled “RecA filament sliding on DNA facilitates Homology search” for consideration at eLife. Your article has been evaluated by a Senior Editor and two reviewers, one of whom is a member of eLife's Board of Reviewing Editors. We are pleased to invite you to submit a revised manuscript, which addresses the comments from the reviewers.

Both reviewers find your work of significant interest and high quality. The results represent a significant advance of the current understanding of homology search during recombination and are suitable for publication in eLife. The following two major comments need to be addressed during revision:

1) All data interpretation and analysis are based on the notion that the RecA-ssDNA filament and the duplex DNA act as two stiff, coaxial rods, and that all relative motion involves linear displacement (sliding) of these two rods relative to one another. This interpretation ignores the possibility of off-axis motion (swinging or other fluctuations), which may be relevant to homology search. While the reviewers find the evidence for sliding, especially data shown in Figures 2 and 3, strong, they think that potential contributions from the off-axis motion, and how they may affect data interpretations and simulations, should be discussed. The authors should discuss this effect in the revised manuscript and, if necessary, provide additional data or analysis to support their arguments. The assay shown in Figure 4 may help assess off-axis motion. In this assay, one expects ∼50% of the binding events to occur in the correct orientation that permits recognition of HS1/HS2, and the remaining 50% to occur with the duplex DNA in the incorrect orientation, which may not lead to recognition of HS1/HS2. A statistical analysis of the fractions may provide clues for whether reorientation or off-axis motion between the duplex DNA and RecA-ssDNA filament occurs. The authors may also find a modification of the Figure 4 assay helpful, namely change the sequence of the HS2 site, such that it would be complementary to a sequence on the opposing strand of the incoming duplex, and then test whether recognition still occurs for both HS1 and HS2 sites or only one of them in a single binding event.

2. The authors suggest that the transitions between the HS1 and HS2 sites observed in Figure 4 are due to sliding. Indeed, the scenario that the unwound ssDNA first detach from the HS1 site, re-anneal with the complementary strand, slide along the RecA filament and then bind to the HS2 site is a highly likely one. However, give the relatively short distance between the HS1 and HS2 sites (5 nucleotides in between), is it possible that the ssDNA segment can transition between the two homology sites through a different mechanism? The authors should discuss whether the 5-nucleotide distance is sufficiently long to conclude that transitions between the two sites are due to sliding, and if necessary, provide additional evidence.

DOI: http://dx.doi.org/10.7554/eLife.00067.018

Author response