Highly mismatch-tolerant homology testing by RecA could explain how homology length affects recombination.

In E. coli, double strand breaks (DSBs) are resected and loaded with RecA protein. The genome is then rapidly searched for a sequence that is homologous to the DNA flanking the DSB. Mismatches in homologous partners are rare, suggesting that RecA should rapidly reject mismatched recombination products; however, this is not the case. Decades of work have shown that long lasting recombination products can include many mismatches. In this work, we show that in vitro RecA forms readily observable recombination products when 16% of the bases in the product are mismatched. We also consider various theoretical models of mismatch-tolerant homology testing. The models test homology by comparing the sequences of Ltest bases in two single-stranded DNAs (ssDNA) from the same genome. If the two sequences pass the homology test, the pairing between the two ssDNA becomes permanent. Stringency is the fraction of permanent pairings that join ssDNA from the same positions in the genome. We applied the models to both randomly generated genomes and bacterial genomes. For both randomly generated genomes and bacterial genomes, the models show that if no mismatches are accepted stringency is ∼ 99% when Ltest = 14 bp. For randomly generated genomes, stringency decreases with increasing mismatch tolerance, and stringency improves with increasing Ltest. In contrast, in bacterial genomes when Ltest ∼ 75 bp, stringency is ∼ 99% for both mismatch-intolerant and mismatch-tolerant homology testing. Furthermore, increasing Ltest does not improve stringency because most incorrect pairings join different copies of repeats. In sum, for bacterial genomes highly mismatch tolerant homology testing of 75 bp provides the same stringency as homology testing that rejects all mismatches and testing more than ∼75 base pairs is not useful. Interestingly, in vivo commitment to recombination typically requires homology testing of ∼ 75 bp, consistent with highly mismatch intolerant testing.

Building Biological Relevance Into Integrative Modelling of Macromolecular Assemblies.

Recent advances in structural biophysics and integrative modelling methods now allow us to decipher the structures of large macromolecular assemblies. Understanding the dynamics and mechanisms involved in their biological function requires rigorous integration of all available data. We have developed a complete modelling pipeline that includes analyses to extract biologically significant information by consistently combining automated and interactive human-guided steps. We illustrate this idea with two examples. First, we describe the ryanodine receptor, an ion channel that controls ion flux across the cell membrane through transitions between open and closed states. The conformational changes associated with the transitions are small compared to the considerable system size of the receptor; it is challenging to consistently track these states with the available cryo-EM structures. The second example involves homologous recombination, in which long filaments of a recombinase protein and DNA catalyse the exchange of homologous DNA strands to reliably repair DNA double-strand breaks. The nucleoprotein filament reaction intermediates in this process are short-lived and heterogeneous, making their structures particularly elusive. The pipeline we describe, which incorporates experimental and theoretical knowledge combined with state-of-the-art interactive and immersive modelling tools, can help overcome these challenges. In both examples, we point to new insights into biological processes that arise from such interdisciplinary approaches.

Influences of ssDNA-RecA Filament Length on the Fidelity of Homologous Recombination.

Chromosomal double-strand breaks can be accurately repaired by homologous recombination, but genomic rearrangement can result if the repair joins different copies of a repeated sequence. Rearrangement can be advantageous or fatal. During repair, a broken double-stranded DNA (dsDNA) is digested by the RecBCD complex from the 5' end, leaving a sequence gap that separates two 3' single-stranded DNA (ssDNA) tails. RecA binds to the 3' tails forming helical nucleoprotein filaments.A three-strand intermediate is formed when a RecA-bound ssDNA with L nucleotides invades a homologous region of dsDNA and forms a heteroduplex product with a length ≤ L bp. The homology dependent stability of the heteroduplex determines how rapidly and accurately homologous recombination repairs double-strand breaks. If the heteroduplex is sufficiently sequence matched, repair progresses to irreversible DNA synthesis. Otherwise, the heteroduplex should rapidly reverse. In this work, we present in vitro measurements of the L dependent stability of heteroduplex products formed by filaments with 90 ≤ L ≤ 420 nt, which is within the range observedin vivo. We find that without ATP hydrolysis, products are irreversible when L > 50 nt. In contrast, with ATP hydrolysis when L < 160 nt, products reverse in < 30 seconds; however, with ATP hydrolysis when L ≥ 320 nt, some products reverse in < 30 seconds, while others last thousands of seconds. We consider why these two different filament length regimes show such distinct behaviors. We propose that the experimental results combined with theoretical insights suggest that filaments with 250 â‰² L â‰² 8500 nt optimize DSB repair.

