- Open Access
A new paradigm of DNA synthesis: three-metal-ion catalysis
Cell & Biosciencevolume 6, Article number: 51 (2016)
The Erratum to this article has been published in Cell & Bioscience 2017 7:32
Enzyme catalysis has been studied for over a century. How it actually occurs has not been visualized until recently. By combining in crystallo reaction and X-ray diffraction analysis of reaction intermediates, we have obtained unprecedented atomic details of the DNA synthesis process. Contrary to the established theory that enzyme-substrate complexes and transition states have identical atomic composition and catalysis occurs by the two-metal-ion mechanism, we have discovered that an additional divalent cation has to be captured en route to product formation. Unlike the canonical two metal ions, which are coordinated by DNA polymerases, this third metal ion is free of enzyme coordination. Its location between the α- and β-phosphates of dNTP suggests that the third metal ion may drive the phosphoryltransfer from the leaving group opposite to the 3′-OH nucleophile. Experimental data indicate that binding of the third metal ion may be the rate-limiting step in DNA synthesis and the free energy associated with the metal-ion binding can overcome the activation barrier to the DNA synthesis reaction.
In 1926 James B. Sumner purified and crystallized the first enzyme urease from jack bean and found that the crystalline form of urease catalyzed the breakdown of urea to ammonium and carbon dioxide . For this feat Sumner received the Nobel Prize in Chemistry in 1946. As early as 1835, the concept of catalysis was suggested by Jacob Berzelius, and in 1935 Henry Eyring, Meredith Gwynne Evans and Michael Polanyi proposed the transition state theory (TST), in which reactants and activated transition state complexes co-exist in a quasi-equilibrium. Based on the TST, activation energies for enzyme catalysis are routinely calculated from experimentally observed reaction rates . It has been taken for granted that enzymes catalyze reactions by stabilizing the transition state and thus reducing the activation energy and accelerating the reaction rate.
Owing to the transient nature of the transition state and technical difficulties of obtaining precise temporal and spatial details of a dynamic process, the exact details of any chemical reaction and how enzymes reduce energy barriers without consuming anything remain unknown. The questions regarding where the activation energy for catalysis, reduced yet still necessary, comes from and how enzymes stabilize the transition state have remained unanswered even after decades of experiments and theoretical calculations. Moreover, in practice many attempts to create artificial enzymes based on the assumed catalytic role of stabilizing transition states have met little success [3–5].
Kinetic rate analyses, particularly at pre-steady state, have generated a wealth of information about enzyme catalysis, intermediate steps and rates of their occurrence [6, 7]. However, various tricks have to be applied to alter reaction processes so that individual intermediate steps can be analyzed and their occurrence rate measured [8–13]. Similar to kinetic studies, to obtain three-dimensional structures of enzyme-substrate complexes, chemical reactions have to be manipulated and stopped by using non-reactive substrate mimics, enzymes inactivated by mutations, or non-permissible cofactors [14–19]. Assuming identical chemical compositions throughout the reaction process as suggested by the transition-state theory, transition states and reaction processes are then reconstructed and modeled based on experimentally measured rates and trapped enzyme-substrate or enzyme-product complexes by rearranging atoms, protons and electrons [20–23].
In recent years, we have developed a new approach to study the kinetics of DNA synthesis reactions and directly visualize reaction intermediates by in crystallo reaction and time-resolved X-ray crystallography. Although at a relatively low temporal resolution (tens of seconds), the actual reaction process is captured at atomic resolution . In this review article we present the surprising finding of a third metal ion that is absent in the enzyme-substrate complex and captured by substrates in the transition state . Binding of this additional third metal ion may provide sufficient activation energy to overcome the barrier to product formation.
Limitation of existing methods of studying catalysis
Enzymes and associated catalytic reactions have been studied in the pre-steady and steady states to yield catalytic requirements, intermediate steps and rates at high temporal resolutions. DNA synthesis reactions are no exception. DNA polymerases catalyze phosphoryltransfer reactions that incorporate dNTPs (A, G, T and C) according to a template sequence one at a time into DNA primer (Fig. 1a; Additional file 1: Movie S1). The reaction is of SN2 type and is pH dependent. All DNA polymerases depend on the metal ions Mg2+ or Mn2+ for catalysis, but differ dramatically in their catalytic rates. Many DNA polymerases undergo dNTP-dependent large conformational changes, while some translesion polymerases do not [12, 26–30]. Extensive kinetic and FRET studies have shown that the large conformational changes in DNA polymerases are faster than the chemical reaction itself, and therefore the phosphoryltransfer reaction is the most critical step in all DNA polymerases [9, 31–33]. After decades of investigation, it remains uncertain how a 3′-OH nucleophile is activated and whether the rate-limiting step of DNA synthesis is chemical or conformational [12, 32, 34–37].
For high spatial resolution, a large number of DNA polymerases have been crystallized in complex with a DNA template and primer pair and incoming dNTP. Catalysis was invariably circumvented by removal of the nucleophilic 3′-OH (replaced by 3′-deoxyribose) at the primer end, substitution of Mg2+ or Mn2+ by Ca2+, which does not support the catalysis, or substitution of dNTPs by non-reactive ddNTP or dNMPNPP analogs [14–19, 38]. Despite the diverse tertiary structures, catalytic rates, and degrees of fidelity in DNA polymerases, the core of the active site contains two Mg2+ or Mn2+, which are coordinated by the conserved carboxylates, DNA and dNTP substrate in the same configuration [19, 38]. How catalysis takes place was not visualized except for by QM/MM (hybrid quantum mechanics/molecular mechanics) modeling [20–23], which never fails to show that products can form by rearranging protons, electrons and atoms.
Kinetic studies of DNA synthesis have also observed two different binding constants of Mg2+ . But if more than one metal ion has a similar binding affinity, such metal-ion titration assays would lead to an under estimation. Because kinetic titration results agree with the number of metal ions found in the active site by X-ray crystallography [14–16, 26], two-metal-ion catalysis has become the widely accepted mechanism for DNA synthesis reactions. In the two-metal-ion catalysis model, one metal ion (B site) is associated with incoming dNTP and the second metal ion (A site) bridges the 3′-OH nucleophile and the incoming dNTP and is suggested to help deprotonating 3′-OH for the nucleophilic attack (Additional file 1: Movie S1). If incorrect or damaged, an incoming dNTP would be rejected before the A-site Mg2+ is recruited. Therefore, binding of two Mg2+ ions, which are notorious for the stringent coordination requirement, helps polymerases to achieve exceedingly high fidelity [39, 40]. As initially proposed by Steitz and Steitz , the DNA synthesis reaction is presumed to be promoted by the two Mg2+ ions, which align the substrates and facilitate the acid–base catalysis and the pentacovalent intermediate formation (Additional file 1: Movie S1).
In crystallo catalysis of DNA synthesis
Because DNA synthesis is pH and metal ion dependent, by reducing pH to 6.0 and using the non-catalytic metal-ion cofactor Ca2+, crystals of native DNA Pol η in complex with DNA substrate and correct incoming dNTP can be grown over a couple of weeks without reaction taking place . After fully grown, crystals were transferred to stabilization buffer to remove free reaction components and raise the pH to 7.0, which is permissible for the phosphoryltransfer reaction. To initiate DNA synthesis, crystals were exposed to Mg2+ or Mn2+, and after a certain reaction time they were cryo-cooled in liquid nitrogen to stop the reaction at the set time point. Ensuing analyses by X-ray diffraction have yielded unprecedented details following the reaction process of DNA synthesis . Binding of the A and B site Mg2+ is necessary to fully align the 3′-OH of DNA primer and incoming dNTP. But only after a 40 s delay after two metal ion binding, does the chemical reaction take place .
Interestingly, in crystallo DNA synthesis has also been achieved with DNA polymerase β. Crystals of DNA polymerase β complexed with DNA substrate (binary complexes) have been grown; however, unlike our Pol η studies, the incoming dNTP and Ca2+ were soaked in later . Similar to the studies of Pol η, catalysis was initiated by exposing Pol β crystals to Mg2+ or Mn2+ and terminated by cryo-cooling in liquid nitrogen. This approach has been applied successfully to study how Pol β incorporates an incorrect incoming nucleotide and oxidatively damaged 8-oxo-dG [42, 43]. Because the reaction rate of Pol β is much faster than Pol η, at the first time point of the data acquisition half of the substrate is already turned into product. Substrate alignment by two metal ions and the delayed chemical reaction were not observed with Pol β [42, 43].
A third metal ion is essential for DNA synthesis
In both in crystallo catalysis by Pol η and β, a third divalent cation coordinated by reaction products is observed (Fig. 1b–d). In the Pol β case, the third metal ion appears at the first time point (20–30 s) when approximately half of the substrate was converted to products. With Pol η, whose reaction rate is much slower, the third Mg2+ was detected ~60 s after the initial appearance of products, when the products were accumulated to 40 % . Because the third Mg2+ is coordinated by the DNA product and pyrophosphate leaving group and its binding is incompatible with the substrate state, it is dubbed “the product metal ion” [42, 43].
We were puzzled about the role of the third metal ion and wondered whether it is involved in catalysis or merely stabilizes the products and acts as a general acid to facilitate protonation of pyrophosphate . Because of the low electron number of Mg2+, it is not possible to detect Mg2+ by X-ray diffraction if its occupancy is lower than 30 %. Therefore, even if the third Mg2+ is required for the catalysis and product formation, it is not detected until the amount of product reaches more than 30 %.
To detect the third metal ion at a low occupancy, we switched the reaction buffer from Mg2+ to electron-rich Mn2+. To our delight, Mn2+ not only allowed detection of the third metal ion at low occupancy, but the affinities for Mn2+ at the two canonical (A and B) and the third (C) sites are significantly different . When in crystallo reaction took place at 10 mM Mn2+, the third Mn2+ appears at the same time as the reaction products and the two are perfectly correlated in appearance, time and occupancy. When the reaction was conducted in 1 mM Mn2+, however, the A and B sites were readily occupied, but the C-site was devoid of Mn2+. As a result, no products could be detected! Our in crystallo metal ion titration reveals the canonical A and B sites both have high-affinity for Mg2+ and Mn2+, and the C-site is of low affinity and thus determines the overall metal-ion concentration requirement for the DNA synthesis reaction [25, 37].
Binding of the third metal ion is rate limiting
A protein sidechain, R61 of Pol η, forms bifurcated salt bridges with the α and β phosphates of the incoming dNTP in the reactant state and overlaps with the C-site metal ion in the polymerase-product complexes if R61 does not change its rotamer conformation . One may imagine that removal of R61 would allow easier and faster binding of the third metal ion and lead to a faster reaction rate. When R61A mutation is introduced, however, the mutant Pol η catalyzes DNA synthesis at a slower rate than WT polymerase [19, 44, 45]. When the R61A mutant Pol η was examined in crystallo, we found that binding of the A and B site metal ions occurred within 40 s as for WT, but without R61 the incoming dNTP is misaligned relative to the 3′-OH by 0.3 Å. Consequently binding of the third metal ion was delayed by 120 s, thus leading to the slower catalytic rate than WT .
All DNA polymerases exhibited reduced catalytic rates when incoming dNTP is substituted with sulfur in the pro-Sp position (Sp-dNTPαS) [8, 34–36]. Before our discovery of the third metal ion, the Sp-sulfur was thought not to be involved in metal ion binding, and the degree of rate reduction has been interpreted as perturbing either chemistry or conformational change necessary for catalysis . But the pro-Sp and the α,β bridging oxygen atoms of dNTP are the only two potential non-water ligands of the third metal ion (Fig. 1). We predicted that sulfur substitution of Sp-dNTPαS therefore would retard the third metal ion binding and reduce the reaction rate.
Experimentally we find that the concentrations of Mg2+ and Mn2+ needed for incorporating Sp-dNTP αS by Pol η in solution are increased by ten and threefold, respectively . This is not surprising because the amount of Mg2+ or Mn2+ needed for the DNA synthesis reaction is determined by the C-site metal ion . The most dramatic changes in the in crystallo reaction are slow product formation and absence of detectable third metal ion even in 20 mM Mg2+ or Mn2+. In addition, although not directly involved in the A and B site metal ion binding, the Sp-sulfur perturbs the occupancy and location of Mg2+ in the A site  probably due to an altered electrostatic environment. Because the reaction-rate reduction caused by Sp-dNTP αS is less than fourfold [35, 36], it was concluded that the rate-limiting step in DNA synthesis is conformational changes. But the high concentrations of Mg2+ (up to 12.5 mM) used in these assays, however, masked the severe defects of the C-site metal ion binding. If 1 mM of Mg2+ had been used, the rate reduction would have been much greater, and thus leading to the opposite conclusion that chemistry is the rate-limiting step!
Our study reveals that binding of the C-site metal ion is essential and that its occurrence determines when the chemical reaction takes place. Binding of the two canonical metal ions at the A and B sites is a prerequisite for capturing the third metal ion but it is not rate limiting. Because delays of the C-site metal ion binding reduce the reaction rate as observed with R61A Pol η and Sp-dNTPαS substitution, we conclude that binding of the third divalent cation is the rate-limiting factor in the DNA synthesis reaction.
Solving the conundrum of when and how the third metal ion binds
The requirement of the third metal ion for the DNA synthesis reaction and its absence in the enzyme-substrate state present a conundrum—when and how the third metal ion binds. With nothing else to turn to, we hypothesized that thermal energy and thermal motion of the well-aligned reactants might provide a transient entrance for the third metal ion. To test the hypothesis, we divided the DNA synthesis reaction into two stages, A- and B-site metal ion binding along with alignment of DNA and dNTP substrate as the first, and capture of the third metal ion and product formation as the second. The first stage was carried out at 1 mM Mn2+. The second stage was conducted in the presence of 5 mM Mn2+ at various temperatures. Variation from 4 to 37 °C did not alter the diffusion rate of Mn2+ as evident in its unchanged binding at the canonical A site, but the increased thermal motions significantly enhanced the third metal ion binding and product formation .
The coordination of the third metal ion by two oxygen atoms of an incoming dNTP is non-ideal in distance and geometry compared to the preferred octahedral geometry of six inner-shell ligands of Mg2+ ion. But its coordination in the product state, when the phosphodiester bond between the α and β phosphates is broken, is nearly perfect octahedral. We suspect that the free- energy gain from binding of the third Mg2+ ion (freeing two inner-shell water ligands and association with dNTP) clears the barrier to the transition state. The stringent preference for octahedral coordination by Mg2+ ion may also drive the phosphotransfer reaction from breaking the existing bond in dNTP to forming a new bond between the 3′-OH and the α phosphate, which is in the reverse direction of the standard textbook version that is starting from the nucleophilic attack to the bond breakage (Fig. 2).
Comparison of the existing DNA polymerases, reverse transcriptases and RNA polymerases demonstrates that the active site compositions are highly conserved. Because both DNA pol η and pol β, which differ in tertiary structures, require the conserved third metal ion in catalysis, we propose that the third metal ion is a general feature and all polymerization reactions of nucleic acids occur by three-metal-ion catalysis (Additional file 2: Movie S2). The different environment surrounding the third metal ion, which is in a stark contrast to the conserved environment surrounding the canonical A and B site metal ions, gives hope that the third metal ion binding can be targeted for species and enzyme specific inhibition of DNA polymerases in treating infectious diseases and cancers.
Extended from kinetic and structural studies, which lay the ground work for understanding DNA synthesis reactions, our in crystallo reaction and time resolved X-ray diffraction analysis has led to the discovery of the third metal ion and its capture being the rate-limiting step of DNA synthesis reaction. This finding reveals the different chemical compositions between the transition state and the reactant state and therefore modifies the transition-state theory of enzyme catalysis. The third metal ion, which is free of enzyme coordination, is probably needed for the phosphoryltransfer reaction with or without enzyme catalysis. The catalytic role of DNA polymerases may not be to reduce the energy barrier (Fig. 3a), which would increase both forward and reverse reaction rates, but rather is to align the substrates perfectly thus increasing the probability for an incoming dNTP to capture the third and the catalytic divalent cation thus overcoming the barrier to product formation (Fig. 3b).
transition state theory
Sumner JB. The isolation and crystallization of the enzyme urease: preliminary paper. JBC. 1926;69:435–41.
Laidler KJ, King MC. Development of transition-state theory. J Phys Chem. 1983;87:2657–64.
Wagner J, Lerner RA, Barbas CF 3rd. Efficient aldolase catalytic antibodies that use the enamine mechanism of natural enzymes. Science. 1995;270(5243):1797–800.
Bolon DN, Mayo SL. Enzyme-like proteins by computational design. Proc Natl Acad Sci USA. 2001;98(25):14274–9.
Jiang L, Althoff EA, Clemente FR, Doyle L, Rothlisberger D, Zanghellini A, Gallaher JL, Betker JL, Tanaka F, Barbas CF 3rd, et al. De novo computational design of retro-aldol enzymes. Science. 2008;319(5868):1387–91.
Schechter AN. Measurement of fast biochemical reactions. Science. 1970;170(3955):273–80.
Johnson KA. Transient-state kinetic analysis of enzyme reaction pathways. Enzymes. 1992;20:1–61.
Patel SS, Wong I, Johnson KA. Pre-steady-state kinetic analysis of processive DNA replication including complete characterization of an exonuclease-deficient mutant. Biochemistry. 1991;30(2):511–25.
Shah AM, Li SX, Anderson KS, Sweasy JB. Y265H mutator mutant of DNA polymerase beta. Proper teometric alignment is critical for fidelity. J Biol Chem. 2001;276(14):10824–31.
Rothwell PJ, Mitaksov V, Waksman G. Motions of the fingers subdomain of klentaq1 are fast and not rate limiting: implications for the molecular basis of fidelity in DNA polymerases. Mol Cell. 2005;19(3):345–55.
Joyce CM, Benkovic SJ. DNA polymerase fidelity: kinetics, structure, and checkpoints. Biochemistry. 2004;43(45):14317–24.
Johnson KA. Role of induced fit in enzyme specificity: a molecular forward/reverse switch. J Biol Chem. 2008;283(39):26297–301.
Tsai MD. How DNA polymerases catalyze DNA replication, repair, and mutation. Biochemistry. 2014;53(17):2749–51.
Pelletier H, Sawaya MR, Kumar A, Wilson SH, Kraut J. Structures of ternary complexes of rat DNA polymerase beta, a DNA template-primer, and ddCTP. Science. 1994;264(5167):1891–903.
Doublie S, Tabor S, Long AM, Richardson CC, Ellenberger T. Crystal structure of a bacteriophage T7 DNA replication complex at 2.2 A resolution. Nature. 1998;391(6664):251–8.
Huang H, Chopra R, Verdine GL, Harrison SC. Structure of a covalently trapped catalytic complex of HIV-1 reverse transcriptase: implications for drug resistance. Science. 1998;282(5394):1669–75.
Ling H, Boudsocq F, Woodgate R, Yang W. Crystal structure of a Y-family DNA polymerase in action: a mechanism for error-prone and lesion-bypass replication. Cell. 2001;107(1):91–102.
Johnson SJ, Taylor JS, Beese LS. Processive DNA synthesis observed in a polymerase crystal suggests a mechanism for the prevention of frameshift mutations. Proc Natl Acad Sci USA. 2003;100(7):3895–900.
Biertumpfel C, Zhao Y, Kondo Y, Ramon-Maiques S, Gregory M, Lee JY, Masutani C, Lehmann AR, Hanaoka F, Yang W. Structure and mechanism of human DNA polymerase eta. Nature. 2010;465(7301):1044–8.
Florian J, Goodman MF, Warshel A. Computer simulations of protein functions: searching for the molecular origin of the replication fidelity of DNA polymerases. Proc Natl Acad Sci USA. 2005;102(19):6819–24.
Lin P, Pedersen LC, Batra VK, Beard WA, Wilson SH, Pedersen LG. Energy analysis of chemistry for correct insertion by DNA polymerase beta. Proc Natl Acad Sci USA. 2006;103(36):13294–9.
Wang L, Broyde S, Zhang Y. Polymerase-tailored variations in the water-mediated and substrate-assisted mechanism for nucleotidyl transfer: insights from a study of T7 DNA polymerase. J Mol Biol. 2009;389(4):787–96.
Lior-Hoffmann L, Wang L, Wang S, Geacintov NE, Broyde S, Zhang Y. Preferred WMSA catalytic mechanism of the nucleotidyl transfer reaction in human DNA polymerase kappa elucidates error-free bypass of a bulky DNA lesion. Nucleic Acids Res. 2012;40(18):9193–205.
Nakamura T, Zhao Y, Yamagata Y, Hua YJ, Yang W. Watching DNA polymerase eta make a phosphodiester bond. Nature. 2012;487(7406):196–201.
Gao Y, Yang W. Capture of a third Mg(2)(+) is essential for catalyzing DNA synthesis. Science. 2016;352(6291):1334–7.
Steitz TA. DNA polymerases: structural diversity and common mechanisms. J Biol Chem. 1999;274(25):17395–8.
Moon AF, Garcia-Diaz M, Batra VK, Beard WA, Bebenek K, Kunkel TA, Wilson SH, Pedersen LC. The X family portrait: structural insights into biological functions of X family polymerases. DNA Repair (Amst). 2007;6(12):1709–25.
Maxwell BA, Suo Z. Recent insight into the kinetic mechanisms and conformational dynamics of Y-Family DNA polymerases. Biochemistry. 2014;53(17):2804–14.
Xia S, Konigsberg WH. RB69 DNA polymerase structure, kinetics, and fidelity. Biochemistry. 2014;53(17):2752–67.
Yang W. An overview of Y-Family DNA polymerases and a case study of human DNA polymerase eta. Biochemistry. 2014;53(17):2793–803.
Showalter AK, Tsai MD. A reexamination of the nucleotide incorporation fidelity of DNA polymerases. Biochemistry. 2002;41(34):10571–6.
Rothwell PJ, Waksman G. Structure and mechanism of DNA polymerases. Adv Protein Chem. 2005;71:401–40.
Zhang H, Cao W, Zakharova E, Konigsberg W, De La Cruz EM. Fluorescence of 2-aminopurine reveals rapid conformational changes in the RB69 DNA polymerase-primer/template complexes upon binding and incorporation of matched deoxynucleoside triphosphates. Nucleic Acids Res. 2007;35(18):6052–62.
Dahlberg ME, Benkovic SJ. Kinetic mechanism of DNA polymerase I (Klenow fragment): identification of a second conformational change and evaluation of the internal equilibrium constant. Biochemistry. 1991;30(20):4835–43.
Washington MT, Prakash L, Prakash S. Yeast DNA polymerase eta utilizes an induced-fit mechanism of nucleotide incorporation. Cell. 2001;107(7):917–27.
Fiala KA, Suo Z. Mechanism of DNA polymerization catalyzed by Sulfolobus solfataricus P2 DNA polymerase IV. Biochemistry. 2004;43(7):2116–25.
Bakhtina M, Lee S, Wang Y, Dunlap C, Lamarche B, Tsai MD. Use of viscogens, dNTPalphaS, and rhodium(III) as probes in stopped-flow experiments to obtain new evidence for the mechanism of catalysis by DNA polymerase beta. Biochemistry. 2005;44(13):5177–87.
Batra VK, Beard WA, Shock DD, Krahn JM, Pedersen LC, Wilson SH. Magnesium-induced assembly of a complete DNA polymerase catalytic complex. Structure. 2006;14(4):757–66.
Yang W, Lee JY, Nowotny M. Making and breaking nucleic acids: two-Mg2+ -ion catalysis and substrate specificity. Mol Cell. 2006;22(1):5–13.
Yang W, Woodgate R. What a difference a decade makes: insights into translesion DNA synthesis. Proc Natl Acad Sci USA. 2007;104(40):15591–8.
Steitz TA, Steitz JA. A general two-metal-ion mechanism for catalytic RNA. Proc Natl Acad Sci USA. 1993;90(14):6498–502.
Freudenthal BD, Beard WA, Shock DD, Wilson SH. Observing a DNA polymerase choose right from wrong. Cell. 2013;154(1):157–68.
Freudenthal BD, Beard WA, Perera L, Shock DD, Kim T, Schlick T, Wilson SH. Uncovering the polymerase-induced cytotoxicity of an oxidized nucleotide. Nature. 2015;517(7536):635–9.
Zhao Y, Gregory MT, Biertumpfel C, Hua YJ, Hanaoka F, Yang W. Mechanism of somatic hypermutation at the WA motif by human DNA polymerase eta. Proc Natl Acad Sci USA. 2013;110(20):8146–51.
Su Y, Patra A, Harp JM, Egli M, Guengerich FP. Roles of Residues Arg-61 and Gln-38 of Human DNA Polymerase eta in Bypass of Deoxyguanosine and 7,8-Dihydro-8-oxo-2′-deoxyguanosine. J Biol Chem. 2015;290(26):15921–33.
Herschlag D, Piccirilli JA, Cech TR. Ribozyme-catalyzed and nonenzymatic reactions of phosphate diesters: rate effects upon substitution of sulfur for a nonbridging phosphoryl oxygen atom. Biochemistry. 1991;30(20):4844–54.
WY drafted this review, YG prepared the figures, and PJW prepared the videos. All authors read and approved the final manuscript.
This work is funded by NIH intramural program (DK036146-08, W.Y.). The authors thank Dr. R. Craigie for editing the manuscript.
No animal or human data is used. The authors declare that they have no competing interests.
This work is funded by NIH intramural program (DK036146-08, W.Y.).
An erratum to this article is available at http://dx.doi.org/10.1186/s13578-017-0159-1.