Skip to main content
  • Research Note
  • Open access
  • Published:

Enzyme kinetics of deoxyuridine triphosphatase from Western corn rootworm



The Western corn rootworm (WCR), Diabrotica virgifera virgifera, is a highly adaptable insect pest that has evolved resistance to a variety of control strategies, including insecticides. Therefore, it is interesting to examine how housekeeping proteins in WCR have been changed under WCR-controlling strategies. In this study, we focused on one of such proteins in WCR, a ubiquitous enzyme 5'-triphosphate nucleotidohydrolase (dUTPase). In the thymidine synthetic pathway, dUTPase hydrolyzes deoxyuridine triphosphate (dUTP) and supplies the substrate, deoxyuridine monophosphate, for the thymidylate synthase (TS). It decreases the cellular content of uracil, reducing uracil misincorporation into DNA. Suppressing the dUTPase activity, therefore, contributes to thymineless death. In this study, we investigated the enzymatic properties of dUTPase.


The WCR dUTPase gene (DUT) was synthesized with the addition of His-tag corresponding DNA sequence and then cloned and expressed in Escherichia coli, and the protein product was purified. The product of WCR DUT hydrolyzed dUTP and was designated as dUTPase. WCR dUTPase did not hydrolyze dATP, dTTP, dCTP, or dGTP. WCR dUTPase was analyzed via size-exclusion chromatography and exhibited a molecular weight corresponding to that of trimer. The present format can be interpreted as nuclear trimer type. Possible isomers will be examined once transcriptome analyses are conducted.

Peer Review reports


Diabrotica virgifera virgifera LeConte, commonly known as Western corn rootworm (WCR), is a major corn pest in North America [1, 2]. WCR is highly adaptable and has evolved resistance to a variety of management strategies, including crop rotation, synthetic insecticides, and genetically modified corn expressing Cry proteins [3]. To understand the mechanisms and effects of acquiring such resistance, it is important to determine whether the functions of any housekeeping proteins have been affected in WCR during the resistance evolution.

The ubiquitous and essential housekeeping enzyme deoxyuridine 5'-triphosphate nucleotidohydrolase (dUTPase) is involved in thymine synthesis (Fig. 1a) [4,5,6,7]. dUTPase removes pyrophosphate from deoxyuridine triphosphate (dUTP) to produce deoxyuridine monophosphate (dUMP), which is the substrate for thymidylate synthase (TS) [8, 9]. TS then methylates dUMP to produce dTMP [10, 11]. dUTPase, therefore, plays two critical roles: production of dUMP, the deoxythymidine triphosphate (dTTP) precursor, and degradation of dUTP to dUMP, preventing uracil misincorporation into DNA. Thymine is required in tissues with active DNA synthesis. Therefore, the thymidylate synthesis pathway, for example, is the point of action for several anticancer agents in humans [12,13,14,15]. However, redundancy of the supply of the substrate dUMP often limits the anticancer effect of pyrimidine antagonists [16, 17]. 5-Fluorouracil (5FU) is a TS inhibitor and used for cancer therapy where a large amount of accumulated dUMP leads to uracil misincorporation [18, 19]. In planarians, administration of 5FU caused death from DNA fragmentation [20, 21].

Fig. 1
figure 1

DUT gene from Western corn rootworm. a Pyrimidine metabolism map adapted from KEGG map00240 [7]. Enzymes supplying dUMP that are identified in WCR and Methanococci are indicated in orange and green, respectively. b Alignment of the dUTPase protein sequences among D. v. virgifera, T. castaneum, D. melanogaster, H. sapiens, S. cerevisiae, and E. coli (see Materials and Methods for the accession numbers). The construct sequence used in this study is shown at the top, where the positions including the removable N-terminal 6-His tag and the thrombin cleavage sites are underlined. The five conserved motifs in dUTPase are highlighted in cyan (M1–M5; see also Additional file 1: Fig. S1a). Substituted amino acids within the conserved motifs are not highlighted. Nuclear localization signals reported in D. melanogaster [25] and H. sapiens [15] are highlighted in gray. dUTPase from D. melanogaster possesses a Drosophila-specific 28-residue segment at the C-terminal [25]. The location of Lys114 in D. v. virgifera and corresponding residues in other species are highlighted in green. c Model for one of the active sites (Additional file 1: Fig. S1b). A substrate and side chains of relevant amino acid were shown in stick model

