Interaction of Rad51 with ATP and Mg2+ induces a conformational change in Rad51

Evolution of meiotic recombination genes in maize and teosinte

Background

Meiotic recombination is a major source of genetic variation in eukaryotes. The role of recombination in evolution is recognized but little is known about how evolutionary forces affect the recombination pathway itself. Although the recombination pathway is fundamentally conserved across different species, genetic variation in recombination components and outcomes has been observed. Theoretical predictions and empirical studies suggest that changes in the recombination pathway are likely to provide adaptive abilities to populations experiencing directional or strong selection pressures, such as those occurring during species domestication. We hypothesized that adaptive changes in recombination may be associated with adaptive evolution patterns of genes involved in meiotic recombination.

Results

To examine how maize evolution and domestication affected meiotic recombination genes, we studied patterns of sequence polymorphism and divergence in eleven genes controlling key steps in the meiotic recombination pathway in a diverse set of maize inbred lines and several accessions of teosinte, the wild ancestor of maize. We discovered that, even though the recombination genes generally exhibited high sequence conservation expected in a pathway controlling a key cellular process, they showed substantial levels and diverse patterns of sequence polymorphism. Among others, we found differences in sequence polymorphism patterns between tropical and temperate maize germplasms. Several recombination genes displayed patterns of polymorphism indicative of adaptive evolution.

Conclusions

Despite their ancient origin and overall sequence conservation, meiotic recombination genes can exhibit extensive and complex patterns of molecular evolution. Changes in these genes could affect the functioning of the recombination pathway, and may have contributed to the successful domestication of maize and its expansion to new cultivation areas.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3486-z) contains supplementary material, which is available to authorized users.

Tyrosine phosphorylation stimulates activity of human RAD51 recombinase through altered nucleoprotein filament dynamics

The DNA strand exchange protein RAD51 facilitates the central step in homologous recombination, a process fundamentally important for accurate repair of damaged chromosomes, restart of collapsed replication forks, and telomere maintenance. The active form of RAD51 is a nucleoprotein filament that assembles on single-stranded DNA (ssDNA) at the sites of DNA damage. The c-Abl tyrosine kinase and its oncogenic counterpart BCR-ABL fusion kinase phosphorylate human RAD51 on tyrosine residues 54 and 315. We combined biochemical reconstitutions of the DNA strand exchange reactions with total internal reflection fluorescence microscopy to determine how the two phosphorylation events affect the biochemical activities of human RAD51 and properties of the RAD51 nucleoprotein filament. By mimicking RAD51 tyrosine phosphorylation with a nonnatural amino acid, p-carboxymethyl-l-phenylalanine (pCMF), we demonstrated that Y54 phosphorylation enhances the RAD51 recombinase activity by at least two different mechanisms, modifies the RAD51 nucleoprotein filament formation, and allows RAD51 to compete efficiently with ssDNA binding protein RPA. In contrast, Y315 phosphorylation has little effect on the RAD51 activities. Based on our work and previous cellular studies, we propose a mechanism underlying RAD51 activation by c-Abl/BCR-ABL kinases.

Without Binding ATP, Human Rad51 Does Not Form Helical Filaments on ssDNA

Construction of the presynaptic filament (PSF) of proper helical structure by Rad51 recombinases is a prerequisite of the progress of homologous recombination repair. We studied the contribution of ATP-binding to this structure of wt human Rad51 (hRad51). We exploited the protein-dissociation effect of high hydrostatic pressure to determine the free energy of dissociation of the protomer interfaces in hRad51 oligomer states and used electron microscopy to obtain topological parameters. Without cofactors ATP and Ca(2+) and template DNA, hRad51 did not exist in monomer form, but it formed rodlike long filaments without helical order. ΔG(diss) indicated a strong inherent tendency of aggregation. Binding solely ssDNA left the filament unstructured with slightly increased ΔG(diss). Adding only ATP and Ca(2+) to the buffer disintegrated the self-associated rods into rings and short helices of further increased ΔG(diss). Rad51 binding to ssDNA only with ATP and Ca bound could lead to ordered helical filament formation of proper pitch size with interface contacts of K(d) ∼ 2 × 10(-11) M, indicating a structure of outstanding stability. ATP/Ca binding increased the ΔG(diss) of protomer contacts in the filament by 16 kJ/mol. The results emphasize that ATP-binding in the PSF of hRad51 has an essential, yet purely structural, role.

Significance of ligand interactions involving Hop2-Mnd1 and the RAD51 and DMC1 recombinases in homologous DNA repair and XX ovarian dysgenesis

