Original Research
Frontiers in Drug Discovery ·
Volume 6 - 2026 | https://doi.org/10.3389/fddsv.2026.1749327
Optimization of a malachite green–pyrophosphatase assay for the Zika virus RNA-dependent RNA polymerase and preliminary identification of natural product modulators
- Vanessa Aitken1
- Anthony F. T. Moore2
- Manuel Terrazas-López1
- Salvador Rodarte3
- Eda Koculi1*corresponding author
- 1Department of Chemistry and Biochemistry, The University of Texas at El Paso, El Paso, TX, United States
- 2Department of Chemistry, University of Central Florida, Orlando, FL, United States
- 3Computational Sciences Program, The University of Texas at El Paso, El Paso, TX, United States
Abstract
The Zika virus (ZIKV) is an Orthoflavivirus linked to several neurological and developmental diseases, and it remains a significant global health concern for which no treatments or vaccines are currently available. In this study, we optimized the malachite green colorimetric assay, previously used to measure the RNA-dependent RNA polymerase (RdRp) activity of the hepatitis C and foot-and-mouth disease viruses, and applied this assay for the preliminary identification of the ZIKV RdRp modulators. RdRp is an enzyme essential for viral replication. Modulator screening was conducted using 204 compounds from the Natural Product Set IV of the National Cancer Institute Developmental Therapeutics Program. We identified two compounds, purpurogallin and digallic acid, as preliminary modulators of the ZIKV RdRp. Future studies will validate these preliminary hits using orthogonal assays, perform dose-response analyses, and elucidate their mechanisms of action.
Keywords
Introduction
The Zika virus (ZIKV) is a positive-sense, single-stranded RNA virus belonging to the Orthoflavivirus genus, which can, in certain cases, cause neurological and developmental disorders (Giraldo et al., 2023; de Sales-Neto et al., 2024; Postler et al., 2023; Godoy et al., 2017). While the primary route of ZIKV transmission is through mosquito bites, it can also be transmitted from mother to fetus and via bodily fluids (Giraldo et al., 2023; de Sales-Neto et al., 2024; Godoy et al., 2017; Rabe et al., 2025). Recently, more than 92 countries have experienced a resurgence of ZIKV outbreaks (Rabe et al., 2025). Moreover, rising global temperatures are expected to increase the risk of ZIKV infections, with an estimated 1.3 billion people projected to be at risk by 2050 (Ryan et al., 2019; Ryan et al., 2021). To date, treatment and vaccines for ZIKV are not available (de Sales-Neto et al., 2024). Therefore, the development of new therapeutic agents to treat ZIKV infections is of significant public health interest (de Sales-Neto et al., 2024). In this study, we optimized a previously used malachite green high-throughput screening (HTS) assay and used it to screen for small-molecule modulators of the ZIKV RNA-dependent RNA polymerase (RdRp), using the Natural Product Set IV from the National Cancer Institute’s Developmental Therapeutics Program (Vardakou et al., 2014; Nguyen et al., 2013; Van Veldhoven and Mannaerts, 1987; Tillotson et al., 2017; Moore et al., 2021; Moore et al., 2025; Sivinski et al., 2022). To the best of our knowledge, this assay has not been used previously to search for ZIKV inhibitors.
The ZIKV genome single-stranded RNA functions as an mRNA (Giraldo et al., 2023; de Sales-Neto et al., 2024). From this mRNA all the proteins required for virus structure, replication and propagation are translated (Giraldo et al., 2023; de Sales-Neto et al., 2024). This RNA is translated by the host cellular machinery into a single-polypeptide chain consisting of ten viral proteins: C, preM, E, NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5 (Giraldo et al., 2023; de Sales-Neto et al., 2024; Sager et al., 2018). The polypeptide is then cleaved by viral and cellular proteases to produce three structural proteins, C, preM, and E, which form the viral shell. The seven nonstructural proteins, NS1, NS2A, NS2B, NS3, NS4A, NS4B, and NS5 are involved in viral replication and immune evasion (Sager et al., 2018; Sirohi and Kuhn, 2017). NS5 contains two domains connected by a linker: an N-terminal methyltransferase domain and a C-terminal RdRp domain (Flory et al., 2021). RdRp is responsible for replicating the ZIKV RNA genome and is essential for ZIKV replication (Flory et al., 2021). RdRp represents an attractive therapeutic target for treating ZIKV infections with potentially minimal side effects because it is encoded exclusively by the ZIKV genome, humans do not possess any homologous protein with RdRp activity, and it plays a critical role in the ZIKV life cycle (Godoy et al., 2017).
Previously, the malachite green HTS assay was employed to measure RdRp activity from hepatitis C and foot-and-mouth disease viruses (Nguyen et al., 2013). The malachite green is a colorimetric assay that detects the complex formed from malachite green, molybdate, and inorganic phosphate (Vardakou et al., 2014; Van Veldhoven and Mannaerts, 1987; Tillotson et al., 2017; Moore et al., 2021; Moore et al., 2025; Sivinski et al., 2022). However, the RdRp from hepatitis C and foot-and-mouth disease viruses, and polymerases in general, produce pyrophosphate during their reactions, instead of inorganic phosphate (Nguyen et al., 2013). Thus, in the previous study, the pyrophosphatase enzyme, which converts pyrophosphate to inorganic phosphate, was included in the RdRp reaction buffers (Nguyen et al., 2013).
Traditional methods for monitoring polymerase activity include the direct detection of radiolabeled nucleic acid products using gel electrophoresis or thin-layer chromatography (Schrier and Wilson, 1973; Shatkin and Sipe, 1968; Sierra et al., 2007; Rodriguez-Wells et al., 2001; Jin et al., 2012). These approaches are highly sensitive. On the other hand, they require strict regulatory handling and careful disposal due to radiation hazards. More importantly, the radiolabeled methods are not compatible with HTS. Thus, while valuable for hit validation, they are impractical for initial large-scale compound screens.
Non-radioactive alternatives, which involve chromophore-labeled nucleotides incorporated into the polymerized RNA molecule, avoid the safety concerns. However, a shortcoming of these assays is less efficient incorporation of modified nucleotides by RdRp than their natural counterparts (Vassiliou et al., 2000; Park et al., 2002). Thus, an inhibitor may appear to work because the polymerase is inefficient in incorporating the chromophore-labeled nucleotides (Vassiliou et al., 2000; Park et al., 2002).
Luminescence-based assays provide another option: pyrophosphate released during nucleotide incorporation is converted to ATP by ATP sulfurylase in the presence of adenosine 5′-phosphosulfate, and the resulting ATP fuels a luciferase reaction that generates detectable photons (Lahser and Malcolm, 2004). However, practical limitations constrain this assay’s utility in HTS. The luminescence signal typically decays rapidly, requiring precisely timed measurements, and background ATP present in reaction components can contribute to the signal, complicating data interpretation (Lahser and Malcolm, 2004).
In contrast, the coupled malachite green/pyrophosphatase colorimetric assay is simple, inexpensive, non-radioactive, and highly adaptable for HTS (Nguyen et al., 2013). By quantifying inorganic phosphate generated during nucleotide incorporation, it provides a robust and reliable readout of polymerase activity (Nguyen et al., 2013). Its proven reliability across diverse enzyme reactions demonstrates its robustness and reproducibility, making it particularly well suited for HTS of compound libraries (Tillotson et al., 2017; Moore et al., 2021; Moore et al., 2025; Park et al., 2002; Wu, 2011).
In order to examine a large library of compounds and accurately determine the modulation of RdRp activity, it is important to design an HTS assay that is both reproducible from day-to-day and offers an expansive dynamic range. While the malachite green HTS assay previously employed for hepatitis C and foot-and-mouth disease virus RdRp was elegantly designed and executed, key performance metrics, such as dynamic range (Z′) and signal window (SW), were not systematically evaluated (Nguyen et al., 2013; Iversen et al., 2006). In this study, we adapted the malachite green HTS assay and performed statistical analysis using Z′ and SW parameters (Iversen et al., 2006). Next, we used the optimized assay to screen for natural product inhibitors of the Zika virus RdRp, and identified two compounds, purpurogallin and digallic acid, that modulate the ZIKV RdRp activity. In its current form, the malachite green assay described here can be used to identify small-molecule modulators of DNA and RNA polymerases in a high-throughput manner.
Materials and methods
Plasmid construction and protein expression
RdRp coding sequence from the ZIKV Puerto Rico strain PRVABC59 was inserted into the pET His6 Sumo TEV LIC (Novagen) by the Protein Production and Analysis Facility at the Sanford Burnham Prebys Medical Discovery Institute in Florida. The final protein construct consisted of a His-tag followed by a SUMO sequence, a TEV protease cut site, eight of the ten amino acid residues of the linker region separating the methyltransferase from RdRp, and the entire sequence of RdRp. The SUMO sequence has been shown to enhance the translation and solubility of proteins expressed in bacteria (Butt et al., 2005). The linker region connecting the methyltransferase and RdRp was retained because this region has been shown to regulate RdRp activity (Flory et al., 2021). The His-tag and SUMO sequence were not cleaved. Thus, the protein assay described here was performed with the RdRp construct containing both the His-tag and SUMO sequence. The inorganic pyrophosphatase from Saccharomyces cerevisiae was purchased from Fisher Scientific.
RdRp expression and purification
Escherichia coli (E. coli) BL21 CodonPlus (DE3) cells were transformed with the pET-His-SUMO-TEV-LIC vector containing the RdRp coding sequence. The cells were grown in Luria Broth medium supplemented with 30 μg/mL kanamycin and 34 μg/mL chloramphenicol at 37 °C until reaching an optical density at 600 nm (OD600) of approximately 0.6. Protein expression was induced with 500 μM isopropyl β-D-1-thiogalactopyranoside, and the cultures were incubated for an additional 16 h at 20 °C. Cells were harvested by centrifugation.
The cell pellets were resuspended in a lysis buffer consisting of a phosphate-buffered saline (PBS: 137 mM NaCl, 2.7 mM KCl, 10 mM Na2HPO4, 1.8 mM KH2PO4), supplemented with 5 mM β-mercaptoethanol (BME) and 1 mM phenylmethylsulfonyl fluoride (PMSF). Cells were lysed by sonication on ice using a Branson Digital Sonifier 450 equipped with a 0.5 mm diameter tip using six 30-s pulses in ice. The lysate was cleared by centrifugation at 11,648 x g for 60 min at 4 °C.
The clarified supernatant was loaded onto an 8 mL Ni-NTA-affinity column pre-equilibrated with PBS containing 5 mM BME, 1 mM PMSF and 20 mM imidazole. The protein was eluted using PBS supplemented with 5 mM BME, 1 mM PMSF and 500 mM imidazole. Fractions containing the RdRp protein were pooled and diluted 1:10 with a buffer consisting of 25 mM Tris-HCl (pH 8.0) and 5 mM BME.
Subsequently, the protein sample was loaded onto a 5 mL heparin affinity column (HiTrap Heparin HP, GE Healthcare Life Sciences). Heparin columns selectively bind nucleic acid–binding proteins and compete with nucleic acids for protein interaction sites (Gao et al., 2021). Therefore, this step facilitated both the further purification of the RdRp protein and the removal of any bound nucleic acids. The column was equilibrated and washed with a buffer consisting of 25 mM Tris-HCl (pH 8.0), 30 mM NaCl, and 5 mM BME. The RdRp was eluted using a linear gradient of the heparin elution buffer consisting of 25 mM Tris-HCl (pH 8.0), 1000 mM NaCl, and 5 mM BME.
To further remove undesired proteins and RdRp aggregates, fractions containing the RdRp were pooled, concentrated, and loaded onto a size-exclusion chromatography column (Superdex 200 10/300 GL, GE Healthcare). The column was equilibrated with the gel-filtration buffer consisting of 25 mM Tris-HCl (pH 8.0), 100 mM KCl, and 2 mM dithiothreitol (DTT). Fractions containing the RdRp were combined and concentrated using 10 kDa molecular weight cutoff centrifugal filters (Amicon Ultra Ultracel-10K, Millipore). Finally, purified RdRp was flash-frozen in small aliquots and stored at −80 °C until use.
Using the malachite green assay for screening inhibitors of the ZIKV RdRp
The RNA polymerase reaction conditions were adapted from a previous study (Vardakou et al., 2014). Pyrophosphate release during the RdRp reaction was measured using the malachite green/molybdate colorimetric assay (Vardakou et al., 2014). Because the malachite green/molybdate colorimetric assay specifically detects inorganic phosphate, but not pyrophosphate, an inorganic pyrophosphatase was added to the reaction mixture to convert a pyrophosphate into an inorganic phosphate (Vardakou et al., 2014; Nguyen et al., 2013; Van Veldhoven and Mannaerts, 1987).
The RNA polymerase reaction initiation complex was formed by annealing poly(C) RNA molecules (200–250 nucleotides in length) with a 10-nucleotide-long poly(G) RNA primer (Lu et al., 2017). The poly(C) RNA was purchased from Sigma-Aldrich and further purified by polyacrylamide gel electrophoresis (Moore et al., 2025; Moore et al., 2017; Woodson and Koculi, 2009; Childs et al., 2016). The HPLC-purified poly(G) oligonucleotides were obtained from Integrated DNA Technologies (IDT). Annealing was performed by mixing the two RNA strands in an annealing buffer containing 10 mM Tris-HCl (pH 7.5), 100 mM NaCl, and 0.2 mM EDTA-KOH (pH 8.0). The mixture was incubated at 95 °C for 2 min and then cooled slowly to room temperature, in order to promote proper duplex formation (Moore et al., 2025).
The RdRp reaction was performed in a 384-well plate at 22 °C, in a 12.5 µL reaction mixture consisting of 5 mM Tris-HCl pH 7.5, 7 mM NaCl, 8 mM MnCl2, 10% glycerol, 0.01% Triton X-100, 10 mM DTT, 0.05 mU/µL inorganic pyrophosphatase, 0.5 mU/µL SUPERase·In, 500 nM annealed RNA, 0.5 mM GTP·Mn, 0–200 nM RdRp, and 1% DMSO (Lu et al., 2017). Reactions for compound screening also contained 10 µM of compound. Because the objective of this study was to optimize, validate, and assess the feasibility of the malachite green–pyrophosphatase assay, only a single polymerase reaction was performed for each small molecule. Replicate confirmation was beyond the scope of this study. The RdRp reactions were initiated by the addition of GTP Mn and allowed to continue for 30–90 min. These reactions were quenched with 10 mM EDTA, which sequesters Mn2+ and halts the RdRp reaction (Moore et al., 2025; Lu et al., 2017).
To detect the orthophosphate generated by pyrophosphatase activity, 75 µL of the malachite green–molybdate reagent was added to the RdRp reaction mixture. After 1 min, 10 µL of 34% sodium citrate was added to stop the reaction and stabilize the color complex. The mixture was then incubated at 22 °C for 50 min before measuring absorbance at 625 nm using a Molecular Devices SpectraMax Plus 384 spectrophotometer and microplate reader. The malachite green–molybdate reagent was prepared as previously described (Moore et al., 2021; Moore et al., 2025). Compounds were classified as modulators if their normalized activity fell below. In Equation 1, X ¯ p is the average absorbance of the positive control, while σp is the standard deviation of the X ¯ p. Equation 1 ensures that selected hits differed from the positive control by at least three standard deviations, corresponding to a 99.7% confidence level under a normal distribution (Pukelsheim, 1994).
X ¯ p − 3 σ p (1)
For the two compounds classified as modulators, the RdRp reaction absorbance value in their presence was normalized to the RdRp positive and negative controls reactions averages. For this normalization, the X ¯ p was defined as 100% enzymatic activity, and the average absorbance of the negative control (X ¯ n) reactions defined as 0% activity. Normalized activity (%) for the two compounds was calculated using the equation below:where C c o m p o u n d is the absorbance value of the RdRp reaction in the presence of the compound.
N o r m a l i z e d A c t i v i t y % = C c o m p o u n d − X ¯ n X ¯ p − X ¯ n × 100 (2)
Results and discussion
Optimization of the malachite green assay for ZIKV RdRp inhibitor screening
The malachite green assay quantifies orthophosphate by detecting the green-colored complex formed from orthophosphate, malachite green, and molybdate (Van Veldhoven and Mannaerts, 1987; Tillotson et al., 2017; Moore et al., 2021; Moore et al., 2025; Sivinski et al., 2022). However, during double-stranded RNA synthesis, the RdRp releases a pyrophosphate as a byproduct (Nguyen et al., 2013). To enable the detection of the pyrophosphate using the assay, an inorganic pyrophosphatase, which converts a pyrophosphate to an orthophosphate, was included in the reaction mixture (Vardakou et al., 2014; Nguyen et al., 2013). We aimed to determine the RdRp reaction and assay conditions that yield reliable, proportional, and accurate measurements. Therefore, we evaluated the linearity of the RdRp reaction, as measured by the malachite green assay, with respect to both time and RdRp concentration. Our results demonstrate that the assay response remains linear for at least 90 min (Figure 1). Based on these findings, we selected a 75-min reaction time with 200 nM RdRp to maximize signal intensity while maintaining assay linearity.