Weaving DNA strands: structural insight on ATP hydrolysis in RecA-induced homologous recombination.

Homologous recombination is a fundamental process in all living organisms that allows the faithful repair of DNA double strand breaks, through the exchange of DNA strands between homologous regions of the genome. Results of three decades of investigation and recent fruitful observations have unveiled key elements of the reaction mechanism, which proceeds along nucleofilaments of recombinase proteins of the RecA family. Yet, one essential aspect of homologous recombination has largely been overlooked when deciphering the mechanism: while ATP is hydrolyzed in large quantity during the process, how exactly hydrolysis influences the DNA strand exchange reaction at the structural level remains to be elucidated. In this study, we build on a previous geometrical approach that studied the RecA filament variability without bound DNA to examine the putative implication of ATP hydrolysis on the structure, position, and interactions of up to three DNA strands within the RecA nucleofilament. Simulation results on modeled intermediates in the ATP cycle bring important clues about how local distortions in the DNA strand geometries resulting from ATP hydrolysis can aid sequence recognition by promoting local melting of already formed DNA heteroduplex and transient reverse strand exchange in a weaving type of mechanism.

Slow extension of the invading DNA strand in a D-loop formed by RecA-mediated homologous recombination may enhance recognition of DNA homology.

DNA recombination resulting from RecA-mediated strand exchange aided by RecBCD proteins often enables accurate repair of DNA double-strand breaks. However, the process of recombinational repair between short DNA regions of accidental similarity can lead to fatal genomic rearrangements. Previous studies have probed how effectively RecA discriminates against interactions involving a short similar sequence that is embedded in otherwise dissimilar sequences but have not yielded fully conclusive results. Here, we present results of in vitro experiments with fluorescent probes strategically located on the interacting DNA fragments used for recombination. Our findings suggest that DNA synthesis increases the stability of the recombination products. Fluorescence measurements can also probe the homology dependence of the extension of invading DNA strands in D-loops formed by RecA-mediated strand exchange. We examined the slow extension of the invading strand in a D-loop by DNA polymerase (Pol) IV and the more rapid extension by DNA polymerase LF-Bsu We found that when DNA Pol IV extends the invading strand in a D-loop formed by RecA-mediated strand exchange, the extension afforded by 82 bp of homology is significantly longer than the extension on 50 bp of homology. In contrast, the extension of the invading strand in D-loops by DNA LF-Bsu Pol is similar for intermediates with ≥50 bp of homology. These results suggest that fatal genomic rearrangements due to the recombination of small regions of accidental homology may be reduced if RecA-mediated strand exchange is immediately followed by DNA synthesis by a slow polymerase.

Residues in the fingers domain of the translesion DNA polymerase DinB enable its unique participation in error-prone double-strand break repair.

