Tension on dsDNA bound to ssDNA-RecA filaments may play an important role in driving efficient and accurate homology recognition and strand exchange

Modeling the Homologous Recombination Process: Methods, Successes and Challenges

Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.

Mechanisms of distinctive mismatch tolerance between Rad51 and Dmc1 in homologous recombination

Homologous recombination (HR) is a primary DNA double-strand breaks (DSBs) repair mechanism. The recombinases Rad51 and Dmc1 are highly conserved in the RecA family; Rad51 is mainly responsible for DNA repair in somatic cells during mitosis while Dmc1 only works during meiosis in germ cells. This spatiotemporal difference is probably due to their distinctive mismatch tolerance during HR: Rad51 does not permit HR in the presence of mismatches, whereas Dmc1 can tolerate certain mismatches. Here, the cryo-EM structures of Rad51-DNA and Dmc1-DNA complexes revealed that the major conformational differences between these two proteins are located in their Loop2 regions, which contain invading single-stranded DNA (ssDNA) binding residues and double-stranded DNA (dsDNA) complementary strand binding residues, stabilizing ssDNA and dsDNA in presynaptic and postsynaptic complexes, respectively. By combining molecular dynamic simulation and single-molecule FRET assays, we identified that V273 and D274 in the Loop2 region of human RAD51 (hRAD51), corresponding to P274 and G275 of human DMC1 (hDMC1), are the key residues regulating mismatch tolerance during strand exchange in HR. This HR accuracy control mechanism provides mechanistic insights into the specific roles of Rad51 and Dmc1 in DNA double-strand break repair and may shed light on the regulatory mechanism of genetic recombination in mitosis and meiosis.

How strand exchange protein function benefits from ATP hydrolysis

Members of the RecA family of strand exchange proteins carry out the central reaction in homologous recombination. These proteins are DNA-dependent ATPases, although their ATPase activity is not required for the key functions of homology search and strand exchange. We review the literature on the role of the intrinsic ATPase activity of strand exchange proteins. We also discuss the role of ATP-hydrolysis-dependent motor proteins that serve as strand exchange accessory factors, with an emphasis on the eukaryotic Rad54 family of double strand DNA-specific translocases. The energy from ATP allows recombination events to progress from the strand exchange stage to subsequent stages. ATP hydrolysis also functions to corrects DNA binding errors, including particularly detrimental binding to double strand DNA.

Mismatch sensing by nucleofilament deciphers mechanism of RecA-mediated homologous recombination

Recombinases polymerize along single-stranded DNA (ssDNA) at the end of a broken DNA to form a helical nucleofilament with a periodicity of ∼18 bases. The filament catalyzes the search and checking for homologous sequences and promotes strand exchange with a donor duplex during homologous recombination (HR), the mechanism of which has remained mysterious since its discovery. Here, by inserting mismatched segments into donor duplexes and using single-molecule techniques to catch transient intermediates in HR, we found that, even though 3 base pairs (bp) is still the basic unit, both the homology checking and the strand exchange may proceed in multiple steps at a time, resulting in ∼9-bp large steps on average. More interestingly, the strand exchange is blocked remotely by the mismatched segment, terminating at positions ∼9 bp before the match-mismatch joint. The homology checking and the strand exchange are thus separated in space, with the strand exchange lagging behind. Our data suggest that the strand exchange progresses like a traveling wave in which the donor DNA is incorporated successively into the ssDNA-RecA filament to check homology in ∼9-bp steps in the frontier, followed by a hypothetical transitional segment and then the post-strand-exchanged duplex.

Mechanisms of fast and stringent search in homologous pairing of double-stranded DNA