Knockout of the dUTPase gene (DUT) in mice [22] and knockdown of DUT mRNA in planarian [20, 21] are shown to be fatal. In insects, DUT silencing efficiently killed fruit fly, Drosophila melanogaster, at the early pupal stage [23,24,25]. Furthermore, knockdown of DUT expression resulted in a 100% mortality in red flour beetle, Tribolium castaneum, at the larval stage [26]. For WCR, feeding of DUT dsRNA to neonates for 9 days killed 54% of the larvae and inhibited the growth of 80% of the survivors [26].

In this study, the WCR DUT gene was synthesized and dUTPase was produced in E. coli and purified. Furthermore, the quaternary structure and enzymatic kinetics of WCR dUTPase were analyzed.

Materials and methods

This study was conducted under the oversight of the Institutional Biosafety Committee at the University of Nebraska–Lincoln, Protocol Number 174.

dUTPase genes and proteins

Genomic sequence (GenBank, NW_021039130.1) of Diabrotica virgifera virgifera LeConte (WCR) contained dUTPase gene (DUT; mRNA, XM_028280744.1; protein, XP_028136545). The dUTPase protein sequences used in this study are Tribolium castaneum (EFA05862.1), Drosophila melanogaster (Q9V3I1), Saccharomyces cerevisiae (P33317), Homo sapiens (P33316), and E. coli (strain K12; P06968).

WCR DUT gene construct

For the construct design, SWISS-MODEL [27,28,29,30] and ProtParam [31] were used for molecular modeling and protein parameter calculation, respectively (Additional file 1: Fig. S1). A DNA fragment encoding WCR dUTPase was synthesized as follows. An N-terminal 6-His tag was introduced to enable metal affinity chromatography purification (Fig. 1b, shown as “construct”). The construct started at the 25th residue of the WCR dUTPase to maintain the ability of the subunits to assemble in trimers based on sequence comparisons and 3D structural modeling (Additional file 1: Fig. S1ab and 1c). To prevent thrombin from cleaving other than the 6-His tag (secondary cleavage), the 89th arginine (R) residue, which is located in a helix and slightly away from the surface, was substituted to a lysine (K) with a G/A mutation (Fig. 1b and Additional file 1: Fig. S1b). Although the Arg89 plays a role in holding substrates, metals, and waters, lysine89 should be able to retain this function as it also has a positive charge. Arg89Lys mutation has also been shown to maintain the activity in the planarian dUTPase [20]. After optimizing the codon utilization for E. coli genes, a DNA fragment encoding the WCR DUT was synthesized by Gene Script (Piscataway, NJ, USA; Additional file 2: Fig. S2). The DUT gene was cloned into pET-15 (Novagen, Madison, WI, USA) using the NcoI and XhoI sites.

Production of dUTPase and size analysis

The WCR dUTPase was prepared as previously described [32,33,34] (Fig. 2). After the chromatographies using Ni–NTA Resin (New England BioLabs, Ipswich, MA, USA) and then HiTrap Q (GE HealthCare, Chicago, IL, USA), approximately 10 mg of WCR dUTPase was purified from 2.5 g of cells. To remove the His-tag, tagged dUTPase was treated with human α-thrombin (Hematologic Technologies, Essex Junction, VT, USA) after dialysis against 20-mM Tris–HCl and 100-mM NaCl (pH 8). His-tagged proteins were removed by the Ni–NTA column, and protein was purified using Benzamidine Sepharose (GE HealthCare). Proteinase inhibitors were used, including 1-mM PMSF and 0.1-mg/mL benzamidine. Tryptic digestion and mass spectrometry confirmed the produced molecules [35,36,37,38] (Additional file 3: Fig. S3). To estimate the molecular weight of WCR dUTPase, size-exclusion chromatography was performed using Superdex Matrix (GE HealthCare). The molecular weight was calibrated using the SEC standard (BioRad, Hercules, CA, USA).