The evolutionarily conserved Hop2-Mnd1 complex is a key cofactor for the meiosis-specific recombinase Dmc1. However, emerging evidence has revealed that Hop2-Mnd1 is expressed in somatic tissues, primary human fibroblasts and cell lines, and that it functions in conjunction with the Rad51 recombinase to repair damaged telomeres via the alternate lengthening of telomeres mechanism. Here, we reveal how distinct DNA-binding activities of Hop2-Mnd1 mediate the stabilization of the RAD51-ssDNA presynaptic filament or stimulate the homologous DNA pairing reaction. We have also endeavored to define the interface that governs the assembly of the higher order complex of Hop2-Mnd1 with RAD51. Unexpectedly, we find that ATP enhances the interaction between Hop2-Mnd1 and RAD51, and that both Hop2 and Mnd1 are involved in RAD51 interaction via their C-terminal regions. Importantly, mutations introduced into these Hop2 and Mnd1 domains, including the HOP2 p.del201Glu mutation present in a patient of XX ovarian dysgenesis, diminish the association and functional synergy of Hop2-Mnd1 with both RAD51 and DMC1. Our findings help delineate the intricate manner in which Hop2-Mnd1 engages and functions with RAD51 and DMC1 in mammalian cells and speak to the possible cause of XX ovarian dysgenesis.

Rad51 protein controls Rad52-mediated DNA annealing

In Saccharomyces cerevisiae, Rad52 protein plays an essential role in the repair of DNA double-stranded breaks (DSBs). Rad52 and its orthologs possess the unique capacity to anneal single-stranded DNA (ssDNA) complexed with its cognate ssDNA-binding protein, RPA. This annealing activity is used in multiple mechanisms of DSB repair: single-stranded annealing, synthesis-dependent strand annealing, and cross-over formation. Here we report that the S. cerevisiae DNA strand exchange protein, Rad51, prevents Rad52-mediated annealing of complementary ssDNA. Efficient inhibition is ATP-dependent and involves a specific interaction between Rad51 and Rad52. Free Rad51 can limit DNA annealing by Rad52, but the Rad51 nucleoprotein filament is even more effective. We also discovered that the budding yeast Rad52 paralog, Rad59 protein, partially restores Rad52-dependent DNA annealing in the presence of Rad51, suggesting that Rad52 and Rad59 function coordinately to enhance recombinational DNA repair either by directing the processed DSBs to repair by DNA strand annealing or by promoting second end capture to form a double Holliday junction. This regulation of Rad52-mediated annealing suggests a control function for Rad51 in deciding the recombination path taken for a processed DNA break; the ssDNA can be directed to either Rad51-mediated DNA strand invasion or to Rad52-mediated DNA annealing. This channeling determines the nature of the subsequent repair process and is consistent with the observed competition between these pathways in vivo.

Real-time measurements of the nucleation, growth and dissociation of single Rad51-DNA nucleoprotein filaments

Human Rad51 (hRad51), the protein central to DNA pairing and strand exchange during homologous recombination, polymerizes on DNA to form nucleoprotein filaments. By making use of magnetic tweezers to manipulate individual DNA molecules, we measured the nucleation and growth of hRad51 nucleoprotein filaments, and their subsequent disassembly in real time. The dependence of the initial polymerization rate upon the concentration of hRad51 suggests that the rate-limiting step is the formation of a nucleus involving 5.5 +/- 1.5 hRad51 monomers, corresponding to one helical turn of the hRad51 nucleoprotein filament. Polymerization is highly cooperative (i.e. a nucleation-limited reaction) at low concentrations and less cooperative (a growth-limited reaction) at high concentrations of the protein. We show that the observed preference of hRad51 to form nucleoprotein filaments on double-stranded DNA rather than on single-stranded DNA is due to the fact that it depolymerizes much faster from ssDNA than from dsDNA: indeed, hRad51 polymerizes faster on ssDNA than on dsDNA. Hydrolysis of ATP by hRad51 does not correlate with its dissociation from dsDNA. This suggests that hRad51 does not depolymerize rapidly from dsDNA after strand exchange but stays bound to the heteroduplex, highlighting the importance of partner proteins to facilitate hRad51 depolymerization from dsDNA.

Magnesium influences the discrimination and release of ADP by human RAD51