Self-organization in the cell relies on the rapid and specific binding of molecules to their cognate targets. Correct bindings must be stable enough to promote the desired function even in the crowded and fluctuating cellular environment. In systems with many nearly matched targets, rapid and stringent formation of stable products is challenging. Mechanisms that overcome this challenge have been previously proposed, including separating the process into multiple stages; however, how particular in vivo systems overcome the challenge remains unclear. Here we consider a kinetic system, inspired by homology dependent pairing between double stranded DNA in bacteria. By considering a simplified tractable model, we identify different homology testing stages that naturally occur in the system. In particular, we first model dsDNA molecules as short rigid rods containing periodically spaced binding sites. The interaction begins when the centers of two rods collide at a random angle. For most collision angles, the interaction energy is weak because only a few binding sites near the collision point contribute significantly to the binding energy. We show that most incorrect pairings are rapidly rejected at this stage. In rare cases, the two rods enter a second stage by rotating into parallel alignment. While rotation increases the stability of matched and nearly matched pairings, subsequent rotational fluctuations reduce kinetic trapping. Finally, in vivo chromosome are much longer than the persistence length of dsDNA, so we extended the model to include multiple parallel collisions between long dsDNA molecules, and find that those additional interactions can greatly accelerate the searching.

Protein folding and the binding of sequence dependent proteins to DNA are examples of self-assembling systems in which the binding energy varies continuously throughout the interaction. Previous theoretical work has highlighted the importance of dividing the interaction into separate stages characterized by interaction times and binding energies that vary by orders of magnitude. Insight into how such a division might naturally arise and promote accurate and efficient self-assembly is provided by our study of a simple tractable model inspired by the homology dependent pairing of double stranded DNA molecules in vivo. In the model, the binding energy is controlled by one single continuously tunable variable whose natural evolution creates stages that efficiently and accurately form stable products.

DNA with Different Local Torsional States Affects RecA-Mediated Recombination Progression

DNA topology is thought to affect DNA enzyme activity. The helical structure of duplex DNA dictates the change of topological states during strand separation when DNA is constrained. During the repair of DNA double-stranded breaks, the RecA nucleoprotein filament invades DNA and carries out consecutive strand exchange reactions coupled with duplex DNA strand separation. It has been suggested that torsional strain could be generated and its accumulation could inhibit strand exchange. We used hairpin and nicked DNA substrates to test how torsional strain alters the RecA-mediated strand exchange efficiency. Single-molecule tethered particle motion (TPM) experiments showed that torsionally constrained hairpin DNA substrates returned nearly no successful strand exchange events catalyzed by RecA. Surprisingly, the strand exchange efficiencies increase in the presence of DNA nicks or loop disruption. The dwell time of transient RecA events in hairpin is shorter compared to those found in nicked or fork DNA substrates, which suggests a limited strand exchange progression in hairpin substrates. Our observation shows that RecA generates local torsional strain during strand exchange, and the inability to dissipate this torsional strain inhibits homologous recombination progression. DNA topological states are thus important regulation measures of DNA recombination.

Enhancement of RecA-mediated self-assembly in DNA nanostructures through basepair mismatches and single-strand nicks

The use of DNA as a structural material for nanometre-scale construction has grown extensively over the last decades. The development of more advanced DNA-based materials would benefit from a modular approach enabling the direct assembly of additional elements onto nanostructures after fabrication. RecA-based nucleoprotein filaments encapsulating short ssDNA have been demonstrated as a tool for highly efficient and fully programmable post-hoc patterning of duplex DNA scaffold. However, the underlying assembly process is not fully understood, in particular when patterning complex DNA topologies. Here, we report the effect of basepair-mismatched regions and single-strand nicks in the double-stranded DNA scaffold on the yield of RecA-based assembly. Significant increases in assembly yield are observed upon the introduction of unpaired basepairs directly adjacent to the assembly region. However, when the unpaired regions were introduced further from the assembly site the assembly yield initially decreased as the length of the unpaired region was increased. These results suggest that an unpaired region acts as a kinetic trap for RecA-based nucleoprotein filaments, impeding the assembly mechanism. Conversely, when the unpaired region is located directly adjacent to the assembly site, it leads to an increase in efficiency of RecA patterning owing to increased breathing of the assembly site.