Robustness, reproducibility, and dynamic range of the malachite green assay employed in this study
HTS robustness testing ensures that minor variations in experimental conditions across different days do not compromise the accuracy of hit detection. The hit detection can be compromised by misidentifying false positives as true hits, or by failing to detect true positives. These evaluations are essential for maintaining consistency and reliability in assay performance and ensuring that inter-day variability remains minimal. To assess the robustness of the malachite green HTS assay used for evaluating ZIKV RdRp activity, absorbance signals were measured for both positive and negative controls of the RdRp reactions. Negative controls lacked the RdRp enzyme, while positive controls contained both the RNA substrate and the enzyme. To evaluate reproducibility, the assays were conducted on four separate days. Our results showed minimal variation in absorbance values across both positive and negative RdRp control reactions. Most of the measurements fell within two standard deviations of X ¯ p and X ¯ n, and all of the values fell within three standard deviations of the mean absorbance value (Figure 2).

The Z′ factor is a dimensionless statistical parameter used to assess the signal dynamic range and variability between the positive and negative controls in an HTS assay (Equation 3) (Iversen et al., 2006). A high Z′ value indicates strong separation between positive and negative control signals and low data noise. Both strong separation between positive and negative controls and low data noise are essential for distinguishing true hits from experimental variability. An acceptable Z′ value ranges from 0.5 to 1 (Iversen et al., 2006). All Z′ values determined from experiments conducted over four separate days fell within the above range, indicating that malachite green assay as optimized by us is robust, exhibits a high signal-to-noise ratio, and is suitable for HTS (Table 1).
Z ′ = 1 − 3 σ p + 3 σ n X ¯ p − X ¯ n (3)
TABLE 1. The malachite green HTS assay is robust, has a large dynamic range, and a high signal-to-noise ratio.
| Plate/Day a | Z’ b | SW c |
|---|---|---|
| 1 | 0.731 | 9.115 |
| 2 | 0.655 | 7.796 |
| 3 | 0.666 | 8.790 |
| 4 | 0.807 | 13.911 |
The malachite-green assay, with the negative and positive controls, was performed on different plates on four different days. The positive control reactions contained both the RNA, substrate and the RdRp, while the negative control reactions lacked the RdRp enzyme.
Z’ values were calculated from experiments performed on four different days using Equation 3. The Z′ values indicate that the separation between the negative and positive controls is sufficient to reliably detect inhibitors. The Z′ values obtained across all 4 days of experimentation fall within the acceptable range of 0.5 < Z′ < 1, demonstrating that the assay exhibits high quality and is suitable for HTS (Iversen et al., 2006).
SW was calculated using Equation 4 for four different plates across multiple days. The SW values are significantly larger than 2 (the accepted SW, value for an HTS-suitable assay) (Iversen et al., 2006).
In Equation 3, X ¯ p, X ¯ n, and σp, are the parameters defined in Equations 1, 2, while σn is the standard deviation of X ¯ n (Iversen et al., 2006).
SW is a dimensionless parameter used to assess the difference in signal between negative and positive controls (Iversen et al., 2006). SW was calculated using Equation 4 across four different days and plates (Table 1) (Iversen et al., 2006). An SW value greater than two is considered adequate for detecting inhibitors in an HTS assay. The SW values obtained on different days were significantly greater than 2, demonstrating once again that our malachite green HTS assay is both robust and sensitive (Iversen et al., 2006).
S W = X ¯ p − X ¯ n − 3 σ p + σ n σ p (4)
Purpurogallin and digallic acid inhibit RdRp activity
The optimized green/molybdate colorimetric assay was used to screen for RdRp inhibitors. A total of 204 out of 419 natural products from the NCI Natural Products Set IV were screened. The screening was performed in the presence of 10 μM of each compound (Supplementary Table S1). The present study is focused on assay optimization, validation and demonstration of feasibility (Figures 1, 2; Table 1). Therefore, only a representative feasibility screening of 204 compounds was performed. An expanded library will be assessed in future work.
We found that purpurogallin (National Service Center (NSC) 35,676) and digallic acid (NSC 59263) modulated the polymerase activity of the ZIKV RdRp based on our threshold of 3σ cutoff (Figure 3; Equation 1; Table 2). While any compound can potentially interfere with a colorimetric assay, such as the malachite green, phytochemicals such as digallic acid and purpurogallin, due to their structural complexity and functional groups, could be more prone to affecting colorimetric assays (Figure 3B). Therefore, these digallic acid and purpurogallin are presented as preliminary RdRp modulator candidates, and the potential limitations of our screening have been discussed in detail further in the manuscript.