hRAD51 lacks cooperative DNA-dependent ATPase activity and appears to function with 5-10-fold less Mg2+ compared to RecA. We have further explored the effect of Mg2+ on adenosine nucleotide binding, ATPase, and DNA strand exchange activities. hRAD51 was saturated with the poorly hydrolyzable analog of ATP, ATPgammaS, at approximately 0.08 mM Mg2+. In contrast, > 0.5 mM Mg2+ was required to saturate hRAD51 with ADP. We found ADP to be a significantly less effective competitive inhibitor of the hRAD51 ATPase at low Mg2+ concentrations (0.08 mM). Mg2+ did not appear to affect the ability of ATPgammaS to competitively inhibit the hRAD51 ATPase. Low Mg2+ (0.08-0.12 mM) enhanced the steady-state ATPase of hRAD51 while higher Mg2+ concentration (> 0.3 mM) was inhibitory. At low Mg2+, hRAD51 appeared capable of nearly complete hydrolysis of available ATP, suggesting a lack of ADP product inhibition. There was a strong correlation between the amount of Mg2+ required for stable ADP binding and the inhibition of hRad51 strand exchange activity. Simultaneous inclusion of exogenous ATP and chelation of Mg2+ with EDTA significantly enhanced ADP-->ATP exchange by hRAD51. These studies are consistent with the hypothesis that Mg2+ influences the discrimination and release of ADP, which may sequentially impose an important regulatory step in the hRAD51 ATPase cycle.

The ATPase motif in RAD51D is required for resistance to DNA interstrand crosslinking agents and interaction with RAD51C

Homologous recombination (HR) is a mechanism for repairing DNA interstrand crosslinks and double-strand breaks. In mammals, HR requires the activities of the RAD51 family (RAD51, RAD51B, RAD51C, RAD51D, XRCC2, XRCC3 and DMC1), each of which contains conserved ATP binding sequences (Walker Motifs A and B). RAD51D is a DNA-stimulated ATPase that interacts directly with RAD51C and XRCC2. To test the hypothesis that ATP binding and hydrolysis by RAD51D are required for the repair of interstrand crosslinks, site-directed mutations in Walker Motif A were generated, and complementation studies were performed in Rad51d-deficient mouse embryonic fibroblasts. The K113R and K113A mutants demonstrated a respective 96 and 83% decrease in repair capacity relative to wild-type. Further examination of these mutants, by yeast two-hybrid analyses, revealed an 8-fold reduction in the ability to associate with RAD51C whereas interaction with XRCC2 was retained at a level similar to the S111T control. These cell-based studies are the first evidence that ATP binding and hydrolysis by RAD51D are required for efficient HR repair of DNA interstrand crosslinks.

Rad51 protein from the thermotolerant yeast Pichia angusta as a typical but thermodependent member of the Rad51 family

The Rad51 protein from the methylotrophic yeast Pichia angusta (Rad51(Pa)) of the taxonomic complex Hansenula polymorpha is a homolog of the RecA-RadA-Rad51 protein superfamily, which promotes homologous recombination and recombination repair in prokaryotes and eukaryotes. We cloned the RAD51 gene from the cDNA library of the thermotolerant P. angusta strain BKM Y1397. Induction of this gene in a rad51-deficient Saccharomyces cerevisiae strain partially complemented the survival rate after ionizing radiation. Purified Rad51(Pa) protein exhibited properties typical of the superfamily, including the stoichiometry of binding to single-stranded DNA (ssDNA) (one protomer of Rad51(Pa) per 3 nucleotides) and DNA specificity for ssDNA-dependent ATP hydrolysis [poly(dC) > poly(dT) > phiX174 ssDNA > poly(dA) > double-stranded M13 DNA]. An inefficient ATPase and very low cooperativity for ATP interaction position Rad51(Pa) closer to Rad51 than to RecA. Judging by thermoinactivation, Rad51(Pa) alone was 20-fold more thermostable at 37 degrees C than its S. cerevisiae homolog (Rad51(Sc)). Moreover, it maintained ssDNA-dependent ATPase and DNA transferase activities up to 52 to 54 degrees C, whereas Rad51(Sc) was completely inactive at 47 degrees C. A quick nucleation and an efficient final-product formation in the strand exchange reaction promoted by Rad51(Pa) occurred only at temperatures above 42 degrees C. These reaction characteristics suggest that Rad51(Pa) is dependent on high temperatures for activity.

Biochemical characterization of the human RAD51 protein. III. Modulation of DNA binding by adenosine nucleotides

