Enzyme kinetics of dUTPase from the planarian Dugesia ryukyuensis
BMC Research Notes volume 12, Article number: 163 (2019)
Planarians including Dugesia ryukyuensis (Dr) have strong regenerative abilities that require enhanced DNA replication. Knockdown of the DUT gene in Dr, which encodes deoxyuridine 5′-triphosphate pyrophosphatase (dUTPase), promotes DNA fragmentation, inhibits regeneration, and eventually leads to death. dUTPase catalyzes the hydrolysis of dUTP to dUMP and pyrophosphate. dUTPase is known to prevent uracil misincorporation in DNA by balancing the intracellular ratio between dUTP and dTTP, and contributes to genome stability. Nevertheless, the catalytic performance of Dr-dUTPase has not been reported.
To confirm the catalytic activity of Dr-dUTPase, we cloned and expressed Dr-DUT in E. coli. Then, we purified Dr-dUTPase using His-tag and removed the tag with thrombin. The resulting Dr-dUTPase had the leading peptide Gly–Ser–His– originating from the vector at the amino terminus, and a mutation, Arg66Lys, to remove the internal thrombin site. We observed the hydrolysis of dUTP by Dr-dUTPase using Cresol Red as a proton sensor. The Km for dUTP was determined to be 4.0 µM, which is similar to that for human dUTPase. Dr-dUTPase exhibited a preference for dUTP over the other nucleotides. We conclude the Dr-dUTPase has catalytic activity.
The planarian Dugesia ryukyuensis is a model system for studying reproductive strategies for survival [1,2,3]. We identified key genes in its sexualization and regeneration, including a DNA repair gene, Dr-Rad51 , and a germ cell regulator gene, Dr-nanos . Since these genes are also involved in genome stability, we hypothesized that the Dr-dut gene involved in nucleoside triphosphate biosynthesis  affects Dr-Rad51. The results of Dr-dut knockdown supported this hypothesis, upon which we observed alteration of the expression levels of Dr-rad51, Dr-rad51c, and the DNA damage response gene Dr-atm . Nevertheless, the catalytic performance of Dr-dUTPase has not been reported.
The DUT gene encodes the enzyme deoxyuridine 5′-triphosphate pyrophosphatase (dUTPase; EC 126.96.36.199) . Human (Homo sapiens, Hs) dUTPase is one of the most intensively studied enzymes  and is involved in pyrimidine metabolism. Specifically, Hs-dUTPase provides dUMP to thymidylate synthase  in the pathway of thymidine nucleotide biosynthesis . By hydrolyzing dUTP, dUTPase decreases the intracellular concentration of dUTP, and hence lowers the probability of uracil being incorporated into DNA [9, 12]. In the biomedical field, Hs-dUTPase has been identified as a new target of cancer chemotherapy, and a dUTPase inhibitor (TAS-114) is under clinical study [13, 14].
Two major isoforms of Hs-dUTPase have been reported. Isoform 2 (164 amino acids; Fig. 1A) is localized in the nucleus, while isoform 3 (252 amino acids) is localized in the mitochondria [15, 16]. The 3D structure of isoform 2 was solved, revealing that the enzyme adopts a homo-trimer form . However, in Dr-dUTPase, RNA sequencing revealed only one form of dUTPase .
Trimer dUTPases have three active sites at the subunit boundaries, consisting of five conserved motifs (Fig. 1A) [12, 19]. The Km of Hs-dUTPase for dUTP was reported to be 1.1 µM [20, 21]. To confirm the catalytic activity of Dr-dUTPase, we cloned, expressed, and purified the enzyme, and then measured its activity.
Materials and methods
Sequence analysis Sequences of trimer dUTPases from Dugesia ryukyuensis (Dr; DDBJ ID, LC421836), Homo sapiens (Hs; UniProt ID, P33316-2), Saccharomyces cerevisiae (Sc; P33317), Escherichia coli (Ec; P06968), and Plasmodium falciparum (Pf; Q8II92) were used because of the availability of kinetic data and 3D structures. Sequences were aligned by MAFFT  using the L-INS-I method  and BLOSUM62 scoring matrix .
Preparation of Dr-dUTPase The deduced amino acid sequence of Dr-dUTPase contained an internal thrombin cleavage site at Arg66 (Fig. 1A). To remove the thrombin cleavage site, the Dr-dut gene was mutated from G to A at position 197. Codon utilization of Dr-DUT was optimized to the expression host E. coli (Codon Optimization Tool, https://www.idtdna.com/codonopt; Additional file 1: Fig. S1, Table S1) and synthesized (Integrated DNA Technologies, Coralville, IA, USA). Dr-dut was cloned into the pET-15b vector (Novagen, San Diego, CA, USA) using NdeI and BamHI sites. Dr-dUTPase was prepared as described previously . Briefly, E. coli BL21 (DE3) was transformed by the construct and cultivated in LB medium. Expression of the Dr-dut gene was induced by IPTG. His-tagged Dr-dUTPase was purified using a Ni–NTA column (GE Healthcare, Chicago, IL, USA). Produced Dr-dUTPase contained the leading peptide MGSSHHHHHHSSGLVPRGSH at the amino terminus. The purified fraction (Additional file 1: Fig. S2) was subjected to desalting by dialysis against 1× thrombin buffer (150 mM NaCl, 1.5 mM CaCl2, 20 mM Tris–HCl, pH 8.0) before the thrombin cleavage. After the thrombin cleavage, the resulting Dr-dUTPase still contained the leading peptide GSH at the amino terminus.
Kinetic measurements The enzymatic activity was monitored by the color change of Cresol Red due to the production of protons by the hydrolysis of nucleotides  using a DU-640 UV/Vis spectrophotometer (Beckman Coulter, Brea, CA, USA). A total of 990 µl of base buffer (100 mM KCl, 5 mM MgCl2, 0.25 mM Bicine, pH 8.2) containing final concentrations of 10 µM dUTP + 25 µM Cresol Red was injected into a cuvette (optical path of 1 cm) containing 10 µl of 0.02 mM enzyme in the base buffer. Absorbance at 573 nm was recorded at intervals of 1.1 s (Additional file 1: Fig. S3). The recorded color change as a function of time was converted to the reaction product, and the Km values were calculated by the integrated Michaelis–Menten method, as described by Larsson et al.  and Inoguchi et al. .
Structural bioinformatics The 3D structures of dUTPase from Hs (PDB ID, 3ehw) and Pf (1vyq) were used. Homology modeling of intact Dr-dUTPase was performed on the SWISS-MODEL server . The system ranked Hs-dUTPase (3ehw) as the template at the top, and the established model obtained a QMEAN value of − 0.21 with the residue range between 7 and 145. The root-mean-square distance (RMSD) between the template and model was 0.1 Å in both monomer and trimer forms. Taking these findings together, the model of D. ryukyuensis dUTPase was judged to be usable for structural mining. Modeling for Dr-dUTPase Arg66Lys resulted in the same folding as with the native Dr-dUTPase, with RMSD against Hs-dUTPase in trimer of 0.1 Å. Structural mining and graphics preparation were performed using the PyMOL Molecular Graphics System, Version v188.8.131.52 (Schrödinger, LLC, New York City, NY, USA).
Amino acid sequence of Dr-dUTPase. The planarian Dr-dut gene  contained an ORF of 462 bp encoding a polypeptide of 153 amino acid residues (Fig. 1A; Additional file 1: Fig. S1). The deduced molecular mass and pI of Dr-dUTPase were 16,554 and 5.7, respectively. The amino acid sequence alignment revealed that Dr-dUTPase possesses the five conserved motifs unique to homotrimeric dUTPases, including Hs-dUTPase . Dr-dUTPase and Hs-dUTPase share sequence identity of 70%.
Production of Dr-dUTPase After the Dr-dut gene had been mutated to exclude the internal thrombin site, its codon utilization was optimized to E. coli (Additional file 1: Fig. S1). The mutation Arg66Lys in motif 2 was chosen to retain the positive charge of the side chain. The pI of mutated dUTPase was thus predicted to remain at 5.7. Alteration of the enzyme may be unavoidable because the residue is located next to Ser67, which is a key residue for the active site. The optimization gave evenly distributed usage of codon utilization (Additional file 1: Table S1), and increased the GC content of the Dr-dut gene from 37 to 50%. The genomic GC content of E. coli was reported to be 50% . The gene was cloned into the pET-15 vector and expressed in E. coli. After affinity purification (Additional file 1: Fig. S2A), the typical yield was about 20 mg of purified tagged protein from 250 ml of culture. Thrombin cleavage was performed to remove the His-tag. The resulting Dr-dUTPase contained the leading GSH peptide and was 156 amino acids in length. The deduced mass and pI were 16,807 and 5.9, respectively. The enzyme was purified uniformly showing the single band in the SDS-PAGE (Additional file 1: Fig. S2B). In the concentration by centrifugation using the 30 kDa-cut membrane (YM-30, MilliporeSigma, Burlington, MA, USA), Dr-dUTPase was retained. Therefore, Dr-dUTPase was potentially existed in oligomer form, including the trimer.
Enzymatic activity of Dr-dUTPase To assess the catalytic activity of Dr-dUTPase, hydrolysis of dUTP was monitored by the color change of Cresol Red. The change in absorbance upon multiple-turnover hydrolysis of dUTP as a function of time was converted to the amount of product at a certain concentration of substrate (Additional file 1: Fig. S3) . To obtain Km and Vmax, the reaction rate against the substrate concentration was fitted to the integrated Michaelis–Menten equation (Additional file 1: Fig. S3 inset). Dr-dUTPase exhibits enzymatic activity, with estimated Km and Vmax values of 4.0 μM and 20.2 μM s−1, respectively, at pH 8.2 and 25 °C (Table 1). kcat was estimated to be 3.4 × 10 s−1, assuming three active sites per trimer of dUTPase. The specificity constant, kcat/Km, for dUTP was thus 8.5 × 106 M−1 s−1. Dr-dUTPase exhibited catalytic activity for dUTP.
To analyze the substrate specificity of Dr-dUTPase, enzyme assays were performed using dATP, dTTP, dCTP, and dGTP as substrates at a final concentration of 100 µM and 4–6 µM dUTPase (Table 1). The Km obtained for the tested nucleotides ranged between 94 and 210 (Table 1). However, the performance of Dr-dUTPase for nucleotides other than dUTP was 2 orders of magnitude lower. Dr-dUTPase in the present format prefers dUTP as its substrate, while allowing the hydrolysis of dATP, dTTP, dCTP, and dGTP at lower levels.
To test the metal dependence of Dr-dUTPase, an enzyme assay was performed after the dialysis, without adding Mg2+. In the dialyzed enzyme, Dr-dUTPase activity was not detected, but the activity was regained as Mg2+ was added to the system. However, the Km value remained unchanged at MgCl2 concentrations from 5 to 25 mM. When the magnesium was replaced with 5 mM MnCl2, Dr-dUTPase maintained its activity (Table 1). When magnesium was replaced with 5 mM CaCl2, enzymatic activity was not observed. These results suggest that the activity of Dr-dUTPase occurs in a manner dependent on the metal ions magnesium and manganese.
Feeding based knockdown of Dr-dut invited planarian death. To obtain further knowledge about the mechanism of the cell death, we analyzed the catalytic performance of Dr-dUTPase.
Dr-dUTPase exhibited catalytic activity, as other trimer dUTPases do (Table 2). Interestingly, its tolerance for other nucleotides was unexpectedly high. For instance, the difference in the specificity constant between dUTP and other nucleotides, and Dr-dUTPase and other dUTPases were 100-fold and about 1000-fold, respectively (Table 2).
Potential causes of the differences in biochemistry included the difference in the amino acid sequence in key residues in conserved motifs (M1-5; Fig. 1A). In the case of Hs-dUTPase, the active site for the A chain consists of Arg85, Ser86, and Gly87; for the B chain, Gly99, Asp102, and Gly110; and for the C chain Phe158  (Fig. 1A, B). In the present form of Dr-dUTPase, substitutions were found as follows: Arg85Lys in the A chain and Phe158Tyr in the C chain. Among them, Phe158 in Hs-dUTPase was suggested to promote stacking enabling uracil recognition . Phe158 in Hs-dUTPase is conserved in other dUTPases, including Sc-dUTPase, Ec-dUTPase, and Pf-dUTPase. By the substitution of Phe158Tyr, Dr-dUTPase potentially obtained alternative stacking between the aromatic ring in the residue and the uracil base. The alternative staking could have allowed the larger ligand tolerance.
In motif 2 of trimer dUTPase, there is only one short α-helix, which holds two amino acid facing active site (Fig. 1A, B). One is the nucleophile Ser, and another is the four residues after it. Combination of the amino acid between the Ser and another residue seemed affect the affinity of the enzyme in terms of the initial capture of the substrate–water complex . In hs-dUTPase and Sc-dUTPase have small hydrophobic residues, and they show a lower affinity then that of Ec-dUTPase which carry the positively charged residues at the corresponding site. In Dr-dUTPase has Ser and Leu on the α-helix. Our results of the activity assay showed Dr-dUTPase has lower affinity to dUTP than that of Hs-dUTPase (Table 2). Nevertheless, Pf-dUTPase had positively charge amino acid at the second position in the helix, this enzyme had a large insertion in the leading sequence of M2. Hence Pf-dUTPase is an exception case.
The present study shows that Dr-dUTPase has catalytic activity that preferentially favors dUTP, but also possesses high tolerance regarding substrate recognition. This result for Dr-dUTPase was unexpected because Hs-dUTPase can slightly recognize dTTP and dCTP, but not dATP and dGTP (Table 2). If this is also the case in vivo, the wide substrate acceptance of Dr-dUTPase would be notable as a planarian-specific phenomenon. Either way, the results of the present study confirm the catalytic activity of dUTPase, increasing our basic understanding of the roles of dUTPase in planarians.
Dr-dUTPase has conserved motif of trimer dUTPase (M1-5), and has conserved residues for the inter-subunit interactions. Then we assumed that Dr-dUTPase keeps a trimer structure. To conform the oligomer state of Dr-dUTPase, we performed the size exclusion chromatography, but unfortunately it could not succeed.
The 3D structures of dUTPase showed that Dr-Tyr139 (B-chain) and Hs-Phe158 (B-chain) in motif 5 (M5), are involved in stacking between the U-base and aromatic ring of the side chain. The Dr-Phe29 (C-chain) and Hs-Tyr48 (C-chain) are located in structural motif 1 (M1), which is a part of the active sites. The aspartate in M1, next to the Dr-Phe29 (C-chain) and Hs-Tyr48, interacts with active site water molecules to stabilize the Mg ion. Therefore, Dr-Phe29 (C-chain) and Hs-Tyr48 are still involved in the substrate binding. Nevertheless, we do not have the crystal structure yet to support this idea.
deoxyuridine 5′-triphosphate pyrophosphatase
Kawakatsu M, Oki I, Tamura S. Taxonomy and geographical-distribution of Dugesia-Japonica and D ryukyuensis in the Far-East. Hydrobiologia. 1995;305(1–3):55–61.
Hoshi M, Kobayashi K, Arioka S, Hase S, Matsumoto M. Switch from asexual to sexual reproduction in the planarian Dugesia ryukyuensis. Integr Comp Biol. 2003;43(2):242–6.
Chinone A, Nodono H, Matsumoto M. Triploid planarian reproduces truly bisexually with euploid gametes produced through a different meiotic system between sex. Chromosoma. 2014;123(3):265–72.
Chinone A, Matsumoto M. DrRad51 is required for chiasmata formation in meiosis in planarian Dugesia ryukyuensis. Mol Reprod Dev. 2014;81(5):409–21.
Nakagawa H, Ishizu H, Chinone A, Kobayashi K, Matsumoto M. The Dr-nanos gene is essential for germ cell specification in the planarian Dugesia ryukyuensis. Int J Dev Biol. 2012;56(1–3):165–71.
Hirmondo R, Lopata A, Suranyi EV, Vertessy BG, Toth J. Differential control of dNTP biosynthesis and genome integrity maintenance by the dUTPase superfamily enzymes. Sci Rep. 2017;7(1):6043.
Alam MS, Moriyama H, Matsumoto M. Inhibition of Dr-dut gene causes DNA damage in planarian. Mol Reprod Dev. 2018;85(3):188–96.
Bertani LE, Haggmark A, Reichard P. Synthesis of pyrimidine deoxyribonucleoside diphosphates with enzymes from Escherichia coli. J Biol Chem. 1961;236:PC67–8.
Mol CD, Harris JM, McIntosh EM, Tainer JA. Human dUTP pyrophosphatase: uracil recognition by a beta hairpin and active sites formed by three separate subunits. Structure. 1996;4(9):1077–92.
Anderson DD, Quintero CM, Stover PJ. Identification of a de novo thymidylate biosynthesis pathway in mammalian mitochondria. Proc Natl Acad Sci USA. 2011;108(37):15163–8.
Hardy LW, Finer-Moore JS, Montfort WR, Jones MO, Santi DV, Stroud RM. Atomic structure of thymidylate synthase: target for rational drug design. Science. 1987;235(4787):448–55.
Vertessy BG, Toth J. Keeping uracil out of DNA: physiological role, structure and catalytic mechanism of dUTPases. Acc Chem Res. 2009;42(1):97–106.
Saito K, Nagashima H, Noguchi K, Yoshisue K, Yokogawa T, Matsushima E, Tahara T, Takagi S. First-in-human, phase I dose-escalation study of single and multiple doses of a first-in-class enhancer of fluoropyrimidines, a dUTPase inhibitor (TAS-114) in healthy male volunteers. Cancer Chemother Pharmacol. 2014;73(3):577–83.
Yano W, Yokogawa T, Wakasa T, Yamamura K, Fujioka A, Yoshisue K, Matsushima E, Miyahara S, Miyakoshi H, Taguchi J, et al. TAS-114, a first-in-class dual dUTPase/DPD inhibitor, demonstrates potential to improve therapeutic efficacy of fluoropyrimidine-based chemotherapy. Mol Cancer Ther. 2018;17(8):1683–93.
Ladner RD, McNulty DE, Carr SA, Roberts GD, Caradonna SJ. Characterization of distinct nuclear and mitochondrial forms of human deoxyuridine triphosphate nucleotidohydrolase. J Biol Chem. 1996;271(13):7745–51.
Ladner RD, Carr SA, Huddleston MJ, McNulty DE, Caradonna SJ. Identification of a consensus cyclin-dependent kinase phosphorylation site unique to the nuclear form of human deoxyuridine triphosphate nucleotidohydrolase. J Biol Chem. 1996;271(13):7752–7.
Varga B, Barabas O, Kovari J, Toth J, Hunyadi-Gulyas E, Klement E, Medzihradszky KF, Tolgyesi F, Fidy J, Vertessy BG. Active site closure facilitates juxtaposition of reactant atoms for initiation of catalysis by human dUTPase. FEBS Lett. 2007;581(24):4783–8.
Ishizuka H, Maezawa T, Kawauchi J, Nodono H, Hirao Y, Nishimura O, Nakagawa H, Sekii K, Tasaka K, Tarui H, et al. The Dugesia ryukyuensis database as a molecular resource for studying switching of the reproductive system. Zool Sci. 2007;24(1):31–7.
Inoguchi N, Chaiseeda K, Yamanishi M, Kim MK, Jang Y, Bajaj M, Chia CP, Becker DF, Moriyama H. Structural insights into the mechanism defining substrate affinity in Arabidopsis thaliana dUTPase: the role of tryptophan 93 in ligand orientation. BMC Res Notes. 2015;8:784.
Larsson G, Nyman PO, Kvassman JO. Kinetic characterization of dUTPase from Escherichia coli. J Biol Chem. 1996;271(39):24010–6.
Quesada-Soriano I, Casas-Solvas JM, Recio E, Ruiz-Perez LM, Vargas-Berenguel A, Gonzalez-Pacanowska D, Garcia-Fuentes L. Kinetic properties and specificity of trimeric Plasmodium falciparum and human dUTPases. Biochimie. 2010;92(2):178–86.
Katoh K, Rozewicki J, Yamada KD. MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization. Brief Bioinform. 2017. https://doi.org/10.1093/bib/bbx108.
Katoh K, Kuma K, Toh H, Miyata T. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005;33(2):511–8.
Henikoff S, Henikoff JF. Amino acid substitution matrices from protein blocks. Proc Natl Acad Sci USA. 1992;89(22):10915–9.
Waterhouse A, Bertoni M, Bienert S, Studer G, Tauriello G, Gumienny R, Heer FT, de Beer TAP, Rempfer C, Bordoli L, et al. SWISS-MODEL: homology modelling of protein structures and complexes. Nucleic Acids Res. 2018;46(W1):W296–303.
Riley M, Abe T, Arnaud MB, Berlyn MK, Blattner FR, Chaudhuri RR, Glasner JD, Horiuchi T, Keseler IM, Kosuge T, et al. Escherichia coli K-12: a cooperatively developed annotation snapshot–2005. Nucleic Acids Res. 2006;34(1):1–9.
Pecsi I, Leveles I, Harmat V, Vertessy BG, Toth J. Aromatic stacking between nucleobase and enzyme promotes phosphate ester hydrolysis in dUTPase. Nucleic Acids Res. 2010;38(20):7179–86.
Tchigvintsev A, Singer AU, Flick R, Petit P, Brown G, Evdokimova E, Savchenko A, Yakunin AF. Structure and activity of the Saccharomyces cerevisiae dUTP pyrophosphatase DUT1, an essential housekeeping enzyme. Biochem J. 2011;437(2):243–53.
Bozoky Z, Rona G, Klement E, Medzihradszky KF, Merenyi G, Vertessy BG, Friedrich P. Calpain-catalyzed proteolysis of human dUTPase specifically removes the nuclear localization signal peptide. PLoS ONE. 2011;6(5):e19546.
MSA, HM and MM designed the study. MSA performed all experiments and wrote the paper. All authors read and approved the final manuscript.
The authors thank to Dr. Sumitaka Hase for their technical assistance.
Availability of data and materials
The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.
Consent to publish
Ethics approval and consent to participate
Keio University Research Grant for Young Researchers (Grant No. 10,2015) and Keio Leading-edge Laboratory of Science & Technology Research Grant for Ph.D. Program (Grant No. 116,2017), Keio University, Japan.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Nucleotide sequence of Dr-DUT gene. Shown sequences are RNA sequencing (Rseq), mutated (G197A), and codon optimized for E. coli. Fig. S2. Production of Dr-dUTPase analyzed by SDS-PAGE. Fig. S3. Hydrolysis of dUTP by Dr-dUTPase. Table S1. Codon utilization in Dr-DUT Rseq and COTEC.
About this article
Cite this article
Alam, M.S., Moriyama, H. & Matsumoto, M. Enzyme kinetics of dUTPase from the planarian Dugesia ryukyuensis. BMC Res Notes 12, 163 (2019). https://doi.org/10.1186/s13104-019-4191-6