The evolutionarily conserved Escherichia coli translesion DNA polymerase IV (DinB) is one of three enzymes that can bypass potentially deadly DNA lesions on the template strand during DNA replication. Remarkably, however, DinB is the only known translesion DNA polymerase active in RecA-mediated strand exchange during error-prone double-strand break repair. In this process, a single-stranded DNA (ssDNA)-RecA nucleoprotein filament invades homologous dsDNA, pairing the ssDNA with the complementary strand in the dsDNA. When exchange reaches the 3' end of the ssDNA, a DNA polymerase can add nucleotides onto the end, using one strand of dsDNA as a template and displacing the other. It is unknown what makes DinB uniquely capable of participating in this reaction. To explore this topic, we performed molecular modeling of DinB's interactions with the RecA filament during strand exchange, identifying key contacts made with residues in the DinB fingers domain. These residues are highly conserved in DinB, but not in other translesion DNA polymerases. Using a novel FRET-based assay, we found that DinB variants with mutations in these conserved residues are less effective at stabilizing RecA-mediated strand exchange than native DinB. Furthermore, these variants are specifically deficient in strand displacement in the absence of RecA filament. We propose that the amino acid patch of highly conserved residues in DinB-like proteins provides a mechanistic explanation for DinB's function in strand exchange and improves our understanding of recombination by providing evidence that RecA plays a role in facilitating DinB's activity during strand exchange.

The positioning of Chi sites allows the RecBCD pathway to suppress some genomic rearrangements.

Bacterial recombinational repair of double-strand breaks often begins with creation of initiating 3' single-stranded DNA (ssDNA) tails on each side of a double-strand break (DSB). Importantly, if the RecBCD pathway is followed, RecBCD creates a gap between the sequences at 3' ends of the initiating strands. The gap flanks the DSB and extends at least to the nearest Chi site on each strand. Once the initiating strands form ssDNA-RecA filaments, each ssDNA-RecA filament searches for homologous double-stranded DNA (dsDNA) to use as a template for the DNA synthesis needed to fill the gap created by RecBCD. Our experimental results show that the DNA synthesis requires formation of a heteroduplex dsDNA that pairs >20 contiguous bases in the initiating strand with sequence matched bases in a strand from the original dsDNA. To trigger synthesis, the heteroduplex must be near the 3' end of the initiating strand. Those experimentally determined requirements for synthesis combined with the Chi site dependence of the function of RecBCD and the distribution of Chi sites in bacterial genomes could allow the RecBCD pathway to avoid some genomic rearrangements arising from directly induced DSBs; however, the same three factors could promote other rearrangements.

ATP hydrolysis provides functions that promote rejection of pairings between different copies of long repeated sequences.

During DNA recombination and repair, RecA family proteins must promote rapid joining of homologous DNA. Repeated sequences with >100 base pair lengths occupy more than 1% of bacterial genomes; however, commitment to strand exchange was believed to occur after testing ∼20-30 bp. If that were true, pairings between different copies of long repeated sequences would usually become irreversible. Our experiments reveal that in the presence of ATP hydrolysis even 75 bp sequence-matched strand exchange products remain quite reversible. Experiments also indicate that when ATP hydrolysis is present, flanking heterologous dsDNA regions increase the reversibility of sequence matched strand exchange products with lengths up to ∼75 bp. Results of molecular dynamics simulations provide insight into how ATP hydrolysis destabilizes strand exchange products. These results inspired a model that shows how pairings between long repeated sequences could be efficiently rejected even though most homologous pairings form irreversible products.

Single molecule identification of homology-dependent interactions between long ssRNA and dsDNA.

Long non-coding RNAs (lncRNAs) are prominently associated with chromosomes in an ever-increasing diversity of roles. To provide further insight into the potential nature of these associations, we have explored, for the first time, the interaction of long single-stranded (ss) RNAs with cognate homologous double-stranded (ds) DNA in vitro Using magnetic tweezers, we measured the effects of ssRNA on force extension curves for dsDNA. We observe that the presence of ssRNA impedes the extension of dsDNA, specifically at low forces, dependent on homology between the RNA and DNA species, and dependent on ssRNA lengths (≥1 kb). The observed effect also depends on the concentration of ssRNA and is abolished by overstretching of the dsDNA. These findings show that significant homologous contacts can occur between long ssRNA and dsDNA in the absence of protein and that these contacts alter the mechanical properties of the dsDNA. We propose that long ssRNA interacts paranemically with long dsDNA via periodic short homologous interactions, e.g. mediated by RNA/DNA triplex-formation, and that dsDNA extension is impeded by formation of RNA secondary structure in the intervening unbound regions. Analogous interactions in vivo would permit lncRNAs to mediate the juxtaposition of two or more DNA regions on the same or different chromosomes.