Adenosine nucleotides affect the ability of RecA small middle dotsingle-stranded DNA (ssDNA) nucleoprotein filaments to cooperatively assume and maintain an extended structure that facilitates DNA pairing during recombination. Here we have determined that ADP and ATP/ATPgammaS affect the DNA binding and aggregation properties of the human RecA homolog human RAD51 protein (hRAD51). These studies have revealed significant differences between hRAD51 and RecA. In the presence of ATPgammaS, RecA forms a stable complex with ssDNA, while the hRAD51 ssDNA complex is destabilized. Conversely, in the presence of ADP and ATP, the RecA ssDNA complex is unstable, while the hRAD51 ssDNA complex is stabilized. We identified two hRAD51 small middle dotssDNA binding forms by gel shift analysis, which were distinct from a well defined RecA small middle dotssDNA binding form. The available evidence suggests that a low molecular weight hRAD51 small middle dotssDNA binding form (hRAD51 small middle dotssDNA(low)) correlates with active ADP and ATP processing. A high molecular weight hRAD51 small middle dotssDNA aggregate (hRAD51 small middle dotssDNA(high)) appears to correlate with a form that fails to process ADP and ATP. Our data are consistent with the notion that hRAD51 is unable to appropriately coordinate ssDNA binding with adenosine nucleotide processing. These observations suggest that other factors may assist hRAD51 in order to mirror RecA recombinational function.

Biochemical characterization of the human RAD51 protein. II. Adenosine nucleotide binding and competition

RecA mediated homologous recombination requires cooperative ATP binding and hydrolysis to assume and maintain an active, extended DNA-protein (nucleoprotein) filament. Human RAD51 protein (hRAD51) lacks the magnitude of ATP-induced cooperativity and catalytic efficiency displayed by RecA. Here, we examined hRAD51 binding and ATPase inhibition pattern by ADP and ATP/adenosine 5'-O-(thiotriphosphate) (ATPgammaS). hRAD51 fully saturates with ATP/ATPgammaS regardless of DNA cofactor (K(D) approximately 5 microm; 1 ATP/1 hRAD51). The binding of ADP to hRAD51 appeared bimodal. The first mode was identical to ATP/ATPgammaS binding (K(app1) approximately 3 microm; 1 ADP/1 hRAD51), while a second mode occurred at elevated ADP concentrations (K(app2) > or = 125 microm; >1 ADP/1 hRAD51). We could detect ADP --> ATP exchange in the high affinity ADP binding mode (K(app1)) but not the low affinity binding mode (K(app2)). At low ATP concentrations (<0.3 mm), ADP and ATPgammaS competitively inhibit the hRAD51 ATPase (K(m)((app)) > K(m)). However, at high ATP (>0.3 mm), the hRAD51 ATPase was stimulated by concentrations of ATPgammaS that were 20-fold above the K(D). Ammonium sulfate plus spermidine decreased the affinity of hRAD51 for ADP substantially ( approximately 10-fold) and ATP modestly ( approximately 3-fold). Our results suggest that ATP binding is not rate-limiting but that the inability to sustain an active nucleoprotein filament probably restricts the hRAD51 ATPase.

Domain structure and dynamics in the helical filaments formed by RecA and Rad51 on DNA

Both the bacterial RecA protein and the eukaryotic Rad51 protein form helical nucleoprotein filaments on DNA that catalyze strand transfer between two homologous DNA molecules. However, only the ATP-binding cores of these proteins have been conserved, and this same core is also found within helicases and the F1-ATPase. The C-terminal domain of the RecA protein forms lobes within the helical RecA filament. However, the Rad51 proteins do not have the C-terminal domain found in RecA, but have an N-terminal extension that is absent in the RecA protein. Both the RecA C-terminal domain and the Rad51 N-terminal domain bind DNA. We have used electron microscopy to show that the lobes of the yeast and human Rad51 filaments appear to be formed by N-terminal domains. These lobes are conformationally flexible in both RecA and Rad51. Within RecA filaments, the change between the "active" and "inactive" states appears to mainly involve a large movement of the C-terminal lobe. The N-terminal domain of Rad51 and the C-terminal domain of RecA may have arisen from convergent evolution to play similar roles in the filaments.

Role of interaction between two silkworm RecA homologs in homologous DNA pairing

Recombinant BmRad51 and BmDmc1, silkworm homologs of the Escherichia coli RecA proteins catalyzing the homologous DNA pairing, were purified from E. coli cells carrying expression vectors. These possessed different enzymatic properties in the joint molecule formation between single-stranded circular DNA and homologous linear double-stranded DNA. The requirement of single-stranded circular DNA for the efficient reaction was twofold higher in BmRad51 than in BmDmc1. Although able to mediate the joint molecule formation independently, a complex of the two enzymes formed prior to single-stranded DNA binding was found to have augmented efficiency of the pairing reaction.

Homologous recombinational repair of DNA ensures mammalian chromosome stability