Fig. 2
figure 2

Purification, quaternary structure, and enzyme activities of WCR dUTPase. a 18% SDS-PAGE. The molecular weights of 6-histidine-tagged dUTPase and dUTPase after being cleaved by thrombin are 17.5 and 15.6 kDa, respectively. b Size-exclusion chromatography. The estimated molecular weight of tag-free dUTPase was 48 kDa. c Hydrolysis of dUTP by dUTPase. Gray dots indicate the observed drop in absorbance and the predicted regression line obtained from a corresponding scan. The inset indicates linear transformation of the data between the arrows according to the integrated Michaelis–Menten equation and the corresponding regression line. d–i. Substrate specificity of WCR dUTPase. In each plot, the orange and blue lines show when a particular substrate is added and no enzyme is added, respectively. The substrates used are dUTP (d and e), dATP (f). dTTP (g), dCTP (h), and dGTP (i). The enzyme is preincubated with 5-mM EDTA in e

Enzyme activity assay

Kinetic assay of WCR dUTPase was conducted using the cresol red method according to Larsson et al. [39] (Figs. 2c–i) on a high-speed spectrophotometer (Hi-Tech SF-61DX2, TgK Scientific, Bradford-on-Avon, UK). Absorbance data from each run was fitted to the function:


where x and y denote time and absorbance, respectively. The predicted points were used to calculate the concentration of product formation at time t, [P]t, using the following equation:

$${[P]}_{t}= \frac{{A}_{0}- {A}_{t}}{{A}_{0}- {A}_{\infty }}{[S]}_{0}$$

where [S]0 denotes the substrate concentration, and A0, At, and A denote absorbance at the start, time t, and end of the reaction, respectively. The time course reaction was calculated using the integrated Michaelis–Menten equation.

$$\frac{{[P]}_{t}}{t}=V- \frac{K}{t}\mathrm{ln}\frac{{[S]}_{0}}{{[S]}_{0}- {[P]}_{t}}$$