Evidence of protein-free homology recognition in magnetic bead force-extension experiments

Earlier theoretical studies have proposed that the homology-dependent pairing of large tracts of dsDNA may be due to physical interactions between homologous regions. Such interactions could contribute to the sequence-dependent pairing of chromosome regions that may occur in the presence or the absence of double-strand breaks. Several experiments have indicated the recognition of homologous sequences in pure electrolytic solutions without proteins. Here, we report single-molecule force experiments with a designed 60 kb long dsDNA construct; one end attached to a solid surface and the other end to a magnetic bead. The 60 kb constructs contain two 10 kb long homologous tracts oriented head to head, so that their sequences match if the two tracts fold on each other. The distance between the bead and the surface is measured as a function of the force applied to the bead. At low forces, the construct molecules extend substantially less than normal, control dsDNA, indicating the existence of preferential interaction between the homologous regions. The force increase causes no abrupt but continuous unfolding of the paired homologous regions. Simple semi-phenomenological models of the unfolding mechanics are proposed, and their predictions are compared with the data.

Binding energies of nucleobase complexes: Relevance to homology recognition of DNA

The binding energies of complexes of DNA nucleobase pairs are evaluated using quantum mechanical calculations at the level of dispersion corrected density functional theory. We begin with Watson-Crick base pairs of singlets, duplets, and triplets and calculate their binding energies. At a second step, mismatches are incorporated into the Watson-Crick complexes in order to evaluate the variation in the binding energy with respect to the canonical Watson-Crick pairs. A linear variation of this binding energy with the degree of mismatching is observed. The binding energies for the duplets and triplets containing mismatches are further compared to the energies of the respective singlets in order to assess the degree of collectivity in these complexes. This study also suggests that mismatches do not considerably affect the energetics of canonical base pairs. Our work is highly relevant to the recognition process in DNA promoted through the RecA protein and suggests a clear distinction between recognition in singlets, and recognition in duplets or triplets. Our work assesses the importance of collectivity in the homology recognition of DNA.

DNA Sequence Alignment during Homologous Recombination

Homologous recombination allows for the regulated exchange of genetic information between two different DNA molecules of identical or nearly identical sequence composition, and is a major pathway for the repair of double-stranded DNA breaks. A key facet of homologous recombination is the ability of recombination proteins to perfectly align the damaged DNA with homologous sequence located elsewhere in the genome. This reaction is referred to as the homology search and is akin to the target searches conducted by many different DNA-binding proteins. Here I briefly highlight early investigations into the homology search mechanism, and then describe more recent research. Based on these studies, I summarize a model that includes a combination of intersegmental transfer, short-distance one-dimensional sliding, and length-specific microhomology recognition to efficiently align DNA sequences during the homology search. I also suggest some future directions to help further our understanding of the homology search. Where appropriate, I direct the reader to other recent reviews describing various issues related to homologous recombination.

Integrating multi-scale data on homologous recombination into a new recognition mechanism based on simulations of the RecA-ssDNA/dsDNA structure

RecA protein is the prototypical recombinase. Members of the recombinase family can accurately repair double strand breaks in DNA. They also provide crucial links between pairs of sister chromatids in eukaryotic meiosis. A very broad outline of how these proteins align homologous sequences and promote DNA strand exchange has long been known, as are the crystal structures of the RecA-DNA pre- and postsynaptic complexes; however, little is known about the homology searching conformations and the details of how DNA in bacterial genomes is rapidly searched until homologous alignment is achieved. By integrating a physical model of recognition to new modeling work based on docking exploration and molecular dynamics simulation, we present a detailed structure/function model of homology recognition that reconciles extremely quick searching with the efficient and stringent formation of stable strand exchange products and which is consistent with a vast body of previously unexplained experimental results.

The poor homology stringency in the heteroduplex allows strand exchange to incorporate desirable mismatches without sacrificing recognition in vivo