The process of homologous recombinational repair (HRR) is a major DNA repair pathway that acts on double-strand breaks and interstrand crosslinks, and probably to a lesser extent on other kinds of DNA damage. HRR provides a mechanism for the error-free removal of damage present in DNA that has replicated (S and G2 phases). Thus, HRR acts in a critical way, in coordination with the S and G2 checkpoint machinery, to eliminate chromosomal breaks before the cell division occurs. Many of the human HRR genes, including five Rad51 paralogs, have been identified, and knockout mutants for most of these genes are available in chicken DT40 cells. In the mouse, most of the knockout mutations cause embryonic lethality. The Brca1 and Brca2 breast cancer susceptibility genes appear to be intimately involved in HRR, but the mechanistic basis is unknown. Biochemical studies with purified proteins and cell extracts, combined with cytological studies of nuclear foci, have begun to establish an outline of the steps in mammalian HRR. This pathway is subject to complex regulatory controls from the checkpoint machinery and other processes, and there is increasing evidence that loss of HRR gene function can contribute to tumor development. This review article is meant to be an update of our previous review [Biochimie 81 (1999) 87].

Role of heteroduplex joints in the functional interactions between human Rad51 and wild-type p53

Our previous work (Dudenhöffer et al., 1999) unveiled a link between the capacity of p53 to regulate homologous recombination processes and to specifically bind to heteroduplex junction DNAs. Here, we show that p53 participates in ternary complex formation after preassembly of nucleoproteins, consisting of the human recombinase hRad51 and junction DNA. The cancer-related mutant p53(273H), which is defective in inhibiting recombination processes, displays a reduced capacity to associate with hRad51-DNA complexes, even under conditions which support DNA-binding. This suggests that hRad51-p53 contacts play a role in targeting p53 to heteroduplex joints and indicates an involvement in recombination immediately following hRad51-mediated strand transfer. To study the initial phase of strand exchange, when heteroduplex joints arise, we applied oligonucleotide based strand transfer assays. We observed that hRad51 stimulates exonucleolytic DNA degradation by p53, when it generates strand transfer intermediates. In agreement with this observation, artificial 3-stranded junction DNAs, designed to mimic nascent recombination intermediates, were found to represent preferred exonuclease substrates, especially when comprising a mismatch within the heteroduplex part. From our data, we propose a model according to which, p53-dependent correction of DNA exchange events is triggered by high-affinity binding to joint molecules and by stabilizing contacts with hRad51 oligomers. Oncogene (2000) 19, 4500 - 4512.

Rad51 uses one mechanism to drive DNA strand exchange in both directions

The Rad51 protein of Saccharomyces cerevisiae, like its bacterial counterpart RecA, promotes strand exchange between circular single-stranded DNA (ssDNA) and linear double-stranded DNA (dsDNA) in vitro. However, the two proteins differ in the requirement for initiating joint molecules and in the polarity of branch migration. Whereas RecA initiates joint molecules from any type of ends on the dsDNA and branch migration proceeds exclusively in the 5'- to 3'-direction with respect to the single strand DNA substrate, initiation mediated by Rad51 requires a complementary 3' or 5' overhanging end of the linear dsDNA and branch migration proceeds in either direction. Here we report that the rates of Rad51-mediated branch migration in either the 5'- to 3'- or 3'- to 5'-directions are affected to the same extent by temperature and MgCl(2). Furthermore, branch migration in both directions is equally impeded by insertions of non-homologous sequences in the dsDNA, inserts of 6 base pairs or more being completely inhibitory. We have also found that the preference of strand exchange in the 5'- to 3'-direction does not change if RPA is replaced by Escherichia coli SSB or T4 gene 32 proteins, suggesting that the preference for the direction of strand exchange is intrinsic to Rad51. Based on these results, we conclude that Rad51-promoted branch migration in either direction occurs fundamentally by the same mechanism, quite probably by stabilizing successively formed heteroduplex base pair.

Spectroscopic methods for analysis of protein secondary structure

Several methods for determination of the secondary structure of proteins by spectroscopic measurements are reviewed. Circular dichroism (CD) spectroscopy provides rapid determinations of protein secondary structure with dilute solutions and a way to rapidly assess conformational changes resulting from addition of ligands. Both CD and Raman spectroscopies are particularly useful for measurements over a range of temperatures. Infrared (IR) and Raman spectroscopy require only small volumes of protein solution. The frequencies of amide bands are analyzed to determine the distribution of secondary structures in proteins. NMR chemical shifts may also be used to determine the positions of secondary structure within the primary sequence of a protein. However, the chemical shifts must first be assigned to particular residues, making the technique considerably slower than the optical methods. These data, together with sophisticated molecular modeling techniques, allow for refinement of protein structural models as well as rapid assessment of conformational changes resulting from ligand binding or macromolecular interactions. A selected number of examples are given to illustrate the power of the techniques in applications of biological interest.