Structure/function relationships in RecA protein-mediated homology recognition and strand exchange.

RecA family proteins include RecA, Rad51, and Dmc1. These recombinases are responsible for homology search and strand exchange. Homology search and strand exchange occur during double-strand break repair and in eukaryotes during meiotic recombination. In bacteria, homology search begins when RecA binds an initiating single-stranded DNA (ssDNA) in the primary DNA-binding site to form the presynaptic filament. The filament is a right-handed helix, where the initiating strand is bound deep within the filament. Once the presynaptic filament is formed, it interrogates nearby double-stranded DNA (dsDNA) to find a homologous sequence; therefore, we provide a detailed discussion of structural features of the presynaptic filament that play important functional roles. The discussion includes many diagrams showing multiple filament turns. These diagrams illustrate interactions that are not evident in single turn structures. The first dsDNA interactions with the presynaptic filament are insensitive to mismatches. The mismatch insensitive interactions lead to dsDNA deformation that triggers a homology testing process governed by kinetics. The first homology test involves ∼8 bases. Almost all interactions are rejected by this initial rapid test, leading to a new cycle of homology testing. Interactions that pass the initial rapid test proceed to a slower testing stage. That slower stage induces nonhomologous dsDNA to reverse strand exchange and begin a new cycle of homology testing. In contrast, homologous dsDNA continues to extend the heteroduplex strand-exchange product until ATP hydrolysis makes strand exchange irreversible.

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.

ssDNA Pairing Accuracy Increases When Abasic Sites Divide Nucleotides into Small Groups.

Accurate sequence dependent pairing of single-stranded DNA (ssDNA) molecules plays an important role in gene chips, DNA origami, and polymerase chain reactions. In many assays accurate pairing depends on mismatched sequences melting at lower temperatures than matched sequences; however, for sequences longer than ~10 nucleotides, single mismatches and correct matches have melting temperature differences of less than 3°C. We demonstrate that appropriately grouping of 35 bases in ssDNA using abasic sites increases the difference between the melting temperature of correct bases and the melting temperature of mismatched base pairings. Importantly, in the presence of appropriately spaced abasic sites mismatches near one end of a long dsDNA destabilize the annealing at the other end much more effectively than in systems without the abasic sites, suggesting that the dsDNA melts more uniformly in the presence of appropriately spaced abasic sites. In sum, the presence of appropriately spaced abasic sites allows temperature to more accurately discriminate correct base pairings from incorrect ones.

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.

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.

RecA homology search is promoted by mechanical stress along the scanned duplex DNA.

A RecA-single-stranded DNA (RecA-ssDNA) filament searches a genome for sequence homology by rapidly binding and unbinding double-stranded DNA (dsDNA) until homology is found. We demonstrate that pulling on the opposite termini (3' and 5') of one of the two DNA strands in a dsDNA molecule stabilizes the normally unstable binding of that dsDNA to non-homologous RecA-ssDNA filaments, whereas pulling on the two 3', the two 5', or all four termini does not. We propose that the 'outgoing' strand in the dsDNA is extended by strong DNA-protein contacts, whereas the 'complementary' strand is extended by the tension on the base pairs that connect the 'complementary' strand to the 'outgoing' strand. The stress resulting from different levels of tension on its constitutive strands causes rapid dsDNA unbinding unless sufficient homology is present.

Changes in the tension in dsDNA alter the conformation of RecA bound to dsDNA-RecA filaments.

