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.

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.

Finding the infectious dose for COVID-19 by applying an airborne-transmission model to superspreader events.

We probed the transmission of COVID-19 by applying an airborne transmission model to five well-documented case studies-a Washington state church choir, a Korean call center, a Korean exercise class, and two different Chinese bus trips. For all events the likely index patients were pre-symptomatic or mildly symptomatic, which is when infective patients are most likely to interact with large groups of people. Applying the model to those events yields results that suggest the following: (1) transmission was airborne; (2) superspreading events do not require an index patient with an unusually high viral load; (3) the viral loads for all of the index patients were of the same order of magnitude and consistent with experimentally measured values for patients at the onset of symptoms, even though viral loads across the population vary by a factor of >108. In particular we used a Wells-Riley exposure model to calculate q, the total average number of infectious quanta inhaled by a person at the event. Given the q value for each event, the simple airborne transmission model was used to determined Sq, the rate at which the index patient exhaled infectious quanta and N0, the characteristic number of COVID-19 virions needed to induce infection. Despite the uncertainties in the values of some parameters of the superspreading events, all five events yielded (N0∼300-2,000 virions), which is similar to published values for influenza. Finally, this work describes the conditions under which similar methods can provide actionable information on the transmission of other viruses.

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.

Interactions between identical DNA double helices.

The molecular mechanism of specific interactions between double stranded DNA molecules has been investigated for many years. Problems remain in how confinement, ions, and condensing agents change the interactions. We consider how the orientational alignment of DNAs contributes to the interactions via free energy simulations. Here we report on the effective interactions between two parallel DNA double helices in 150-mM NaCl solution using all atom models. We calculate the potential of mean force (PMF) of DNA-DNA interactions as a function of two coordinates, interhelical separation of parallel double helices and relative rotation of a DNA molecule with respect to the other about the helical axis. We generate the two-dimensional PMF to better understand the effective interactions when a DNA molecule is in juxtaposition with another. The analysis of the ion and solvent distributions around the DNA and particularly in the interface region shows that certain alignments of the DNA pair enhance the interactions. At local free energy minima in distance and alignment, water molecules and Na^{+} ions form a hydrogen bonded network with the phosphates from each DNA. This network contributes an attractive energy component to the DNA-DNA interactions. Our results provide a molecular mechanism whereby local DNA-DNA interactions, depending on the helical orientation, give a potential mechanism for stabilizing pairing of much larger lengths of homologous DNA that have been seen experimentally. The study suggests an atomically detailed local picture of relevance to certain aspects of DNA condensation or aggregation.

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.

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.

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.

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.

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.

Chromosomes Progress to Metaphase in Multiple Discrete Steps via Global Compaction/Expansion Cycles.

Mammalian mitotic chromosome morphogenesis was analyzed by 4D live-cell and snapshot deconvolution fluorescence imaging. Prophase chromosomes, whose organization was previously unknown, are revealed to comprise co-oriented sister linear loop arrays displayed along a single, peripheral, regularly kinked topoisomerase II/cohesin/condensin II axis. Thereafter, rather than smooth, progressive compaction as generally envisioned, progression to metaphase is a discontinuous process involving chromosome expansion as well as compaction. At late prophase, dependent on topoisomerase II and with concomitant cohesin release, chromosomes expand, axes split and straighten, and chromatin loops transit to a radial disposition around now-central axes. Finally, chromosomes globally compact, giving the metaphase state. These patterns are consistent with the hypothesis that the molecular events of chromosome morphogenesis are governed by accumulation and release of chromosome stress, created by chromatin compaction and expansion. Chromosome state could evolve analogously throughout the cell cycle.

RecA-mediated sequence homology recognition as an example of how searching speed in self-assembly systems can be optimized by balancing entropic and enthalpic barriers.

Ideally, self-assembly should rapidly and efficiently produce stable correctly assembled structures. We study the tradeoff between enthalpic and entropic cost in self-assembling systems using RecA-mediated homology search as an example. Earlier work suggested that RecA searches could produce stable final structures with high stringency using a slow testing process that follows an initial rapid search of ∼9-15 bases. In this work, we will show that as a result of entropic and enthalpic barriers, simultaneously testing all ∼9-15 bases as separate individual units results in a longer overall searching time than testing them in groups and stages.