Watching DNA polymerase [eegr] make a phosphodiester bond
Teruya Nakamura,Ye Zhao,Wei Yang,Teruya Nakamura,Yuriko Yamagata,Ye Zhao,Yue-jin Hua
DOI: https://doi.org/10.1038/nature11181
IF: 64.8
2012-01-01
Nature
Abstract:DNA synthesis has been extensively studied, but the chemical reaction itself has not been visualized. Here we follow the course of phosphodiester bond formation using time-resolved X-ray crystallography. Native human DNA polymerase η, DNA and dATP were co-crystallized at pH 6.0 without Mg2+. The polymerization reaction was initiated by exposing crystals to 1 mM Mg2+ at pH 7.0, and stopped by freezing at desired time points for structural analysis. The substrates and two Mg2+ ions are aligned within 40 s, but the bond formation is not evident until 80 s. From 80 to 300 s structures show a mixture of decreasing substrate and increasing product of the nucleotidyl-transfer reaction. Transient electron densities indicate that deprotonation and an accompanying C2′-endo to C3′-endo conversion of the nucleophile 3′-OH are rate limiting. A third Mg2+ ion, which arrives with the new bond and stabilizes the intermediate state, may be an unappreciated feature of the two-metal-ion mechanism. Formation and breakage of chemical bonds underlie all life processes. DNA replication, which is essential for cell proliferation, is but one example1, 2. In each reaction cycle, a dNTP complementary to the templating base is incorporated into DNA in a nucleotidyl-transfer reaction catalysed by a polymerase, during which a new bond is formed between the 3′-OH of the primer strand and the α-phosphate of the dNTP, and the phosphodiester bond between the α- and β-phosphates of dNTP is broken (Fig. 1a). As a result the primer strand is extended by one nucleotide and a pyrophosphate is released. This reaction has been shown to require two Mg2+ ions (A and B) and is inferred to be SN2-type, forming a pentacovalent phosphate intermediate3, 4. A similar two-metal-ion-dependent mechanism is believed to be shared by all DNA and RNA polymerases and many nucleases5. DNA synthesis has been analysed by kinetic measurements, dynamic simulations and structural studies2, 6, 7. Crystal structures of a number of DNA polymerases in substrate-bound forms have been determined using non-reactive substrate analogues or with Ca2+ instead of Mg2+ to prevent the nucleotidyl-transfer reaction5, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18. These structures show active-site configurations in an array of pre-reaction ground states. Pre-steady-state kinetic studies of several DNA polymerases indicate that, after large conformational changes induced by substrate binding, unspecified subtle changes in the active site are the rate-limiting step19, 20, 21, 22, 23, 24, 25, 26, 27. The actual process of DNA synthesis has, to our knowledge, never been visualized. Human DNA polymerase η (Pol η) is specialized in lesion bypass and has a preformed catalytic centre that undergoes limited conformational change during DNA synthesis28. Crystal structures of the catalytic domain of Pol η complexed with DNA and a non-reactive dNTP analogue poised for catalysis have been determined at 1.8 Å resolution29, 30. Taking advantage of the slow reaction rate and hindered post-reaction DNA translocation in crystallo, we follow the course of DNA synthesis by X-ray crystallography and report here the structural changes and transient elements that are associated with the reaction. The Pol η activity in a reaction buffer containing 0.9 mM Mg2+ and 0.1 mM Ca2+ was only one-third of that with 1 mM Mg2+ (Fig. 1b), indicating that both active-site metal ions have to be Mg2+ for the catalysis. In accordance with the bell-shaped pH-dependence curve characteristic for the acid–base catalysis17, 31, the catalytic rate of Pol η is extremely low at pH 6.0 and rises with increasing pH from 6.0 to 8.0 (Fig. 1c). To prevent the nucleotidyl-transfer reaction, crystals of Pol η–substrate complexes were grown at pH 6.0 with only one Ca2+ per protein–DNA–dATP complex (Methods). After soaking in MES (pH 6.0–7.2) or HEPES buffer (pH 7.0–7.5) with Na+ or K+ but no divalent cation for 5 to 30 min, these crystals maintained excellent diffraction up to pH 7.0, but decayed at pH 7.5 with reduced resolution and increased diffuse scattering. The structure of Pol η crystal equilibrated at pH 6.8 was refined at 1.50 Å. A Ca2+ ion clearly occupies the B metal-ion-binding site and is coordinated in the octahedral geometry by oxygen atoms of the dATP and the active-site residues (Fig. 1d). The A metal-ion-binding site has low occupancy of a monovalent cation. The occupancy of Na+ or K+ in the A site increases with the pH (Supplementary Fig. 1a, b), which is probably correlated with deprotonation of the active-site carboxylates. Pol η binds Na+ more readily than K+, particularly below pH 7.0, probably because Na+ is similar to Mg2+ in size and is smaller than K+. As expected, there is no nucleotidyl-transfer reaction, and the structure is termed the ground-state ternary complex. The ground-state structure is similar to the Pol η ternary complex with a non-reactive dATP analogue dAMPNPP (2′-deoxyadenosine-5′-[(α, β)-imido]triphosophate) and two Mg2+ ions (Protein Data Bank ID 3MR2)29. However, without an A-site divalent cation, the 3′-OH of the primer strand shifts away from the dATP and forms hydrogen bonds with the side chains of the S113 and D115 residues, and the active-site carboxylates D13 and E116 adopt different rotamer conformations (Fig. 1d). Replacement of dAMPNPP with dATP leads to stabilization of R61, which forms bidentate hydrogen bonds with the α-phosphate. The nucleotidyl transfer was initiated by transferring Pol η crystals to a pH 6.8 or 7.0 reaction buffer containing 1 mM Mg2+ but no dATP (Methods). After incubation at 293 K for 40 to 300 s, the reaction was terminated in ~40-s intervals by freezing crystals in liquid nitrogen at 77 K. Diffraction data were collected to 1.50–1.95 Å Bragg spacings (Supplementary Table 1). The reaction process was monitored by the electron density corresponding to the new chemical bond in the Fo − Fc map compared with the refined ground-state structure (Fig. 2a and Supplementary Fig. 2). The reaction time courses at pH 6.8 and 7.0 in crystallo (Fig. 2b) are approximately 20–100-fold slower than in solution (Fig. 1c), probably owing to the reduced thermal motion. By 40 s, the A site was fully occupied with Mg2+ (Fig. 2a and Supplementary Fig. 1c, d). About 50% of Ca2+ in the B site was replaced by Mg2+, and exchange of the remaining Ca2+ took place slowly (Supplementary Fig. 3). Binding of two Mg2+ ions leads to the alignment of the 3′-OH and dATP18, 32. The refined structure is nearly identical to that with dAMPNPP, except for tighter dATP coordination by R61 and closer proximity of the Mg2+ ions (3.4 Å apart versus 3.6 Å). As there is no sign of bond formation, the structure is termed the reactant state. Electron density corresponding to a new bond between the 3′-OH and α-phosphorus of dATP begun to emerge at 80 s, increased quickly over the next 60 s and reached a maximum after 200 s, when the reaction was 60–70% complete (Fig. 2b). The seemingly reduced rate after 140 s is probably due to the reverse reaction in crystallo. The slight decline of product after 250 s (Fig. 2b) is due to a sideway product translocation, which is clear with longer incubation (Supplementary Fig. 4a). To alleviate an impediment to proper DNA translocation by the crystal lattice, we replaced an AT base pair with a mismatch at the DNA end that formed lattice contacts (Methods and Supplementary Table 2). The DNA with a TG mismatch led to isomorphous crystals and proper translocation of the DNA product (Supplementary Fig. 4b). Interestingly, the time courses of nucleotidyl transfer in the AT and TG crystals are nearly identical and are unaffected by the lattice contacts (Fig. 2b). At the peak of chemical-bond formation between 200 and 250 s, the scissile phosphate can be refined in a pentacovalent transition state without restraints (Fig. 2c); however, the bond distances between the phosphorus and the attacking or leaving oxygen atoms are 2.2–2.5 Å (Fig. 2c), much longer than the expected 2.0 Å observed with transition-state mimics such as AlF4 or MgF3 (refs 33–35). Moreover, in the Fo − Fc map residual electron densities are observed around the new and scissile phosphodiester bonds (Fig. 2c). The same diffraction data, however, can be well fitted as a mixture of the reactant state and the product state immediately before and after the nucleotidyl transfer (Fig. 2d). Because the transition state is transient and unstable, the structures obtained between 80 and 300 s are refined as a mixture of the reactant state and product state at different ratios (Fig. 2b and Supplementary Table 1a). For instance, the 1.52 Å structure at 230 s consists of 40% substrate and 60% product. Between the reactant and product state, the protein, DNA and dATP are superimposable except for atoms in the reaction centre. Most notably, the α-phosphorus of dATP moves 1.4 Å along a straight line between the attacking and leaving oxygen atoms (separated by 4.6 Å; Fig. 2d). The shift of the α-phosphate is accompanied by alteration of residue R61 of Pol η, which flips away from the scissile phosphate and is replaced by a new metal ion and water molecules (Fig. 2a, d; see details below). On the primer strand, changes are confined to the 3′ nucleotide. The 3′-OH, together with the deoxyribose, moves towards the α-phosphate by 0.5 Å, and the sugar pucker changes from C2′-endo in the reactant state to C3′-endo in the product state (Fig. 2d). With the loss of the nucleophile and α-phosphate as ligands, the A-site Mg2+ dissociates in the product state, as evidenced by the declining occupancy (Supplementary Fig. 1c, d). Concomitantly, D13 assumes a second conformation and forms a hydrogen bond with K224. The C3′-endo (A form) conformation at the 3′ primer end was observed in ternary complexes with the A-, B- and X-family DNA polymerases and is thought to be important for forming a shallow minor groove for dNTP selection9, 11, 12, 14, 15, 17, 18, 32. Among the Y-family polymerases, the primer end has always been observed as C2′-endo28. In the Pol η ground-state and reactant-state complexes, the dATP and the nucleotide 5′ to the primer end have the A-form conformation, but only in the product-state structure does the primer end adopt the C3′-endo conformation to avoid clashes between its C2′ atom and the non-bridging oxygen of dATP during and immediately after the nucleophilic attack (Figs 2d and 3a). Because the electron density is weak for the sugar moiety at the 3′-primer end, we tested the effect of the A-form conformation by using a primer with a ribonucleotide at its 3′ end. The catalytic efficiency (kcat/Km) of Pol η is comparable whether a ribonucleotide or deoxyribonucleotide is at the primer end (Fig. 3b), as observed for DNA pol β (ref. 36), indicating that the A-form conformation is probably necessary for DNA synthesis in general. Two unexpected spheres of electron density appear in the Fo − Fc maps in the course of new bond formation (Fig. 2a). The first is within hydrogen-bonding distance of the 3′-OH. Its electron density peaks at 80 s and declines considerably after 140 s when the product state becomes prominent. Although absent in the ground state and the 40-second reactant-state structure, the electron density is superimposable with a water molecule observed in the Pol η–dAMPNPP ternary complex. This water molecule is hydrogen-bonded with the 3′-OH, the O4′ of dATP and another water molecule in the 80-second reactant-state structure, but it is incompatible with the C3′-endo conformation in the product state (Fig. 3a). The 3′-OH is also hydrogen-bonded to a water molecule in an X-family DNA polymerase crystallized with dUMPNPP (2′-deoxyuridine-5′-[(α, β)-imido]triphosophate)32. A water-mediated and substrate-assisted catalytic mechanism has previously been proposed37, 38. In addition to the water molecule, S113 is hydrogen-bonded with the 3′-OH in the ground state and may pick up the proton and pass it to E116 in the reactant state (Fig. 1d). However, the S113A-mutant Pol η retains 95% of kcat with a threefold increase of Km (Fig. 3b and Supplementary Table 3). Interestingly, mutations of the S113 equivalent in the A- and B-family DNA polymerases (highly conserved histidine and threonine, respectively) also have limited impact on catalytic efficiency39, 40, 41. We therefore propose that 3′-OH is prone to deprotonation when coordinated by Mg2+ and aligned with the incoming nucleotide4. The proton can be passed on to the transient water molecule and then to bulk solvent. When the 3′-end of the primer is a ribonucleotide, the O2′ may relay the proton out instead. An incoming dNTP most probably participates in deprotonation of the nucleophile because, with the nitrogen substitution in dAMPNPP, the water molecule is bound stably to the 3′-OH29, 30. The second emerging sphere of electron density is very close to the α-phosphate on the opposite side of the A- and B-site Mg2+ ions (Fig. 2a). It starts to appear at 140 s and intensifies with the reaction time. By 230 s with the structure refined at 1.52 Å resolution, the octahedral geometry and the short coordination distances indicate that this is a divalent cation and most probably Mg2+ (Fig. 4). This third metal ion is liganded by four water molecules, two of which occupy the space of the departed R61 (Fig. 2d). Two additional ligands are the leaving oxygen (bridging between the α- and β-phosphate) and the non-bridging oxygen of the α-phosphate (Fig. 4). This Mg2+ thus bridges the two reaction products destined to separate and prevents DNA translocation. In the TG crystals with 1 mM Mg2+ in the reaction buffer, the mixed reactant-state and product-state intermediate stays at equilibrium for up to 15 min, but with reduced Mg2+ (5 µM) and addition of dATP (5 µM), product release occurs much faster (Supplementary Information). The physiological role of the third metal ion is probably to neutralize the negative charge built up in the transition state and may also facilitate protonation of the pyrophosphate. Three metal ions have been structurally observed in a number of enzymes catalysing phosphoryl-transfer reactions, including alkaline phosphatase, P1 nuclease and endonuclease IV33, 42, 43, 44. In the case of endonuclease IV, the three Zn2+ ions observed in both the substrate and product state44 are reminiscent of the Mg2+ with Pol η (Supplementary Fig. 5). The two metal ions essential for catalysis flank the scissile phosphate on one side, and the third metal ion bridges the reaction products on the other. A similar arrangement of three metal ions was proposed for group-I introns based on chemical probing45. However, the third metal ion was absent in a crystal structure of a splicing intermediate, probably owing to its transient nature46. Previously, flash–freeze was used to trap covalent intermediates and conformational changes associated with a chemical reaction47, 48, 49. Here, we extend the technology to record DNA synthesis in real time and at atomic resolution (Fig. 5 and Supplementary Movie 1). The arrival of a water molecule to deprotonate the 3′-OH occurs only after two-Mg2+-ion-induced substrate alignment and in the presence of a correct dNTP but not a non-reactive dAMPNPP. Mg2+-dependent and water-mediated deprotonation may occur in all DNA polymerases. The third metal ion that replaces the side chain of R61 during the nucleotidyl transfer probably stabilizes the transition state and facilitates product release. In the one-metal-ion-dependent TraI-like nucleases and topoisomerases a conserved Lys or Arg also occupies the position equivalent to R61 in Pol η (ref. 50). The involvement of an additional transient metal ion may be a general feature of both the one- and two-metal-ion mechanism. S113A-mutant Pol η was made using QuikChange (Stratagene). Wild-type Pol η (residues 1–432) and S113A-mutant proteins were expressed and purified as described29. The ternary complexes of Pol η, DNA and dATP were mixed with Ca2+ at 1:1 molar ratio and concentrated to a protein concentration of ~3 mg ml−1. Oligonucleotides used for original (AT pair at the DNA end) and three ‘loosely packed’ crystals (TG or GT mispaired and an unpaired A) are shown in Supplementary Table 2. An unpaired nucleotide at the primer 5′ end was disordered in the original Pol η structures29, and was therefore removed to improve the diffraction quality. All crystals were obtained using the hanging-drop vapour-diffusion method against a reservoir solution containing 0.1 M MES (pH 6.0) and 9–17% (w/v) PEG2K-MME and streak seeding. Based on diffraction qualities, the complexes containing the original AT and TG-end DNA were selected for chemical reaction in crystallo. These crystals were first transferred and incubated in a pre-reaction buffer containing 0.1 M MES (pH 6.8–7.2), 5 μM dATP, 20% (w/v) PEG2K-MME and 1 mM dithiothreitol (DTT) for ~30 min. For the pH 6.8 time course (Fig. 2b, Supplementary Table 1b) and translocation in the TG crystal, the pre-reaction buffer also contained 5 µM MgCl2. The different pre-reaction buffer might cause the larger variations of the pH 6.8 reactions. Chemical reaction was initiated by transferring the crystals with a nylon loop into a reaction buffer containing 0.1 M MES (pH 6.8–7.2), 1 mM MgCl2, 20% (w/v) PEG2K-MME and 1 mM DTT. After incubation for a desired time point (40 s and intervals of ~40 s), the crystals were quickly dipped in a cryo-solution supplemented with 20% (w/v) glycerol and flash-cooled in liquid nitrogen. For in crystallo reactions, MES buffer was titrated by KOH and not NaOH to avoid binding of Na+ ion at the A site. To observe product release, TG crystals were soaked in 5 µM dATP and 5 µM MgCl2 for 30 min and chased with addition of 1 mM MgCl2 for 1–8 min. Diffraction data were collected at 100 K on beam lines 22ID and 22BM of the advanced photon source. Data were processed with HKL200051 or XDS52, and converted to structure factors by TRUNCATE53. All data were isomorphic in the P61 space group and processed with identical unit-cell parameters a = b = 98.8 Å and c = 82.4 Å for structure and electron-density comparison. All structural figures were drawn using PyMOL (http://www.pymol.org). To monitor the new bond formation (Fig. 2a), Fo − Fc maps are calculated with PHENIX54 using the Fo of 40 s, 80 s, 140 s and 230 s after addition of 1 mM Mg2+ (pH 7.0) and the Fc of the refined reactant-state structure (40 s) with the 3′-OH of the primer end, the A-site Mg2+ and the α-phosphate of dATP omitted. To determine the reaction time courses (Fig. 2b and Supplementary Figs 1–3), all data were scaled against and compared to the refined ground-state structure, and electron densities for the new bond (between the 3′-OH of the primer strand and the α-phosphorus of dATP) and the metal ions were measured in the Fo − Fc map (with the B-site metal ion omitted) at each time point using programs in CCP453. After refinement of the ground-state structure, the intermediate structures were refined as a mixture of reactant state and product state using PHENIX54 and Coot55. The ratio of reactant state and product state at each time point was determined according to the time course of the new bond formation (Fig. 2b) and difference (Fo − Fc) in electron-density maps. For high-resolution refinement (~1.5 Å resolution), the translation–libration–screw refinement was applied. Data collection and refinement statistics are summarized in Supplementary Table 1. Kinetic measurements including Km and kcat were done using the template/primer shown in Supplementary Table 3. The basic reaction mixture contained 2.5 nM Pol η, 5 μM 5′-fluorescein-labelled primer and template, and 0–80 µM dNTP in 40 mM Tris–HCl (pH 7.5), 5 mM MgCl2, 10 mM dithiothreitol, 100 mM KCl, 0.1 mg ml−1 bovine serum albumin and 5% glycerol. For the metal-ion-competition assay, Mg2+ was increased from 0 to 1.0 mM in 0.1-mM steps and Ca2+ was decreased to keep the combined concentration of 1.0 mM. The effects of pH were screened from pH 6.0 to pH 8.0. All reactions took place at room temperature for 4 or 8 min and were stopped by addition of 10× formamide loading buffer to the final concentrations of 8% formamide, 1 mM EDTA (pH 8.0) and 0.1 mg ml−1 xylene cyanol. After heating to 90 °C for 3 min and immediately placing on ice, products were resolved on 20% polyacrylamide sequencing gels containing 5.5 M urea. Quantification and curve fitting were carried out as described20. Km, kcat and (kcat/Km) of wild-type and S113A-mutant Pol η catalysing single-nucleotide incorporation using deoxyribonucleotide versus ribonucleotide primer are summarized in Supplementary Table 3 and plotted in Fig. 3b. Download references We thank D. Leahy, M. Gellert and R. Craigie for editing the manuscript. The research was supported by the intramural research program of the National Institute of Diabetes and Digestive and Kidney Diseases, National Institutes of Health (W.Y., T.N. and Y.Z.); the Japan Society for the Promotion of Science Institutional Program for Young Researcher Overseas visits, Kumamoto University, and the Kumayaku Alumni Research Fund (T.N.); Chinese Ministry of Education scholarship (Y.Z.); National Natural Science Foundation of China (Y.-J.H.); and Grant-in-Aid for Scientific Research from the Ministry of Education, Culture, Sports, Science, and Technology of Japan (Y.Y.). Atomic coordinates and structure factors for the reported crystal structures have been deposited with the Protein Data Bank with accession codes from 4ECQ to 4ECZ, 4ED0 to 4ED3, and 4ED6 to 4ED8.