MCM8 and MCM9 form a functional helicase complex (MCM8/9) that plays an essential role in DNA homologous recombination repair for DNA double-strand break. However, the structural characterization of MCM8/9 for DNA binding/unwinding remains unclear. Here, we report structures of the MCM8/9 complex using cryo-electron microscopy single particle analysis. The structures reveal that MCM8/9 is arranged into a heterohexamer through a three-fold symmetry axis, creating a central channel that accommodates DNA. Multiple characteristic hairpins from the N-terminal oligosaccharide/oligonucleotide (OB) domains of MCM8/9 protrude into the central channel and serve to unwind the duplex DNA. When activated by HROB, the structure of MCM8/9’s N-tier ring converts its symmetry from C3 to C1 with a conformational change that expands the MCM8/9’s trimer interface. Moreover, our structural dynamic analyses revealed that the flexible C-tier ring exhibited rotary motions relative to the N-tier ring, which is required for the unwinding ability of MCM8/9. In summary, our structural and biochemistry study provide a basis for understanding the DNA unwinding mechanism of MCM8/9 helicase in homologous recombination.
This paper presents important findings on the hexametric structure of MCM8/9, which potentially explain its role as a DNA helicase in homologous recombination. This solid work will be of interest to biologists studying DNA transactions.
DNA helicases play critical roles in multiple cellular processes including DNA replication, transcription, recombination and repair. Dysfunction of multiple DNA helicases have been correlated with various human diseases including cancers, Bloom Syndrome, Werner Syndrome, Fanconi Anemia, reproductive deficiencies and infertility, etc (1,2). One of the most extensively studied subfamilies of DNA helicases is the minichromosome maintenance (MCM) proteins. In the human genome, the MCM family contains two subgroups of AAA+ ATPase/helicase complexes, consisting of eight MCM members (MCM2-9) (3). Structural studies of these complexes are critical for elucidating their mechanisms.
The hexameric MCM2-7 complex, also known as a replicative helicase, plays a vital role in the initiation and elongation of eukaryotic chromosome replication. Previous structural and biochemical studies have revealed that the MCM2-7 complex exhibits high structural flexibility during DNA replication (4). It can adopt multiple conformations by interacting with its loaders and activators to form different intermediates, including the helicase-loading intermediate Orc-Cdc6-Cdt1-MCM2-7 (OCCM), the inactive MCM2-7 Double Hexamer (DH) and the active replicative helicase complex Cdc45-MCM2-7-GINS (CMG) (5–8). For loading onto DNA, the apo MCM2-7 hexamer first binds to the licensing factor Cdt1. In the presence of ATP, the initiators ORC and Cdc6 assemble on the origin DNA to form an active loading platform and load the Cdt1-bound MCM2-7 onto DNA to form the OCCM complex (4). The origin DNA is encircled into the central channel of the MCM2-7 hexamer via the MCM2-MCM5 entry “gate” (9). While ATP hydrolysis occurs in OCCM, it seems that Cdt1, Cdc6 and ORC are released in sequential disassembly steps, and then the first loaded MCM2-7 encircling DNA, recruits a second ORC-Cdc6 complex, which in turn loads the second Cdt1-bound MCM2-7 to form a head-to-head MCM2-7 DH (4). During double-hexamer assembly, the N-tier ring of MCM2-7 rotates by nearly 30° when aligned with the MCM2-7 hexamer in its loading intermediate OCCM (10). As the MCM2-7 double hexamer on DNA is inactive, lots of regulators along with the helicase activators Cdc45 and GINS complex, are required to convert the double hexamer into two activated CMG replicative helicases that translocate along the leading strands with a bi-directional replication mechanism (11–13). The aforementioned structural studies have also revealed that layers of characteristic hairpin loops inside the chamber of MCM2-7 hexamer, including the helix-2 insertion loops (H2I), the presensor-1 (PS1) as well as β-turn motifs of oligosaccharide/oligonucleotide (OB) domains, contribute to DNA binding and separation (14). Additionally, another allosteric communication loop (ACL) conserved in all MCM helicases was also found to be important for helicase activity. The ACL mediates inter-subunit interactions between the N-terminal DNA processing domain and the C-terminal AAA+ motor domain (15,16).
MCM8 and MCM9 make up the other helicase complex, which is the homolog of MCM2-7 (17). In contrast to the replicative helicase MCM2-7, MCM8 and MCM9 function as a DNA helicase complex (MCM8/9) in HR-mediated DNA repair for DNA double-strand breaks (DSBs) and DNA interstrand crosslinks (ICLs) (18). Several studies have suggested that cells lacking MCM8 and MCM9 are highly sensitive to DNA cross-linking reagents (19,20). The MCM8 and MCM9 knock-out mice were found to be sterile and revealed gametogenesis deficiency as well as chromosomal damage due to impaired HR (18). During HR repair, MCM8/9 was rapidly recruited to the DNA damage sites and colocalized with the recombinase Rad51 (21). It also interacted with the nuclease complex MRN (MRE11-RAD50-NBS1), which was required for DNA resection at DSBs to facilitate the HR repair (22). Recently, HROB (also known as C17orf53 or MCM8IP) has been identified as an essential factor in loading the MCM8/9 complex to the sites of DNA damage and stimulating its helicase activity to promote replication fork progression during DNA recombination and replication. Loss of HROB led to HR defects because it is necessary for the recruitment and activation of MCM8/9 (23,24). The HROB harbors a flexible N-terminal domain (NTD), a central proline-rich region (PRR) and an OB fold domain at the C-terminus. It was reported that the MCM8/9 binding motif (MBM) is embedded in the PRR region and binds to the NTD domains of MCM8/9 (24). Although progress has been made in understanding the role of MCM8/9 HR during the past decade, the structure/function relationship, especially the structural characterization of MCM8/9 helicase activity remains unknown. Notably, MCM8/9 is composed of two subunits, MCM8 and MCM9. It is challenging to extrapolate its role through a simple comparison with MCM2-7, which consists of six subunits.
MCM8 and MCM9 proteins are structurally similar, with a MCM domain at the N-terminus (NTD) and an AAA+ ATPase domain at the C-terminus (CTD). The NTD and CTD are connected by flexible linkers (N-C linkers), which allow the movement of the AAA+ domain during unwinding (25). The NTDs of MCM8 and MCM9 share an analogous domain organization and can be divided into three domains: the zinc finger domain (ZF), the helical domain (HD) and the OB domain. In addition to the AAA+ ATPase domain, MCM8 possesses an additional WHD domain in the CTD while MCM9 contains a putative helix-turn-helix (HTH) domain and an extensional C-tail at the C-terminus (26).
In the previous study, we reported a 6.6 Å cryo-EM structure of human MCM8/9 (hMCM8/9) NTD ring. We also analyzed the crystal structures of MCM8 and MCM9 NTDs separately, which exhibited conformational changes when assembled into the MCM8/9 hexamer (27). However, the low-resolution structure of the MCM8/9 NTD ring is insufficient to fully illustrate the assembly and activation mechanisms of the MCM8/9 hexamer. Here, we present a near-atomic resolution structure of the chicken MCM8/9 (gMCM8/9) complex and a 3.95 Å cryo-EM structure of the human MCM8/9 NTD ring under HROB induction. Based on these structures, we conducted further investigation on the structural features of the DNA-binding region of MCM8/9 and the conformational changes required for its activation by HROB. Our structural and biochemical studies shed light on the structural and functional relationship of the MCM8/9 helicase.
Overall structure of the gMCM8/9 helicase complex
As our previous study shows the hMCM8/9 helicases complex is unstable for high-resolution structure determination, we first solved the cryo-EM structure of MCM8/9 complex from Gallus gallus. (Gallus gallus, gMCM8/9; amino acid, gMCM8: 50-830; gMCM9:1-691) (Figure 1A, Table 1, Materials and Methods). The gMCM8/9 shares up to 80% sequence homology with hMCM8/9. Based on the electron density map, the structures of the NTD and CTD were reconstructed at the resolutions of 3.66 Å and 5.21 <, respectively, and then aligned into a composite structure (Figure 1B, 1C and S1).
The structure reveals that gMCM8 and gMCM9 alternately assemble into a pinwheel-like hexamer with a central channel. The trimerization of the gMCM8/9 heterodimer is in line with our previously reported hMCM8/9 structure at 6.6 Å resolution (27). When we performed structural superposition, the NTD ring of gMCM8/9 superimposes well onto that of hMCM8/9, indicating that eukaryotic MCM8 and MCM9 form a conserved hexamer complex during evolution (Figure S2, S3 and S4). The diameter of the NTD ring is ∼132 < and that of the inner channel is ∼ 28 < which is large enough to accommodate dsDNA to pass through (Figure 1C and S5). The assembly pattern of gMCM8/9 complex with a three-fold (C3) symmetry axis is unique among MCM helicases. It arranges the three MCM8 and three MCM9 in the same planes, enclosing a channel with multiple symmetrical hairpins. This assembly pattern differs from that of the eukaryotic MCM2-7 heterohexamer with non-symmetry and the archaeal MCM homohexamer with 6-fold symmetry (7,28). In addition, a distinguishing feature of the gMCM8/9 is that the helical domain of gMCM8 contains a long α5 helix that protrudes around the hexameric ring (Figure 1B).
The domain structures of gMCM8/9’s, including the MCM domains, the AAA+ motor domains and the WH domain, fit well into their electron density maps (Figure S5). For each gMCM8/9 heterodimer, the positions of characteristic hairpin loops of MCM helicase, such as H2I, PS1, OB hairpins (referred to OB-hps) and the ACL were well defined. Among them, the OB-hps extend to the central channel, while the ACLs are close to the H2I and PS1 hairpins within the AAA+ domains (Figure 1B).
The assembly of gMCM8 and gMCM9 in the NTD ring
As our structural analysis revealed that the gMCM8/9 hexamer is the trimerization of the MCM8/9 heterodimer, this leads to the formation of three dimer interfaces within each MCM8/9 heterodimer and three putative trimer interfaces between each neighboring MCM8/9 heterodimer (Figure 2A). The high-resolution structure of the gMCM8/9 NTD ring allows us to analyze these two types of interfaces in detail, which serve as the structural basis of MCM8/9 assembly.
The dimer interface between NTDs of gMCM8 and gMCM9 is largely composed of the β-strands from the OB domain of gMCM9 and the flanking ZF and OB domain of gMCM8. Both the polar and hydrophobic interactions between respective β-strands are responsible for the binding of the MCM8/9 dimer. For example, the side chains of R220gMCM9 and Q248gMCM9 form hydrogen bonds with the carbonyl groups in the main chains of I172gMCM8 and T226gMCM8, respectively. In addition, the main chain of T119gMCM9 also forms hydrogen bond with the neighboring carbonyl group in the side chain of D228gMCM8. Furthermore, the benzene ring of F253gMCM9 inserts into a hydrophobic pocket surrounded by residues L179, P217, F181, L201 and Y199 from gMCM8 and W250 from gMCM9, which further stabilize the interaction between gMCM8 and gMCM9 (Figure 2B, S6A and S6B).
Compared to the strong interaction between the dimer interface, the interaction within trimer interface is weaker. The trimer interface is dominated by electrostatic interaction and hydrogen-bonding interactions. The side chain of R152gMCM8 forms salt bridges with both the side chains of D200gMCM9 and D232gMCM9. In addition, the side chain of K158gMCM8 forms a hydrogen bond with the side chain of E37gMCM9 (Figure 2C, S6C and S6D). Based on the preceding analysis, we assume that opening a “gate” on the trimer interface for DNA loading will be easier than on the dimer interface. Taken together, these interactions facilitate and stabilize the assembly of the NTD hexameric ring of gMCM8/9.
Structural comparison of MCM8/9 and MCM2-7
To explore the functional state of gMCM8/9, we performed structural comparisons between gMCM8/9 and various MCM2-7-containing complexes. We first compared the gMCM8/9 structure with MCM2-7 in the OCCM complex by aligning their respective C-tier rings. As shown in Figure 3A, the NTDs of gMCM8/9 fit well to those of MCM7, MCM3 and MCM5, but rotated slightly in the following MCM2, MCM6 and MCM4. This minor variation is due to the opening “gate” of the MCM5-MCM2 interface in the OCCM complex. When superimposing the hexameric gMCM8/9 complex onto one hexamer of the MCM2-7 DH, we found that the two C-tier rings aligned well while their N-tier rings had a ∼30° clockwise rotation (Figure 3B and 3C). Such rotation could be caused by the two hexameric complex’s different assemblypatterns. A more plausible explanation is that the assembly of the MCM2-7 DH induces the conformational change in its N-tier ring (10). We also compared the gMCM8/9 structure with that of MCM2-7 in the apo CMG complex and found that the NTD ring has a similar translation and rotation but slighter than comparison with MCM2-7 DH (data not shown). As the NTD ring of MCM2-7 that transits from the DH to the CMG just requires a small domain translation (29). Overall, our structure of gMCM8/9 was conformationally similar to the structure of MCM2-7 complex in its loading intermediate OCCM rather than that in the DH. This is also consistent with previous finding that the recombinant MCM8/9 complex alone exhibited limited helicase activity in vitro (24). Therefore, we presume the current structure of gMCM8/9 may represent a loading state of MCM8/9 helicase.
HROB activates MCM8/9 probably by inducing conformational change of its NTD ring
To explore the activation mechanism of hMCM8/9, we examined its activity in the presence of OB-fold containing protein HROB (Figure 4A). Using a helicase assay, we found that the HROB-MBM motif was unable to stimulate the helicase activity of hMCM8/9. However, it’s CTD, which includes the OB fold domain, could significantly activate hMCM8/9 (Figure 4B). The HROB-CTD not only directly interacts with hMCM8/9 NTD, but also forms a complex with it, as confirmed by the GST pulldown assay (Figure 4C).
We collected the complex formed by the HROB-CTD and hMCM8/9 NTD proteins described above and reconstructed their structure using cryo-electron microscopy (Figure 4D, Table 2, Materials and Methods). Unfortunately, the density map of HROB-CTD was not observed in the reconstruction. And the 3.95 Å structure of hMCM8/9 NTD ring was reconstructed with C1 symmetry, in contrast to the structure of gMCM8/9, which has C3 symmetry (Figure S7). We also found that the structure of hMCM8/9 NTD ring underwent a large conformational change when compared to the previously described hMCM8/9 NTD (27). As a result, the two conformations were designated as Conformation I (hMCM8/9 NTD at 6.6 Å) and Conformation II (hMCM8/9 NTD at 3.95 Å) (Figure 4E and 4F). In line with our observation, when superimposed these two conformations of hMCM8/9 NTD, we found that one of the trimer interfaces in Conformation II was expanded by ∼5 Å (Figure 4G). This expansion further pushed its neighboring MCM8/9 heterodimers to shift to either side. The left dimer translates downward by ∼8 Å and tilts by ∼6° while the right one shifts by ∼7 Å and tilts by ∼8° (Figure 4I and 4J). Interestingly, the ZF of MCM8 shifts by ∼3 Å to the left and the ZF of MCM9 shifts by ∼4.5 Å to the right, resulting in an overall expansion of ∼6 Å (Figure 4G). These changes are consistent with our analysis for the trimer interface of gMCM8/9, which is easier to open for accommodating DNA.
Notably, the N-tier ring of MCM2-7 from the CMG-DNA complex translates laterally by ∼12 Å and tilts by ∼12° when compared to the structure of MCM2-7 DH (10). Without DNA, the C-tier ring of MCM2-7 from the apo CMG also adopts two different conformations including tilted and untilted, indicating a conformational change of MCM2-7 during activation (29). Based on the above observations and the stimulation effect of HROB on the helicase activity of MCM8/9, we propose that HROB activates MCM8/9 by inducing a conformational change of its NTD ring.
Spatial distribution and functional role of the OB fold hairpins
As the hairpin loops of OB fold domains are essential for DNA binding and/or unwinding in the reported hexameric MCM helicase complexes (14), we analyzed the two types of hairpin loops from OB domains in the NTD rings of gMCM8/9 and hMCM8/9. The first type of OB hairpins of gMCM8 and gMCM9 (refer to OB-hp) protrude into the central chamber in a three-fold symmetry manner and narrow down the channel (Figure 5A). In particular, the OB-hp of gMCM8/9 form two layers. The three OB-hp of gMCM8 form the upper layer and the three gMCM9 OB-hp comprise the lower layer (Figure 5C). Importantly, the tip of the gMCM8 OB-hp contains several positively charged residues (K337, K339, K342) while the gMCM9 OB-hp contains continuous polar residues (H254, Q255, D256) (Figure 5D), which are conserved in hMCM8/9 (R345, K347, K350 in hMCM8 and Q252, Q253, D254 in hMCM9). As seen in the recently reported cryo-EM structure of the S. cerevisiae CMG on a forked DNA, certain conserved residues, such as K364 and K367 from MCM7 OB hairpin, R449 and R451 from MCM4 OB hairpin loop and polar residues Q308 and N307 from MCM3 OB fold domain, are involved in the interaction with forked DNA and may contribute to its separation (4). We artificially docked a forked DNA into the central channel to generate a gMCM8/9-DNA model and found that the OB-hps of gMCM8 are capable to closely contact with it and insert their highly positively charged terminal loops into the major or minor grooves of the DNA strand, implying that they could be involved in substrate DNA processing and/or unwinding (Figure 5C).
The OB domains of gMCM8 and gMCM9 also harbor another type of interesting hairpin loop known as the allosteric communication loops (ACLs) (Figure 5B). These loops are located on the opposite side of the OB-hp but extend to the CTD ring of gMCM8/9. Several studies have found that the ACLs of M. thermautotrophicus MCM (mtMCM) and Sulfolobus solfataricus MCM (SsoMCM) helicases facilitate interactions between NTDs and AAA+ catalytic domains of hexameric helicases as well as modulate the positioning of the OB-hp and then regulate the helicase activities (15,16). As shown in Figure 5E and 5F, the ACLs show a high degree of conservation when performing superposition of gMCM8 NTD onto the structure of mtMCM(15). The sequence alignment of the ACLs also highlights some key residues including Q235, E236 and P249 in gMCM8 (Q294, E295, P308 in hMCM8) that are conserved in mtMCM and important for the mtMCM helicase activity (Figure 5E and 5F) (15). Thus, we propose that the ACLs are required for the helicase activity of MCM8/9.
To investigate the role of OB domain hairpins in regulating the helicase activities of MCM8/9, we deleted the tip residues (345-RNKNDK-350 in hMCM8) to shorten the OB-hp or mutated the polar residues (Q252, Q253, D254 in hMCM9) to change the charge property of OB-hp. Then we performed an in vitro helicase assay based on plasmid-based substrate (30,31), using the recombinant hMCM8/9 complex and corresponding mutants (Figure S8). In the helicase assay, we observed that both the wild type (WT) hMCM8/9 complex and its mutants only exhibited limited DNA unwinding activity. In the presence of HROB, the WT hMCM8/9 exhibited ∼6-fold stimulation of the helicase activity, while the mutants deleting the tip residues of the hMCM8 OB-hp or mutating that of hMCM9 OB-hp both showed decreased helicase activities. However, the hMCM8 ACL mutant (hMCM8_E295R) corresponding to that in mtMCM had no significant change even though in the presence of HROB, implying that the MCM8’s NTD doesn’t interact with its CTD directly through ACL hairpins (Figure 5G and 5H). Collectively, the above structural analyses and biochemical experiments suggest that the OB-hps of MCM8/9 play an important role in DNA unwinding.
Structural dynamic study of the gMCM8/9 CTD ring
Since the C-tier of MCM2-7 undergoes conformational change coupling with ATP hydrolysis (32), we propose that the CTD AAA+ domain of MCM8/9 is also mobile for DNA unwinding, serving as a motor domain. Notably, the electron density of gMCM8/9 CTD ring is a little blurry and its structure was reconstructed at 5.21 Å (Figure S1 and S9). On the basis of our cryo-EM dataset, we performed structural flexibility analyses using multi-body refinement to investigate the dynamics of the gMCM8/9 complex. In this method, we used separate focused refinements and generated movies to describe the most important motions on the relative orientations of the rigid bodies, which allows us to gain insight into the molecular motions of gMCM8/9, particularly regarding the structural flexibility of its CTD ring (33–35). As expected, twelve components were generated, each with a discrete number of independently moving bodies, providing a characterization of the motions in gMCM8/9 CTD ring. As shown in the movie 1, the CTD AAA+ motor ring reciprocally rotates right and left relative to the NTD ring. It would be very interesting to see similar rotary movement along DNA, which may underlie the translocation/unwinding mode for MCM8/9 helicase.
The N-C linkers is required for MCM8/9’s helicase activity
Unlike MCM8/9, such rotary motion between NTD and CTD rings has yet to be observed in the structures of MCM2-7. We analyzed the linkers connecting the NTD and CTD (N-C linkers) of MCM8/9 and found that these N-C linkers, which are 30-40 residues long, are longer than those in MCM2-7. In MCM2-7, the N-C linker of MCM 3 is the shortest, which may restrict the rotary motion of MCM2-7 (Figure S10). To test whether the N-C linkers are crucial for the relative motion between N-tier and C-tier rings of MCM8/9, we shortened the N-C linkers by deleting several amino acids in loop regions of hMCM8 (MCM8Δ369-377) or hMCM9 (MCM9Δ283-287) (Figure S11). In the helicase assay in vitro, we found that both single and combined deletions have significantly lower unwinding ability than the WT MCM8/9 complex. The presence or absence of HROB had no effect on this result, indicating that the N-C linkers are required for the MCM8/9 helicase activity (Figure 6A and 6B).
We next studied whether the N-C linkers are important for the in vivo functions of MCM8 or MCM9. To do so, we employed a chemoresistance assay because MCM8 and MCM9 had been reported to mediate cellular resistance to the DNA cross-linking reagent cisplatin, which induces DNA interstrand crosslinks that are mostly repaired by HR (17,19,36). As shown in the cell-viability assay, we found a notable reduction in survival of MCM8- or MCM9- deficient DT40 cells following cisplatin treatment compared to the WT DT40, indicating that the KO cells exhibit significant sensitivity to cisplatin. Expression of WT MCM8 or MCM9 cDNAs in their respective KO cells almost completely complemented chemoresistance. In contrast, expression of their N-C linker deletion mutants had minimal complementary effects (Figure 6C, 6D, 6E and 6F). These results suggest that the N-C linkers of MCM8/9 are required for its cellular resistance to DNA-damaging reagents.
In this study, we reconstructed two structures of the MCM8/9 helicase complex at near-atomic resolution. We found that MCM8 and MCM9 form a single heterohexameric ring with a centric channel that allows the DNA to pass through. The OB hairpin loops within the centric channel may be involved in interacting with and separating substrate DNA. The N-C linkers of MCM8/9 play important roles in helicase activity. Notably, the regulator HROB induces conformational change of MCM8/9 NTD ring, which may underlie the helicase activation.
MCM8/9 and MCM2-7 are both MCM proteins with highly conserved domain organization and sequence similarity (3,37). Compared with the heterohexameric structure of MCM2-7, the MCM8/9 hexamer is arranged in an alternate mode with the three-fold symmetry axis. Interestingly, it was reported that the Escherichia coli helicase DnaB also exists in a quaternary state with C3 or C6 symmetry and can change the symmetry states reversibly depending on the pH value of the buffers or the presence of nucleotides (38). It has been proposed that the symmetry transition plays a role in the loading of regulators on DNA and is directly linked to DNA translocation (39,40). Similarly, we doubt that MCM8/9 changes its symmetry from C3 to C1 symmetry in the presence of HROB. In other words, helicases may orchestrate the inner channels by controlling the relative positions of subunits, assisting in DNA loading and unwinding.
Taking inspiration from the structural and unwinding mechanisms of MCM2-7, we also investigated the key elements of MCM8/9 involved in DNA interaction or separation. In the central channel of the MCM2-7 complex, several structural elements involved in DNA binding and strand separation were illustrated by a series of CMG structures complexed with different types of DNA (4). According to recent research, the OB hairpin loops of MCM3, MCM4, MCM6 and MCM7 form “dam-and-diversion tunnel” to block and divert the lagging strand from the leading strand (4). For the MCM8/9 complex, the OB-hps of MCM8 and MCM9 are essential for DNA unwinding. They form a two-layered structure with three-fold symmetry. In which, three OB-hps of MCM9 form the upper layer, while three OB-hps of MCM8 form the lower layer, suggesting a synergistic mechanism for MCM8/9 OB hairpins to interact with DNA. When modelled a fork DNA into the channel, the OB-hp loops make very close contact with the DNA strand that passes through the channel, capable of inserting their loops into the major or minor grooves of the DNA in a sequential manner. While in the Conformation II of hMCM8/9 NTD ring, the OB-hps of hMCM8 and hMCM9 lose their C3 symmetry without significant conformational change. This symmetry transition could be caused by HROB while inducing trimer interface expansion. In addition, the highly conserved residue F363 from OB hairpin loop of MCM7 serves as a strand separation pin to unwind forked DNA via making a π-π interaction with DNA (41). The highly conserved phenylalanine was also found on the OB-hps of MCM8 and MCM9 (F353 of hMCM8, F251 of hMCM9) (Figure 5D). They will most likely perform functions similar to those found in MCM2-7.
We also present a new conformation (Conformation II) of the MCM8/9 NTD ring induced by HROB, which may enhance the DNA binding ability of MCM8/9 complex. We found that when MCM8/9 is complexed with the CTD domain of HROB, it becomes more stable and exhibits significantly higher helicase activity. The HROB CTD domain is about 20 kDa and is composed of a flexible loop and an OB fold domain. The flexible loop containing ∼20 amino acids is sufficient to interact with the MCM8/9 complex, but it alone is unable to stimulate the helicase activity of MCM8/9 (24). Our helicase assay demonstrated that the OB fold domain of HROB is required for MCM8/9 activation. The OB fold domain is a nucleic acid-binding motif that likely help to load the DNA into the central channel of MCM8/9. Induced by HROB, the structure of MCM8/9 NTD ring convert its C3 symmetry to C1 and two of the MCM8/9 heterodimers translate and tilt following the expansion of the trimer interface. When the MCM2-7 double hexamer encircles the double-stranded origin DNA, the MCM5-MCM2 interface characterizes as a “gate” that undergoes open and closed conformation change to engage the DNA entry (9). The expansion of the MCM8/9 trimer interface may involve a transient conformation state to facilitate the DNA entry into the inner channel of the MCM8/9 hexamer. Thus, HROB not only aids in the loading of DNA into MCM8/9, but also acts as an activator, collaborating with MCM8/9 to promote DNA unwinding.
Benefit from the rapid development of cryo-electron microscopy technology over the past decade, a dozen structures of MCM2-7-containing intermediates have been resolved including the OCCM, DH and CMG with/ without DNA (4). Based on these structures, the hand-over-hand rotary mode and the “dam-and-diversion tunnel” model had been proposed for MCM2-7 helicase in DNA translocation and unwinding (4,12). For MCM8/9 helicase, there is currently no established model to describe its unwinding mechanism. Here, our structural and biochemical studies have provided some hints in this regard. First of all, the structural dynamics analysis reveals that the MCM8/9 CTD ring rotates relative to the NTD ring. Secondly, the length of the MCM8 and MCM9 N-C linkers are surprisingly much longer than those of MCM2-7 complex, which may permit the rotations of MCM8/9. The fact that single-site structure-guided mutations on the ACL loops has no effect on the helicase activity of MCM8/9 further supports the idea that the N-C linkers, rather than the ACL loop, mediated the interaction between the N-tier ring and the C-tier ring of MCM8/9. Last but not least, HROB can convert inactive MCM8/9 into an active helicase by recruiting and inducing conformational change in the MCM8/9 NTD ring.
Although a rotary model is in sight and our findings also hint at a distinct unwinding mechanism of MCM8/9, many questions remain unsolved. The most pressing issue is that the type of DNA substrate unwound by MCM8/9 is unknown. In vivo sequencing and in vitro DNA binding results are required to screen appropriate substrate. Furthermore, as the HROB only stimulate ∼60% of the MCM8/9 helicase in vitro, we believe that other regulators remain to be discover for the full activation of MCM8/9 helicase. Based on these results, the high-resolution structures of MCM8/9 with DNA in the presence of ADP or ATP are necessary to fully clarify the unwinding mechanism of MCM8/9. Despite these limitations, the MCM8/9 complex structures provided in this study will be a rich source of information for further research into the mechanism of DNA helicase in homologous recombination, as well as for the analysis of multiple disease mutations in MCM8/9 associated with premature ovarian failure as well as cancers.
Material and methods
The human MCM8/9 genes were obtained by reverse transcription of mRNA from HeLa cells. The gene fragments coding MCM8/9 of Gallus gallus were also from corresponding reverse transcripts and the total RNA was extracted from the chicken DT40 cell line. Human MCM8_61-840 and MCM9_1-684 were inserted into the pFastbac-Dual vector (Invitrogen). Gallus gallus MCM8_50-830 and MCM9_1-691 were inserted into the pFastbac-1 vector, respectively. All the MCM8 were fused with 6× His tag at their C-terminus for subsequent affinity purification. The genes coding hMCM8/9 NTDs (MCM8_61-376, MCM9_1-276) and relative mutants were cloned into multiple cloning site I (MCS I) of pRSFDuet-1 vector (Novagen). All of the mutations were introduced by the standard PCR-based mutagenesis method described before (42). Human HROB (NM_001171251.3) was obtained by reverse transcription of mRNA from 293T cells. HROB-CTD_391-580 was inserted into pET28S-SUMO vector which expresses the N-terminal SUMO tag. All constructs were confirmed by DNA sequencing.
Protein expression and purification
To express the core hMCM8/9 and gMCM8/9 complex proteins, recombinant baculoviruses were prepared using the Bac-to-Bac expression system (Invitrogen). The proteins were expressed in Trichoplusia ni (BTI-Tn5B1–4, Hi5) insect cells for 60 hours at 27℃ before harvesting. The cell pellet was collected by centrifugation at 800× g and then freshly frozen by liquid nitrogen and stored at -80 ℃ before use. For purification, the cell pellet was resuspended with lysis buffer (50 mM HEPES, pH 8.0, 350 mM KCl, 50 mM sodium glutamate, 10% glycerol and 30 mM imidazole) and lysed by cell homogenizer (Avestin Emulsiflex C3). Then the lysate was clarified by centrifugation at 16,000× rpm for 30 min. The supernatant was mixed with nickel-NTA resin (Novagen) and was continuously stirred for 1 hour at 4 ℃ for adequate binding. The beads were then collected into a column and extensively washed with lysis buffer to remove undesired proteins. Target proteins were then eluted down with elution buffer (50 mM HEPES, pH 7.8, 80 mM KCl, 20 mM sodium glutamate, 5% glycerol and 300 mM imidazole). The proteins were further purified by HiTrap Heparin HP 5ml column (GE Healthcare) and size exclusion column (Superose 6 Increase 10/300 GL, GE Healthcare) in buffer (20 mM HEPES, pH 7.8, 150 mM KCl, 1 mM DTT and 5% glycerol). The strategies of purification of MCM8/9 mutants were as the same as that used for the core hMCM8/9 complex.
For expression of HROB-CTD, 10 mL BL21 codon plus strain containing the recombinant pET28a-HROB-CTD plasmid was inoculated into 1 L LB medium. The bacteria solution was firstly cultured at 37℃ until the OD600 ≈ 0.8 and then cooled down to 16℃. The expression of recombinant protein was induced by IPTG (isopropyl-β-d-thiogalactoside) at a concentration of 0.2 mmol/L and kept expressing for another 20 hours. Subsequently, the cells were harvested and disrupted by a high-pressure cell homogenizer in lysing buffer. After centrifugation, the supernatant was mixed with nickel-NTA beads and retained for 1 h at 4 ℃ with rolling. Before elution, the SUMO tag was cleaved by Ulp Protease and finally the HROB-CTD protein was eluted in the buffer containing 50 mM HEPES, pH 7.8, 80 mM KCl, 20 mM sodium glutamate, 5% glycerol. Following Ni-NTA affinity purification, the eluted protein was further purified by HiTrap Heparin HP 5ml column (GE Healthcare) and size exclusion chromatography (HiLoadTM16/600 Superdex TM pg, GE Healthcare) in the same buffer using for hMCM8/9. The GST tag HROB-CTD was also expressed in BL21 (DE3) E. coli cells and purified using GSH-Sepharose affinity chromatography (GE Healthcare) following by the size-exclusion chromatography. And the HROB-MBM peptides (HROB_391-413) were commercial synthesized by China Peptide Co., Ltd.
Cryo-EM sample preparation
The fresh gMCM8/9 proteins were used to prepare the cryo-EM sample. Briefly, 0.8 µL crosslinker BS(PEG)9 (Thermo scientific) and 10 μM protein were incubated for 1 hour at 4℃ and filtered by Superose 6 in buffer (20 mM HEPES, pH 7.8, 150 mM KCl) after centrifugation. The fractions were identified by SDS-PAGE and the peak fraction was diluted to about 0.3 mg/mL to prepare cryo-EM samples. 4 µL fresh sample was applied on glow-discharged Cu holey carbon grids (Quantifoil R 1.2/1.3) and incubated for 1 min at room temperature. Grids were then blotted for 3 s in 100% humidity at 4℃ and plunged frozen in liquid ethane cooled by liquid nitrogen using Vitrobot (Thermo).
For the hMCM8/9-HROB complex, the hMCM8/9 NTDs and HROB proteins were mixed by the ratio 1:3. 4 µL fresh sample (0.48 mg/mL) was applied on glow-discharged Au holey carbon grids (Quantifoil R 0.6/1) and the following operation were the same as describe above for preparing the sample of gMCM8/9 complex.
Cryo-EM data acquisition
The grids of gMCM8/9 were loaded onto a transmission electron microscope operated at 200 kV (FEI Talos Arctica) in liquid nitrogen temperatures. Images were recorded on a Gatan K2 camera in counting mode, at a nominal magnification of 29,000×. The total dose rate was about 50 e-/ Å2 for each micrograph stack. The total dose was fractionalized to 32 frames equally, corresponding to a physical pixel size of 1.0 Å (0.5 Å super-resolution pixel size). The defocus range was set from -1.7 to -2.2 μm. 2,891 raw movie micrographs were collected using the serialEM software (43).
The grids of hMCM8/9 NTD complex were loaded onto a transmission electron microscope operated at 300kV (FEI Titan Krios) equipped with a K3 direct electron detector. The parameter settings of total dose and defocus range were as same as that of gMCM8/9. While the images’ physical pixel size is 0.827 Å and the super-resolution pixel size is 0.4135 Å. 4078 raw movie micrographs were collected using the Thermo Fisher Scientific EPU software.
Single particle data processing
The data set of gMCM8/9 was processed by Cryosparc-3.2 (44). All the movies were processed using patch motion correction and patch CTF estimation. A total of 2,134,191 particles were automatically picked by template picker and 1,574,821 particles were produced after the inspect particle picks treatment. Three parallel runs of 2D classification (K = 100) were performed using the data (binned 4) with a box size at 256 pixels. A small subset of the “good” particles with clear secondary-structure features was selected from 2D classification results and was used to generate 4 initial models. The 434,320 particles in the best initial model class were re-extracted with original pixel size and further refined to a 6.7 Å low-resolution reconstruction using non-uniform refinement (New) (45). Four resolution gradient templates generated from non-uniform refinement were used for the multi-reference heterogeneous refinement, resulting in a 4.31 Å global reconstruction (290,080 particles) renamed consensus map. The NTD map was improved to 3.73 Å with local refinement followed by applying a tight N-terminal mask with C3 symmetry. Local CTF refinement further improved the resolution to 3.66 Å. For the CTD map, a reference was morphed from a predicted gMCM8/9 CTD model using AlphaFold (Figure S12) (46). Before processing CTD local refinement (applying C3 symmetry), the projection of NTD particles from the 3.66 Å NTD map were subtracted from the 4.31 Å consensus map, and generated a 6.57 Å CTD hexamer map. Finally, 3D classification and local refinement resulted in a 5.21 Å (mean resolution, validation from Phenix (47)) CTD map using 72,948 particles. The 3.66 Å NTD map and 5.21 Å CTD map were aligned to the consensus map, respectively, and then merged together to generatea final map.
The EMDB accession code is EMD-32346 and the PDB entry ID is 7W7P. The image processing and 3D reconstruction steps were illustrated in Figure S13.
RELION-3.1.1 (48) was also used to process the same data, resulting in a 4.2 Å NTD map and a 5.4 Å CTD map. All the movie frames were motion-corrected by MotionCorr2 (49) with a patch alignment of 5×5. CTF parameters were estimated using Gctf (50). 965,175 particles were auto-picked and extracted to 4×binned particles, which rescales the initial box size at 256 pixels to the final box size at 64 pixels. The particles were subjected to 2 rounds of 2D classifications (K = 100). 388,378 particles with clear secondary-structure features were selected for 3D initial model generation. Three parallel runs of multi-reference 3D classification (K = 9) were performed and the best class was merged. 144,202 particles were re-extracted at original pixel size and a final map was imported as reference and were submitted to the final 3D auto-refinement. Finally, we yielded a 6.6 Å global 3D density map estimated by the gold-standard Fourier shell correlation at a correlation cutoff value of 0.143.
As the density of the gMCM8/9 CTD was a little blurred, we then used the multi-body refinement program (51) to process the dataset. 6.6 Å gMCM8/9 global map’s NTD and CTD were subtracted as the first and the second subtracted body respectively using Chimera (52). The corresponding masks were generated from the final map using mask creation in RELION. The multi-body refinement resulted in a 4.2Å 3D density map of gMCM8/9 NTD and a 5.4 Å 3D density map of gMCM8/9 CTD. It generated 12 components, with each including a discrete number of independently moving bodies. The first five components explained 86.77% of the variance in the data. The image processing and 3D reconstruction steps were illustrated in Figure S14.
The data set of hMCM8/9 NTD complex was processed by Cryosparc-3.3. In brief, a total of 1,855,585 particles were automatically picked by template picker and 677,960 particles were produced after the inspect particle picks treatment. After the processing of 2D classification, Ab-Initio Reconstruction and Heterogeneous refinement (K=3), the 95,257 particles in the best 3D class were re-extracted with original pixel size and further refined to a 3.95 Å resolution reconstruction (C1 symmetry) using Non-uniform refinement (New). The image processing and 3D reconstruction steps were illustrated in Figure S15.
Model building and refinement
For the gMCM8/9 model building, the gMCM8 NTD and gMCM9 NTD monomer models were predicted through SWISS-MODEL by providing the hMCM8 NTD (PDB: 7DP3) and hMCM9 NTD (PDB: 7DPD) structural template, respectively. Firstly, the predicted structures were used as the initial model directly docking into the EM map. As the gMCM8/9 NTD has a C3 symmetry axis, another two copies were generated by horizontal rotating for subsequent rigid body-fitting into the 3D density map with program Chimera. Secondly, using the automatic refinement plugin function in COOT (53) to refine the initial model. The final model was refined against the corresponding map using PHENIX (54) in real space with secondary structure and geometry restraints. The gMCM8/9 CTD hexamer model was predicted as above described and the processing of docking and refinement was the same as the gMCM8/9 NTD model.
For the hMCM8/9 NTD complex model building, in brief, hMCM8 NTD (PDB: 7DP3) and hMCM9 NTD (PDB: 7DPD) models were used as the initial model directly docking into the EM map. The other processing procedure was similar to that of the gMCM8/9 model building. All the structures of the gMCM8/9 complex were validated through examination of the Molprobity scores and statistics of the Ramachandran plots. Molprobity scores were calculated as described in Table 2.
To prepare the substrate, the oligonucleotide (5’-(dT)40GTTTTCCCAGTCACGACG-TTGTAAAACGACGGCCAGTGCC-3’) containing a 40 nt region complementary to the M13mp18(+) stand and a 40 nt oligo-dT at the 5′ end was labeled at the 3′ terminus with [α-32P] dCTP (Perkin Elmer) and annealed to the single-stranded DNA M13mp18(24). 0.1 nM (in molecules) DNA substrates were mixed with 5 µg recombinant MCM8/9 complex as well as its mutants as indicated within each 15 µl volume reaction in the helicase buffer (25 mM HEPES, pH 7.5, 1 mM magnesium acetate, 25 mM sodium acetate, pH 5.2, 4 mM ATP, 0.1 mg/ml BSA, 1 mM DTT). 2.5 µg HROB was used as an activator. To avoid re-annealing, the reaction was supplemented with a 100-fold unlabeled oligonucleotide. The reactions were then incubated at 37 °C for 60 min and stopped by adding 1 µl of stop buffer (0.4% SDS, 30 mM EDTA, and 6% glycerol) and 1µl of proteinase K (20 mg/ml, Sigma) into the reaction for another 10 min incubation at 37 °C. The products were separated by 15% polyacrylamide gel electrophoresis in 1× TBE buffer and analyzed by the Amersham typhoon (Cytiva).
Cell-viability assays were performed as described previously (17). Briefly, the MCM8 or MCM9 KO DT40 cell lines were transfected with plasmids containing wildtype MCM8, MCM9 or the indicated mutants by electroporation. 1.0×103 cells were seeded into each well of a 96-well plate and treated with a range of concentrations of cisplatin. After 48_h of incubation, the cell viability was measured by the Cell counting kit 8 (CCK8) according to the standard procedures. The absorbance readings at 405nm were conducted using an Epoch Microplate Spectrophotometer (BioTek Instrumentals, Inc., Winooski, VT, USA). All experiments were performed at least 3 times.
Pull Down Assay
The HEK293T cells transfected with Flag-hMCM8/9-FL or Flag-hMCM8/9-NTD were cultured overnight and washed twice with cold phosphate-buffered saline (PBS). Cell pellets were resuspended with lysis buffer (20 mM Tris, pH7.5, 150 mM NaCl, 5mM EDTA, 0.5% NP-40, 10% glycerol, protease inhibitor cocktail (Roche, 04693132001)). After incubation for 45 min at 4°C with gentle agitation, the whole-cell lysates were collected by centrifugation (12,000 × g for 15 min, at 4 °C). GST beads coupled with 2 μg GST-HROB or GST alone were then incubated with an equal volume of above HEK293T cell lysates at 4°C for 4h. The beads were washed four times with lysis buffer. Proteins bound to the beads were separated by SDS– PAGE and subsequently immunoblotted with anti-Flag antibody (Cytiva).
The statistical analysis was performed using Prism (GraphPad Software Inc.) on the data from at least three independent experiments. Statistical significance was assessed by two-tailed unpaired Student’s t-test. P ≤ 0.05 was considered significant.
We thank the Cryo-EM Centre, Southern University of Science and Technology for assistance in cryo-EM data collection. We are grateful to Hongjie Zhang (Institute of Biophysics, CAS) for technical assistance in helicase assays.
The 3D cryo-EM map of the gMCM8/9 complex and the hMCM8/9 NTD ring have been deposited to the EMDB database with accession number EMD-32346 and EMD-33989, respectively. The corresponding atomic models have been deposited at the RSCB PDB under the accession code 7W7P and 7YOX, respectively.
The work was supported by the grants from Shenzhen Science and technology planning project (project No. JCYJ20200109142412265, ZDSYS20220606100803007 to Y.L., project No.RCBS20200714114922284 to Z.W.); Pearl River Scholar Team Project of Guangdong Province (project No.2019ZT08Y464 to Y.L.) and National Natural Science Foundation of China (project No.32000860, 32271321 to Z.W.)
Declaration of interests
The authors declare no competing interests.
- 1.History of DNA HelicasesGenes (Basel 11
- 2.Regulation of homologous recombination in eukaryotesAnnu Rev Genet 44:113–139
- 3.MCM proteins and DNA replicationCurrent opinion in cell biology 18:130–136
- 4.Molecular mechanisms of eukaryotic origin initiation, replication fork progression, and chromatin maintenanceThe Biochemical journal 477:3499–3525
- 5.Structural basis of Mcm2-7 replicative helicase loading by ORC-Cdc6 and Cdt1Nature structural & molecular biology 24:316–324
- 6.Mechanism of head-to-head MCM double-hexamer formation revealed by cryo-EMNature 575:704–710
- 7.Structure of the eukaryotic MCM complex at 3.8 ÅNature 524:186–191
- 8.Cryo-EM structures of the eukaryotic replicative helicase bound to a translocation substrateNat Commun 7
- 9.A unique DNA entry gate serves for regulated loading of the eukaryotic replicative helicase MCM2-7 onto DNAGenes Dev 28:1653–1666
- 10.Cryo-EM structure of Mcm2-7 double hexamer on DNA suggests a lagging-strand DNA extrusion modelProc Natl Acad Sci U S A 114:E9529–e9538
- 11.Structure of eukaryotic CMG helicase at a replication fork and implications to replisome architecture and origin initiationProc Natl Acad Sci U S A 114:E697–e706
- 12.Molecular Basis for ATP-Hydrolysis-Driven DNA Translocation by the CMG Helicase of the Eukaryotic ReplisomeCell Rep 28:2673–2688
- 13.The mechanism of eukaryotic CMG helicase activationNature 555:265–268
- 14.MCM structure and mechanics: what we have learned from archaeal MCMSub-cellular biochemistry 62:89–111
- 15.Coupling of DNA binding and helicase activity is mediated by a conserved loop in the MCM proteinNucleic Acids Res 36:1309–1320
- 16.Intersubunit allosteric communication mediated by a conserved loop in the MCM helicaseProc Natl Acad Sci U S A 106:1051–1056
- 17.Mcm8 and Mcm9 form a complex that functions in homologous recombination repair induced by DNA interstrand crosslinksMol Cell 47:511–522
- 18.MCM8- and MCM9-deficient mice reveal gametogenesis defects and genome instability due to impaired homologous recombinationMol Cell 47:523–534
- 19.Inhibiting the MCM8-9 complex selectively sensitizes cancer cells to cisplatin and olaparibCancer Sci 110:1044–1053
- 20.Motifs of the C-terminal domain of MCM9 direct localization to sites of mitomycin-C damage for RAD51 recruitmentJ Biol Chem 296
- 21.The MCM8-MCM9 complex promotes RAD51 recruitment at DNA damage sites to facilitate homologous recombinationMol Cell Biol 33:1632–1644
- 22.MCM8-9 complex promotes resection of double-strand break ends by MRE11-RAD50-NBS1 complexNat Commun 6
- 23.Control of homologous recombination by the HROB-MCM8-MCM9 pathwayGenes Dev 33:1397–1415
- 24.MCM8IP activates the MCM8-9 helicase to promote DNA synthesis and homologous recombination upon DNA damageNat Commun 11
- 25.Crystal structure of a near-full-length archaeal MCM: functional insights for an AAA+ hexameric helicaseProc Natl Acad Sci U S A 105:20191–20196
- 26.The MCM8/9 complex: A recent recruit to the roster of helicases involved in genome maintenanceDNA Repair (Amst 76:1–10
- 27.Structural study of the N-terminal domain of human MCM8/9 complexStructure 29:1171–1181
- 28.The structure and function of MCM from archaeal M. ThermoautotrophicumNature structural biology 10:160–167
- 29.Structure of the eukaryotic replicative CMG helicase suggests a pumpjack motion for translocationNature structural & molecular biology 23:217–224
- 30.Processive DNA helicase activity of the minichromosome maintenance proteins 4, 6, and 7 complex requires forked DNA structuresProc Natl Acad Sci U S A 98:54–59
- 31.MCM9 Is Required for Mammalian DNA Mismatch RepairMol Cell 59:831–839
- 32.The ring-shaped hexameric helicases that function at DNA replication forksNature structural & molecular biology 25:122–130
- 33.Multi-body Refinement of Cryo-EM Images in RELIONMethods in molecular biology (Clifton, N.J.) 2215:145–160
- 34.A 3.8 Å resolution cryo-EM structure of a small protein bound to an imaging scaffoldNat Commun 10
- 35.Structure and activation mechanism of the BBSome membrane protein trafficking complexeLife 9
- 36.The dimeric Mcm8-9 complex of Xenopus laevis likely has a conserved function for resistance to DNA damageCell Cycle 12:1338–1339
- 37.MCM8 is an MCM2-7-related protein that functions as a DNA helicase during replication elongation and not initiationCell 120:315–328
- 38.pH-controlled quaternary states of hexameric DnaB helicaseJ Mol Biol 303:383–393
- 39.Flexibility of the rings: structural asymmetry in the DnaB hexameric helicaseJ Mol Biol 321:839–849
- 40.Quaternary polymorphism of replicative helicase G40P: structural mapping and domain rearrangementJ Mol Biol 357:1063–1076
- 41.Cryo-EM Structure of the Fork Protection Complex Bound to CMG at a Replication ForkMol Cell 78:926–940
- 42.Structure of BAI1/ELMO2 complex reveals an action mechanism of adhesion GPCRs via ELMO family scaffoldsNat Commun 10
- 43.Automated electron microscope tomography using robust prediction of specimen movementsJ Struct Biol 152:36–51
- 44.cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determinationNat Methods 14:290–296
- 45.Non-uniform refinement: adaptive regularization improves single-particle cryo-EM reconstructionNat Methods 17:1214–1221
- 46.AlphaFold Protein Structure Database: massively expanding the structural coverage of protein-sequence space with high-accuracy modelsNucleic Acids Res 50:D439–d444
- 47.Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in Phenix. Acta crystallographica. Section DStructural biology 75:861–877
- 48.RELION: implementation of a Bayesian approach to cryo-EM structure determinationJ Struct Biol 180:519–530
- 49.MotionCor2: anisotropic correction of beam-induced motion for improved cryo-electron microscopyNat Methods 14:331–332
- 50.Gctf: Real-time CTF determination and correctionJ Struct Biol 193:1–12
- 51.Characterisation of molecular motions in cryo-EM single-particle data by multi-body refinement in RELIONeLife 7
- 52.UCSF Chimera--a visualization system for exploratory research and analysisJ Comput Chem 25:1605–1612
- 53.Features and development of CootActa Crystallogr D Biol Crystallogr 66:486–501
- 54.Towards automated crystallographic structure refinement with phenix.refineActa Crystallogr D Biol Crystallogr 68:352–367