RecA family proteins are responsible for homology search and strand exchange. In bacteria, homology search begins after RecA binds an initiating single-stranded DNA (ssDNA) in the primary DNA-binding site, forming the presynaptic filament. Once the filament is formed, it interrogates double-stranded DNA (dsDNA). During the interrogation, bases in the dsDNA attempt to form Watson–Crick bonds with the corresponding bases in the initiating strand. Mismatch dependent instability in the base pairing in the heteroduplex strand exchange product could provide stringent recognition; however, we present experimental and theoretical results suggesting that the heteroduplex stability is insensitive to mismatches. We also present data suggesting that an initial homology test of 8 contiguous bases rejects most interactions containing more than 1/8 mismatches without forming a detectable 20 bp product. We propose that, in vivo, the sparsity of accidental sequence matches allows an initial 8 bp test to rapidly reject almost all non-homologous sequences. We speculate that once the initial test is passed, the mismatch insensitive binding in the heteroduplex allows short mismatched regions to be incorporated in otherwise homologous strand exchange products even though sequences with less homology are eventually rejected.

The differential extension in dsDNA bound to Rad51 filaments may play important roles in homology recognition and strand exchange

RecA and Rad51 proteins play an important role in DNA repair and homologous recombination. For RecA, X-ray structure information and single molecule force experiments have indicated that the differential extension between the complementary strand and its Watson–Crick pairing partners promotes the rapid unbinding of non-homologous dsDNA and drives strand exchange forward for homologous dsDNA. In this work we find that both effects are also present in Rad51 protein. In particular, pulling on the opposite termini (3′ and 5′) of one of the two DNA strands in a dsDNA molecule allows dsDNA to extend along non-homologous Rad51-ssDNA filaments and remain stably bound in the extended state, but pulling on the 3′5′ ends of the complementary strand reduces the strand-exchange rate for homologous filaments. Thus, the results suggest that differential extension is also present in dsDNA bound to Rad51. The differential extension promotes rapid recognition by driving the swift unbinding of dsDNA from non-homologous Rad51-ssDNA filaments, while at the same time, reducing base pair tension due to the transfer of the Watson–Crick pairing of the complementary strand bases from the highly extended outgoing strand to the slightly less extended incoming strand, which drives strand exchange forward.

Simplified biased random walk model for RecA-protein-mediated homology recognition offers rapid and accurate self-assembly of long linear arrays of binding sites

Inspired by RecA-protein-based homology recognition, we consider the pairing of two long linear arrays of binding sites. We propose a fully reversible, physically realizable biased random walk model for rapid and accurate self-assembly due to the spontaneous pairing of matching binding sites, where the statistics of the searched sample are included. In the model, there are two bound conformations, and the free energy for each conformation is a weakly nonlinear function of the number of contiguous matched bound sites.

Complementary strand relocation may play vital roles in RecA-based homology recognition

RecA-family proteins mediate homologous recombination and recombinational DNA repair through homology search and strand exchange. Initially, the protein forms a filament with the incoming single-stranded DNA (ssDNA) bound in site I. The RecA–ssDNA filament then binds double-stranded DNA (dsDNA) in site II. Non-homologous dsDNA rapidly unbinds, whereas homologous dsDNA undergoes strand exchange yielding heteroduplex dsDNA in site I and the leftover outgoing strand in site II. We show that applying force to the ends of the complementary strand significantly retards strand exchange, whereas applying the same force to the outgoing strand does not. We also show that crystallographically determined binding site locations require an intermediate structure in addition to the initial and final structures. Furthermore, we demonstrate that the characteristic dsDNA extension rates due to strand exchange and free RecA binding are the same, suggesting that relocation of the complementary strand from its position in the intermediate structure to its position in the final structure limits both rates. Finally, we propose that homology recognition is governed by transitions to and from the intermediate structure, where the transitions depend on differential extension in the dsDNA. This differential extension drives strand exchange forward for homologs and increases the free energy penalty for strand exchange of non-homologs.