TABLE 2. Summary of compounds exhibiting significant inhibition of ZIKV RdRp in the malachite green assay.
| Compound a | Raw signal b | X ¯ p a | X ¯ n c | NormalizedActivity %d | σp a | 3σp cutoff a |
|---|---|---|---|---|---|---|
| Digallic acid | 0.861 | 1.213 | 0.406 | 56% | 0.054 | 1.015 |
| Purpurogallin | 0.777 | 1.213 | 0.406 | 46% | 0.054 | 1.015 |
Modulators were defined as compounds with raw signals below X ¯ *p* − 3 σ p, as described in Equation 1. X ¯ *p* is the average signal of the positive controls, and σ p is the standard deviation of X ¯ *p*.
Raw signal corresponds to the RdRp reaction malachite green absorbance measured in the presence of the compound.
X ¯ n represents the average absorbance signal of the negative controls.
Normalized activity was calculated as described in Equation 2.
Purpurogallin is a benzotropolone derivative (Bentley, 2008). This compound occurs naturally in the nut galls of the Quercus species or can be synthesized via pyrogallol oxidation (Bentley, 2008). Purpurogallin has been shown to protect cells from hydrogen peroxide-induced DNA damage and apoptosis, both of which are associated with oxidative stress (Bentley, 2008). Importantly, ZIKV infection induces oxidative stress, which facilitates viral replication (Kayesh et al., 2025). Therefore, purpurogallin may serve as a dual-function inhibitor by targeting ZIKV RdRp and by mitigating oxidative stress associated with infection.
Moreover, purpurogallin has demonstrated anti-inflammatory effects in human endothelial cells (Kim et al., 2012; Park et al., 2013; Ku et al., 2013). Acute inflammation is considered a central mechanism underlying ZIKV-associated pathology. Hence, modulation of the cellular inflammatory response has been proposed as a promising therapeutic strategy for ZIKV infection (de Sales-Neto et al., 2024). Importantly, the anti-inflammatory compound MK-591 has emerged as a potential candidate for the treatment of ZIKV infection (Paul et al., 2024). Accordingly, purpurogallin may function as a multifunctional inhibitor by attenuating the inflammatory response, inhibiting RdRp activity, and reducing oxidative stress.
The other compound shown to modulate the activity of ZIKV RdRp is digallic acid (Figure 3; Table 2). Digallic acid is derived from Pistacia lentiscus fruits or synthesized in laboratory via partial hydrolysis of tannic acid (Bhouri et al., 2010; Verzele et al., 1983). A previous study also reported that digallic acid inhibits the reverse transcriptase activity of the murine leukemia virus (MLV) and the human immunodeficiency virus (HIV) (Nakane et al., 1990). In that study, digallic acid inhibited the activity of reverse transcriptase by competing for the template RNA binding site (Nakane et al., 1990). In addition to inhibiting MLV and HIV reverse transcriptase, digallic acid was shown to moderately inhibit the activity of human DNA polymerases alpha and beta, but did not inhibit human DNA polymerase gamma or E. coli DNA polymerase I (Nakane et al., 1990). These cross-polymerase inhibitions, combined with mechanistic studies, suggest that digallic acid may interact with conserved structural motifs involved in DNA or RNA template binding. Future structural studies of the ZIKV RdRp and HIV reverse transcriptase in the presence of digallic acid could aid in the design of inhibitors structurally similar to digallic acid. These inhibitors could be highly specific for ZIKV and HIV polymerases and spare the human DNA polymerases. These dual-acting inhibitors may hold therapeutic potential for treating co-infections caused by ZIKV and HIV.
To further evaluate the mechanism of inhibition, we considered whether antioxidant properties might contribute to the observed decrease on the RdRp activity in the presence of purpurogallin and digallic acid. Thus, we compared the inhibitory activity of digallic acid and purpurogallin to that of other natural products tested that possessed antioxidant activity (Supplementary Table S2). Of the 204 compounds investigated, 14 were identified as having antioxidant properties (Supplementary Table S1). With the exception of digallic acid and purpurogallin, none of these antioxidants modulated ZIKV RdRp activity (Supplementary Table S1; Figure 3; Table 2). Therefore, the observed modulation by digallic acid and purpurogallin is unlikely to result from their general antioxidant properties and instead indicates a more specific interaction with the RdRp enzyme.
Despite their promising inhibitory activity, both digallic acid and purpurogallin present limitations in terms of stability and selectivity. Purpurogallin is generally stable under physiological conditions but may undergo degradation when exposed to elevated oxidative or radiolytic stress (Prasad et al., 1994; Ratti et al., 2023). In contrast, digallic acid is not very stable in cellular environments and is susceptible to hydrolysis into gallic acid monomers (Zheng et al., 2024; Rodriguez-Duran et al., 2011; Erdogan and Karadag, 2006). Furthermore, because of their potential for non-specific interactions and promiscuous binding, both compounds present limitations in their specificity for ZIKV RdRp (Bhouri et al., 2010; Ku and Bae, 2014; Veser, 1987; Souto et al., 2020; Zhen et al., 2019; Richard, 2019; Glanzer et al., 2016; Alfonso et al., 2022; Chatterjee et al., 2015; Namkung et al., 2011; Suzuki et al., 2014; Kumar et al., 2024). Nevertheless, these molecules may serve as preliminary chemical scaffolds for further inhibitor developments. Future computational modeling and structural studies of RdRp in complex with purpurogallin or digallic acid may inform rational scaffold optimization, which could improve both stability of these compounds and their selectivity for ZIKV RdRp. The findings presented here provide an early foundation for the development of two new classes of inhibitors.
Conclusions
The optimized malachite green–pyrophosphatase assay we employed here is robust, reproducible, and sensitive across different experimental days (Table 1). This assay, as optimized by our team, could be employed for the HTS of DNA and RNA polymerases, as well as other enzymes that produce pyrophosphate during their catalytic cycles. Through the screening of NCI natural product libraries, purpurogallin and digallic acid were found as preliminary inhibitors of the activity of the RdRp (Figure 3; Supplementary Table S1).
Limitations
Because the primary objectives of this work were to optimize, validate, and assess the feasibility of the malachite green–pyrophosphatase assay for HTS rather than to identify or characterize lead inhibitors for clinical development, orthogonal assays were not included, and all compounds were tested only once at a single screening concentration (10 μM). Thus, assay interference cannot be excluded.
In addition, because the objectives of this study were the optimization, validation, and feasibility assessment of the malachite green–pyrophosphatase assay for screening ZIKV RdRp inhibitors, the preliminary inhibitory activity of purpurogallin and digallic acid was evaluated only under in vitro conditions. Lastly, the determinations of the purpurogallin and digallic acid concentrations required to inhibit 50% of ZIKV RdRp activity (IC50), the measurements of these chemicals cellular antiviral efficacy, and cytotoxicity, remain to be performed.
Future studies will validate the preliminary hits, purpurogallin and digallic acid (Figure 3; Table 2), using assays that directly monitor RdRp-mediated RNA synthesis. In addition, these studies will characterize the IC50 values, evaluate antiviral activity in ZIKV-infected cell cultures, and determine the cytotoxicity of purpurogallin and digallic acid in various human cell lines.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fddsv.2026.1749327/full#supplementary-material
Supplementary material Supplementary material at the source
Footnotes
Edited by: Octavio Luiz Franco, Catholic University of Brasilia (UCB), Brazil
Reviewed by: Maria Letizia Barreca, University of Perugia, Italy
Hrithika Pandey, Sharda University, India
Statements
- Conflict of interest
The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
- Funding
The author(s) declared that financial support was received for this work and/or its publication. This work was supported in part by the National Institute of General Medical Sciences grant R01-GM131062, the University of Texas System Rising STARs Program, and the start-up from the Chemistry and Biochemistry Department at the University of Texas at El Paso (to EK).
- Data availability
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author/s.
- Author contributions
VA: Data curation, Formal Analysis, Writing – original draft, Writing – review and editing. AM: Data curation, Formal Analysis, Writing – original draft, Writing – review and editing, Conceptualization, Investigation, Methodology. MT: Data curation, Formal Analysis, Writing – original draft, Writing – review and editing. SR: Data curation, Formal Analysis, Writing – original draft, Writing – review and editing. EK: Data curation, Formal Analysis, Writing – original draft, Writing – review and editing, Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Resources, Supervision, Validation, Visualization.
- Acknowledgements
We are grateful to Luis Gracia Mazuca for the discussion regarding this manuscript.
References
- Alfonso E. E., Troche R., Deng Z., Annamalai T., Chapagain P., Tse-Dinh Y. C. (2022). Potent inhibition of bacterial DNA gyrase by digallic acid and other gallate derivatives. ChemMedChem 17 (23), e202200301. doi:10.1002/cmdc.202200301
- Bentley R. (2008). A fresh look at natural tropolonoids. Nat. Prod. Rep. 25 (1), 118 – 138. doi:10.1039/b711474e
- Bhouri W., Derbel S., Skandrani I., Boubaker J., Bouhlel I., Sghaier M. B. (2010). Study of genotoxic, antigenotoxic and antioxidant activities of the digallic acid isolated from Pistacia lentiscus fruits. Toxicol Vitro 24 (2), 509 – 515. doi:10.1016/j.tiv.2009.06.024
- Butt T. R., Edavettal S. C., Hall J. P., Mattern M. R. (2005). SUMO fusion technology for difficult-to-express proteins. Protein Expr. Purif. 43 (1), 1 – 9. doi:10.1016/j.pep.2005.03.016
- Chatterjee T., Sheikh I. A., Chakravarty D., Chakrabarti P., Sarkar P., Saha T. (2015). Effects of small molecule calcium-activated chloride channel inhibitors on structure and function of accessory cholera enterotoxin (ace) of Vibrio cholerae. PLoS One 10 (11), e0141283. doi:10.1371/journal.pone.0141283
- Childs J. J., Gentry R. C., Moore A. F., Koculi E. (2016). The DbpA catalytic core unwinds double-helix substrates by directly loading on them. RNA 22 (3), 408 – 415. doi:10.1261/rna.052928.115
- de Sales-Neto J. M., Madruga Carvalho D. C., Arruda Magalhaes D. W., Araujo Medeiros A. B., Soares M. M., Rodrigues-Mascarenhas S. (2024). Zika virus: antiviral immune response, inflammation, and cardiotonic steroids as antiviral agents. Int. Immunopharmacol. 127, 111368. doi:10.1016/j.intimp.2023.111368
- Erdogan G., Karadag R. (2006). Potentiometric and spectrophotometric studies of the mixed-ligand complexes of quercetin and digallic acid with Iron(II) and Magnesium(II). Rev. Anal. Chem. 25 (4), 257 – 270. doi:10.1515/revac.2006.25.4.257
- Flory C., Chan K. W. K., Bifani A. M., Choy M. M. J., Lescar J., Ooi E. E. (2021). Optimal flexibility of the linker region of zika virus NS5 methyltransferase-polymerase is critical for virus replication. Antivir. Res. 195, 105194. doi:10.1016/j.antiviral.2021.105194
- Gao Y., Cao D., Pawnikar S., Akhter S., Miao Y., Liang B. (2021). Efficient purification and assembly of ribonucleoprotein complex for interaction analysis by MST assay coupled with GaMD simulations. Star. Protoc. 2 (1), 100315. doi:10.1016/j.xpro.2021.100315
- Giraldo M. I., Gonzalez-Orozco M., Rajsbaum R. (2023). Pathogenesis of zika virus infection. Annu. Rev. Pathol. 18, 181 – 203. doi:10.1146/annurev-pathmechdis-031521-034739
- Glanzer J. G., Endres J. L., Byrne B. M., Liu S., Bayles K. W., Oakley G. G. (2016). Identification of inhibitors for single-stranded DNA-Binding proteins in eubacteria. J. Antimicrob. Chemother. 71 (12), 3432 – 3440. doi:10.1093/jac/dkw340
- Godoy A. S., Lima G. M., Oliveira K. I., Torres N. U., Maluf F. V., Guido R. V. (2017). Crystal structure of zika virus NS5 RNA-Dependent RNA polymerase. Nat. Commun. 8, 14764. doi:10.1038/ncomms14764
- Iversen P. W., Eastwood B. J., Sittampalam G. S., Cox K. L. (2006). A comparison of assay performance measures in screening assays: signal window, Z' factor, and assay variability ratio. J. Biomol. Screen 11 (3), 247 – 252. doi:10.1177/1087057105285610
- Jin Z., Leveque V., Ma H., Johnson K. A., Klumpp K. (2012). Assembly, purification, and pre-steady-state kinetic analysis of active RNA-Dependent RNA polymerase elongation complex. J. Biol. Chem. 287 (13), 10674 – 10683. doi:10.1074/jbc.M111.325530
- Kayesh M. E. H., Kohara M., Tsukiyama-Kohara K. (2025). Effects of oxidative stress on viral infections: an overview. Npj Viruses 3 (1), 27. doi:10.1038/s44298-025-00110-3
- Kim T. H., Ku S. K., Lee I. C., Bae J. S. (2012). Anti-inflammatory functions of purpurogallin in LPS-Activated human endothelial cells. BMB Rep. 45 (3), 200 – 205. doi:10.5483/BMBRep.2012.45.3.200
- Ku S. K., Bae J. S. (2014). Antiplatelet and antithrombotic activities of purpurogallin in vitro and in vivo. BMB Rep. 47 (7), 376 – 381. doi:10.5483/bmbrep.2014.47.7.195
- Ku S. K., Lee I.C., Bae J. S. (2013). Inhibitory effects of purpurogallin on the endothelial protein C receptor shedding in vitro and in vivo. J. Korean Soc. Appl. Biol. Chem. 56, 519 – 524. doi:10.1007/s13765-013-3169-7
- Kumar N., Goel N., Gandara J. S. (2024). Advancement of phenolic acids in drug discovery: fundamentals and applications.
- Lahser F. C., Malcolm B. A. (2004). A continuous nonradioactive assay for RNA-Dependent RNA polymerase activity. Anal. Biochem. 325 (2), 247 – 254. doi:10.1016/j.ab.2003.10.034
- Lu G., Bluemling G. R., Collop P., Hager M., Kuiper D., Gurale B. P. (2017). Analysis of ribonucleotide 5'-Triphosphate analogs as potential inhibitors of zika virus RNA-dependent RNA polymerase by using nonradioactive polymerase assays. Antimicrob. Agents Chemother. 61 (3). doi:10.1128/AAC.01967-16
- Moore A. F., Gentry R. C., Koculi E. (2017). DbpA is a region-specific RNA helicase. Biopolymers 107 (3), e23001. doi:10.1002/bip.23001
- Moore A. F. T., de Victoria A. L., Koculi E. (2021). Interactions of the C-Terminal truncated DEAD-box protein DDX3X with RNA and nucleotide substrates. ACS Omega 6 (19), 12640 – 12646. doi:10.1021/acsomega.1c00700
- Moore A. F. T., Berhie Y., Weislow S. I., Koculi E. (2025). Substrate specificities of DDX1: a human DEAD-box protein. ACS Omega. 10 (3), 2598 – 2607. doi:10.1021/acsomega.4c07522
- Nakane H., Fukushima M., Ono K. (1990). Differential inhibition of reverse transcriptase and various DNA polymerases by digallic acid and its derivatives. J. Nat. Prod. 53 (5), 1234 – 1240. doi:10.1021/np50071a015
- Namkung W., Phuan P. W., Verkman A. S. (2011). TMEM16A inhibitors reveal TMEM16A as a minor component of calcium-activated chloride channel conductance in airway and intestinal epithelial cells. J. Biol. Chem. 286 (3), 2365 – 2374. doi:10.1074/jbc.M110.175109
- Nguyen H. T., Chong Y., Oh D. K., Heo Y. S., Viet P. T., Kang L. W. (2013). An efficient colorimetric assay for RNA synthesis by viral RNA-Dependent RNA polymerases, using thermostable pyrophosphatase. Anal. Biochem. 434 (2), 284 – 286. doi:10.1016/j.ab.2012.11.019
- Park C., Kee Y., Park J., Myung H. (2002). A nonisotopic assay method for hepatitis C virus NS5B polymerase. J. Virol. Methods 101 (1-2), 211 – 214. doi:10.1016/s0166-0934(01)00443-8
- Park H. Y., Kim T. H., Kim C. G., Kim G. Y., Kim C. M., Kim N. D. (2013). Purpurogallin exerts anti-inflammatory effects in lipopolysaccharide-stimulated BV2 microglial cells through the inactivation of the NF-kappaB and MAPK signaling pathways. Int. J. Mol. Med. 32 (5), 1171 – 1178. doi:10.3892/ijmm.2013.1478
- Paul S., Patra S., Mondal A., Adhikari G.S., Ghosh P., Bysack M. (2024). The role of microglia in zika virus pathogenesis: possible diagnostic and therapeutic strategies. Microbe 5, 100180. doi:10.1016/j.microb.2024.100180
- Postler T. S., Beer M., Blitvich B. J., Bukh J., de Lamballerie X., Drexler J. F. (2023). Renaming of the genus flavivirus to orthoflavivirus and extension of binomial species names within the family flaviviridae. Arch. Virol. 168 (9), 224. doi:10.1007/s00705-023-05835-1
- Prasad K., Kapoor R., Lee P. (1994). Purpurogallin, a scavenger of polymorphonuclear leukocyte-derived oxyradicals. Mol. Cell Biochem. 139 (1), 27 – 32. doi:10.1007/BF00944200
- Pukelsheim F. (1994). The three sigma rule. Am. Statistician 48 (48), 88 – 91. doi:10.1080/00031305.1994.10476030
- Rabe I. B., Hills S. L., Haussig J. M., Walker A. T., Dos Santos T., San Martin J. L. (2025). A review of the recent epidemiology of zika virus infection. Am. J. Trop. Med. Hyg. 112 (5), 1026 – 1035. doi:10.4269/ajtmh.24-0420
- Ratti A., Fassi E. M. A., Forlani F., Zangrossi M., Mori M., Cappitelli F. (2023). Unlocking the antibiofilm potential of natural compounds by targeting the NADH:quinone oxidoreductase WrbA. Antioxidants (Basel) 12 (8), 1612. doi:10.3390/antiox12081612
- Richard S. A. (2019). Elucidating the novel biomarker and therapeutic potentials of high-mobility group box 1 in subarachnoid hemorrhage: a review. AIMS Neurosci. 6 (4), 316 – 332. doi:10.3934/Neuroscience.2019.4.316
- Rodriguez-Duran L. V., Valdivia-Urdiales B., Contreras-Esquivel J. C., Rodriguez-Herrera R., Aguilar C. N. (2011). Novel strategies for upstream and downstream processing of tannin acyl hydrolase. Enzyme Res. 2011, 823619. doi:10.4061/2011/823619
- Rodriguez-Wells V., Plotch S. J., DeStefano J. J. (2001). Primer-dependent synthesis by poliovirus RNA-Dependent RNA polymerase (3D(pol)). Nucleic Acids Res. 29 (13), 2715 – 2724. doi:10.1093/nar/29.13.2715
- Ryan S. J., Carlson C. J., Mordecai E. A., Johnson L. R. (2019). Global expansion and redistribution of Aedes-borne virus transmission risk with climate change. PLoS Negl. Trop. Dis. 13 (3), e0007213. doi:10.1371/journal.pntd.0007213
- Ryan S. J., Carlson C. J., Tesla B., Bonds M. H., Ngonghala C. N., Mordecai E. A. (2021). Warming temperatures could expose more than 1.3 billion new people to zika virus risk by 2050. Glob. Chang. Biol. 27 (1), 84 – 93. doi:10.1111/gcb.15384
- Sager G., Gabaglio S., Sztul E., Belov G. A. (2018). Role of host cell secretory machinery in zika virus life cycle. Viruses 10 (10), 559. doi:10.3390/v10100559
- Schrier B. K., Wilson S. H. (1973). Investigation of methods for measurement of radioactivity in tritiated DNA and applications to assays for DNA polymerase activity. Anal. Biochem. 56 (1), 196 – 207. doi:10.1016/0003-2697(73)90184-x
- Shatkin A. J., Sipe J. D. (1968). RNA polymerase activity in purified reoviruses. Proc. Natl. Acad. Sci. U. S. A. 61 (4), 1462 – 1469. doi:10.1073/pnas.61.4.1462
- Sierra M., Airaksinen A., Gonzalez-Lopez C., Agudo R., Arias A., Domingo E. (2007). Foot-and-mouth disease virus mutant with decreased sensitivity to ribavirin: implications for error catastrophe. J. Virol. 81 (4), 2012 – 2024. doi:10.1128/JVI.01606-06
- Sirohi D., Kuhn R. J. (2017). Zika virus structure, maturation, and receptors. J. Infect. Dis. 216 (Suppl. l_10), S935 – S944. doi:10.1093/infdis/jix515
- Sivinski J., Ngo D., Zerio C. J., Ambrose A. J., Watson E. R., Kaneko L. K. (2022). Allosteric differences dictate GroEL complementation of E. coli. FASEB J. 36 (3), e22198. doi:10.1096/fj.202101708RR
- Souto J. A., Sarno F., Nebbioso A., Papulino C., Alvarez R., Lombino J. (2020). A new family of jumonji C domain-containing KDM inhibitors inspired by natural product purpurogallin. Front. Chem. 8, 312. doi:10.3389/fchem.2020.00312
- Suzuki T., Suzuki J., Nagata S. (2014). Functional swapping between transmembrane proteins TMEM16A and TMEM16F. J. Biol. Chem. 289 (11), 7438 – 7447. doi:10.1074/jbc.M113.542324
- Tillotson J., Kedzior M., Guimaraes L., Ross A. B., Peters T. L., Ambrose A. J. (2017). ATP-competitive, marine derived natural products that target the DEAD box helicase, eIF4A. Bioorg Med. Chem. Lett. 27 (17), 4082 – 4085. doi:10.1016/j.bmcl.2017.07.045
- Van Veldhoven P. P., Mannaerts G. P. (1987). Inorganic and organic phosphate measurements in the nanomolar range. Anal. Biochem. 161 (1), 45 – 48. doi:10.1016/0003-2697(87)90649-x
- Vardakou M., Salmon M., Faraldos J. A., O'Maille P. E. (2014). Comparative analysis and validation of the malachite green assay for the high throughput biochemical characterization of terpene synthases. MethodsX 1, 187 – 196. doi:10.1016/j.mex.2014.08.007
- Vassiliou W., Epp J. B., Wang B. B., Del Vecchio A. M., Widlanski T., Kao C. C. (2000). Exploiting polymerase promiscuity: a simple colorimetric RNA polymerase assay. Virology 274 (2), 429 – 437. doi:10.1006/viro.2000.0492
- Verzele M., Delahaye P., van Dijck J. (1983). Digallic acid. Bull. Des. Sociétés Chim. Belg. 92 (2), 181 – 186. doi:10.1002/bscb.19830920212
- Veser J. (1987). Kinetics and inhibition studies of catechol O-methyltransferase from the yeast Candida tropicalis. J. Bacteriol. 169 (8), 3696 – 3700. doi:10.1128/jb.169.8.3696-3700.1987
- Woodson S. A., Koculi E. (2009). Analysis of RNA folding by native polyacrylamide gel electrophoresis. Methods Enzymol. 469, 189 – 208. doi:10.1016/S0076-6879(09)69009-1
- Wu Z. L. (2011). Phosphatase-coupled universal kinase assay and kinetics for first-order-rate coupling reaction. PLoS One 6 (8), e23172. doi:10.1371/journal.pone.0023172
- Zhen A. X., Piao M. J., Hyun Y. J., Kang K. A., Ryu Y. S., Cho S. J. (2019). Purpurogallin protects keratinocytes from damage and apoptosis induced by ultraviolet B radiation and particulate matter 2.5. Biomol. Ther. Seoul. 27 (4), 395 – 403. doi:10.4062/biomolther.2018.151
- Zheng F., Mai S., Cen X., Zhao P., Ye W., Ke J. (2024). Discovery of digallic acid as XOD/URAT1 dual target inhibitor for the treatment of hyperuricemia. Bioorg Chem. 147, 107381. doi:10.1016/j.bioorg.2024.107381
Citation
Aitken, V., et al. (2026). Optimization of a malachite green–pyrophosphatase assay for the Zika virus RNA-dependent RNA polymerase and preliminary identification of natural product modulators. Frontiers in Drug Discovery, Volume 6. https://doi.org/
- Received
- Accepted
- Published
- Volume
- 6 - 2026
Copyright
Copyright © 2026 Aitken, Moore, Terrazas-López, Rodarte and Koculi.
License: https://creativecommons.org/licenses/by/4.0/
*Correspondence: Eda Koculi
Article metrics
- Views
- 0
- Downloads
- 0
- Citations
- 0
- Shares
- 0