The RecA protein is an ATPase that mediates recombination via strand exchange. In strand exchange a single-stranded DNA (ssDNA) bound to RecA binding site I in a RecA/ssDNA filament pairs with one strand of a double-stranded DNA (dsDNA) and forms heteroduplex dsDNA in site I if homology is encountered. Long sequences are exchanged in a dynamic process in which initially unbound dsDNA binds to the leading end of a RecA/ssDNA filament, while heteroduplex dsDNA unbinds from the lagging end via ATP hydrolysis. ATP hydrolysis is required to convert the active RecA conformation, which cannot unbind, to the inactive conformation, which can unbind. If dsDNA extension due to RecA binding increases the dsDNA tension, then RecA unbinding must decrease tension. We show that in the presence of ATP hydrolysis decreases in tension induce decreases in length whereas in the absence of hydrolysis, changes in tension have no systematic effect. These results suggest that decreases in force enhance dissociation by promoting transitions from the active to the inactive RecA conformation. In contrast, increases in tension reduce dissociation. Thus, the changes in tension inherent to strand exchange may couple with ATP hydrolysis to increase the directionality and stringency of strand exchange.

Single-molecule studies of the stringency factors and rates governing the polymerization of RecA on double-stranded DNA.

RecA is a key protein in homologous recombination. During recombination, one single-stranded DNA (ssDNA) bound to site I in RecA exchanges Watson-Crick pairing with a sequence-matched ssDNA that was part of a double-stranded DNA molecule (dsDNA) bound to site II in RecA. After strand exchange, heteroduplex dsDNA is bound to site I. In vivo, direct polymerization of RecA on dsDNA through site I does not occur, though it does in vitro. The mechanisms underlying the difference have been unclear. We use single-molecule experiments to decouple the two steps involved in polymerization: nucleation and elongation. We find that elongation is governed by a fundamental clock that is insensitive to force and RecA concentration from 0.2 and 6 µM, though rates depend on ionic conditions. Thus, we can probe nucleation site stability by creating nucleation sites at high force and then measuring elongation as a function of applied force. We find that in the presence of ATP hydrolysis a minimum force is required for polymerization. The minimum force decreases with increasing RecA or ATP concentrations. We propose that force reduces the off-rate for nucleation site binding and that nucleation site stability is the stringency factor that prevents in vivo polymerization.

Study of force induced melting of dsDNA as a function of length and conformation.

We measure the constant force required to melt double-stranded (ds) DNA as a function of length for lengths from 12 to 100,000 base pairs, where the force is applied to the 3'3' or 5'5' ends of the dsDNA. Molecules with 32 base pairs or fewer melt before overstretching. For these short molecules, the melting force is independent of the ends to which the force is applied and the shear force as a function of length is well described by de Gennes theory with a de Gennes length of less than 10 bp. Molecules with lengths of 500 base pairs or more overstretch before melting. For these long molecules, the melting force depends on the ends to which the force is applied. The melting force as a function of length increases even when the length exceeds 1000 bp, where the length dependence is inconsistent with de Gennes theory. Finally, we expand de Gennes melting theory to 3'5' pulling and compare the predictions with experimental results.

The structure of DNA overstretched from the 5'5' ends differs from the structure of DNA overstretched from the 3'3' ends.

It has been suggested that the structure that results when double-stranded DNA (dsDNA) is pulled from the 3'3' ends differs from that which results when it is pulled from the 5'5' ends. In this work, we demonstrate, using lambda phage dsDNA, that the overstretched states do indeed show different properties, suggesting that they correspond to different structures. For 3'3' pulling versus 5'5' pulling, the following differences are observed: (i) the forces at which half of the molecules in the ensemble have made a complete force-induced transition to single stranded DNA are 141 +/- 3 pN and 122 +/- 4 pN, respectively; (ii) the extension vs. force curve for overstretched DNA has a marked change in slope at 127 +/- 3 pN for 3'3' and 110 +/- 3 pN for 5'5'; (iii) the hysteresis (H) in the extension vs. force curves at 150 mM NaCl is 0.3 +/- 0.8 pN microm for 3'3' versus 13 +/- 8 pN for 5'5'; and (iv) 3'3' and 5'5' molecules show different changes in hysteresis due to interactions with beta-cyclodextrin, a molecule that is known to form stable host-guest complexes with rotated base pairs, and glyoxal that is known to bind stably to unpaired bases. These differences and additional findings are well-accommodated by the corresponding structures predicted on theoretical grounds.