where K and V denote KM and Vmax, respectively. Equations 2 and 3 were used to plot the graph of \(\frac{{[P]}_{t}}{t}\) against \(\frac{\mathrm{ln}\frac{{[S]}_{0}}{{[S]}_{0}-[{S]}_{t}}}{t}\). Linear regression was used to calculate the slope that gives the KM for each experiment. Kcat and Kcat/KM were calculated using the propagation of error approach. Three experiments were performed under the same conditions and the average value was reported. All calculations were performed using R v3.5.3 (R Foundation, Vienna, Austria; /


dUTPase production

There are monomeric, dimeric, and trimeric dUTPases [40]. The WCR dUTPase protein sequence was compared with the trimeric dUTPases (Fig. 1b). Five conserved motifs characteristic to the trimeric dUTPases (M1–M5) [40] were also found in WCR dUTPase (Figs. 1b and Additional file 1: Fig. S1a). We synthesized and cloned a DNA fragment encoding WCR dUTPase (shown as “construct” in Figs. 1b and Additional file 2: Fig. S2). After purification, SDS-PAGE revealed that the purified protein has a molecular weight of 15.6 kDa (Fig. 2a). Tryptic mass spectrometry confirmed the identity of the purified protein (Additional file 3: Fig. S3). A total of 32 fragments were exclusive to the amino acid sequence of WCR dUTPase. Size-exclusion chromatography revealed that the molecular weight of WCR dUTPase was 48 kDa (Fig. 2b). These observations indicate that WCR dUTPase forms a trimer, which is consistent with that indicated by sequence similarity and the five conserved motifs as described above.

dUTPase activity

Cresol red assay [41] revealed that the purified dUTPase has enzymatic activity (Fig. 2c–d). The addition of EDTA at the final concentration of 0.5 mM prevented the color change of cresol red (Fig. 2e). We also observed no significant color change when the substrate was changed to either dATP, dTTP, dCTP, or dGTP (Fig. 2f–i). We recorded the reaction trace with enzyme and dUTP concentrations of 50 nM and 1 mM, respectively (Fig. 1c), and estimated KM as 0.7 ± 0.1 μM (p < 0.01 for curve fitting and linear regression; Table 1).

Table 1 Kinetics of dUTPase of different species against dUTP

Comparison of dUTPase activity between organisms

The KM value for WCR dUTPase was estimated to be 0.7 μM, whereas those reported for other organisms were between 0.4 and 3.6 μM (Table 1). WCR dUTPase had a higher specificity constant (kcat/KM) than other eukaryotes. WCR dUTPase exhibited a strict preference for dUTP. However, human [8] and D. melanogaster dUTPases have exhibited slight activities against dTTP and dCTP [42, 45]. The multiple alignment and 3D structural modeling of WCR dUTPase showed that one of the amino acids within the M3 conserved motif, Lys114, is different from the corresponding amino acids in D. melanogaster and human dUTPases (Gly110 and Gly106, respectively; Fig. 1b). These conserved motifs are involved in the catalysis and interaction with the deoxyribose ring [43]. The presence of a Lys residue in the M3 motif can reduce the flexibility of the loop at this position due to possible additional hydrogen bonds. Thus, it potentially has a smaller chance of accepting a nucleotide other than dUTP.

The effect of sugar moiety on substrate specificity is evaluated as follows. In the dUTPase models from human (3ehw) and WCR, the C3–OH of 2'-deoxyribose forms a hydrogen bond with Asp102 (106 in WCR dUTPase) and C2 facing Tyr105 (109 in WCR dUTPase; Fig. 1c). This configuration seemed too narrow to accommodate C2–OH; repulsion due to pi and anion arrangement is also expected. In yeast dUTPase, the Tyr88Ala mutant enzyme with reduced steric hindrance has been reported to have equivalent reactivity toward both dUTP and UTP [46]. Side chains corresponding to Arg86, Asp87, Thr89, Glu91 in yeast dUTPase are toward outward from the active site. E. coli dUTPase has been reported to be active against UTP [47]. In M3, Val101 of human dUTPase (105Val in WCR dUTPase) is changed to Ile89 E. coli dUTPase. WCR dUTPase seemed to have difficulty accepting DUT as a substrate.


Previous knockdown of DUT expression was reported to result in 89%–100% mortality in T. castaneum larvae (via injection of 150 nL of 0.01–250-ng/μL dsDUT) [26]. In the same study, feeding assay with 500 ng/cm2 of dsDUT revealed only 54% mortality in WCR neonates at 9 days [26].

As illustrated in the metabolic pathway presented in Fig. 1a, four enzymes besides dUTPase (UniProt ID, A0A6P7FUX9_DIAVI) can produce dUMP (cf. KEGG map 00240). Among the four enzymes, potential WCR homologs were identified for dCMP deaminase (A0A6P7GGP3_DIAVI), dTMP kinase (A0A6P7FXJ0_DIAVI), and deoxyribonucleotide kinase (dNK, A0A6P7FA76_DIAVI). Those three enzymes were also observed in both D. melanogaster and T. castaneum. dUMP-Forming dCTP deaminase is observed in Methanococci but not in the three aforementioned insects. It is noteworthy that while mammals have two types of thymidine kinases (TK1-like and non-TK1-like [48]), insects have only one multisubstrate enzyme, dNK. Phylogenetic analyses revealed that insect dNK might have evolved from a more specialized TK2 (non-TK1-like) enzyme [49]. Thus, higher TK activity does not interfere with the redundant supply of dUMP described below.

The redundant dUMP supply may explain why suppression of dUTPase expression had limited effect on WCR mortality. Independent of the dUMP availability, administration of 5FU (KEGG, map00983), a TS inhibitor commonly used to treat cancer in humans [50], is expected to suppress thymine production (Fig. 1a). More investigations on dUTPase function, gene repertoire, and 5FU metabolism during the development of WCR and other coleopterans are warranted.


This study found that WCR DUT encodes a uracil-specific pyrophosphatase, which could explain the physiological effects of WCR DUT knockdown [26]. Future research is warranted to fully characterize WCR DUT.

The identification, confirmation, and characterization of possible isoforms of dUTPase are necessary to understand its structure and function. The levels of glycosylation and phosphorylation should also be considered, as these can affect the activity of the protein. Research on dUTPase in humans is ongoing, with the goal of developing new cancer chemotherapeutic agents and malaria treatments. Rácz et al. [51] have identified two additional isoforms of dUTPase in humans, DUT-N and DUT-M, which are localized to the nucleus and mitochondria, respectively. Future studies are warranted to investigate a molecular species that remains in the cytoplasm, which may represent a novel isoform of WCR dUTPase.

A study of dUTPase levels post-knockdown is key to understanding dUTPase inhibitors and pest control potential. TAS-114 [52] is a potent dUTPase inhibitor with chemotherapy applications. Although its inhibitory activity against WCR dUTPase has not been tested, it is hoped that similar environmentally friendly compounds will be developed.

Availability of data and materials

The datasets and plasmid DNA constructs used in this study are available from the corresponding author upon reasonable request.


  1. Bazok R, Lemic D, Chiarini F, Furlan L. Western corn rootworm (Diabrotica virgifera virgifera LeConte) in Europe: current status and sustainable pest management. Insects. 2021;12:195.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Gray ME, Sappington TW, Miller NJ, Moeser J, Bohn MO. Adaptation and invasiveness of western corn rootworm: intensifying research on a worsening pest. Annu Rev Entomol. 2009;54:303–21.

    Article  CAS  PubMed  Google Scholar 

  3. Kuwar SS, Mishra R, Banerjee R, Milligan J, Rydel T, Du Z, Xie Z, Ivashuta S, Kouadio JL, Meyer JM, Bonning BC. Engineering of Cry3Bb1 provides mechanistic insights toward countering western corn rootworm resistance. Curr Res Insect Sci. 2022;2: 100033.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Vertessy BG. Flexible glycine rich motif of Escherichia coli deoxyuridine triphosphate nucleotidohydrolase is important for functional but not for structural integrity of the enzyme. Proteins. 1997;28:568–79.

    Article  CAS  PubMed  Google Scholar 

  5. 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:4783–8.

    Article  CAS  PubMed  Google Scholar 

  6. Zhang Y, Moriyama H, Homma K, Van Etten JL. Chlorella virus-encoded deoxyuridine triphosphatases exhibit different temperature optima. J Virol. 2005;79:9945–53.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Kanehisa M, Furumichi M, Sato Y, Ishiguro-Watanabe M, Tanabe M. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49:D545–51.

    Article  CAS  PubMed  Google Scholar 

  8. Toth J, Varga B, Kovacs M, Malnasi-Csizmadia A, Vertessy BG. Kinetic mechanism of human dUTPase, an essential nucleotide pyrophosphatase enzyme. J Biol Chem. 2007;282:33572–82.

    Article  CAS  PubMed  Google Scholar 

  9. Islam Z, Gurevic I, Strutzenberg TS, Ghosh AK, Iqbal T, Kohen A. Bacterial versus human thymidylate synthase: kinetics and functionality. PLoS ONE. 2018;13: e0196506.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Anderson DD, Quintero CM, Stover PJ. Identification of a de novo thymidylate biosynthesis pathway in mammalian mitochondria. Proc Natl Acad Sci USA. 2011;108:15163–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Reyes CL, Sage CR, Rutenber EE, Nissen RM, Finer-Moore JS, Stroud RM. Inactivity of N229A thymidylate synthase due to water-mediated effects: isolating a late stage in methyl transfer. J Mol Biol. 1998;284:699–712.

    Article  CAS  PubMed  Google Scholar 

  12. Chon J, Stover PJ, Field MS. Targeting nuclear thymidylate biosynthesis. Mol Aspects Med. 2017;53:48–56.

    Article  CAS  PubMed  Google Scholar 

  13. Martin DS, Stolfi RL, Sawyer RC, Nayak R, Spiegelman S, Young CW, Woodcock T. An overview of thymidine. Cancer. 1980;45:1117–28.

    Article  CAS  PubMed  Google Scholar 

  14. Vertessy BG, Toth J. Keeping uracil out of DNA: physiological role, structure and catalytic mechanism of dUTPases. Acc Chem Res. 2009;42:97–106.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. 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: e19546.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  16. Wilson PM, Danenberg PV, Johnston PG, Lenz HJ, Ladner RD. Standing the test of time: targeting thymidylate biosynthesis in cancer therapy. Nat Rev Clin Oncol. 2014;11:282–98.

    Article  CAS  PubMed  Google Scholar 

  17. Bommera RK, Kethireddy S, Govindapur RR, Eppakayala L. Synthesis, biological evaluation and docking studies of 1,2,4-oxadiazole linked 5-fluorouracil derivatives as anticancer agents. BMC Chem. 2021;15:30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Yokogawa T, Yano W, Tsukioka S, Osada A, Wakasa T, Ueno H, Hoshino T, Yamamura K, Fujioka A, Fukuoka M, et al. dUTPase inhibition confers susceptibility to a thymidylate synthase inhibitor in DNA-repair-defective human cancer cells. Cancer Sci. 2021;112:422–32.

    Article  CAS  PubMed  Google Scholar 

  19. Koehler SE, Ladner RD. Small interfering RNA-mediated suppression of dUTPase sensitizes cancer cell lines to thymidylate synthase inhibition. Mol Pharmacol. 2004;66:620–6.

    CAS  PubMed  Google Scholar 

  20. Alam MS, Moriyama H, Matsumoto M. Enzyme kinetics of dUTPase from the planarian Dugesia ryukyuensis. BMC Res Notes. 2019;12:163.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Alam MS, Moriyama H, Matsumoto M. Inhibition of Dr-dut gene causes DNA damage in planarian. Mol Reprod Dev. 2018;85:188–96.

    Article  CAS  PubMed  Google Scholar 

  22. Palinkas HL, Racz GA, Gal Z, Hoffmann OI, Tihanyi G, Rona G, Gocza E, Hiripi L, Vertessy BG. CRISPR/Cas9-mediated knock-out of dUTPase in mice leads to early embryonic lethality. Biomolecules. 2019.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Muha V, Horvath A, Bekesi A, Pukancsik M, Hodoscsek B, Merenyi G, Rona G, Batki J, Kiss I, Jankovics F, et al. Uracil-containing DNA in Drosophila: stability, stage-specific accumulation, and developmental involvement. PLoS Genet. 2012;8: e1002738.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  24. Horvath A, Bekesi A, Muha V, Erdelyi M, Vertessy BG. Expanding the DNA alphabet in the fruit fly: uracil enrichment in genomic DNA. Fly. 2013;7:23–7.

    Article  CAS  PubMed  Google Scholar 

  25. Muha V, Zagyva I, Venkei Z, Szabad J, Vertessy BG. Nuclear localization signal-dependent and -independent movements of Drosophila melanogaster dUTPase isoforms during nuclear cleavage. Biochem Biophys Res Commun. 2009;381:271–5.

    Article  CAS  PubMed  Google Scholar 

  26. Knorr E, Fishilevich E, Tenbusch L, Frey MLF, Rangasamy M, Billion A, Worden SE, Gandra P, Arora K, Lo W, et al. Gene silencing in Tribolium castaneum as a tool for the targeted identification of candidate RNAi targets in crop pests. Sci Rep. 2018;8:2061.

    Article  PubMed  PubMed Central  Google Scholar 

  27. 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:W296–303.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Bordogna A, Pandini A, Bonati L. Predicting the accuracy of protein-ligand docking on homology models. J Comput Chem. 2011;32:81–98.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Forrest LR, Tang CL, Honig B. On the accuracy of homology modeling and sequence alignment methods applied to membrane proteins. Biophys J. 2006;91:508–17.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Benkert P, Biasini M, Schwede T. Toward the estimation of the absolute quality of individual protein structure models. Bioinformatics. 2011;27:343–50.

    Article  CAS  PubMed  Google Scholar 

  31. Wilkins MR, Gasteiger E, Bairoch A, Sanchez JC, Williams KL, Appel RD, Hochstrasser DF. Protein identification and analysis tools in the ExPASy server. Methods Mol Biol. 1999;112:531–52.

    CAS  PubMed  Google Scholar 

  32. Bajaj M, Moriyama H. Purification, crystallization and preliminary crystallographic analysis of deoxyuridine triphosphate nucleotidohydrolase from Arabidopsis thaliana. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2007;63:409–11.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Homma K, Moriyama H. Crystallization and crystal-packing studies of Chlorella virus deoxyuridine triphosphatase. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2009;65:1030–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Badalucco L, Poudel I, Yamanishi M, Natarajan C, Moriyama H. Crystallization of Chlorella deoxyuridine triphosphatase. Acta Crystallogr Sect F Struct Biol Cryst Commun. 2011;67:1599–602.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Kesinger E, Liu J, Jensen A, Chia CP, Demers A, Moriyama H. Influenza D virus M2 protein exhibits ion channel activity in Xenopus laevis oocytes. PLoS ONE. 2018;13: e0199227.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Nakachi M, Matsumoto M, Terry PM, Cerny RL, Moriyama H. Identification of guanylate cyclases and related signaling proteins in sperm tail from sea stars by mass spectrometry. Mar Biotechnol. 2008;10:564–71.

    Article  CAS  Google Scholar 

  37. Keller A, Nesvizhskii AI, Kolker E, Aebersold R. Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal Chem. 2002;74:5383–92.

    Article  CAS  PubMed  Google Scholar 

  38. Nesvizhskii AI, Keller A, Kolker E, Aebersold R. A statistical model for identifying proteins by tandem mass spectrometry. Anal Chem. 2003;75:4646–58.

    Article  CAS  PubMed  Google Scholar 

  39. Larsson G, Nyman PO, Kvassman JO. Kinetic characterization of dUTPase from Escherichia coli. J Biol Chem. 1996;271:24010–6.

    Article  CAS  PubMed  Google Scholar 

  40. Larsson G, Svensson LA, Nyman PO. Crystal structure of the Escherichia coli dUTPase in complex with a substrate analogue (dUDP). Nat Struct Biol. 1996;3:532–8.

    Article  CAS  PubMed  Google Scholar 

  41. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Kovari J, Barabas O, Takacs E, Bekesi A, Dubrovay Z, Pongracz V, Zagyva I, Imre T, Szabo P, Vertessy BG. Altered active site flexibility and a structural metal-binding site in eukaryotic dUTPase: kinetic characterization, folding, and crystallographic studies of the homotrimeric drosophila enzyme. J Biol Chem. 2004;279:17932–44.

    Article  CAS  PubMed  Google Scholar 

  43. 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:243–53.

    Article  CAS  PubMed  Google Scholar 

  44. Barabas O, Pongracz V, Kovari J, Wilmanns M, Vertessy BG. Structural insights into the catalytic mechanism of phosphate ester hydrolysis by dUTPase. J Biol Chem. 2004;279:42907–15.

    Article  CAS  PubMed  Google Scholar 

  45. 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:178–86.

    Article  CAS  PubMed  Google Scholar 

  46. Tongaonkar HB, Kulkarni JN, Kamat MR. Solitary metastases from renal cell carcinoma: a review. J Surg Oncol. 1992;49:45–8.

    Article  CAS  PubMed  Google Scholar 

  47. Bjornberg O, Nyman PO. The dUTPases from herpes simplex virus type 1 and mouse mammary tumour virus are less specific than the Escherichia coli enzyme. J Gen Virol. 1996;77(Pt 12):3107–11.

    Article  PubMed  Google Scholar 

  48. Slot Christiansen L, Munch-Petersen B, Knecht W. Non-Viral deoxyribonucleoside kinases-diversity and practical use. J Genet Genomics. 2015;42:235–48.

    Article  PubMed  Google Scholar 

  49. Piskur J, Sandrini MP, Knecht W, Munch-Petersen B. Animal deoxyribonucleoside kinases: ‘forward’ and ‘retrograde’ evolution of their substrate specificity. FEBS Lett. 2004;560:3–6.

    Article  CAS  PubMed  Google Scholar 

  50. Ladner RD, Lynch FJ, Groshen S, Xiong YP, Sherrod A, Caradonna SJ, Stoehlmacher J, Lenz HJ. dUTP nucleotidohydrolase isoform expression in normal and neoplastic tissues: association with survival and response to 5-fluorouracil in colorectal cancer. Cancer Res. 2000;60:3493–503.

    CAS  PubMed  Google Scholar 

  51. Racz GA, Nagy N, Varady G, Tovari J, Apati A, Vertessy BG. Discovery of two new isoforms of the human DUT gene. Sci Rep. 2023;13:7760.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Doi T, Yoh K, Shitara K, Takahashi H, Ueno M, Kobayashi S, Morimoto M, Okusaka T, Ueno H, Morizane C, et al. First-in-human phase 1 study of novel dUTPase inhibitor TAS-114 in combination with S-1 in Japanese patients with advanced solid tumors. Invest New Drugs. 2019;37:507–18.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank Drs. Etsuko Moriyama, Blair D. Siegfried, and Ana Maria Velez Arango for their collaboration and support. We thank Dr. Javier Seravalli and Professor Donald Becker in the Department of Biochemistry at the University of Nebraska—Lincoln (UNL), USA, for their technical support in the stopped-flow activity assays. We thank Drs. Mike Naldrett and Sophie Alvarez at Proteomics & Metabolomics Facility, Center for Biotechnology, UNL, USA, for their technical support in protein identification.


Research Council Interdisciplinary grant, Office of Research & Economic Development, UNL.

Author information

Authors and Affiliations



CRR conducted the protein production, enzyme kinetics, and structural mining and analyzed the data. HM led the study and designed the experiments, and both authors wrote and approved the final manuscript.

Corresponding author

Correspondence to Hideaki Moriyama.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Fig. S1.

Structural modeling. a The monomer model of WCR dUTPase (subunit A). The five conserved motifs (M1-M5), Arg89, N- and C-termini are indicated. b The dUTPase trimer model. Chains B and C are added to the monomer view (a). The stick model represents the key residues in the active site. The boxed area is enlarged in Fig 1c. Method: Structural models were built using SWISS-MODEL [27] with a human dUTPase (PDB ID: 3ehw) as the template. Models for dUTPase from WCR in the native and mutated version, Arg89Lys, had QMEAN values of −0.3 and −0.35, respectively. Models are usable because the root mean square deviation between the model and template was less than 2 Å [28, 29] and the QMEAN was less than 1 [30]. Structural mining was performed using PyMOL (Version 2.0, Schrödinger, New York, NY, USA) and the ProtParam server was used to calculate protein parameters [31].

Additional file 2: Fig. S2.

The WCR DUT construct sequence after being optimized for codon utilization of E. coli. The DUT construct sequence was cloned into pET-15 using NcoI and XhoI sites. One internal mutation (G → A at the 253th nucleotide) was placed to produce Arg89Lys. After thrombin cleavage (indicated by //), the N-terminal sequence, including the His tag, was removed.

Additional file 3: Fig. S3.

Peptides identified by mass spectrophotometry. a. Mass spectrum. b. Identified peptides. The full-length of the dUTPase protein construct sequence was identified by MS/MS. Sixty-seven exclusive unique spectra were identified out of the 968 total spectra. Likely deaminated glutamines are highlighted in cyan. All 33 unique peptide fragments are aligned against the dUTPase construct sequence. The numbers after each aligned peptide fragment represent the absolute abundance out of the 968 spectra. Method: The SDS-band was cut [35, 36] and submitted for mass spectrometry at the University of Nebraska–Lincoln Proteomics and Metabolomic Research Core Facility (Lincoln, NE, USA) [37, 38].

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Riera-Ruiz, C., Moriyama, H. Enzyme kinetics of deoxyuridine triphosphatase from Western corn rootworm. BMC Res Notes 16, 336 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: